Skip to content

Labeling points with xYplot

10 messages · justin bem, ratna ghosal, Frank E Harrell Jr +3 more

Original
John Poulsen
# 
Hello R-Users,

I am trying to label points on an xYplot.  I can label the points, but 
then I lose the error bars.  I suspect that I should be using 
panel.xYplot rather than panel.xyplot, but that didn't work either...

Can anybody help me out?  Example code is below.  First xYplot has error 
bars, second tries to label the points and the bars disappear.

Thanks,
John


a=c(1:4)
b=rnorm(4,15,5)
blab=b+1

xYplot(Cbind(b, b-1.5, b+1.5)~a,
        ylim=c(0,max(b)+3))

xYplot(Cbind(b, b-1.5, b+1.5)~a, ylim=c(0,max(b)+3),
        panel=function(x,y,...){
          panel.xyplot(x,y,...)
          ltext(x=list(1,2,3,4),
                y=list(13,14,15,16),
                labels=list("A","B","C","D"), pos=2)})
Frank E Harrell Jr
#
John Poulsen wrote:
One of the purposes of xYplot is to use panel.xYplot
Frank
John Poulsen
# 
Yes, using panel.xYplot is the best thing to do.  I couldn't get it to 
work, which probably means I was specifying it incorrectly.  I would be 
very happy if someone had an example of how to do it.

Thanks again -- John
Frank E Harrell Jr wrote:
John Poulsen
# 
That did it!  Thanks for the help!!! -- John
David Winsemius wrote:
David Winsemius
# 
See if this does what you are expecting.

  xYplot(Cbind(b, b-1.5, b+1.5)~a, ylim=c(0,max(b)+3),
       panel=function(x,y,...){
         panel.xYplot(x,y,subscripts=x,...)
         panel.superpose(x,y, subscripts=x, groups=x,...)
         ltext(x=list(1,2,3,4),
               y=list(13,14,15,16),
               labels=list("A","B","C","D"), pos=2)})


---

David Winsemius, MD
Heritage Labs
comes with the usual warrantee... if it breaks you get to keep the parts
On Nov 10, 2008, at 11:26 AM, John Poulsen wrote:

            

      
      
      
    

        

      
      
    

  
  


  


      

    

    

      
      

        


  

        

      Dieter Menne
    

    

      
      #
    

  


  

      
ratna ghosal <ratnaghosal <at> yahoo.com> writes:

            

      
      
      
    

        
When lost, use str:

library(nlme)
fm1Orth.lme = lme(distance ~ I(age-11), data = Orthodont,
                     random = ~ I(age-11) | Subject)
ints = intervals(fm1Orth.lme)
str(ints)
# $ fixed   : num [1:2, 1:3] 23.168 0.518 24.023 0.66 24.878 ...
#  ..- attr(*, "dimnames")=List of 2
# Estimate of slope 
ints$fixed[2,2]


Dieter