[R-br] Matrix Modelo

[Eder David Borges da Silva]

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
>
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110604/5fb4bab9/attachment.html>