Skip to content

scatterplot to boxplot translation?

7 messages · Vining, Kelly, Uwe Ligges, Tom Fletcher +1 more

Vining, Kelly
# 
My apologies if anyone is seeing this twice...looks like my previous message didn't come through...

Dear UseRs,
I have a feeling this is a relatively simple question, but I'm having a hard time getting my head around it. I have a simple x-y scatterplot with many points, as shown below(attached). I'd like to make a boxplot of this by interval, such that there is one box representing the points in the 0-100 interval, one for the 101-200 interval, and so on. How do I structure my R data frame to be able to generate such a boxplot?


From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org] On Behalf Of Vining, Kelly
Sent: Friday, December 09, 2011 11:01 AM
To: r-help at r-project.org
Subject: [R] scatterplot to boxplot translation?


-------------- next part --------------
A non-text attachment was scrubbed...
Name: C_count_vs_RPKM.png
Type: image/png
Size: 16970 bytes
Desc: C_count_vs_RPKM.png
URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20111209/9bc8354f/attachment.png>
David Winsemius
#
On Dec 9, 2011, at 2:11 PM, Vining, Kelly wrote:

            

      
      
      
    

        
?cut

            

      
      
      
    

        

      
      
    

        
David Winsemius, MD
West Hartford, CT

  
  


  


      

    

    

      
      

        


  

        

      Vining, Kelly
    

    

      
      #
    

  


  

      
Thanks for the tip on "cut," seems like it should work. I must still be missing something, though. Here, I'm cutting on the y variable, then attempting the boxplot:

cutRPKM <- cut(count$RPKM, breaks=4)

head(cutRPKM)
[1] (-0.0995,24.8] (-0.0995,24.8] (-0.0995,24.8] (-0.0995,24.8] (-0.0995,24.8]
[6] (-0.0995,24.8]
Levels: (-0.0995,24.8] (24.8,49.8] (49.8,74.7] (74.7,99.6]

boxplot(as.numeric(cutRPKM))

This gives me a single box instead of five boxes. ??

Thanks again,
--Kelly V.

        

      
      
    

  
  


  


      

    

    

      
      

        


  

        

      Uwe Ligges
    

    

      
      #
    

  


  

      
On 09.12.2011 20:41, Vining, Kelly wrote:

            

      
      
      
    

        
You obviously want:

boxplot(count$RPKM ~ cut(count$RPKM, breaks=seq(0, max(count$RPKM), 
by=100)))


Uwe Ligges

            

      
      
      
    

        

      
      
    

  
  


  


      

    

    

      
      

        


  

        

      Tom Fletcher
    

    

      
      #
    

  


  

      
You need to create some grouping for your cut points (0-100, etc). 

See ?cut

Then, you can use boxplot and formula (y ~ NEWVARIABLE from cut) 
boxplot(y ~ cut(x)) 

There may be other ways to do this, but the above should work.

TF

-----Original Message-----
From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org] On Behalf Of Vining, Kelly
Sent: Friday, December 09, 2011 1:11 PM
To: Vining, Kelly; r-help at r-project.org
Subject: Re: [R] scatterplot to boxplot translation?

My apologies if anyone is seeing this twice...looks like my previous message didn't come through...

Dear UseRs,
I have a feeling this is a relatively simple question, but I'm having a hard time getting my head around it. I have a simple x-y scatterplot with many points, as shown below(attached). I'd like to make a boxplot of this by interval, such that there is one box representing the points in the 0-100 interval, one for the 101-200 interval, and so on. How do I structure my R data frame to be able to generate such a boxplot?


From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org] On Behalf Of Vining, Kelly
Sent: Friday, December 09, 2011 11:01 AM
To: r-help at r-project.org
Subject: [R] scatterplot to boxplot translation?

  
  


  


      

    

    

      
      

        


  

        

      David Winsemius
    

    

      
      #
    

  


  

      
On Dec 9, 2011, at 2:50 PM, Uwe Ligges wrote:

        

            

      
      
      
    

        
In that context (having defined a cut-variable with single-integer  
break argument),  would have thought this should work:

  boxplot(count$RPKM ~ cutRPKM)