include nsnps and their weights from qq plots

precimed · Nov 12, 2018 · 23bb52d · 23bb52d
1 parent 32efaf8
commit 23bb52d
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 0 deletions.
diff --git a/BGMG_cpp_qq_plot.m b/BGMG_cpp_qq_plot.m
@@ -83,6 +83,7 @@
     qq_options.lamGC_data = BGMG_util.lamGCfromQQ(data_logpvec, hv_logp);
     qq_options.lamGC_model = BGMG_util.lamGCfromQQ(model_logpvec, hv_logp);
     qq_options.n_snps = sum(options.mask);
+    qq_options.sum_data_weights = sum(weights_bgmg(options.mask))
     qq_options.qqlimy = hv_logp(find(isfinite(data_logpvec), 1, 'last' ))*1.05;
 
     plot_data.qq_options = qq_options;

diff --git a/BGMG_cpp_stratified_qq_plot.m b/BGMG_cpp_stratified_qq_plot.m
@@ -78,6 +78,9 @@
         plot_data{conditional_trait, zthresh_index}.model_logpvec = model_logpvec;
         plot_data{conditional_trait, zthresh_index}.params = params;
         plot_data{conditional_trait, zthresh_index}.zthresh = zthresh;
+        plot_data{conditional_trait, zthresh_index}.n_snps = sum(abs(zvec2)>=zthresh);
+        plot_data{conditional_trait, zthresh_index}.sum_data_weights = sum(weights_bgmg(abs(zvec2)>=zthresh))
+
         if (conditional_trait==1) && (zthresh_index == 1)
             plot_data{conditional_trait, zthresh_index}.pdf = pdf;
             plot_data{conditional_trait, zthresh_index}.pdf_zgrid = zgrid;