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@yahoo.com.br alexandre.santos@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@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