[R-br] Matrix Modelo

[Benilton Carvalho]

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 em 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 em 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 em 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)