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.