SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
What is the meaning density in CummeRbund plots? Carlos Borroto Bioinformatics 16 12-29-2016 09:05 PM
cummeRbund get isoforms paolo.kunder Bioinformatics 6 04-23-2014 07:59 AM
Biological replicates with cuffdiff, cummeRbund turnersd Bioinformatics 15 11-19-2012 05:59 AM
CDS vs cDNA s13ep MGISEQ (FKA Complete Genomics) 4 04-20-2012 02:19 PM
How can I look up CDS from NCBI litc Bioinformatics 1 12-06-2010 01:03 PM

Reply
 
Thread Tools
Old 03-05-2012, 03:15 AM   #1
dvanic
Member
 
Location: Sydney, Australia

Join Date: Jan 2012
Posts: 61
Default CummeRbund error without CDS

Hi! I'm trying to run cummeRbund on a cufflinks dataset without a reference annotation of any sort (i.e. just the transcripts assembled from bowtie/tophat results), which leads to my cds.diff being empty.
Hence, I get the following error when trying to create a cummeRbund database:

Quote:
Writing CDS table
Error in sqliteExecStatement(con, statement, bind.data) :
bind.data must have non-zero dimensions
How do I get around this (other than by using a reference annotation)? I'd still like to be able to explore the gene/tss/isoform information?

Thanks in advance!
dvanic is offline   Reply With Quote
Old 03-06-2012, 01:58 PM   #2
lgoff
Member
 
Location: Cambridge, MA

Join Date: Feb 2008
Posts: 82
Default

Hi dvanic,
This has been addressed in the development version of cummeRbund. (now v1.1.3). Please try updating to the new version and see if this resolves this issue.

compbio.mit.edu/cummeRbund

Cheers,
Loyal
lgoff is offline   Reply With Quote
Old 03-11-2012, 10:45 PM   #3
christophe
Junior Member
 
Location: Phoenix

Join Date: May 2010
Posts: 7
Default

Hi loyal Goff,
I have been trying running CummeRbund. As suggested I installed the latest version of cummeRbund (v1.1.3) within R 2.15.0 . The Error << Error in data.frame(index = c(1:length(samples)), sample_name = samples) arguments imply differing number of rows: 2, 0
Calls: readCufflinks -> loadGenes -> populateSampleTable -> data.frame
>>
disappeared , but a new one showed up as copied below:
NOTE: my working directory is the ouput directory of cuffdiff. I have to mention that CDS, promoter and TSS files are empty in my case. Only Isoforms and Genes files contain results.



Loading required package: RSQLite
Loading required package: methods
Loading required package: DBI
Loading required package: ggplot2
Loading required package: reshape2
Creating database /cuffdiff_outputs/cuffData.db
Reading /cuffdiff_outputs/genes.fpkm_tracking
Checking samples table...
Populating samples table...
Writing genes table
Reshaping geneData table
Recasting
Error in as.data.frame(cast(genemelt, ... ~ measurement)) :
could not find function "cast"
Calls: readCufflinks -> loadGenes -> as.data.frame



Thanks for your help.
Best,
Christophe
christophe is offline   Reply With Quote
Old 03-11-2012, 10:47 PM   #4
christophe
Junior Member
 
Location: Phoenix

Join Date: May 2010
Posts: 7
Default

Addendum to my previous post:
I have tried to run cummeRbund both on linux and Windows and got the same error:

Error in as.data.frame(cast(genemelt, ... ~ measurement)) :
could not find function "cast"
Calls: readCufflinks -> loadGenes -> as.data.frame
Execution halted
christophe is offline   Reply With Quote
Old 03-12-2012, 06:54 AM   #5
lgoff
Member
 
Location: Cambridge, MA

Join Date: Feb 2008
Posts: 82
Default Temp fix

Hi cristophe,
This was a bug inadvertently created when I switched to reshape2. I have fixed it for future releases, but for now, the workaround is to import the original 'reshape' library before calling 'readCufflinks':

Code:
>library(cummeRbund)
>library(reshape)
>cuff<-readCufflinks()
Let me know if this doesn't work.

Cheers,
Loyal
lgoff is offline   Reply With Quote
Old 03-12-2012, 09:24 AM   #6
christophe
Junior Member
 
Location: Phoenix

Join Date: May 2010
Posts: 7
Default

Hi Loyal,
That works.
Only reshape2 was loading automatically and not reshape.
Thanks
Christophe
christophe 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 04:06 PM.


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