<html>
<head>
<meta content="text/html; charset=windows-1252"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Walmes,<br>
<br>
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:<br>
<br>
## Criando 30 níveis com 10 repetições<br>
trat<-gl(30,10)<br>
<br>
## Se \beta = C\mu então \mu = C^{-1}\beta.<br>
## Número de níveis do fator.<br>
k <- 75<br>
<br>
## Row 1.<br>
r1 <- rep(1/k, k)<br>
<br>
## Matriz diagonal.<br>
D <- diag(k-1)<br>
<br>
## Positive ones.<br>
po <- cbind(D, 0)<br>
<br>
## Negative ones.<br>
no <- cbind(0, -D)<br>
<br>
## Matriz de contrastes "diferenças em sequência".<br>
C <- rbind(r1, po+no)<br>
rownames(C) <- c("mu", paste0(1:(k-1), "-", 2:k))<br>
C<br>
<br>
contrasts(trat)<-t(C[-1,]) ## Retira o col do intercepto, deixa
k-1.<br>
<br>
trat.ortho <- model.matrix(~trat)[, -1]<br>
<br>
head(trat.ortho)<br>
<br>
Que não equivale a uma comparação completa entre todos os níveis,
poderia me ajudar a conseguir a fazer isso?<br>
<br>
Obrigado,<br>
<br>
<pre class="moz-signature" cols="72">--
======================================================================
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)
<a class="moz-txt-link-abbreviated" href="mailto:e-mails:alexandresantosbr@yahoo.com.br">e-mails:alexandresantosbr@yahoo.com.br</a>
<a class="moz-txt-link-abbreviated" href="mailto:alexandre.santos@cas.ifmt.edu.br">alexandre.santos@cas.ifmt.edu.br</a>
Lattes: <a class="moz-txt-link-freetext" href="http://lattes.cnpq.br/1360403201088680">http://lattes.cnpq.br/1360403201088680</a>
======================================================================</pre>
<br>
<br>
<br>
<br>
<br>
<div class="moz-cite-prefix">On 30/03/2015 15:59, walmes . wrote:<br>
</div>
<blockquote
cite="mid:CAFU=EkZgK5T-XgxxTL0za0tOb3hyx8gY-2C1bdTzO-w6Jt9j5g@mail.gmail.com"
type="cite">
<div dir="ltr">
<div class="gmail_default" style="font-family:trebuchet
ms,sans-serif">Alexandre,<br>
<br>
</div>
<div class="gmail_default" style="font-family:trebuchet
ms,sans-serif">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.<br>
<br>
<span style="font-family:monospace,monospace">## Número de
níveis do fator.<br>
k <- 5<br>
<br>
## Row 1.<br>
r1 <- rep(1/k, k)<br>
<br>
## Matriz diagonal.<br>
D <- diag(k-1)<br>
<br>
## Positive ones.<br>
po <- cbind(D, 0)<br>
<br>
## Negative ones.<br>
no <- cbind(0, -D)<br>
<br>
## Matriz de contrastes "diferenças em sequência".<br>
C <- rbind(r1, po+no)<br>
rownames(C) <- c("mu", paste0(1:(k-1), "-", 2:k))<br>
C<br>
</span><br>
À disposição.<br>
</div>
<div class="gmail_default" style="font-family:trebuchet
ms,sans-serif">Walmes.</div>
<br>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
R-br mailing list
<a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
<br>
<pre class="moz-signature" cols="72">
</pre>
</body>
</html>