[R-br] Matrix Modelo

[Eder David Borges da Silva]

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 em 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 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)
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110604/5c29dadd/attachment.html>