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

[Marcelo Claro de Souza]

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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120417/856c00be/attachment.html>