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 <[email protected]>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 <[email protected]>: > > 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 > > > > > > > > -- > Successful people ask better questions, and as a result, they get > better answers. (Tony Robbins) > _______________________________________________ > R-br mailing list > [email protected] > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >
_______________________________________________ R-br mailing list [email protected] https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
