search

joepickrell / gwas-pw

Pariwise analysis of GWAS
GNU General Public License v2.0
36 stars 22 forks source link

.MLE file columns #6

Open ssjj791 opened 5 years ago

ssjj791 commented 5 years ago

Can anyone explain what the 2nd, 3rd and 4th column in the output .MLE file represent? For example,

pi_0 0.288863 0.401738 0.497371 pi_1 0.00605293 0.040412 0.098162 pi_2 0.411552 0.516692 0.622863 pi_3 0.000653681 0.0120106 0.0519648 pi_4 0.00129967 0.0291475 0.0912592

The first column means different models. What do the values mean?

hrafnfaedhir commented 1 year ago

The third column is the regional level prior for each chunk. I am still uncertain as to what parameters the 2nd and 4th fields contain.

hrafnfaedhir commented 1 year ago

Looking at the code, Column 2 is the low estimate of the Prior for each model in the final MCMC convolution. Column 3 is the Prior for each model used in the segbfs.gz file and used for the calculation of the Posteriors. Column 4 is the high estimate of the Priors in the final convolution of the MCMC.