Skip to content

extracting effects estimates from aov model with Error()

2 messages · Peter Ho, Brian Ripley

#
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
#
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: