about summary refs log tree commit diff
diff options
context:
space:
mode:
authorPjotr Prins2017-10-11 12:51:51 +0000
committerPjotr Prins2017-10-13 15:27:24 +0000
commit25ad9756ebecfdb2f01b6f87c35bd731e3a3186d (patch)
tree4ace2989827d9caf75839742cf013248b5cedae5
parent6ae92af3335be3c6abb2588a7c8f14239a6f0a02 (diff)
downloadpangemma-25ad9756ebecfdb2f01b6f87c35bd731e3a3186d.tar.gz
Fix header for not showing se
-rw-r--r--src/lmm.cpp6
-rwxr-xr-xtest/dev_test_suite.sh2
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`
 }