[R-br] Fwd: Enc: model.matrix() (continuação)

walmes . walmeszeviani em gmail.com
Segunda Março 30 16:59:42 BRT 2015


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.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150330/d6b1b3f2/attachment.html>


Mais detalhes sobre a lista de discussão R-br