Muito obrigado Walmes, Era exatamente isto o que eu queria, fazendo o contraste manualmente e extrair a estimativa da média de cada nível sem usar a função glht() do pacote multicomp, Redobrados agradecimentos, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ====================================================================== On 25/02/2015 13:09, walmes . wrote:
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.
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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.
--- Este email foi escaneado pelo Avast antivírus. http://www.avast.com