Benilton,<div>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...)</div>
<div>Att</div><div><br></div><div><br></div><div><br></div><div><br><br><div class="gmail_quote">Em 4 de junho de 2011 17:14, Benilton Carvalho <span dir="ltr"><<a href="mailto:beniltoncarvalho@gmail.com">beniltoncarvalho@gmail.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">a matriz que vc quer pode ser criada... mas vc nao deve faze-lo.<br>
<br>
da forma que vc descreve, a X'X nao sera inversivel e vc nao<br>
conseguira' estimar o modelo.<br>
<br>
use o que o model.matrix() te oferece e defina os contrastes apropriadamente.<br>
<br>
b<br>
<br>
2011/6/4 Eder David Borges da Silva <<a href="mailto:eder@leg.ufpr.br">eder@leg.ufpr.br</a>>:<br>
<div><div></div><div class="h5">> Esqueci, poderia criar esta matrix por partes, porem no problema real os<br>
> niveis são desbalanceados, complicando a confecção na "mão".<br>
><br>
> Em 4 de junho de 2011 16:54, Eder David Borges da Silva <<a href="mailto:eder@leg.ufpr.br">eder@leg.ufpr.br</a>><br>
> escreveu:<br>
>><br>
>> Pessoal,<br>
>> gostaria da ajuda de vocês para criar uma matrix do modelo, como exemplo<br>
>> abaixo:<br>
>> LOCAL <- rep(c('local1','local2'),each=4)<br>
>> Rep <- 1:8<br>
>> CLONE <- rep(1:2,4)<br>
>> BLOCO <- rep(1:4,each=2)<br>
>> dados <- data.frame(LOCAL,BLOCO,CLONE,Rep)<br>
>> dados$Int <- paste(dados$LOCAL,dados$CLONE,sep='')<br>
>> dados$Obs <- 1<br>
>> dados$TCH <- sort(rnorm(nrow(dados),10,1))<br>
>> dados$IND <- as.factor(dados$IND)<br>
>> dados$CLONE <- as.factor(dados$CLONE)<br>
>> dados$BLOCO <- as.factor(dados$BLOCO)<br>
>> summary(dados)<br>
>> op <- options(contrasts=c("contr.treatment", "contr.poly"))<br>
>> MM <- model.matrix(~BLOCO+CLONE+LOCAL:CLONE,dados)<br>
>> MM<br>
>><br>
>>   (Intercept) BLOCO2 BLOCO3 BLOCO4 CLONE2 CLONE1:LOCALlocal2<br>
>> CLONE2:LOCALlocal2<br>
>> 1           1      0      0      0      0                  0<br>
>> 0<br>
>> 2           1      0      0      0      1                  0<br>
>> 0<br>
>> 3           1      1      0      0      0                  0<br>
>> 0<br>
>> 4           1      1      0      0      1                  0<br>
>> 0<br>
>> 5           1      0      1      0      0                  1<br>
>> 0<br>
>> 6           1      0      1      0      1                  0<br>
>> 1<br>
>> 7           1      0      0      1      0                  1<br>
>> 0<br>
>> 8           1      0      0      1      1                  0<br>
>> 1<br>
>><br>
>> esta foi a matrix criado pelo model.matrix, porem gostaria de uma matrix<br>
>> diferente, pois nesta criada pela model.matrix o bloco1 não aparece pois é<br>
>> estimado por direfença, assim como o clone1, eu quero que todos os niveis<br>
>> aparecem tendo como nome das colunas:<br>
>><br>
>>  (Intercept) BLOCO1 BLOCO2 BLOCO3 BLOCO4 CLONE1 CLONE2<br>
>> CLONE1:LOCAL1 CLONE1:LOCAL2 CLONE2:LOCAL1 CLONE2:LOCAL2<br>
>><br>
>> Alguem tem uma dica?<br>
>><br>
>> Att<br>
><br>
</div></div>> _______________________________________________<br>
> R-br mailing list<br>
> <a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
> <a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
><br>
><br>
<font color="#888888"><br>
<br>
<br>
--<br>
Successful people ask better questions, and as a result, they get<br>
better answers. (Tony Robbins)<br>
_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
</font></blockquote></div><br></div>