Peter Ho
2003-Jun-25 21:36 UTC
[R] extracting effects estimates from aov model with Error()
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
Prof Brian Ripley
2003-Jun-26 06:30 UTC
[R] extracting effects estimates from aov model with Error()
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