[R] Fwd: Re: bwplot puts the bars in the wrong place

[James Rome]

Peter, that finally worked.

        hrs = seq(1, 24, 1)                       
        g = bwplot(tt~OnHFact |Runway, data=gdf, ylab="Taxi time (min)",
        main=title, xlab="Hour of day",  par.strip.text=list(cex=0.7),
rot=90,  xlim=c(0, 25),
        drop.unused.levels=FALSE, scales=list(x = list(rot=90, cex=.6
        ,at=hrs, labels=hrlabs
        )),
        panel = function(x, ...) {
            panel.grid(h = -1, v = 24)                               
                           
            panel.bwplot(x, ...)                                       
                                                                       
                                                                       
                                                               
      } )
        print(g)

But, it seems a lot easier and safer to call xyplot() and panel.bwplot()
without factors. The scary thing here is that you must look very closely
at the data and the plot to see that there is an error.

Thanks for the help,
Jim

On 4/19/2010 10:44 PM, Peter Ehlers wrote:James,

It's actually the bars for hour 3 (which don't exist) that are
missing. You still need the 'drop.unused.levels=FALSE' and if
you make 'OnHour' into a factor then you won't need the
'horizontal=FALSE'.

 -Peter Ehlers

On 2010-04-19 8:27, James Rome wrote:
> You were right about the gdf$. Something else must have been giving me
> the errors.
>
> So I tried it with your scale suggestion:
>          hrs = seq(1, 24, 1)
>          hrlabs = as.character(seq(0,23,1))
>              g = bwplot(tt~OnHour |Runway, data=gdf, horizontal=FALSE,
> ylab="Taxi time (min)",
>              main=title, xlab="Hour of day",
> par.strip.text=list(cex=0.7), rot=90,  xlim=c(0, 25),
>              scales=list(x = list(rot=90, cex=.6
>              ,at=hrs, labels=hrlabs
>              )),
>              panel = function(x, ...) {
>                  panel.grid(h = -1, v = -1)
>
>                  panel.bwplot(x, ...)
>
>
>            } )
>          print(g)
> And now the bars for hour 23 are missing
>
> On 4/17/2010 4:16 PM, Peter Ehlers wrote:James,
>
> Your problem with bwplot is with your scales definition.
> bwplot places the plots at locations 1:24; you can then
> supply arbitrary labels for these locations. Try changing
> the following:
>
>   xlim = c(-1, 24) to xlim = c(0, 25)
>   at = hrs to at = hrs + 1 (or redefine hrs)
>
> As to needing to use 'gdf$': this makes no sense to me.
> What error occurs if you leave that off?
>
>   -Peter Ehlers
>
> On 2010-04-17 13:00, James Rome wrote:
>> David,
>>
>> I did post a solved message:
>>           hrs = seq(0, 23, 1)
>>           hrlabs = as.character(seq(0,23,1))
>>           g = xyplot(gdf$tt~gdf$OnHour |gdf$Runway, data=gdf, ylab="Taxi
>> time (min)", main=title, xlab="Hour of day",
>>               xlim=c(-1, 24), scales=list(x = list(rot=90, cex=.6,
>> alternating=c(3,3,3,3),
>>               at=hrs, labels=hrlabs
>>               )),
>>               panel=function(x, ...) {
>>                   panel.grid(h = -1, v = 24)
>>                   panel.bwplot(x, horizontal=FALSE, col="black",...)
>>
>>               }
>>           )
>>           print(g)
>> created the plot I sent you (withxyplot.png), with the 0 data plotted on
>> zero, and the 3AM data blank.
>>
>> I arrived at this by doing xyplot and observing that the points were
>> plotted correctly. Then I added panel.bwplot() and it was still correct,
>> so I removed panel.xyplot() and it worked. But replacing the main call
>> to xyplot with bwplot gets it wrong.
>> For example, I did it with
>>
>>           hrs = seq(0, 23, 1)
>>           hrlabs = as.character(seq(0,23,1))
>>           g = bwplot(gdf$tt~gdf$OnHFact |gdf$Runway, data=gdf,
>> ylab="Taxi
>> time (min)", main=title, xlab="Hour of day",
>>               xlim=c(-1, 24), horizontal=FALSE,
>> drop.unused.levels=FALSE,
>>               scales=list(x = list(rot=90, cex=.6,
>> alternating=c(3,3,3,3),
>>               at=hrs, labels=hrlabs
>>               )),
>>               panel=function(x, ...) {
>>                   panel.grid(h = -1, v = 24)
>>                   panel.bwplot(x,
>> col="black",...)
>>               }
>>           )
>>           print(g)
>> and get the attached withbwplotAndFactors.png, which is wrong. It was
>> also wrong with OnHour, and with drop.unused.levels=TRUE.
>>
>> And I did indeed get errors if I left the gdf$ out of the formulas,
>> because it is defined in many data frames.
>>
>> Thanks,
>> Jim
>>
>> On 4/17/2010 12:35 PM, David Winsemius wrote:
>> On Apr 17, 2010, at 12:29 PM, James Rome wrote:
>>
>>> It has all the data needed, and only that data, which I got dinged on
>>> before.
>>
>> I was seeing a different number of points in some panels in your plot
>> compared to my plot. I have not downloaded the most recent dataset, and
>> from what you are saying here should  I assume it had not changed? But
>> your plot had a title that could not have been created by the posted
>> code , so i just don't know where you are with all of this. Are you
>> going to post a SOLVED message when it is finally answered?
>>
>>
>>
>>
>>
>> ______________________________________________
>> 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.
>
>



-------------- next part --------------
A non-text attachment was scrubbed...
Name: bwplotGood.PNG
Type: image/png
Size: 52549 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20100420/6a165dba/attachment.png>