[R] extracting effects estimates from aov model with Error()

[Prof Brian Ripley]

I think you are applying effects to the summary, but you would need to 
apply it to each stratum in turn, e.g.  effects(fit[[2]])

Consider

example(aov)
effects(npk.aov[[3]])


On Wed, 25 Jun 2003, Peter Ho wrote:

> Hi all,
> 
> Does anyone know how I might be able to extract effects estimates from 
> an unreplicated blocked factorial design?
>  My aov model was the following:
> 
> aov(AM ~ Tf*Tm*Pl*Ps*Mc*M0*M1*M2*M3+Error(Batch), initial)
> 
> Using effects() , I get the following response:
> 
>  > effects(CV.unreplicated.block.aov.summary)
> Error in effects(CV.unreplicated.block.aov.summary) :
>         no applicable method for "effects"
> 
> Thanks
> 
> Peter
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://www.stat.math.ethz.ch/mailman/listinfo/r-help
> 

-- 
Brian D. Ripley,                  ripley at stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272866 (PA)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595