[R-br] Fwd: Enc: model.matrix()

[ASANTOS]

Bom dia Pessoal,

             Continuando a discussão sobre a função model.matrix() eu 
gostaria de tentar utilizar o objeto criado por model.matrix() 
(Treat.ortho <- model.matrix(~-1+dados$avaliacao)) no CRM abaixo para 
criar contrastes entre os três níveis de um tratamento, mas a função 
contrasts() não funciona diretamente e não quero que ter que especificar 
via cbind(c(0,1,0),c(0,0,1)), ou seja, gostaria de criar um método que 
operacionalizasse todos os contrastes em função dos níveis informados sendo:

#Dados artificiais
colony<-as.factor(sort(rep(c("colony1", "colony2","colony3"), 100)))
y1 <- c(rnorm(100,1,0.1),rnorm(100,5,0.1),rnorm(100,3.5,0.1))
y2 <- c(rnorm(100,10,0.3),rnorm(100,11,0.6),rnorm(100,5,0.6))
y3 <- c(rnorm(100,10,2.3),rnorm(100,11,2.6),rnorm(100,11,2.5))
y4 <- c(rnorm(100,5,0.5),rnorm(100,7,0.5),rnorm(100,22,0.5))
y5 <- c(rnorm(100,11,0.1),rnorm(100,11,0.1),rnorm(100,11,0.1))


#Cria o data frame
avaliacao <- as.factor(colony)
espectro <- cbind(y1,y2,y3,y4,y5)
dados <- data.frame(avaliacao = I(as.matrix(avaliacao)), bands = 
I(as.matrix(espectro)))


## Cria uma matriz de contraste
##Niveis
levels(avaliacao)

Treat.ortho <- model.matrix(~-1+dados$avaliacao)


## Cria um fator para cada contraste
col1vs2 <-
col1vs3 <-

    Aqui não sei como fazer isso?

Obrigado,

-- 
======================================================================
Alexandre dos Santos
Proteção Florestal
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso
Campus Cáceres
Caixa Postal 244
Avenida dos Ramires, s/n
Bairro: Distrito Industrial
Cáceres - MT                      CEP: 78.200-000
Fone: (+55) 65 8132-8112 (TIM)   (+55) 65 9686-6970 (VIVO)
e-mails:alexandresantosbr em yahoo.com.br
         alexandre.santos em cas.ifmt.edu.br
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================




On 24/02/2015 15:02, walmes . wrote:
> Até onde sei, a partir de um fator não tem como você dizer para a 
> model.matrix() que você quer remover colunas durante a construção. O 
> mais fácil parece ser construir a matriz e remover o não desejado. 
> Como as colunas tem nomes fica fácil remover a partir deles.
>
> À disposição.
> Walmes.
> ​
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.



---
Este email foi escaneado pelo Avast antivírus.
http://www.avast.com
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150225/b263787a/attachment.html>