click here to see plots + legends
qR.rb -r ./marg-lhood-bf.R
Simulate GWAS summary stats, and analyse under various prior dn or variance, storing PP and PP-rank at the true CV
rm ~/scratch/fm-priors/* ## cleanup
for cvtype in friendly lonely; do
for nn in 2000 5000 10000; do
qR.rb -r -y 0-9 ./sim.R --args cvtype=$cvtype NN=$nn
done
done
## wait
qR.rb -r ./collate.R
Over all simulations with p < 5e-8:
Rscript -e 'rmarkdown::render("Readme.Rmd")'
mv Readme.html docs