diwis / bufet Goto Github PK

Hi there,

I was using BUFET to enrich my DE miRNAs, but found there were no involved genes or miRNAs in the results. Could you help me with that?

Best,

We were trying to understand how do you calculate the Observed.Target.Gene.Overlap.Proportio by doing the operations by hand.
Let me provide a concrete term example:

This is our results:

GO.term.ID GO.term.size Observed.Target.Gene.Overlap.Proportio Mean.Random.Simulated.MicroRNA.Target.Overlap.Proportion
hsa03010~Ribosome           27                             0.00276243                                               0.00445511
   One.sided.empirical.p.value Benjamini.Hochberg.0.05.FDR
                   0.345034                           -

This is involved-genes.txt

>hsa03010
MIR1	
MIR2	GEN1,GEN2
MIR3	
MIR4
MIR5	
MIR6	
MIR7	
MIR8	
MIR9	
MIR10

What we are expecting is
Observed.Target.Gene.Overlap.Proportio = 2 / 27
Being 2 (GEN1,GEN2) the number of target genes and 27 the number of genes annotated in hsa03010.
But 2/27 (0.07407407) is not the reported Observed.Target.Gene.Overlap.Proportio (0.00276243).