
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@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ====================================================================== --- Este email foi escaneado pelo Avast antivírus. http://www.avast.com