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

[ASANTOS]

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 em yahoo.com.br
         alexandre.santos em 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 em 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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150226/20b60094/attachment.html>