search

gkichaev / PAINTOR_V3.0

Fast, integrative fine mapping with functional data
55 stars 21 forks source link

Computation of credible set #12

Closed MichelNivard closed 7 years ago

MichelNivard commented 7 years ago

This may be blindingly obvious, but I could not find the way in which you compute the XX% credible set.

Is it the smallest set with a cumulative sum of posterior probabilites which is XX% of the total sum of posterior probabilities?

Best, Michel

ruthjohnson95 commented 7 years ago

Yes, exactly!

In the code, we rank all of the SNPs by their calculated posterior probability and then add to a set until the cumulative probability of the set is equal to the (XX%)*(total posterior prob), which is what you mentioned above.