SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
CummeRbund questions glados Bioinformatics 15 07-10-2017 02:13 AM
Some cummeRbund questions rvaerle RNA Sequencing 61 02-25-2015 04:01 AM
cummeRbund get isoforms paolo.kunder Bioinformatics 6 04-23-2014 07:59 AM
CummeRbund Volcanoplot error Him26 Bioinformatics 1 04-20-2012 10:50 AM
CummeRbund error without CDS dvanic Bioinformatics 5 03-12-2012 09:24 AM

Reply
 
Thread Tools
Old 06-15-2012, 11:11 AM   #1
newbietonextgen
Member
 
Location: USA

Join Date: Nov 2010
Posts: 56
Default CummeRbund error

Hi all,

I am just running the cummerbund R script and running into a problem when i make the volcano plots. Its only with the volcano plots and other plots look fine, density and scatter. What am i doing wrong? Below is the output...

> library(cummeRbund)
Loading required package: RSQLite
Loading required package: DBI
Loading required package: ggplot2
Loading required package: reshape2
> data <-readCufflinks ("all_samples")
Creating database all_samples/cuffData.db
Reading all_samples/genes.fpkm_tracking
Checking samples table...
Populating samples table...
Writing genes table
Reshaping geneData table
Recasting
Writing geneData table
Reading all_samples/promoters.diff
Writing promoterDiffData table
Reading all_samples/isoforms.fpkm_tracking
Checking samples table...
OK!
Writing isoforms table
Reshaping isoformData table
Recasting
Writing isoformData table
Reading all_samples/isoform_exp.diff
Writing isoformExpDiffData table
Reading all_samples/tss_groups.fpkm_tracking
Checking samples table...
OK!
Writing TSS table
Reshaping TSSData table
Recasting
Writing TSSData table
Reading all_samples/tss_group_exp.diff
Writing TSSExpDiffData table
Reading all_samples/splicing.diff
Writing splicingDiffData table
Reading all_samples/cds.fpkm_tracking
Checking samples table...
OK!
Writing CDS table
Reshaping CDSData table
Recasting
Writing CDSData table
Reading all_samples/cds_exp.diff
Writing CDSExpDiffData table
Reading all_samples/cds.diff
Writing CDSDiffData table
Indexing Tables...
> csDensity(genes(data))
> csScatter(genes(data), 'Infected', 'Uninfected')
> csVolcano(genes(data), 'Infected', 'Uninfected')
Error in if (nrow(layer_data) == 0) return() : argument is of length zero


Also is there a way i can change the X axis of the volcano plots so that they all look similar? When i searched around ...I found the code in another thread..So can i switch the scale_y_continuous to scale_x_continuous? or do i just use xlimits= c(-x, x ).

v<-csVolcano(genes(cuff),'Sample_A','Sample_B')
v+ scale_y_continuous(limits = c(0, 20))

Session info

> sessionInfo()
R version 2.15.0 (2012-03-30)
Platform: x86_64-pc-linux-gnu (64-bit)

locale:
[1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
[3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8
[5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8
[7] LC_PAPER=C LC_NAME=C
[9] LC_ADDRESS=C LC_TELEPHONE=C
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C

attached base packages:
[1] stats graphics grDevices utils datasets methods base

other attached packages:
[1] cummeRbund_1.2.0 reshape2_1.2.1 ggplot2_0.9.1 RSQLite_0.11.1
[5] DBI_0.2-5

loaded via a namespace (and not attached):
[1] colorspace_1.1-1 dichromat_1.2-4 digest_0.5.2 grid_2.15.0
[5] labeling_0.1 MASS_7.3-18 memoise_0.1 munsell_0.3
[9] plyr_1.7.1 proto_0.3-9.2 RColorBrewer_1.0-5 scales_0.2.1
[13] stringr_0.6 tools_2.15.0
>
newbietonextgen is offline   Reply With Quote
Old 06-15-2012, 11:24 AM   #2
newbietonextgen
Member
 
Location: USA

Join Date: Nov 2010
Posts: 56
Default

i did (rebuild=T), did not work. Kindly help
newbietonextgen is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 12:08 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO