Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • 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

Latest Articles

Collapse

  • seqadmin
    Techniques and Challenges in Conservation Genomics
    by seqadmin



    The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

    Avian Conservation
    Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
    03-08-2024, 10:41 AM
  • seqadmin
    The Impact of AI in Genomic Medicine
    by seqadmin



    Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
    02-26-2024, 02:07 PM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 03-14-2024, 06:13 AM
0 responses
34 views
0 likes
Last Post seqadmin  
Started by seqadmin, 03-08-2024, 08:03 AM
0 responses
72 views
0 likes
Last Post seqadmin  
Started by seqadmin, 03-07-2024, 08:13 AM
0 responses
81 views
0 likes
Last Post seqadmin  
Started by seqadmin, 03-06-2024, 09:51 AM
0 responses
68 views
0 likes
Last Post seqadmin  
Working...
X