[R] Split plot analysis problems

Richard M. Heiberger rmh at temple.edu
Wed Jul 22 20:38:35 CEST 2009


Jean-Paul Maalouf wrote:
> Do you have any idea on how can I verify preliminary assumptions in 
> this model (normality of the residuals and variance homogeneity), 
> since R is not able to extract residuals?

Of course, R extracts residuals.  Use the proj() function.  See ?proj 
for the example
to get the projection of an aovlist object onto the terms of a linear model

proj(npk.aovE)


To get the projections into a simple data.frame, use
tmpE <- proj(npk.aovE)
tmpE.df <- do.call("data.frame", tmpE)
tmpE.df

Mark Difford's solution effectively retrieved columns 3 and 10 from tmpE.df

Rich




More information about the R-help mailing list