.prsice should report PRS R2, not full model R2

choishingwan · May 25, 2020 · 3667274 · 3667274
1 parent ba87617
commit 3667274
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/inc/prsice.hpp b/inc/prsice.hpp
@@ -329,7 +329,7 @@ class PRSice
                              std::unique_ptr<std::ostream>& prsice_out)
     {
         (*prsice_out) << pheno_name << "\t" << region_name << "\t"
-                      << cur_threshold << "\t" << res.r2;
+                      << cur_threshold << "\t" << res.r2 - m_null_r2;
         if (has_prevalence && m_binary_trait)
             (*prsice_out) << "\t" << get_adjusted_r2(res.r2, top, bot);
         else if (has_prevalence)