[R-sig-ME] "mixed" MANOVAs

[peter dalgaard]

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