Benilton, Eu não usarei ela para estimar o modelo, eu estimei o modelo com a lmer, e tenho os efeitos, o que querro e representar o valor pelos seus componentes, multiplicando os efeito pela matrix criada, tendo desta forma na soma de linha (efeito de mu+bloco+gen+gen:local...) Att Em 4 de junho de 2011 17:14, Benilton Carvalho <beniltoncarvalho@gmail.com>escreveu:
a matriz que vc quer pode ser criada... mas vc nao deve faze-lo.
da forma que vc descreve, a X'X nao sera inversivel e vc nao conseguira' estimar o modelo.
use o que o model.matrix() te oferece e defina os contrastes apropriadamente.
b
2011/6/4 Eder David Borges da Silva <eder@leg.ufpr.br>:
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 < eder@leg.ufpr.br> 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 R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
-- Successful people ask better questions, and as a result, they get better answers. (Tony Robbins) _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br