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.