Here's a lattice solution using some faked data. library(lattice) testdat <- data.frame(time = rep(rep(1:10, each = 4), 8), y = rnorm(320), gender = factor(rep(rep(c('F', 'M'), each = 80), 2)), grade = factor(rep(c('1-3', '4-6'), each = 160)), gp = factor(rep(rep(c('Control', 'Intervention'), each = 40), 4))) testmn <- aggregate(y ~ time + gp + gender + grade, data = testdat, FUN = mean) xyplot(y ~ time | gender * grade, groups = gp, data = testmn, type = c('p', 'l'), auto.key = list(space = 'top', lines = TRUE, points = FALSE, columns = 2))
I have four replicate observations at each of ten times for eight combinations of gender, grade level and treatment group. aggregate() is used to average the four observations at each time point, and then point/line graphs of average response vs. time in each treatment is plotted for each gender and grade level combination. I hope I caught the essential characteristics of what you wanted. Dennis On Fri, Aug 13, 2010 at 2:58 PM, Robin Jeffries <rjeffr...@ucla.edu> wrote: > Hello, > > I need to plot the means of some outcome for two groups (control vs > intervention) over time (discrete) on the same plot, for various subsets > such as gender and grade level. What I have been doing is creating all > possible subsets first, using the aggregate function to create the means > over time, then plotting the means over time (as a simple line plot with > both control & intervention on one plot) for one subset. I then use par() > and repeat this plot for each gender x grade level subset so they all > appear > on one page. > > > This appears to me to be very similar to an xyplot, something like > mean(outcome) ~ gender + gradelevel. However, I can't figure out how I > could get both control and intervention lines in the same plot. > > Any suggestions? What i'm doing now -works-, but just seems to be the long > way around. > > -Robin > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help@r-project.org mailing list > 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. > [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list 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.