edsgard / geneiase

Gene-based test of allele-specific expression without requiring haplotype information
GNU General Public License v3.0
14 stars 2 forks source link

Do zero p.nom and fdr values indicate significant ASE? #4

Open jaavedm opened 3 years ago

jaavedm commented 3 years ago

Hello,

I have a simple question. After running geneiase in static mode on my dataset, I see lots of genes reporting p.nom == 0 and fdr == 0. I also see these zero p.nom and fdr values in the output from the test datasets provided in the package. Here is an example from the test dataset:

    feat n.vars mean.s median.s  sd.s  cv.s liptak.s p.nom   fdr
   <dbl>  <dbl>  <dbl>    <dbl> <dbl> <dbl>    <dbl> <dbl> <dbl>
 1  12.1     11   4.72     4.10  4.58 0.970    15.7      0     0
 2 124.       4   7.81     6.20  6.86 0.878    15.6      0     0
 3 151.       4   7.51     5.58  7.10 0.946    15.0      0     0

Do these genes with zero values have significant allele-specific expression? I'm curious if the zero values means lack of machine precision to represent very small p-values, or if these zero values are a placeholder value to indicate a failure in the statistical test for these genes?

Thank you for your help.

jaavedm commented 3 years ago

Hello @edsgard, Can you help me with my question regarding interpreting p.nom and fdr == 0? Thanks.