Hi Shakeel,
Assuming that you are starting with a bunch of dates:
# make a vector of character strings that can be converted to dates
rep_dates<-paste(sample(1:30,500,TRUE),sample(1:12,500,TRUE),
sample(2013:2017,500,TRUE),sep="/")
# if this isn't your format, change it
date_format<-"%d/%m/%Y"
# create a data frame with a column of dates
rep_df<-data.frame(rep_dates=as.Date(rep_dates,format=date_format))
# add the week of the year
rep_df$rep_week<-format(rep_df$rep_dates,"%V")
# add the year
rep_df$rep_year<-format(rep_df$rep_dates,"%Y")
# get a table of the weekly counts by year
rep_tab<-table(rep_df$rep_week,rep_df$rep_year)
# get the row means (5 year averages)
rep5<-apply(rep_tab,1,mean)
# plot the 5 year weekly averages
plot(rep5,type="b",ylim=c(0,4),xlab="Week",ylab="Reports per week")
# add the 2017 weekly counts
points(rep_tab[,5],type="b",col="red")
legend(1,4,c("5 yr average","2017"),pch=1,lty=1,col=c("black","red"))
Jim
On Wed, May 9, 2018 at 4:37 PM, Shakeel Suleman
<shakeel.sule...@phe.gov.uk> wrote:
> Hi,
>
> I am fairly new to 'R' and would like advice on the following. I want to
> calculate a weekly average number of reports (e.g. of flu, norovirus) based
> on the same weeks for the last five years. I will then use this to plot a
> chart with 52 points for the average based on the last five years; another
> line will then plot the current year, enabling a comparison of current weekly
> counts against a five year average for the same week. I would like some
> advice on how this can be done in 'R' . My data is disaggregated data - with
> dates in the format in 01/01/2018.
>
> Thanks
>
> Shakeel Suleman
>
>
>
> **************************************************************************
> The information contained in the EMail and any attachments is confidential
> and intended solely and for the attention and use of the named addressee(s).
> It may not be disclosed to any other person without the express authority of
> Public Health England, or the intended recipient, or both. If you are not the
> intended recipient, you must not disclose, copy, distribute or retain this
> message or any part of it. This footnote also confirms that this EMail has
> been swept for computer viruses by Symantec.Cloud, but please re-sweep any
> attachments before opening or saving. http://www.gov.uk/PHE
> **************************************************************************
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
