Esqueci, poderia criar esta matrix por partes, porem no problema real os niveis são desbalanceados, complicando a confecção na "mão".
Em 4 de junho de 2011 16:54, Eder David Borges da Silva <[email protected]>escreveu: > 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 > >
_______________________________________________ R-br mailing list [email protected] https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
