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
    Exploring the Dynamics of the Tumor Microenvironment
    by seqadmin




    The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...
    07-08-2024, 03:19 PM
  • seqadmin
    Exploring Human Diversity Through Large-Scale Omics
    by seqadmin


    In 2003, researchers from the Human Genome Project (HGP) announced the most comprehensive genome to date1. Although the genome wasn’t fully completed until nearly 20 years later2, numerous large-scale projects, such as the International HapMap Project and 1000 Genomes Project, continued the HGP's work, capturing extensive variation and genomic diversity within humans. Recently, newer initiatives have significantly increased in scale and expanded beyond genomics, offering a more detailed...
    06-25-2024, 06:43 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 07-10-2024, 07:30 AM
0 responses
25 views
0 likes
Last Post seqadmin  
Started by seqadmin, 07-03-2024, 09:45 AM
0 responses
201 views
0 likes
Last Post seqadmin  
Started by seqadmin, 07-03-2024, 08:54 AM
0 responses
211 views
0 likes
Last Post seqadmin  
Started by seqadmin, 07-02-2024, 03:00 PM
0 responses
193 views
0 likes
Last Post seqadmin  
Working...
X