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
