[R-br] Matrix Modelo
Eder David Borges da Silva
eder em leg.ufpr.br
Sábado Junho 4 16:54:59 BRT 2011
Pessoal,
gostaria da ajuda de vocês para criar uma matrix do modelo, como exemplo
abaixo:
LOCAL <- rep(c('local1','local2'),each=4)
Rep <- 1:8
CLONE <- rep(1:2,4)
BLOCO <- rep(1:4,each=2)
dados <- data.frame(LOCAL,BLOCO,CLONE,Rep)
dados$Int <- paste(dados$LOCAL,dados$CLONE,sep='')
dados$Obs <- 1
dados$TCH <- sort(rnorm(nrow(dados),10,1))
dados$IND <- as.factor(dados$IND)
dados$CLONE <- as.factor(dados$CLONE)
dados$BLOCO <- as.factor(dados$BLOCO)
summary(dados)
op <- options(contrasts=c("contr.treatment", "contr.poly"))
MM <- model.matrix(~BLOCO+CLONE+LOCAL:CLONE,dados)
MM
(Intercept) BLOCO2 BLOCO3 BLOCO4 CLONE2 CLONE1:LOCALlocal2 CLONE2:LOCALlocal2
1 1 0 0 0 0 0 0
2 1 0 0 0 1 0 0
3 1 1 0 0 0 0 0
4 1 1 0 0 1 0 0
5 1 0 1 0 0 1 0
6 1 0 1 0 1 0 1
7 1 0 0 1 0 1 0
8 1 0 0 1 1 0 1
esta foi a matrix criado pelo model.matrix, porem gostaria de uma
matrix diferente, pois nesta criada pela model.matrix o bloco1 não
aparece pois é estimado por direfença, assim como o clone1, eu quero
que todos os niveis aparecem tendo como nome das colunas:
(Intercept) BLOCO1 BLOCO2 BLOCO3 BLOCO4 CLONE1 CLONE2 CLONE1:LOCAL1
CLONE1:LOCAL2 CLONE2:LOCAL1 CLONE2:LOCAL2
Alguem tem uma dica?
Att
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110604/fe7a88c4/attachment.html>
Mais detalhes sobre a lista de discussão R-br