[R] Overlaying a moving average curve on top of a barplot
José Santos Alegria
Jose.Alegria at oni.pt
Mon Feb 3 17:33:02 CET 2003
It makes sense and it works!
Thanks and best regards,
José A. S. Alegria
-----Original Message-----
From: Marc Schwartz [mailto:MSchwartz at MedAnalytics.com]
Sent: segunda-feira, 3 de Fevereiro de 2003 4:18 PM
To: José Santos Alegria; r-help at stat.math.ethz.ch
Subject: RE: [R] Overlaying a moving average curve on top of a barplot
>-----Original Message-----
>From: r-help-admin at stat.math.ethz.ch
>[mailto:r-help-admin at stat.math.ethz.ch] On Behalf Of José
>Santos Alegria
>Sent: Monday, February 03, 2003 5:30 AM
>To: r-help at stat.math.ethz.ch
>Subject: [R] Overlaying a moving average curve on top of a barplot
>
>
>I'm using standard barplot (Windows version 1.6.2 of R) to
>represent a certain weekly metric "v" and I would like to
>properly overlay on top of it its moving average "mean.8"
>(window of 8 weeks). I must be doing something wrong since the
>moving average (using "lines") doesn't overlay properly, i.e.,
>both x-scales do not match!
>
>...
>barplot(v[8:length(v)], col=7)
>lines(mean.8[1:length(mean.8)], lty=1, lwd=2, col=2)
>...
>
>How do I make sure that both graphics are in synch as far as
>the x-scale and y-scales are concerned?
>
>Thanks,
>
>José A. S. Alegria
The problem that you are having is that you need to get the bar midpoints from barplot() in order to use those values as the x coordinate values for lines(). barplot() returns the bar midpoints as an invisible vector (or matrix where appropriate).
Modify your code to:
mp <- barplot(v[8:length(v)], col=7)
lines(mp, mean.8[1:length(mean.8)], lty=1, lwd=2, col=2)
This puts the bar midpoints (x axis values) in mp, which you can then use for lines().
HTH,
Marc
More information about the R-help
mailing list