25 Fev
2015
25 Fev
'15
17:09
Considere o CMR abaixo. da <- data.frame(col=gl(3,100, labels=c("A","B","C"))) da$y <- with(da, rnorm(length(col), mean=as.numeric(col), sd=1)) ## Modelo com a restrição de não haver nível de referência, aqui estima ## a média de cada nível. m0 <- lm(y~0+col, data=da) summary(m0) require(multcomp) ## Matriz de contrastes (lembrar de ter o 0 na formula!). Xcontr <- contrMat(n=1:nlevels(da$col), type="Tukey") str(Xcontr) ## Estimativa e erro padrão. Xcontr%*%coef(m0) sqrt(diag(Xcontr%*%vcov(m0)%*%t(Xcontr))) ## Via glht(). summary(glht(m0, linfct=mcp(col="Tukey"))) À disposição. Walmes.