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
Biological replicates with cuffdiff, cummeRbund turnersd Bioinformatics 15 11-19-2012 05:59 AM
BLAST for spliced isoforms? tboothby Bioinformatics 8 09-07-2011 08:20 AM
looking for isoforms papori De novo discovery 2 06-22-2011 03:23 AM
Cufflinks-Tophat-isoforms ? repinementer Bioinformatics 1 07-19-2010 07:37 AM

Reply
 
Thread Tools
Old 02-16-2012, 02:31 AM   #1
paolo.kunder
Member
 
Location: Milano, Italy

Join Date: Aug 2011
Posts: 93
Default cummeRbund get isoforms

Hi,
I am facing some problems with cummeRbund R package:

library(cummeRbund)
cuff <- readCufflinks()

#get the top 5 diff expr isoforms
iso.diff <- diffData(isoforms(cuff))
iso.diff.top <- iso.diff[order(iso.diff$q_value),][1:5,]

iso.diff.top

isoform_id sample_1 sample_2 status value_1 value_2 ln_fold_change
20524 TCONS_00020591 controlli casi OK 13.0221 1.18723 -3.45528e+00
39253 TCONS_00039380 controlli casi OK 0.0000 7.12244 1.79769e+308
75370 TCONS_00075672 controlli casi OK 31.1991 4.17233 -2.90258e+00
75504 TCONS_00075806 controlli casi OK 18.1982 2.19050 -3.05446e+00
75485 TCONS_00075787 controlli casi OK 13.1554 1.39641 -3.23587e+00


myIsoIds <- iso.diff.top$isoform_id
myIsoIds
TCONS_00039380.....etc

how I do extract Isoforms values and plot for instance an expression bar plot?

how I interact with Isoforms database and extract isoforms values like I do with genes values (myGenes <- getGenes(cuff, myGeneIds))

Any suggestions will be appreciated,

CIAO!
Paolo
paolo.kunder is offline   Reply With Quote
Old 02-16-2012, 07:36 AM   #2
lgoff
Member
 
Location: Cambridge, MA

Join Date: Feb 2008
Posts: 82
Default

Hi Paolo,
This is something that I am working to build into the next public release of cummeRbund. I recognize the need for a 'getFeatures()' method similar to getGenes() but not tied to the gene-based hierarchy.

In the meantime, you should be able to take the results of the call to fpkm(isoforms(cuff)), apply your filtering, and then create a ggplot2 object using that data directly. Everything you need to plot expression values should be output by this method.

Let me know if I can help with anything else.

Cheers,
Loyal
lgoff is offline   Reply With Quote
Old 02-21-2012, 05:47 PM   #3
vyellapa
Member
 
Location: phoenix

Join Date: Oct 2011
Posts: 59
Default cummeRbund Installation

While installling cummeRbund 1.0.0 on R 2.14 , Im getting a "reshape2" package not found. While trying to manually install "respahe2", I have dependency errors to install packages "stringr", "evaluate" and "testhat". Im not sure how to install the dependencies as there is a deadlock while installing packages. For example, when I try to install "stringr" the error is that "evaluate" dependency is not install and vice-versa. Is there an easy way to install reshape2 with all its dependencies?
vyellapa is offline   Reply With Quote
Old 02-21-2012, 08:16 PM   #4
lgoff
Member
 
Location: Cambridge, MA

Join Date: Feb 2008
Posts: 82
Default

You can install reshape2 and all of its dependencies by starting an R session and typing:

Code:
install.packages('reshape2')
Cheers,
Loyal
lgoff is offline   Reply With Quote
Old 06-03-2012, 05:45 PM   #5
beta.gal
Junior Member
 
Location: Australia

Join Date: Jun 2012
Posts: 1
Default

Hi,

Now that the newer beta version of cummeRbund is out, I was wondering if the getFeatures, or something like it, had been introduced? I have not been able to detect anything like this from the manual.

Regards,
Paula
beta.gal is offline   Reply With Quote
Old 06-04-2012, 04:18 AM   #6
lgoff
Member
 
Location: Cambridge, MA

Join Date: Feb 2008
Posts: 82
Default

Quote:
Originally Posted by beta.gal View Post
Hi,

Now that the newer beta version of cummeRbund is out, I was wondering if the getFeatures, or something like it, had been introduced? I have not been able to detect anything like this from the manual.

Regards,
Paula
Hi Paula,
A getFeatures() method has indeed been introduced ( I think it was introduced in v.1.99.2) it works with the same type of arguments as getGenes(), but returns a cuffFeatureSet instead of retrieving all gene information. You can find details by calling

Code:
> ?getFeatures()
Cheers,
Loyal
lgoff is offline   Reply With Quote
Old 04-23-2014, 07:59 AM   #7
mbjohnson
Member
 
Location: Boston, MA

Join Date: Apr 2009
Posts: 15
Default

Hi Loyal -
I successfully added additional annotations at the gene and isoform levels using addFeatures(). However, these extra annotations were not included in the output of getFeatures() using either level="genes" or level="isoforms". The extra annotations were there using getGenes() but I'm mostly interested in getting only the significant isoforms (with all my added annotations). So I was happy to find the getFeatures() function but it looks like this may be a bug in it.
Thanks!
Matt
mbjohnson 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 10:19 AM.


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