Re: [R] Doubt about nested aov output

2005-09-07 Thread John Wilkinson \(pipex\)
Ronaldo , It looks as though you have specified you model incorrectly. In the Rats example ,the Treatment is the only fixed effect,Rat and Liver are random effects In aov testing for sig of 'Means' of Random Effects is pointless and that is why 'p' values are not given.Further more the interacti

Re: [R] Doubt about nested aov output

2005-09-06 Thread Douglas Bates
On 9/6/05, Ronaldo Reis-Jr. <[EMAIL PROTECTED]> wrote: > Hi Spencer, > > Em Dom 04 Set 2005 20:31, Spencer Graves escreveu: > > Others may know the answer to your question, but I don't. However, > > since I have not seen a reply, I will offer a few comments: > > > > 1. What versi

Re: [R] Doubt about nested aov output

2005-09-06 Thread Ronaldo Reis-Jr.
Hi Spencer, Em Dom 04 Set 2005 20:31, Spencer Graves escreveu: > Others may know the answer to your question, but I don't. However, > since I have not seen a reply, I will offer a few comments: > > 1. What version of R are you using? I just tried superficially > similar things w

Re: [R] Doubt about nested aov output

2005-09-05 Thread Spencer Graves
That's great. In your example, the aov fit is following by a warning: > res2 <- aov(y ~ fix + Error(fix/R1/R2), data = dat) Warning message: Error() model is singular in: aov(y ~ fix + Error(fix/R1/R2), data = dat) Moreover, to test whether your change was statistically sign

Re: [R] Doubt about nested aov output

2005-09-05 Thread Christoph Buser
Hi I think the problem is related to specifying "Tratamento" both as a fixed factor and in the error term. I attached a script with a reproducible example that shows a similar output. I do not know the details of the original data and the questions of interest, but maybe a model including "Tratame

Re: [R] Doubt about nested aov output

2005-09-04 Thread Spencer Graves
Others may know the answer to your question, but I don't. However, since I have not seen a reply, I will offer a few comments: 1. What version of R are you using? I just tried superficially similar things with the examples in ?aov in R 2.1.1 patched and consistently got F

[R] Doubt about nested aov output

2005-08-30 Thread Ronaldo Reis-Jr.
Hi, I have two doubts about the nested aov output. 1) I have this: > anova.ratos <- aov(Glicogenio~Tratamento+Error(Tratamento/Rato/Figado)) > summary(anova.ratos) Error: Tratamento Df Sum Sq Mean Sq Tratamento 2 1557.56 778.78 Error: Tratamento:Rato Df Sum Sq Mean Sq F