## Manuel,
## Please look at the maiz example in ?mmc in the HH package.
install.packages("HH") ## if necessary
library(HH)
?mmc
## After the full maiz example, then you need one more command
maiz.proj <- proj(maiz.aov)
maiz.proj
## I think you are looking for
maiz.proj$Within[, "Residuals"]
## Rich
On Wed, Jul 16, 2014 at 9:46 AM, Manuel Lambert
<lambert_man...@hotmail.com> wrote:
>
>
>
> Hello,
> I'm having problems to analyse results from a RCB Design Experiment.I have
> three blocks. For each block: four treatments (factor A), randomized. And for
> each factor A treatment, I have 6 differents treatments (factor C),
> randomized.
> myAOV=aov(response ~ factorA*factorC + Block + Error(field), data)
> The model i use would be response = factorA + factorC + factorA:factorC +
> Block (fixed effect) + Field (subplot, random effect).
> I have some good results but my problem is the following: i would like to
> have a look at the residuals distribution and the command:
> residuals(myAOV)
> returns a NULL object. How to do then to extract the residuals? The function
> works if i take off "Error(field)" from my aov function but this doesn't give
> me relevant residuals right?
> Thanks a lot!
>
>
>
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
