[R] Panel plots for means of cyclical observations
I want to create a panel plot using xyplot of a line graph whose
x-axis is months of the year and y-axis is the average rainfall in a
given month over the 6 years the data spans.
There should be two levels in this panel plot: odd and even months.
Creating this plot without splitting it into levels is quite
straightforward (creating a for loop to compute a vector of averages)
but the approach is not useful if you want to split the plots into
different levels.
Here is the code:
dfmt <- "%d/%m/%Y"
date <- seq(as.Date("01/01/2010", dfmt), as.Date("31/12/2015", dfmt),
"day")
month <- months(date)
rainfall <- runif(2191, 0, 150)
monthsOfYear <- c("January", "February", "March", "April",
"May", "June", "July", "August", "September", "October",
"November", "December")
parity <- match(month, monthsOfYear) %% 2
# even parity = 0, odd parity = 1
z <- data.frame(rainfall, date, month, parity)
The problem with using xyplot( y ~ x | f, data=z .. ) is the x and y I
want to plot are not columns in z but rather some kind of statistical
summary of columns.
__
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.
Re: [R] Panel plots for means of cyclical observations
## Kieran,
## I think the root problem is that you allowed the levels of month to
be alphabetical.
## continuing with your example
##You need to take control of the levels with a statement like
levels(z$month)
z$month <- factor(z$month, levels=monthsOfYear)
levels(z$month)
## now you can write something like
z$month.parity <- factor(z$month,
levels=unlist(matrix(levels(z$month), 6, 2, byrow=TRUE)))
levels(z$month.parity)
## [1] "January" "March" "May" "July" "September" "November"
## [7] "February" "April" "June" "August""October" "December"
xyplot(rainfall ~ date | month, layout=c(2, 6), data=z)
xyplot(rainfall ~ date | month.parity, layout=c(6, 2), data=z)
## I hope the modifications to deal with your summary measures will be
straightforward.
## Rich
On Fri, Jan 15, 2016 at 7:16 AM, Kieran wrote:
> I want to create a panel plot using xyplot of a line graph whose
> x-axis is months of the year and y-axis is the average rainfall in a
> given month over the 6 years the data spans.
>
> There should be two levels in this panel plot: odd and even months.
>
> Creating this plot without splitting it into levels is quite
> straightforward (creating a for loop to compute a vector of averages)
> but the approach is not useful if you want to split the plots into
> different levels.
>
> Here is the code:
>
> dfmt <- "%d/%m/%Y"
> date <- seq(as.Date("01/01/2010", dfmt), as.Date("31/12/2015", dfmt),
> "day")
> month <- months(date)
> rainfall <- runif(2191, 0, 150)
> monthsOfYear <- c("January", "February", "March", "April",
> "May", "June", "July", "August", "September", "October",
> "November", "December")
>
> parity <- match(month, monthsOfYear) %% 2
> # even parity = 0, odd parity = 1
>
> z <- data.frame(rainfall, date, month, parity)
>
> The problem with using xyplot( y ~ x | f, data=z .. ) is the x and y I
> want to plot are not columns in z but rather some kind of statistical
> summary of columns.
>
> __
> 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.
Re: [R] Panel plots for means of cyclical observations
Hi
Continuing on from the data.frame code for z that you supplied
Using the data without summarizing first
z$mth = format(date "%m")
z$mth = format(date, "%m")
xyplot(rainfall ~ as.numeric(mth), z,
groups = yr,
type = "l",
auto.key = T,
scales = list(x = list(at = 1:12,
labels = month.abb,
rot = 60)),
panel = panel.superpose,
panel.groups = function(x,y, ...){
panel.average(x,y, fun = mean, horizontal = FALSE,
...)
}
) ## xyplot
You may want to look into the zoo package as it has several date grouping
functions.
untested
z.tom <- aggregate(rainfall ~ year +month, z, mean, na.rm = T)
xyplot(rainfall ~ month, z.tom, groups = year,
scales = ... ,
panel = panel.superpose)
Regards
Duncan
Duncan Mackay
Department of Agronomy and Soil Science
University of New England
Armidale NSW 2351
Email: home: mac...@northnet.com.au
-Original Message-
From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of Kieran
Sent: Friday, 15 January 2016 22:16
To: R-help@r-project.org
Subject: [R] Panel plots for means of cyclical observations
I want to create a panel plot using xyplot of a line graph whose
x-axis is months of the year and y-axis is the average rainfall in a
given month over the 6 years the data spans.
There should be two levels in this panel plot: odd and even months.
Creating this plot without splitting it into levels is quite
straightforward (creating a for loop to compute a vector of averages)
but the approach is not useful if you want to split the plots into
different levels.
Here is the code:
dfmt <- "%d/%m/%Y"
date <- seq(as.Date("01/01/2010", dfmt), as.Date("31/12/2015", dfmt),
"day")
month <- months(date)
rainfall <- runif(2191, 0, 150)
monthsOfYear <- c("January", "February", "March", "April",
"May", "June", "July", "August", "September", "October",
"November", "December")
parity <- match(month, monthsOfYear) %% 2
# even parity = 0, odd parity = 1
z <- data.frame(rainfall, date, month, parity)
The problem with using xyplot( y ~ x | f, data=z .. ) is the x and y I
want to plot are not columns in z but rather some kind of statistical
summary of columns.
__
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.
