From 86323ccaf26ad0a3b706a67a0014dd04b9965823 Mon Sep 17 00:00:00 2001
From: Pjotr Prins
Date: Thu, 5 Oct 2017 11:55:03 +0000
Subject: LMM: skip beta (again) with LRT

---
 src/lmm.cpp | 14 ++++++++++----
 1 file changed, 10 insertions(+), 4 deletions(-)

(limited to 'src')

diff --git a/src/lmm.cpp b/src/lmm.cpp
index e2f23a2..1193700 100644
--- a/src/lmm.cpp
+++ b/src/lmm.cpp
@@ -107,8 +107,10 @@ void LMM::WriteFiles() {
   }
 
   auto common_header = [&] () {
-    outfile << "beta" << "\t"
-            << "se" << "\t";
+    if (a_mode != 2)
+      outfile << "beta" << "\t";
+
+    outfile << "se" << "\t";
 
     outfile << "logl_H1" << "\t";  // we may make this an option
 
@@ -135,8 +137,12 @@ void LMM::WriteFiles() {
   };
 
   auto sumstats = [&] (SUMSTAT st) {
-    outfile << scientific << setprecision(6) << st.beta << "\t"
-    << st.se << "\t";
+    outfile << scientific << setprecision(6);
+
+    if (a_mode != 2)
+      outfile << st.beta << "\t";
+
+    outfile << st.se << "\t";
 
     outfile << st.logl_H1 << "\t";
 
-- 
cgit v1.2.3