<div dir="ltr"><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Considere o CMR abaixo.<br><br><span style="font-family:monospace,monospace">da <- data.frame(col=gl(3,100, labels=c("A","B","C")))<br>da$y <- with(da, rnorm(length(col), mean=as.numeric(col), sd=1))<br><br>## Modelo com a restrição de não haver nível de referência, aqui estima<br>## a média de cada nível.<br>m0 <- lm(y~0+col, data=da)<br>summary(m0)<br><br>require(multcomp)<br><br>## Matriz de contrastes (lembrar de ter o 0 na formula!).<br>Xcontr <- contrMat(n=1:nlevels(da$col), type="Tukey")<br>str(Xcontr)<br><br>## Estimativa e erro padrão.<br>Xcontr%*%coef(m0)<br>sqrt(diag(Xcontr%*%vcov(m0)%*%t(Xcontr)))<br><br>## Via glht().<br>summary(glht(m0, linfct=mcp(col="Tukey")))<br></span><br>À disposição.<br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Walmes.<br></div>​</div>