From 25ad9756ebecfdb2f01b6f87c35bd731e3a3186d Mon Sep 17 00:00:00 2001
From: Pjotr Prins
Date: Wed, 11 Oct 2017 12:51:51 +0000
Subject: Fix header for not showing se

---
 src/lmm.cpp            | 6 +++---
 test/dev_test_suite.sh | 2 +-
 2 files changed, 4 insertions(+), 4 deletions(-)

diff --git a/src/lmm.cpp b/src/lmm.cpp
index 68d3c4f..a49c8c5 100644
--- a/src/lmm.cpp
+++ b/src/lmm.cpp
@@ -104,10 +104,10 @@ void LMM::WriteFiles() {
   }
 
   auto common_header = [&] () {
-    if (a_mode != 2)
+    if (a_mode != 2) {
       outfile << "beta" << "\t";
-
-    outfile << "se" << "\t";
+      outfile << "se" << "\t";
+    }
 
     outfile << "logl_H1" << "\t";  // we may make this an option
 
diff --git a/test/dev_test_suite.sh b/test/dev_test_suite.sh
index 0fc4423..136ba76 100755
--- a/test/dev_test_suite.sh
+++ b/test/dev_test_suite.sh
@@ -44,7 +44,7 @@ testBXDLMMLikelihoodRatio() {
     assertEquals 0 $?
 
     outfn=output/$outn.assoc.txt
-    assertEquals "80498" `wc -w < $outfn`
+    assertEquals "73180" `wc -w < $outfn`
     assertEquals "3088458212.93" `perl -nle 'foreach $x (split(/\s+/,$_)) { $sum += sprintf("%.2f",(substr($x,,0,6))) } END { printf "%.2f",$sum }' $outfn`
 }
 
-- 
cgit v1.2.3