[R-br] Criar comparações múltiplas entre fatores

ASANTOS alexandresantosbr em yahoo.com.br
Quarta Março 9 12:18:13 BRT 2016


Caros listeiros,


        Gostaria de achar uma maneira de automatizar a criação de 
comparações múltiplas entre fatores, no meu caso:

#Fator
treat <- gl(3, 20, labels = paste("t", 1:3, sep=""))

#Variables
sp1_1 = rnorm(10, 5, 0.25)
sp1_2 = rnorm(50, 2.5, 0.25)
sp1 = c(sp1_1,sp1_2)
sp2 = rnorm(60, 2.5, 0.25)
sp3_1 = rnorm(10, 12, 0.25)
sp3_2 = rnorm(50, 2.5, 0.25)
sp3 = c(sp3_1,sp3_2)
sp4 = rnorm(60, 2.5, 0.25)
sp<-cbind(sp1,sp2,sp3,sp4)

#Criação da matriz de contrastes
Treat_Imp<-model.matrix(~ treat-1)

Treat_Imp1_2<-Treat_Imp[, 1] - Treat_Imp[, 2]
Treat_Imp2_3<-Treat_Imp[, 2] - Treat_Imp[, 3]
Treat_Imp1_3<-Treat_Imp[, 1] - Treat_Imp[, 3]

         Eu acabei criando manualmente as comparações, através dos 
objetos Treat_Imp1_2, Treat_Imp2_3 e Treat_Imp1_3. Mas alguém poderia me 
dar uma ideia de como sistematizar isso, e ter a criação desses objetos 
de uma maneira direta e que não envolva pacotes como o multcomp?

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
OrcID: orcid.org/0000-0001-8232-6722
Researchgate: https://www.researchgate.net/profile/Alexandre_Santos10
LinkedIn: https://br.linkedin.com/in/alexandre-dos-santos-87961635
======================================================================



Mais detalhes sobre a lista de discussão R-br