0

I am trying to get an ANOVA table for my split-split plot design where it will use type III Sum sq instead of type I. This is what I have done so far;

> Attach(Data)
> library(car)
> options(contrasts = c("contr.sum", "contr.poly"))
> mod <- aov(Response ~ A*B*C + Error(Block/A/B/C))
> Anova(mod, type='III')
Error in vcov.default(mod) :
     there is no vcov() method for models of class aovlist, listof

I don't understand why I keep getting this error message, or what to do about it.

Any help trying to get around this error message would be greatly appreciated. Thanks in advance!

Katie
  • 1
  • 1
  • 1
  • As far as I remember, it says an error term isn't allowed. – cuttlefish44 Nov 22 '17 at 20:50
  • Then how would I add in the error term? I need the error term to analyze this as a split-split plot. – Katie Nov 23 '17 at 15:58
  • I think It might be good idea to use `lme` package, such as [https://stats.stackexchange.com/questions/153479/](https://stats.stackexchange.com/questions/153479/) – cuttlefish44 Nov 24 '17 at 02:42

0 Answers0