Hi all: I've got it... it appears that total.density was also defined in two separate data frames (se.predict.data and dc.predict.data) with levels order 16, 32, 8. Using relevel(), I moved 8 to the first position and it's solved the plotting problem.
Ista's 'minimal' reproducible code request prompted me to discover my error; thanks all. - Justin Montemarano Graduate Student Kent State University - Biological Sciences http://www.montegraphia.com On Wed, Mar 21, 2012 at 12:42 PM, R. Michael Weylandt < michael.weyla...@gmail.com> wrote: > You'll also want to use dput() to send us an exact encoding of your > data when making that reproducible example: there might be something > subtle at play here that print methods won't show. > > Michael > > On Wed, Mar 21, 2012 at 12:28 PM, Ista Zahn <istaz...@gmail.com> wrote: > > On Wed, Mar 21, 2012 at 12:00 PM, Justin Montemarano <jmont...@kent.edu> > wrote: > >> Ista: > >> > >> Your attached code did work for me; moreover, the facets were presented > in > >> the desired order with facet_wrap() and facet_grid(), which is what I'm > >> using because I have a second factor used in facet_grid(). > >> > >> Still, my plots with total.density as a facet are coming out in 16, 32, > 8, > >> and I'm not seeing why. Below is my plot code - > >> > >>> ggplot(ag.tab[ag.tab$plant.sp == 'EC',], aes(x = days.out, y = > >>> per.remain)) + facet_grid(total.density ~ prop.ec) + > >>> #add point and error bar data > >>> theme_set(theme_bw()) + > >>> geom_point() + geom_errorbar(aes(ymin = per.remain - se, ymax = > >>> per.remain + se), width = 3) + > >>> #add predicted model data > >>> geom_line(data = se.predict.data[se.predict.data$plant.sp == > 'EC',], > >>> aes(x = x.values, y = predicted.values), colour = c('red')) + > >>> geom_line(data = dc.predict.data[dc.predict.data$plant.sp == > 'EC',], > >>> aes(x = x.values, y = predicted.values), colour = c('blue'), linetype = > >>> c('dashed')) + > >>> > >>> xlab('Day') + ylab('Percent Mass Remaining') + > opts(panel.grid.major = > >>> theme_blank(), panel.grid.minor = theme_blank()) > >> > >> Is there anything odd about it that might be producing the odd ordering > >> problem? FYI, avoiding subsetting ag.tab doesn't do the trick. > > > > I don't know. Please create a minimal example that isolates the > > problem. You can start with > > > > levels(ag.tab$total.density) > > > > ggplot(ag.tab[ag.tab$plant.sp == 'EC',], aes(x = days.out, y = > per.remain)) + > > facet_grid(total.density ~ prop.ec) + > > geom_point() > > > > Best, > > Ista > > > >> - > >> Justin Montemarano > >> Graduate Student > >> Kent State University - Biological Sciences > >> > >> http://www.montegraphia.com > >> > >> > >> On Wed, Mar 21, 2012 at 11:42 AM, Ista Zahn <istaz...@gmail.com> wrote: > >>> > >>> Hi Justin, > >>> > >>> this gives the correct order (8, 16, 32) on my machine: > >>> > >>> total.density <- > >>> > >>> > c(8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,16,16,16,16,16,16,16,16,16,16,16,16,16,16,16,16,16,16,16,16,16,16,16,16,32,32,32,32,32,32,32,32,32,32,32,32,32,32,32,32,32,32,32,32,32,32,32,32) > >>> total.density <- factor(total.density, levels=c(8, 16, 32), > ordered=TRUE) > >>> str(total.density) > >>> > >>> order(levels(total.density)) > >>> > >>> dat <- data.frame(td = total.density, v1 = > rnorm(1:length(total.density))) > >>> > >>> ggplot(dat, aes(x = v1)) + > >>> geom_density() + > >>> facet_wrap(~td) > >>> > >>> Does it work for you? If yes, then you need to tell us what you're > >>> doing that is different from this example. If no, please give use the > >>> output of sessionInfo(). > >>> > >>> best, > >>> Ista > >>> > >>> On Wed, Mar 21, 2012 at 11:16 AM, Justin Montemarano < > jmont...@kent.edu> > >>> wrote: > >>> > I think I understand, but I believe my original interest is in the > order > >>> > of > >>> > levels(total.density), since ggplot appears to be using that to order > >>> > the > >>> > facets. Thus, I'm still getting three graphs, ordered (and displayed > >>> > as) > >>> > 16 to 32 to 8, rather than the more intuitive, 8 to 16 to 32. I'm > sorry > >>> > if > >>> > I wasn't clear and/or I've missed your message. > >>> > - > >>> > Justin Montemarano > >>> > Graduate Student > >>> > Kent State University - Biological Sciences > >>> > > >>> > http://www.montegraphia.com > >>> > > >>> > [[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. > >> > >> > > > > ______________________________________________ > > 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.
