Walmes, Tenho mais uma questão sobre o problema ainda, pois no caso de eu ter 30 níveis com 10 repetições e emprego a estrutura diagonal proposta por você, eu sempre estou comparando 1com 2, 2 com 3 ate a enésima comparação, teria como eu fazer algo para comparar o nível 1 com todos os outros, o nível 2 com todos os outros assim sucessivamente, pois no meu exemplo proposto tenho: ## Criando 30 níveis com 10 repetições trat<-gl(30,10) ## Se \beta = C\mu então \mu = C^{-1}\beta. ## Número de níveis do fator. k <- 75 ## Row 1. r1 <- rep(1/k, k) ## Matriz diagonal. D <- diag(k-1) ## Positive ones. po <- cbind(D, 0) ## Negative ones. no <- cbind(0, -D) ## Matriz de contrastes "diferenças em sequência". C <- rbind(r1, po+no) rownames(C) <- c("mu", paste0(1:(k-1), "-", 2:k)) C contrasts(trat)<-t(C[-1,]) ## Retira o col do intercepto, deixa k-1. trat.ortho <- model.matrix(~trat)[, -1] head(trat.ortho) Que não equivale a uma comparação completa entre todos os níveis, poderia me ajudar a conseguir a 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 30/03/2015 15:59, walmes . wrote:
Alexandre,
Daria para fazer com um for(), mas aí percebi a estrutura diagonal de 1's positivo e 1's negativos e achei mais fácil somar matrizes diagonais. O CMR esclarece melhor o ponto.
## Número de níveis do fator. k <- 5
## Row 1. r1 <- rep(1/k, k)
## Matriz diagonal. D <- diag(k-1)
## Positive ones. po <- cbind(D, 0)
## Negative ones. no <- cbind(0, -D)
## Matriz de contrastes "diferenças em sequência". C <- rbind(r1, po+no) rownames(C) <- c("mu", paste0(1:(k-1), "-", 2:k)) C
À 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.