Multivariate Bayesian variable selection regression

SuSiE paper results based on estimated prior

This notebook displays another version of this page -- in the SuSiE paper analysis we focused on fixed scaled prior variance (PVE = 0.1). Here we show the results of using estimated prior variance.

Conclusion

Results are largely robust to setting prior to 0.1 vs estimating it.

  • Results on the ROC and PIP direct comparisons now appear closer to that of DAP-G.
  • CS power results remained mostly the same. Estimated prior did not help at improving coverage (reducing false positives).
In [2]:
%cd ../dsc
/project/mstephens/SuSiE/mvarbvs/dsc

PIP calibration

In [3]:
%preview susie_comparison/PIP_comparison_1008_calibrated_prior_0p0_null_0p0.png
%preview susie_comparison/PIP_comparison_1008_calibrated_prior_0p0_null_0p0.png
> susie_comparison/PIP_comparison_1008_calibrated_prior_0p0_null_0p0.png (38.4 KiB):

Variable-level power vs false discovery

In [5]:
%preview susie_comparison/ROC_1008_prior_0p0.png
%preview susie_comparison/ROC_1008_prior_0p0.png
> susie_comparison/ROC_1008_prior_0p0.png (47.4 KiB):

CS power comparison

In [6]:
%preview susie_comparison/cs_eval/cs_comparison_prior_0p0_null_0p0.png
%preview susie_comparison/cs_eval/cs_comparison_prior_0p0_null_0p0.png
> susie_comparison/cs_eval/cs_comparison_prior_0p0_null_0p0.png (15.1 KiB):

PIP direct comparison

Agreement with DAP is now tighter for the one-causal case.

In [7]:
%preview susie_comparison/PIP_comparison_1008_prior_0p0_null_0p0.susie_vs_dap.png
%preview susie_comparison/PIP_comparison_1008_prior_0p0_null_0p0.susie_vs_dap.png
> susie_comparison/PIP_comparison_1008_prior_0p0_null_0p0.susie_vs_dap.png (85.2 KiB):

Reliability of reported CS coverage

In [9]:
%preview susie_comparison/Coverage_1008_prior_0p0.png
%preview susie_comparison/Coverage_1008_prior_0p0.png
> susie_comparison/Coverage_1008_prior_0p0.png (45.9 KiB):

Copyright © 2016-2020 Gao Wang et al at Stephens Lab, University of Chicago