<html>
<head>
<meta content="text/html; charset=UTF-8" http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<span style="font-family:monospace,monospace">Walmes,<br>
<br>
Conhece alguma função para automatizar a formação dos
contrastes abaixo:<br>
</span><br>
<br>
<span style="font-family:monospace,monospace">## Se \beta = C\mu
então \mu = C^{-1}\beta.<br>
C <- rbind("mu"=c(1,1,1,1)/4,<br>
"1vs2"=c(1,-1,0,0),<br>
"2vs3"=c(0,1,-1,0),<br>
"3vs4"=c(0,0,1,-1))<br>
<br>
</span><br>
<span style="font-family:monospace,monospace">Pois em um fator com 4
níveis, tudo bem é bem razoável fazer "na mão", porém e se eu
tivesse <span style="font-family:monospace,monospace">um fator
com 40 níveis como em trat <-gl(40, 30)</span></span><span
style="font-family:monospace,monospace"><span
style="font-family:monospace,monospace">, fica um pouco
complicado, então eu gostaria de achar um jeito de criar os
vetores das comparações com maior facilidade,<br>
<br>
<br>
Obrigado,<br>
<br>
</span></span><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>
<span style="font-family:monospace,monospace"><span
style="font-family:monospace,monospace"></span> </span><br>
<br>
<div class="moz-cite-prefix">On 22/03/2015 11:48, walmes . wrote:<br>
</div>
<blockquote
cite="mid:CAFU=EkYZjM-v6a6OqGobKgkvt68ikbyjhjtv8xfnsMPwF4DxrQ@mail.gmail.com"
type="cite">
<div dir="ltr">
<div class="gmail_default" style="font-family:trebuchet
ms,sans-serif">Uma matriz de contraste tem que ter no máximo
k-1 funções e essas devem ser linearmente independentes. Os
contrastes de Tukey formam uma matriz com k*(k-1)/2 hipóteses
(ou funções) lineares. Além do mais, as funções não são
linearmente independentes, o posto é sempre k-1. É por isso
que você não consegue usar os contrastes de Tukey como
definição de restrições na estimação. Abaixo segue CMR para
esclarecer e como montar seus próprios contrastes.<br>
</div>
<div class="gmail_default" style="font-family:trebuchet
ms,sans-serif"><br>
<span style="font-family:monospace,monospace">## Fator de 4
níveis.<br>
trat <- gl(4, 3)<br>
levels(trat)<br>
<br>
## k níveis => k-1 funções/combinações lineares
independentes.<br>
contrasts(trat)<br>
contrasts(C(trat, contr="contr.sum"))<br>
contrasts(C(trat, contr="contr.helmert"))<br>
<br>
Xcontr <- multcomp::contrMat(n=1:nlevels(trat),
type="Tukey")<br>
dim(Xcontr)<br>
<br>
## Posto das matrizes (número de autovalores não nulos).<br>
eigen(crossprod(Xcontr))$values<br>
eigen(crossprod(contrasts(trat)))$values<br>
<br>
## Como implementar contrastes próprios?<br>
## Ex: Contrastes entre pares na sequência (1-2, 2-3, 3-4,
etc).<br>
## Se \beta = C\mu então \mu = C^{-1}\beta.<br>
C <- rbind("mu"=c(1,1,1,1)/4,<br>
"1vs2"=c(1,-1,0,0),<br>
"2vs3"=c(0,1,-1,0),<br>
"3vs4"=c(0,0,1,-1))<br>
K <- solve(C) ## K = C^{-1}.<br>
<br>
contrasts(trat) <- K[,-1] ## Retira o col do intercepto,
deixa k-1.<br>
<br>
y <- rnorm(trat, as.numeric(trat), sd=0.01)<br>
tapply(y, trat, mean)<br>
mean(y)<br>
diff(tapply(y, trat, mean))<br>
<br>
m0 <- lm(y~trat)<br>
coef(m0)<br>
anova(m0)<br>
summary(m0)<br>
</span><br>
</div>
À
<div class="gmail_default" style="font-family:trebuchet
ms,sans-serif;display:inline"> disposição.<br>
</div>
<div class="gmail_default" style="font-family:trebuchet
ms,sans-serif;display:inline">Walmes.</div>
</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>
<br /><br />
<hr style='border:none; color:#909090; background-color:#B0B0B0; height: 1px; width: 99%;' />
<table style='border-collapse:collapse;border:none;'>
<tr>
<td style='border:none;padding:0px 15px 0px 8px'>
<a href="http://www.avast.com/">
<img border=0 src="http://static.avast.com/emails/avast-mail-stamp.png" alt="Avast logo" />
</a>
</td>
<td>
<p style='color:#3d4d5a; font-family:"Calibri","Verdana","Arial","Helvetica"; font-size:12pt;'>
Este email foi escaneado pelo Avast antivírus.
<br><a href="http://www.avast.com/">www.avast.com</a>
</p>
</td>
</tr>
</table>
<br />
</body>
</html>