[R-br] modelo e gráfico

[Rodrigo Coster]

O que tu considera "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.

# Simula os dados
set.seed(31415)
dose <- rep(seq(.1,1,.1),each=20)
grupo <- gl(2,1,200,labels=c('Grupo C','Grupo H'))
pthi <- .7 + dose*12 + as.real(grupo)*8 + rnorm(200)

# Ajusta o modelo
(modCH<-lm(pthi~dose+grupo))
summary(modCH)

# Gráfico
plot(dose,pthi)
abline(coef(modCH)[1],coef(modCH)[2],col=1) # Para o grupo C
abline(coef(modCH)[1]+coef(modCH)[3],coef(modCH)[2],col=2) # Para o grupo H
legend('topleft',c('Grupo C','Grupo H'),col=1:2,lty=1)



2012/11/30 <alanarocha em sapo.pt>

> Presados,
> tenho este modelo como os dois níveis
> como sei o valor estimado do nível default criado pelo R?
> modCH<-lm(pthi~dose+grupo)
> modCH
> summary(modCH)
> outra questão é que sei que abline() é para modeos do tipo y~x
> eu queria um exemplo de como fazer as rectas para o meu modelo que envolve
> factor
> desde já obrigada
> Ana
>
> ______________________________**_________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/**cgi-bin/mailman/listinfo/r-br<https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br>
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-**guia<http://www.leg.ufpr.br/r-br-guia>)
> e forneça código mínimo reproduzível.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121130/98899902/attachment.html>