[R-sig-ME] "mixed" MANOVAs

peter dalgaard pdalgd at gmail.com
Sun Jul 1 11:39:21 CEST 2012

On Jun 30, 2012, at 21:22 , peter dalgaard wrote:

>> m <-  aov(cbind(yield,foo) ~  N*P*K + Error(block), npk)
> Warning message:
> In aov(cbind(yield, foo) ~ N * P * K + Error(block), npk) :
>  Error() model is singular
> (The two do basically the same thing, manova() does little more than adding an extra class to the result.)
> The warning is a bit obscure; it probably just means that something inside aov() didn't expect a matrix response. 

On closer examination, the warning is just a buglet: Inside aov() we should be looking at NROW(er.fit$coefficients), not length(). So the situation seems to be that the feature (multistratum MANOVA) has indeed implemented itself...

Peter Dalgaard, Professor,
Center for Statistics, Copenhagen Business School
Solbjerg Plads 3, 2000 Frederiksberg, Denmark
Phone: (+45)38153501
Email: pd.mes at cbs.dk  Priv: PDalgd at gmail.com

More information about the R-sig-mixed-models mailing list