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@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@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.