[R] Overlaying a moving average curve on top of a barplot

[José Santos Alegria]

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