[R-br] Dificuldade em montar matriz de contraste com sapply()

[ASANTOS]

Caros Listeiros,

           Estou tentando montar um matriz de contraste entre 3 níveis 
(colony1, colony2 e colony3) dentro de um fator chamado P1 com a função 
sapply, mas não consigo ver onde estou errando, pois os coeficientes 
entre as combinações de níveis não estão ficando corretos, alguém 
poderia dar uma ajuda, segue CRM:

## Cria uma matriz de contraste
# Níveis dentro de um fator P1
P1 <- c("colony1", "colony2", "colony3")
dados <- data.frame(P1)
comb <- cbind(combn(unique(c(P1,P1)), 2)[2:1,])## Combinação entre os níveis
comb
# Comparação entre níveis com sapply()
sapply(1:ncol(comb), function(x) with(dados,
       2*( (P1==comb[1,x] & P1 == comb[2,x]) | (P1==comb[1,x] & P1 == 
comb[2,x])) -
        (P1 == P1 & P1 %in% comb[,x]) ) )
  dados[ , 2:4] <- sapply(1:ncol(comb), function(x) with(dados,
       2*( (P1==comb[1,x] & P1 == comb[2,x]) | (P1==comb[1,x] & P1 == 
comb[2,x])) -
        (P1 == P1 & P1 %in% comb[,x]) ) )
  dados## Os coeficientes não estão corretos

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
======================================================================


---
Este email foi escaneado pelo Avast antivírus.
http://www.avast.com