SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
filter above value Pol8 Bioinformatics 2 03-03-2015 06:10 AM
Filter PE reads together nr23 Illumina/Solexa 3 11-29-2012 04:03 AM
GATK filter zhangtao13039 Bioinformatics 0 11-08-2012 07:16 PM
using filter in Samtools emilyjia2000 Bioinformatics 1 08-01-2011 02:50 PM
DGE - filter or not filter masterpiece Bioinformatics 0 07-11-2011 09:55 PM

Reply
 
Thread Tools
Old 04-09-2015, 06:30 AM   #1
sghignone
Junior Member
 
Location: Italy

Join Date: Nov 2008
Posts: 6
Default cummeRbund :: filter on repFpkm

Dear all,
I'm wondering whether anyone has already implemented a filter based on repFpkm in cummeRbund, so to select all those genes for which the Fpkm is above of a certain threshold for all the replicates belonging to a condition OR to another one.
I mean that, for example, a gene must have in control at least fpkm = 1, for all replicates #1, #2 and #3, even if it doesn't happen for the treatment; and vice versa.

i'm using

gene.diff<-diffData(genes(cuff))

gene.diff.q1.filtered <- gene.diff[gene.diff$value_1 > 1,] #condition q1
gene.diff.q2.filtered <- gene.diff[gene.diff$value_2 > 1,] #condition q2


but it selects only genes having 'global' Fpkm above 1 (takes all reps together).


i need to get rid of cases like this:

>gene.diff.filtered[gene.diff.filtered$gene_id == "comp13340",]
gene_id sample_1 sample_2 status value_1 value_2 log2_fold_change test_stat p_value q_value significant
2102 comp13340 q1 q2 OK 0 1.01699 Inf NA 5e-05 0.00135144 yes



> gene.repFpkm[gene.repFpkm$gene_id == "comp13340",]

gene_id sample_name replicate rep_name raw_frags internal_scaled_frags external_scaled_frags fpkm
2102 comp13340 q1 0 q1_0 0 0.00000 0.00000 0.00000
35771 comp13340 q1 1 q1_1 0 0.00000 0.00000 0.00000
69440 comp13340 q1 2 q1_2 0 0.00000 0.00000 0.00000
103109 comp13340 q2 0 q2_0 2 2.47685 2.47685 1.26962
136778 comp13340 q2 1 q2_1 0 0.00000 0.00000 0.00000
170447 comp13340 q2 2 q2_2 8 8.39216 8.39216 1.78133
effective_length status
2102 NA OK
35771 NA OK
69440 NA OK
103109 NA OK
136778 NA OK
170447 NA OK

it's clear that not all reps have fpkm > 1 ..........

Thanks for any help or hint!

Cheers,

Stefano
sghignone is offline   Reply With Quote
Reply

Tags
cummerbund, tophat

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:55 PM.


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