Marc Carpentier
2009-Mar-12 10:53 UTC
[R] problem with 'abline' in a regression with repeated measures
Hi everybody. I'm just learning to use R. I've got a strange problem and I couldn't find any answer to my problem. So here am I. I want to visualize a model of regression on repeated measures : 5 measures (y) by x-value (lx for log10(x) ; lxf colonn is for as.factor) ; 5 x-values. I fail to use 'abline' command (regression's line don't display) in this particular context, on a 'stripchart' display. Any idea ? Here are the few commands that I execute, the lattest being my problem :> tab <- read.csv2("ex3.csv") > sapply(tab,data.class)sujet x xf lx lxf y "numeric" "numeric" "numeric" "numeric" "numeric" "numeric"> attach(tab) > lxf <- as.factor(lxf)> model1 <- lm(y~lx) > model2 <- lm(y~lxf) > model3 <- lm(y~lx+lxf)> ym <- tapply(y,lxf,mean) > ym-0.9 -0.6 -0.3 0 0.3 43.0 60.6 77.2 92.2 109.4> stripchart(y~lx, method="jitter", jitter=.1, vertical=T, pch=20) > lines (1:5, ym, type="b", pch=4, cex=2, lty=2) > abline(model1) ## I don't understand why it doesn't work ; neither do abline(coef(model1) or abline(92,45)## It usually works on any regression, but stripchart maybe is a special plot ; I suspect an absciss misunderstanding by abline... Thanks for any help. I simulated this data set : lx lxf y -0,9 -0,9 47 -0,9 -0,9 41 -0,9 -0,9 44 -0,9 -0,9 38 -0,9 -0,9 45 -0,6 -0,6 58 -0,6 -0,6 63 -0,6 -0,6 58 -0,6 -0,6 58 -0,6 -0,6 66 -0,3 -0,3 78 -0,3 -0,3 79 -0,3 -0,3 78 -0,3 -0,3 77 -0,3 -0,3 74 0 0 96 0 0 91 0 0 97 0 0 85 0 0 92 0,3 0,3 109 0,3 0,3 114 0,3 0,3 111 0,3 0,3 104 0,3 0,3 109 [[alternative HTML version deleted]]