Seria isso? str(exp) contrasts(exp$trat)
mat.c <- matrix(c(-1,0,1,2,-1,-1), nc=2)
contrasts(exp$trat) <- mat.c
m0 <- aov(y~trat, data=exp)
summary(m0, split=list("trat"=list("C.orto.1"=1, "C.orto.2"=2)))
À disposição.
Walmes.
_______________________________________________ R-br mailing list [email protected] https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
