Hello Dieter and everyone,
Thank you for your advice... but I didn't manage to solve my problem... :-(
I actually like the fact that R tells me which ones of my regressions didn't
achieve convergence and gives me a
warning saying that there was a false or singular convergence...
the problem is t
Hi!
I'm comparing the differences in contaminant concentration between 2
different groups of people ( N=36, N=37). When using a simple linear
regression model I found no differences between groups, but when evaluating
the diagnostic plots of the residuals I found my independent variable to
have
christiaan pauw wrote:
Hi Everybody
I need to create a lot of frequency tables with frequencies and percentages
(and cumilative freq and % as well) for a report. freq() in prettyR give
more or less what I need.
I am trying to export the result of freq() to html but the html doesn't look
look the
My question is: Why would you want a data
structure that is clearly not expressive of the
data involved?
Let me guess. You are coming from statistical
software where data are always rectangular.
Patrick Burns
patr...@burns-stat.com
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of "The R
James Widman wrote:
I am trying to duplicate the example by Spencer Graves in the wiki,
using lmer with the Nozzle data.
http://wiki.r-project.org/rwiki/doku.php?id=guides:lmer-tests
However the Chisq value and the fitAB values that are calculated are
different compared to those in the example.
Hello
I want to make some variables with the ifelse-function, but i don't know how
to do it.
I want to make these five variables;
b2$PRRSvac <- ifelse(b2$status=='A' | b2$status=='Aa',1,0)
b2$PRRSdk <- ifelse(b2$status=='B' | b2$status=='Bb',1,0)
b2$sanVac <- ifelse(b2$status=='C' | b2$sta
Hi,
I have very simple question.
How to open/read/write the file with .rda extension which is usually
available in the data folder of any package?
By
Sukhbir Singh Rattan
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing li
hey Joe,
I had a quick look at your code. In this line:
b2$PRRSvac <- ifelse(b2$status=="MS-X, PRRS-pos"|b2$status=="R�d
SPF+Myc+Ap2+DK+Vac",1,0)
some special characters are used. You must use only plain ascii characters
Hope this help.
Cheers,
Patrizio
2009/2/13 joe1985 :
>
> Hello
>
> I want
See ?library, ?data and for direct access of the file ?load
Uwe Ligges
Sukhbir Rattan wrote:
Hi,
I have very simple question.
How to open/read/write the file with .rda extension which is usually
available in the data folder of any package?
By
Sukhbir Singh Rattan
[[alternative HTM
lme(..., control=lmeControl(opt="optim"))
Uwe
Ron Burns wrote:
I am using R 2.7.0 on a linux platform. I am trying to reproduce a 2002
example using lme from the nlme library.
I want to change the otimizer from the default (nlminb) to optim.
Specifically, this is what I am trying to do:
R>
Remko Duursma wrote:
Dear R-helpers,
making a package (windowsXP), that includes a bunch of functions, but
none have examples (all example code is within \dontrun{} blocks).
I do R CMD check Maeswrap, all bits get OK-ed, except:
"Running examples in 'Maeswrap-Ex.R' failed."
What does the
On Thu, 2009-02-12 at 21:32 -0300, milton ruser wrote:
> Dear all,
> I was thinking how much of R capabilities SAS Institute could incorporate on
> SAS support?
>
> Cheers
>
> miltinho
> brazil
Most likely, as many as they currently provide for Tomcat, Postgres,
Apache or other open source prod
Hi
r-help-boun...@r-project.org napsal dne 13.02.2009 10:38:31:
>
> Hello
>
> I want to make some variables with the ifelse-function, but i don't know
how
> to do it.
>
> I want to make these five variables;
>
> b2$PRRSvac <- ifelse(b2$status=='A' | b2$status=='Aa',1,0)
> b2$PRRSdk <- if
Dear R-Community,
I'm currently trying to find a way to conduct a meta-analysis in R.
I would like to analyze data from mostly-cross-sectional survey-studies. The
effect sizes would be correlations.
The R packages "meta" and "rmeta" are, as far as I can see, set up for
analysis with effect sizes
Remko Duursma wrote:
Dear R-helpers,
making a package (windowsXP), that includes a bunch of functions, but
none have examples (all example code is within \dontrun{} blocks).
I do R CMD check Maeswrap, all bits get OK-ed, except:
"Running examples in 'Maeswrap-Ex.R' failed."
When I run the 'Ma
You can readily create a dynamic display for using qqplot and similar functions
in conjunction with either the playwith or TeachingDemos packages.
For example, to investigate the effect of the shape parameter in the skew
normal distribution on its qqplot relative to the normal distribution:
li
I'm glad you are happy with your solution. Just one inline comment
below, expanding on my private reply to you:
Žroutík wrote:
Resolved.
nchar(unlist(strsplit(as.character(n),"\\.")))[2]
in the function:
NumberPrecision <- function(n) {
PocetCyklu <- 0
if(n != round(n)) {
PocetCyklu <- nchar
Sebastian Stegmann schrieb:
[...]
The R packages "meta" and "rmeta" are, as far as I can see, set up for
analysis with effect sizes for differences (i.e. comparison of the
means/odds-ratios of experimental and control group).
Only the function "metagen" from the "meta"-package looks like it
?try
--
David Winsemius
On Feb 13, 2009, at 3:06 AM, lauramorg...@bluewin.ch wrote:
Hello Dieter and everyone,
Thank you for your advice... but I didn't manage to solve my
problem... :-(
I actually like the fact that R tells me which ones of my
regressions didn't achieve convergence a
Dear Sebastian,
Many researchers may transform the Pearson coefficients into Fisher's
z scores first by using
z <- 0.5*log((1+r)/(1-r)).
The standard errors of the Fisher's z scores are z.SE <- 1/sqrt(n-3)
where n are the sample sizes (see
http://en.wikipedia.org/wiki/Fisher_transformation).
Eit
--
Message: 22
Date: Thu, 12 Feb 2009 09:38:10 -0600
From: Dirk Eddelbuettel
Subject: Re: [R] r and c++
To: JC
Cc: r-help@r-project.org
Message-ID: <18836.16994.244794.92...@ron.nulle.part>
Content-Type: text/plain; charset=us-ascii
Javier,
On 12 February 2009 at
Hi
johan...@dsr.life.ku.dk napsal dne 13.02.2009 13:17:06:
> Hi
>
> I tried what you said, but i gives FALSE all the way down! And I think
its
> because there are other words in phrases as well. So if there is 'San'
> in the phrase it has to be TRUE even thoug there are other words before
a
On Thu, Feb 12, 2009 at 3:41 PM, Michael Bibo <
michael_b...@health.qld.gov.au> wrote:
> Mark Ungrin utoronto.ca> writes:
>
> > > * Installing *source* package 'RGtk2' ...
> > > checking for pkg-config... /usr/bin/pkg-config
> > > checking pkg-config is at least version 0.9.0... yes
> > > checkin
Dear lettice gurus,
the code below works just fine to produce a dotplot. However, I am not
successful changing the color of the lines in the legend (auto.key).
If I add col=..., it only changes the color of the letters in the
legend, but not the color of the lines.
Really appreciate any advice on t
Thank you very much, Duncan! Now it is totally clear for me. Reading the 'R
inferno' tract suggested to me in private by a user (nut one can google it),
and wondering about heaven and hell
Best regards, Z
On Fri, Feb 13, 2009 at 1:05 PM, Duncan Murdoch wrote:
I'm glad you are happy with your
Sundar Dorai-Raj gmail.com> writes:
>
> Sorry, that was my lack of understanding on how coplot works.
> Try the following:
>
> coplot(lbxglu~lbxgh|eth, data = reg.dat.5,
>panel= function(...) {
> panel.smooth(...)
> abline(h = 126, col = "red")
> abline(v = 6.5, col = "blue"
Philippe Grosjean wrote:
James Widman wrote:
I am trying to duplicate the example by Spencer Graves in the wiki,
using lmer with the Nozzle data.
http://wiki.r-project.org/rwiki/doku.php?id=guides:lmer-tests
However the Chisq value and the fitAB values that are calculated are
different compare
Thank you very much. Thank you again regarding the suggestion below. I will
give that a shot and I guess I've got my work counted out for me. I counted 45
different distributions.
Is the best way to get a QQPlot of each, to run through producing a data set
for each distribution and then us
Very quick one please team;
How do I return the number of possible unique subsets in a data set (not
combinations i.e. not (1,2,3,4) and (4,3,2,1)) please ?
e.g. In a list of integers 1 to 12, there is 495 sets with 4 elements in
regards
glenn
[[alternative HTML version de
do you mean:
combn(12, 4)
Best,
Dimitris
glenn wrote:
Very quick one please team;
How do I return the number of possible unique subsets in a data set (not
combinations i.e. not (1,2,3,4) and (4,3,2,1)) please ?
e.g. In a list of integers 1 to 12, there is 495 sets with 4 elements i
You might also want to look at the idealized situation:
library(playwith)
library(sn)
playwith(qqnorm(qsn(1:99/100, shape = shape)),
parameters = list(shape = seq(-3, 3, .1)))
On Fri, Feb 13, 2009 at 6:43 AM, Gabor Grothendieck
wrote:
> You can readily create a dynamic display for using
Hi,
I was trying to copy a directory from one location to another. Is there a
command to do this (similar to the file.copy command)?
thanks!
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.
Hi all,
I was trying to copy a directory from one location to another. Is there a
command (similar to file copy command) that will let me do this?
thanks.
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https
The other solutions offered are perfectly workable. Here is a strategy
that is generalizable to other matrix designs (and on checking the
source of upper.tri and lower.tri, it's not surprising that they use
precisely the same strategy):
n <- 9
dm <- matrix(0, nrow=n, ncol=n)
dm[col(dm) >=
Hello
I have a big dataset where i have pulled out some periods of time acording
to a herd-number (see below):
CHR_NR DATO_TEST SALMO_INDEX didiffdatoperiode
113772003-06-20 7.929
113772003-09-01 4.359
113772004-02-24 5
Hi again,
Thanks a lot for all the suggestions. It will take me a little bit to wrap my
head around to understand what is what, though! This will help me quite a bit.
One difference in the result output between you're solution and Mark's solution
is this:
loc sp tot
L1L1 b 60
L2.5
Thank you a lot!!! That was fantastic! It worked perfectly!!!
Laura
Hi Laura,
try using the functions try() and inherits() to "catch" any errors in the loop:
summaryList <- list()
for (i in 1:35)
{
tempSummary <- try(summary(resultList[[i]]), silent = TRUE)
if (inherits(tempSummary, "tr
Hi all,
For my research I have to use a Multinomial Probit model. I saw that
there are two packages, that include a method to estimate my
parameters. The first one is the MNP-package of Imai and van Dyk. The
second one is part of the bayesm-package of Rossi.
The results for both packages are not t
Step 1: Read the R-FAQ 7.31:
Step 2: Ask yourself whether it is likely that such a small value can
be sensibly distinguished from zero using a program that stores values
in double precision.
Step 3: Consider alternate methods, perhaps logarithms or programs
that are equipped to deal with
?list files
... in particular the pattern argument
--
David Winsemius
On Feb 12, 2009, at 3:38 PM, Stropharia wrote:
Thanks a lot Levi. Your code was much shorter and more elegant. With
a few
minor alterations I got this (see below) to work.
Does anyone know if there is a way to automa
Hi,
Ive written some code that fetches data from an Access Database (2003),
processes the data, then saves the modified data back into a table in the
Access database.
It works if I only pass through the code once, but if I put a loop around
the code so that I fetch data from a different source t
Tim Smith wrote:
Hi all,
I was trying to copy a directory from one location to another. Is there a
command (similar to file copy command) that will let me do this?
I suggest to make a system() (or shell, if you are under Windows) call
to the OS in order to copy recursively.
Uwe Ligges
Hi All,
Would be grateful,if anyone can answer my queries.
I need to share code. For example, if I am working in C/C++, I would put
some function declarations in .h files that you would include. In PHP, I
would create files with the common functions in them and then "include()"
them. So far, I
On Fri, 13 Feb 2009, Tim Smith wrote:
I was trying to copy a directory from one location to another. Is
there a command to do this (similar to the file.copy command)?
file.copy() does this in the development version of R. Most people
use
system(paste("cp -r", from, to))
in earlier versi
This is probably the right time to issue a warning about the error of
making transformations on the dependent variable before doing your
analysis. The classic error that newcomers to statistics commit is to
decide that they want to "make their data normal". The assumptions of
most regressio
Hi,
I was wondering if anybody can help me. In the small data set below I would
like to select the index which doesn't contain the letter 'N' in the variable
'letters'. How can I discard these rows when the letter has a different
position everytime (but the same letter for the whole column)?
i
guRus:
I have a variable "beta" as an argument to R's beta function. So
essentially I have a case of beta(alpha, beta). What surprises me is that
R doesn't barf on this stupid programming practice. R gets the right
answer. How does R know "beta the variable" from "beta the function"?
Josep
Hello,
I have to run a general linear mixed model which looks at 2 dependent
variables at the same time (var1 divided by var2). I have tryed to search
for such a kind of model structure but since I just started using R my
search was not successful. Especielly since I only have an old SAS GENMOD
co
Hi Tim,
This is what i usually try and do (see ?file.copy)
## R Start...
> src.dir <- 'PFO-SBS001\\Redirected\\tonyb\\Desktop\\folderA\\'
> dest.dir <- 'PFO-SBS001\\Redirected\\tonyb\\Desktop\\folderB\\'
>
> file.names <- dir(src.dir)
> sapply(file.names, function(x) {
+
one way is:
index <- c(1, 2, 3)
let <- c("CCTTGGAA", "NNTTGGAAT", "AACCTTNN")
z <- data.frame(index, let)
index[-grep("N", let)]
# or
z[-grep("N", let), ]
I hope it helps.
Best,
Dimitris
naomi.duijveste...@ipg.nl wrote:
Hi,
I was wondering if anybody can help me. In the small data set bel
Hello
Currently i'm working on a paper for my study and, on advice of my tutor,
i'm using R.
But now i'm stuck on something.
I've used the bootstrap method to determine if a found value was
significant. Lets say this found value is x, the collection of random
samples is y (n=1e+06)
I've compare
Hello,
I've been trying to use odfWeave and prettyR packages to create documents with
both text and graphs, but so far I
haven't been very lucky...
With the function R2html () in prettyR package, when I try to use a source
file, which works perfectly if I run it
directly form R using "source("f
Moumita Das wrote:
Hi All,
Would be grateful,if anyone can answer my queries.
I need to share code. For example, if I am working in C/C++, I would put
some function declarations in .h files that you would include. In PHP, I
would create files with the common functions in them and then "inc
Dear R experts,
I have a problem with a function I wrote. The fuction looks like this:
series<-function(x,s){
foo<-list();
ind3<-integer();
for (j in diff){
for (i in 1:(n-12)){
if
(!(x[i,j]==0)&!(x[i+1,j]==0)&!(x[i+2,j]==0)&!(x[i+3,j]==0)&!(x[i+4,j]==0)&!(x[i+5,j]==0)&!(x[i+6,j]==0)
&!(x[i+
See ?get and try:
sin <- 1
get("sin", mode = "numeric")
get("sin", mode = "function")
On Fri, Feb 13, 2009 at 10:20 AM, wrote:
> guRus:
>
> I have a variable "beta" as an argument to R's beta function. So
> essentially I have a case of beta(alpha, beta). What surprises me is that
> R doesn't
Thank you everyone for the help.
Dieter's advice resolved the issue.
charis
Dieter Menne wrote:
>
> charis kaskiris.com> writes:
>
>>
>>
>> Thank you for helping.
>>
>> I am running the 32-bit version of R on a 64-bit Windows XP machine.
>> After
>> reinstalling the rJava package I starte
You really should not assign a new value to "letters" since it is a
very useful constant vector that allows working with lowercase
letters. Notice that Domitris did not copy your code exactly.
Next time you start R type letters at the R prompt. You cannot do so
now, since you overwrote that
On Fri, Feb 13, 2009 at 10:47 AM, Gabor Grothendieck
wrote:
> See ?get and try:
Interesting. I hadn't paid attention to the 'mode' argument before.
Where would it be advisable to use anything but mode='any' or mode='function'?
-s
__
R-help@r
Why do you care what distribution your data comes from?
That is a serious question, the more we know about what your actual
question/goal is, the more we can help. It is a common mistake for people who
know enough statistics to be dangerous to focus on the distribution of the data
rather than
Hi ,When I try to load the 'mgcv" package, often, but not always, get this
error message. What am I doing wrong?
__
This is mgcv 1.4-1.1
Error in runif(1) :
.Random.seed is not an integer vector but of type 'list'
Error : .onAttach failed in 'attachNamespace'
Error: package/namesp
No example and I find the specification unclear, but perhaps you can
follow the steps outlined in this earlier posting to get closer what
you want. or, as always, provide a base example and
description of what is needed.
http://finzi.psych.upenn.edu/R/Rhelp08/archive/154176.html
-
One possibility is something along the lines of:
plot(density(bootstraps))
abline(v=original.value)
Patrick Burns
patr...@burns-stat.com
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of "The R Inferno" and "A Guide for the Unwilling S User")
TomPoes wrote:
Hello
Currently i'm working o
Stavros Macrakis wrote:
> On Fri, Feb 13, 2009 at 10:47 AM, Gabor Grothendieck
> wrote:
>
>> See ?get and try:
>>
>
> Interesting. I hadn't paid attention to the 'mode' argument before.
>
> Where would it be advisable to use anything but mode='any' or mode='function'?
>
>
don't know i
Have you created an object called .Random.seed? The error looks to be
triggered by the mgcv startup code calling the RNG via `runif', and the RNG
finding something unexpected in .Random.seed. (See ?.Random.seed for more
information).
best,
Simon
On Friday 13 February 2009 16:25, Veerappa Ch
Hello! I have encountered a really weird problem. Maybe you've
encountered it before?
I have a large data frame "importances". It has one factor ($A) with 3
levels: 3, 9, and 15. $B is a regular numeric variable.
Below I am picking a really small sub-frame (just 3 rows) based on
"indices". "indices
lauramorgana bluewin.ch bluewin.ch> writes:
>
> Hello,
> I've been trying to use odfWeave and prettyR packages to create documents with
both text and graphs, but so far I haven't been very lucky...
...
> library(odfWeave) I get this message:
>
> Loading: XML
> Error: package 'XML' not loaded
Hi
I'm trying to fit a model in betareg and I'm getting errors, but have no
idea what they mean or how to solve them. Does anyone have experience with
this?
> model <- betareg(ACT ~ ST*SoilT, data = actDL_F)
Warning messages:
1: In sqrt(W) : NaNs produced
2: In sqrt(W) : NaNs produced
3: In sqrt
Think of the levels as a table you are going to index into. The
factors that you see (2, 1, 3) are the indices into the levels so you
get 9, 3, 15 as the result.
What were you expecting? It is working as it is supposed to.
On Fri, Feb 13, 2009 at 12:09 PM, Dimitri Liakhovitski wrote:
> Hello!
on 02/13/2009 11:09 AM Dimitri Liakhovitski wrote:
> Hello! I have encountered a really weird problem. Maybe you've
> encountered it before?
> I have a large data frame "importances". It has one factor ($A) with 3
> levels: 3, 9, and 15. $B is a regular numeric variable.
> Below I am picking a real
G'day Stavros,
On Fri, 13 Feb 2009 11:11:28 -0500
Stavros Macrakis wrote:
> On Fri, Feb 13, 2009 at 10:47 AM, Gabor Grothendieck
> wrote:
> > See ?get and try:
>
> Interesting. I hadn't paid attention to the 'mode' argument before.
>
> Where would it be advisable to use anything but mode='an
On Fri, Feb 13, 2009 at 12:24 PM, Marc Schwartz
wrote:
> on 02/13/2009 11:09 AM Dimitri Liakhovitski wrote:
>> Hello! I have encountered a really weird problem. Maybe you've
>> encountered it before?
>> I have a large data frame "importances". It has one factor ($A) with 3
>> levels: 3, 9, and 15.
Dear friends,
i've a little problem with a study about mediation. I will call my dependent
variable "Y", the independent "X" and "Z" will be my mediator variable...
Using a linear regression, i can say that there isn't a direct effect of X on Y
because the coefficient of X isn't significative
Combining the various approaches on the list, here's a simple
one-liner that puts the NAs at the end:
t(apply(mat,1,function(r) { dr<-duplicated(r); c( r[!dr],
rep(NA,sum(dr)) ) ))
If you don't care where the NAs are, the following is a tad shorter
and perhaps clearer:
mat[ t(apply(mat
On Fri, Feb 13, 2009 at 12:25 PM, Berwin A Turlach
wrote:
> On Fri, 13 Feb 2009 11:11:28 -0500
> Stavros Macrakis wrote:
>> Where would it be advisable to use anything but mode='any' or
>> mode='function'?
>
> I guess the answer to this question is more often than not in the source. :)
...
> The
Sorry - one clarification:
When I run:
> test$xx - the what I am currently seeing is:
[1] 9 3 15
Levels: 3 9 15
But what I am expecting to be seeing is:
[1] 9 3 15
Levels: 9 3 15
Or maybe: Levels: 2 1 3
On Fri, Feb 13, 2009 at 12:38 PM, Dimitri Liakhovitski wrote:
> On Fri, Feb 13, 2009
Stavros Macrakis wrote:
> On Fri, Feb 13, 2009 at 12:25 PM, Berwin A Turlach
> wrote:
>
>> On Fri, 13 Feb 2009 11:11:28 -0500
>> Stavros Macrakis wrote:
>>
>>> Where would it be advisable to use anything but mode='any' or
>>> mode='function'?
>>>
>> I guess the answer to this quest
Perhaps you could look at the coin package's symmetry tests and see if
they satisfy:
http://finzi.psych.upenn.edu/R/library/coin/html/SymmetryTests.html
--
David Winsemius
On Feb 11, 2009, at 11:26 AM, joseph.g.bo...@gsk.com wrote:
Does anyone know if any R package has a function that will
It comes down to 2 simple rules:
1. If you don't care about the order of the factor levels, then it doesn't
matter how R codes the relationship
2. If you do care about the order, then tell R what order you want.
Consider the following:
> x <- c(9,3,15,9,15,9,3)
> factor(x)
[1] 9 3 15 9 15
Berwin A Turlach wrote:
> G'day Stavros,
>
> On Fri, 13 Feb 2009 11:11:28 -0500
> Stavros Macrakis wrote:
>
>
>> On Fri, Feb 13, 2009 at 10:47 AM, Gabor Grothendieck
>> wrote:
>>
>>> See ?get and try:
>>>
>> Interesting. I hadn't paid attention to the 'mode' argument before.
>>
>>
(typos corrected)
Combining the various approaches on the list, here's a simple
one-liner that puts the NAs at the end:
t(apply(mat,1,function(r) { dr<-duplicated(r); c( r[!dr],
rep(NA,sum(dr)) ) }))
If you don't care where the NAs are, the following is a tad shorter
and perhaps clearer:
Dear R-users,
I am very enthusiastic about the capacities offered by the np package and
pretty fond of it.
Nevertheless, trying to estimate Ichimura and Klein and Spady models on my
data, I would like to be able to provide the npindex function with my
guesses for the bandwidth (eventually compu
on 02/13/2009 11:38 AM Dimitri Liakhovitski wrote:
> On Fri, Feb 13, 2009 at 12:24 PM, Marc Schwartz
> wrote:
>> on 02/13/2009 11:09 AM Dimitri Liakhovitski wrote:
>>> Hello! I have encountered a really weird problem. Maybe you've
>>> encountered it before?
>>> I have a large data frame "importanc
Jason,
Just to answer your direct question, there is Mathowrld.wolfram.com,
where there are 87 continuous distributions listed.
I have also used the book Statistical Distributions, 2nd ed, Merran Evans, et
al.
which has most of the usual distributions with pictures and relationships.
Of course al
Hi all,
I'm tryint to turn my ANOVA result from anova(model) in to p values. The
default output is:
term ; Df ; Deviance ; resid Df; Resid. Dev
Is there a way to add p values to this output so the significance of the
terms can be easily determined?
Thanks in advance.
Richie
[[alternativ
That is very nice. Maybe just one slight improvement so
to express it in a non-destructive form:
replace(mat, t(apply(mat,1,duplicated)), NA)
On Fri, Feb 13, 2009 at 12:58 PM, Stavros Macrakis
wrote:
> (typos corrected)
>
> Combining the various approaches on the list, here's a simple
> on
It seems that SAS is willing to offer this interface if you are
willing to purchase one of their expensive add-on packages. Not surprising.
Gerard Smits
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.e
Both Greg and Marc - thank you so much!
It helped a lot. What I just discovered also works (similar to Greg's
suggestions) is to make it first a character and THEN to do:
as.factor(as.numeric(original character vector))).
Wow! R never stops surprizing one - and I am just in the beginning of
the j
Richie -
There is a test= argument that can be set to "Chisq", "F" or
"Cp", for various different tests. See the help file for
anova.glm for details (or look at anova.xxx if your model is
of class "xxx").
- Phil
On Fri, 13 Feb 2009, R User R
It depends on what type of model "model" is.
In some cases you can do something like
anova(model, test='Chisq')
See the appropriate help to find out the options, but note that if the p-values
are not provided by default, it could be because there is doubt about the
accuracy of the approximation
Hi,
We are trying to calculate ML-parameter estimates of a mixed effects models
where the observations are weighted. The "weights" option in lmer() with
RELM=FALSE seems attractive. Does anyone know the mechanism it uses to
calculate weighted ML estimates? (Is there a paper?). Also, our model
i
hi Im having some problems reshaping
Ive managed to apply it but have some problems
the attached document will explain
any help is appreciated
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the po
MCI
A1
A2
A13
A14
A23
A24
A33
A34
Grouped together
56766
N/A
N/A
N/A
N/A
N/A
N/A
N/A
N/A
N/A
6459
N/A
N
Hello fellows: I've een trying to set up a function that performs 100
loops producing the coresponding 100 series. I want to save all those
datasets in a dataframe, so I wrote this...
prep <- function()
# Clase[1]/Categoria[2]/Phi[3]/Rf[4]
peso <- c(.0,.03,.3,.6)
# Extension del calculo
Dimitri Liakhovitski gmail.com> writes:
> the code below works just fine to produce a dotplot. However, I am not
> successful changing the color of the lines in the legend (auto.key).
> If I add col=..., it only changes the color of the letters in the
> legend, but not the color of the lines.
I
Hello, everybody.
A student asked me a howto question I can't answer. We want the
length of the drawn axes to fill the full width and height of the
plot, like so:
|
| *
| *
| *
---|-
However, when we use plot with axes=F and then use the axis co
Try defining the argument 'pos' in the axis command-line, like:
x <- rnorm(100)
z <- gl(2,50)
y <- rnorm(100, mean= 1.8*as.numeric(z))
plot(x,y,type="n", axes=F)
points(x,y, pch="$",cex=0.7, col=z)
axis(1, col="green", col.axis="green", pos=0)
axis(2, col="red", col.axis="red", pos=-2)
_
on 02/13/2009 01:25 PM Paul Johnson wrote:
> Hello, everybody.
>
> A student asked me a howto question I can't answer. We want the
> length of the drawn axes to fill the full width and height of the
> plot, like so:
>
>|
>| *
>| *
>| *
> ---|
Something like this might work where you are using 'lapply' to create
a list of values from the equivalent of the 'for' loop that you had.
You can then 'cbind' them to create a return matrix, or you can
convert it to a dataframe:
prep <- function()
{
# Clase[1]/Categoria[2]/Phi[3]/Rf[4]
peso <-
Nicole Schneider wrote:
Hello,
I have to run a general linear mixed model which looks at 2 dependent
variables at the same time (var1 divided by var2). I have tryed to search
for such a kind of model structure but since I just started using R my
search was not successful. Especielly since I only
On Fri, Feb 13, 2009 at 1:51 PM, Marc Schwartz
wrote:
> on 02/13/2009 01:25 PM Paul Johnson wrote:
>> Hello, everybody.
>>
>> A student asked me a howto question I can't answer. We want the
>> length of the drawn axes to fill the full width and height of the
>> plot, like so:
> Paul,
>
> I am gu
1 - 100 of 130 matches
Mail list logo
