[R-br] Comparar duas regressões no mesmo plot

[Benilton Carvalho]

Testar diferencas entre "reg1" e "reg2" resume-se a ajustar um modelo
unico e testar interacoes e efeitos diferenciais...

Usando o seu exemplo, eu comecaria por determinar:

x = c(x1, x2)
y = c(y1, y2)
z = factor(c(rep(1, length(x1)), rep(2, length(x2))))
fit = lm(y~x*z)

se reg1 e reg2 forem "iguais", entao teriamos a equivalencia com o modelo

fit0 = lm(y~x)

e avaliar a equivalencia seria

anova(fit0, fit)

(note: nao testei o codigo... e o walmes eh a referencia no assunto)

sobre "inserir legenda com significado da reta", eu nao entendi o q vc
quer... mas talvez seja apenas:

legend('bottomright', c('D', 'W'), pch=c(19, 2), lty=2:1)

?

b

2012/4/17 Marcelo Claro de Souza <marcelo_claro em yahoo.com.br>:
> Plotei a seguinte regressão hipotética e gostaria de saber se existe alguma
> forma se verificar se as reg 1 e 2 são significativamente diferentes, talvez
> comparando ângulo de inclinação das retas?
> Tb gostaria de saber como faço p inserir na legenda o que significa cada
> reta, pois só consegui colocar o que significam os pontos.
>
> x1=c(1,2,3,4,5,6,7,8,9,10)
> x2=c(13,14,15,16,17,18,19)
> y1=c(30,40,50,60,70,80,90,100,110,120)
> y2=c(20,22,23,24,25,26,27)
> x=c(1,2,3,4,5,6,7,8,9,10,13,14,15,16,17,18,19)
> y=c(30,40,50,60,70,80,90,100,110,120,20,22,23,24,25,26,27)
> reg1=lm(x1~y1)
> reg2=lm(x2~y2)
> simbolos=c(rep(19,10),rep(2,7))
> plot(x~y,pch=simbolos)
> abline(reg1,lty=2)
> abline(reg2,lty=1)
> legend("bottomright",c("D","W"),pch=c(19,2))
>
> Muito obrigado pela ajuda.
>
> Marcelo Claro de Souza
> Biologist, PhD student in Plant Biology
> Institute of Bioscience - UNESP, Brazil
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.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) e forneça código
> mínimo reproduzível.