30 Mar
2015
30 Mar
'15
19:59
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.