[R] Unable to specify order of a factor

[David Winsemius]

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.weylandt at 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 <istazahn at gmail.com>  
>> wrote:
>>> On Wed, Mar 21, 2012 at 12:00 PM, Justin Montemarano <jmontema at 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 <istazahn at 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 <
>> jmontema at 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 at 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 at 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 at 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