[R-br] Fwd: Enc: model.matrix()

[walmes .]

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.
​
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150225/c7eddcf1/attachment.html>