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
Problem with cummeRbund and RSQLite noplacetoh1de Bioinformatics 9 11-03-2014 01:33 PM
cummeRbund get isoforms paolo.kunder Bioinformatics 6 04-23-2014 07:59 AM
cummeRbund installation m_elena_bioinfo Bioinformatics 3 08-02-2012 01:42 AM

Reply
 
Thread Tools
Old 04-21-2012, 05:30 PM   #1
ekimmike
Member
 
Location: EU

Join Date: Apr 2012
Posts: 14
Default Galaxy -> cummeRbund How?

Hi,

I process my data using Galaxy due to the lack of good server, everything work fine Tophat->Cufflink->Cuffcompare->Cuffdiff (these files look fine)

Problem starts when try to open them from R and I cannot force cummeRbund to load the files to db that I import from Galaxy, should I combine them somehow, or simply place in one directory.

I would really appreciate your help.

Best!
ekimmike is offline   Reply With Quote
Old 04-21-2012, 05:41 PM   #2
ekimmike
Member
 
Location: EU

Join Date: Apr 2012
Posts: 14
Default

some more details...

> cuff<-readCufflinks("~XXX/RNAseq/proteinA")
Creating database ~XXX/RNAseq/proteinA/cuffData.db
Reading ~~XXX/RNAseq/proteinA/genes.fpkm_tracking
Checking samples table...
Populating samples table...
Writing genes table
Reshaping geneData table
Recasting
Writing geneData table
Error in sqliteExecStatement(con, statement, bind.data) :
RS-DBI driver: (unable to bind data for parameter ':status')
In addition: Warning messages:
1: In `[<-.factor`(`*tmp*`, ri, value = c(0, 127890, 493903, 358.358, :
invalid factor level, NAs generated
2: In `[<-.factor`(`*tmp*`, ri, value = c(0, 127470, 481005, 262.821, :
invalid factor level, NAs generated
3: In `[<-.factor`(`*tmp*`, ri, value = c(0, 128310, 506801, 453.896, :
invalid factor level, NAs generated
4: In `[<-.factor`(`*tmp*`, ri, value = c(0, 886.049, 8123.86, 2.0214, :
invalid factor level, NAs generated
5: In `[<-.factor`(`*tmp*`, ri, value = c(0, 619.614, 0, 0, 11.9375, :
invalid factor level, NAs generated
6: In `[<-.factor`(`*tmp*`, ri, value = c(0, 1152.48, 4298850, 4.56339, :
invalid factor level, NAs generated
ekimmike is offline   Reply With Quote
Old 04-24-2012, 07:35 AM   #3
Carlos Borroto
Member
 
Location: Baltimore, MD

Join Date: Mar 2011
Posts: 19
Default

Hi,

If you only tell readCufflinks the directory with the files from Galaxy, you will need to first make sure you rename all the files to the default names cuffdiff usually uses. My recommendation would be to instead specify each argument manually.

Something like this
Code:
 
cuff <- readCufflinks(dir = "path/to/galaxy/output/",
                         dbFile = "output_database",
                         geneFPKM = "genes_fpkm_tracking",
                         geneDiff = "genes_exp",
                         isoformFPKM = "isoforms_fpkm_tracking",
                         isoformDiff = "isoforms_exp",
                         TSSFPKM = "tss_groups_fpkm_tracking",
                         TSSDiff = "tss_groups_exp",
                         CDSFPKM = "cds_fpkm_tracking",
                         CDSExpDiff = "cds_exp_diff",
                         CDSDiff = "cds_diff",
                         promoterFile = "promoters_diff",
                         splicingFile = "splicing_diff",
                         rebuild = T)
Hope it helps,
Carlos
Carlos Borroto is offline   Reply With Quote
Old 10-03-2013, 07:57 PM   #4
saran.ela123
Junior Member
 
Location: Singapore

Join Date: Aug 2013
Posts: 5
Default

hi,

Doubt clarification.

what is that output_database ?
saran.ela123 is offline   Reply With Quote
Old 10-30-2013, 12:19 AM   #5
IBseq
Member
 
Location: uk

Join Date: Jul 2012
Posts: 56
Default

hi all,
can anyone tell em how to fix this?

head(features(myGene))
Error in function (classes, fdef, mtable) :
unable to find an inherited method for function "features", for signature "CuffGene"


thanks,
ibseq
IBseq is offline   Reply With Quote
Old 07-25-2019, 05:03 AM   #6
EdRyder
Junior Member
 
Location: UK

Join Date: Jun 2019
Posts: 1
Default

Just found this thread as I had the same problem as above for:
...
head(features(myGene))
Error in function (classes, fdef, mtable) :
unable to find an inherited method for function "features", for signature "CuffGene"
...
and managed to find the cause

The problem is if you load the library("GenomicFeatures") explicitly, to make a txdb object for example, it somehow stops the features(myGene) from working. The trick is to load CummeRbund, load GenomicFeatures, do what you need with the latter, then unload GenomicFeatures again. It will complain that it can't, but seems to fix the issue nonetheless.

Using R 3.6.0 and BioConductor 3.9, in windows and CentOS Linux.

--Ed
EdRyder 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:59 PM.


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