<html>
<head>
<meta content="text/html; charset=UTF-8" http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Caros Listeiros,<br>
<br>
Gostaria de criar uma matriz de contraste com vários
níveis dentro de um fator, conforme o CRM abaixo, quando utilizo 3
níveis, funciona corretamente a função contrasts(), porém quando
utilizo 4 ou mais níveis não funciona, alguém poderia me sugerir uma
solução,<br>
<br>
Obrigado,<br>
<br>
<br>
#Dados artificiais<br>
colony<-as.factor(sort(rep(c("colony1", "colony2","colony3",
"colony4"), 100)))<br>
y1 <- c(rnorm(100,1,0.1),rnorm(100,5,0.1),rnorm(100,3.5,0.1))<br>
y2 <- c(rnorm(100,10,0.3),rnorm(100,11,0.6),rnorm(100,5,0.6))<br>
y3 <- c(rnorm(100,10,2.3),rnorm(100,11,2.6),rnorm(100,11,2.5))<br>
y4 <- c(rnorm(100,5,0.5),rnorm(100,7,0.5),rnorm(100,22,0.5))<br>
y5 <- c(rnorm(100,11,0.1),rnorm(100,11,0.1),rnorm(100,11,0.1))<br>
<br>
#Cria o data frame<br>
avaliacao <- as.factor(colony)<br>
espectro <- cbind(y1,y2,y3,y4,y5)<br>
dados <- data.frame(avaliacao = I(as.matrix(avaliacao)), bands =
I(as.matrix(espectro)))<br>
<br>
## Cria uma matriz de contraste<br>
##Niveis<br>
levels(avaliacao)<br>
Xcontr <- contrMat(n=1:nlevels(avaliacao), type="Tukey")<br>
contrasts(avaliacao) <- Xcontr<br>
<br>
<br>
<br>
<br>
<div class="moz-cite-prefix">On 25/02/2015 13:09, walmes . wrote:<br>
</div>
<blockquote
cite="mid:CAFU=EkYRXfG6qjU+HPbzFgbqkUTiY427ipO3wmNtTHrVkBGxAA@mail.gmail.com"
type="cite">
<div dir="ltr">
<div class="gmail_default" style="font-family:trebuchet
ms,sans-serif">Considere o CMR abaixo.<br>
<br>
<span style="font-family:monospace,monospace">da <-
data.frame(col=gl(3,100, labels=c("A","B","C")))<br>
da$y <- with(da, rnorm(length(col), mean=as.numeric(col),
sd=1))<br>
<br>
## Modelo com a restrição de não haver nível de referência,
aqui estima<br>
## a média de cada nível.<br>
m0 <- lm(y~0+col, data=da)<br>
summary(m0)<br>
<br>
require(multcomp)<br>
<br>
## Matriz de contrastes (lembrar de ter o 0 na formula!).<br>
Xcontr <- contrMat(n=1:nlevels(da$col), type="Tukey")<br>
str(Xcontr)<br>
<br>
## Estimativa e erro padrão.<br>
Xcontr%*%coef(m0)<br>
sqrt(diag(Xcontr%*%vcov(m0)%*%t(Xcontr)))<br>
<br>
## Via glht().<br>
summary(glht(m0, linfct=mcp(col="Tukey")))<br>
</span><br>
À disposição.<br>
</div>
<div class="gmail_default" style="font-family:trebuchet
ms,sans-serif">Walmes.<br>
</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">--
======================================================================
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 />
<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" />
</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>