O que tu considera "<span style="font-family:arial,sans-serif;font-size:13px">valor estimado do nível default criado pelo R"? Quando a gente trabalha com variaveis categóricas, sempre devemos 'sacrificar' uma categoria e trabalhar com as outras em relação a ela. No caso abaixo temos que o grupoH tem uma média 8,08 maior que o grupoC. A estimativa da média do grupoC (condicionada a uma dose) pode ser obtida vendo o intercepto e a inclinação. Segue também um exemplo de gráfico.</span><div>
<span style="font-family:arial,sans-serif;font-size:13px"><br></span></div><div><div><font face="arial, sans-serif"># Simula os dados</font></div><div><font face="arial, sans-serif">set.seed(31415)</font></div><div><font face="arial, sans-serif">dose <- rep(seq(.1,1,.1),each=20)</font></div>
<div><font face="arial, sans-serif">grupo <- gl(2,1,200,labels=c('Grupo C','Grupo H'))</font></div><div><font face="arial, sans-serif">pthi <- .7 + dose*12 + as.real(grupo)*8 + rnorm(200)</font></div>
<div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif"># Ajusta o modelo</font></div><div><font face="arial, sans-serif">(modCH<-lm(pthi~dose+grupo))</font></div><div><font face="arial, sans-serif">summary(modCH)</font></div>
<div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif"># Gráfico</font></div><div><font face="arial, sans-serif">plot(dose,pthi)</font></div><div><font face="arial, sans-serif">abline(coef(modCH)[1],coef(modCH)[2],col=1) # Para o grupo C</font></div>
<div><font face="arial, sans-serif">abline(coef(modCH)[1]+coef(modCH)[3],coef(modCH)[2],col=2) # Para o grupo H</font></div><div><font face="arial, sans-serif">legend('topleft',c('Grupo C','Grupo H'),col=1:2,lty=1)</font></div>
<div style="font-family:arial,sans-serif;font-size:13px"><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">2012/11/30  <span dir="ltr"><<a href="mailto:alanarocha@sapo.pt" target="_blank">alanarocha@sapo.pt</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Presados,<br>
tenho este modelo como os dois níveis<br>
como sei o valor estimado do nível default criado pelo R?<br>
modCH<-lm(pthi~dose+grupo)<br>
modCH<br>
summary(modCH)<br>
outra questão é que sei que abline() é para modeos do tipo y~x<br>
eu queria um exemplo de como fazer as rectas para o meu modelo que envolve factor<br>
desde já obrigada<br>
Ana<br>
<br>
______________________________<u></u>_________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/<u></u>cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-<u></u>guia</a>) e forneça código mínimo reproduzível.<br>
</blockquote></div><br></div>