<html>
<head>
<meta content="text/html; charset=UTF-8" http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Bom dia Pessoal,<br>
<br>
Continuando a discussão sobre a função model.matrix() eu
gostaria de tentar utilizar o objeto criado por model.matrix()
(Treat.ortho <- model.matrix(~-1+dados$avaliacao)) no CRM abaixo
para criar contrastes entre os três níveis de um tratamento, mas a
função contrasts() não funciona diretamente e não quero que ter que
especificar via cbind(c(0,1,0),c(0,0,1)), ou seja, gostaria de criar
um método que operacionalizasse todos os contrastes em função dos
níveis informados sendo:<br>
<br>
#Dados artificiais<br>
colony<-as.factor(sort(rep(c("colony1", "colony2","colony3"),
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>
<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>
<br>
## Cria uma matriz de contraste<br>
##Niveis<br>
levels(avaliacao)<br>
<br>
Treat.ortho <- model.matrix(~-1+dados$avaliacao)<br>
<br>
<br>
## Cria um fator para cada contraste<br>
col1vs2 <- <br>
col1vs3 <- <br>
<br>
Aqui não sei como 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>
<div class="moz-cite-prefix">On 24/02/2015 15:02, walmes . wrote:<br>
</div>
<blockquote
cite="mid:CAFU=EkbSoXJPvugg2kahtG+oHD6gd-104rP_BV6mbw=q1Y_B+g@mail.gmail.com"
type="cite">
<div dir="ltr">
<div class="gmail_default" style="font-family:trebuchet
ms,sans-serif">Até onde sei, a partir de um fator não tem como
você dizer para a model.matrix() que você quer remover colunas
durante a construção. O mais fácil parece ser construir a
matriz e remover o não desejado. Como as colunas tem nomes
fica fácil remover a partir deles.<br>
<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">
</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>