[R] ANOVA for RCBD, two factors: how to plot residuals?..

[Richard M. Heiberger]

## 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_manuel at 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 at 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.