On Mar 21, 2012, at 1:07 PM, Justin Montemarano wrote:
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.
I've had the experience in the last few years that almost all of my
questions to Rhelp have needed to be peacefully euthanized after being
subjected to the rack of hammering into a "reproducible" condition.
--
David.
-
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.
David Winsemius, MD
West Hartford, CT
______________________________________________
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.
