SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
extract subset of fastq based on length sequence?? dapizarro Ion Torrent 12 08-26-2014 01:01 PM
Limit findSimilar() to a subset of genes in the whole Cufflinks data set [cummeRbund] Gus Bioinformatics 1 03-23-2014 11:24 AM
cummeRbund 'features' error shengzizhang Bioinformatics 2 12-04-2013 05:43 AM
DESeq2-test a subset of genes quinne5 RNA Sequencing 2 10-16-2013 05:19 AM
cuffdiff on subset of alignments EGrassi Bioinformatics 0 06-20-2012 08:05 AM

Reply
 
Thread Tools
Old 10-01-2014, 10:11 PM   #1
beeman
Member
 
Location: Australia

Join Date: May 2012
Posts: 20
Default CummeRbund: subset Cuffdiff genes based >1 isoform or other features

I using cummeRbund to analyse some cuffdiff output. I'd like to be able to create a geneset that contains genes with >1 isoform (or other defined sig diff feature like >! tss or promoter etc).

I can create an object containing only sig DE genes as per the instructions in the manual using;
> mySigGeneIds<-getSig(cuff,alpha=0.05,level='genes')
> mySigGenes<-getGenes(cuff,mySigGeneIds)

But I was wondering if it's possible to filter these gids further to get those items with differential isoform expression between samples. At the moment I only want to get gids with >1 isoform, however I'd also like to apply this strategy/filter to other features such as TSS, promoter, etc etc..


Additionally, once such a list is created what would be the correct way of looping through the genes and plotting data for each gene?

eg;
1. get sig genes with >1 isoform
> foobar
2. For each gene in list, plot expression
> for gene_id in x_list:
> myGeneId<-gene_id
> myGene<-getGene(cuff,myGeneId)
> gl.iso.rep<-expressionPlot(isoforms(myGene),replicates=T


I'm not very experienced with R, but I've read all cummeRbund manual and some other documentation and I can't seem to figure it out. Any assistance would be warmly welcomed!!

Thanks,
beeman is offline   Reply With Quote
Reply

Tags
cuffdiff, cufflinks, cummerbund, getsig

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 01:15 AM.


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