Displaying 1 result from an estimated 1 matches for "rhoag".
Did you mean:
hoag
2012 Oct 07
1
Why do I get different results for type III anova using the drop1 or Anova command?
...ts can I rely on? I based my script on:
http://mcfromnz.wordpress.com/2011/03/0 ... explained/
<http://mcfromnz.wordpress.com/2011/03/02/anova-type-iiiiii-ss-explained/>
but I surely made mistakes.
Thanks for help.
Kind regards,
Christiane
#0
library(car)
P<-read.table("rhoPgegenrhoAgewichtet-mit Befruchtungssystem0912.txt",
sep="\t",dec=",", header=T,encoding="latin1")
summary(P)
str(P)
#1
modD3 <- lm(P$rhoPug ~ P$rhoAg*P$dprime*P$Befruchtungssystem, data=P)
Anova(modD3,contrasts=list(rhoAg=contr.sum, dprime=contr.poly),
type="III&...