LMM: skip beta (again) with LRT

genenetwork · Oct 5, 2017 · 86323cc · 86323cc
1 parent d672c81
commit 86323cc
Show file tree

Hide file tree

Showing 3 changed files with 12 additions and 5 deletions.
diff --git 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";
 

diff --git 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 "87816" `wc -w < $outfn`
+    assertEquals "80498" `wc -w < $outfn`
     assertEquals "3088458212.93" `perl -nle 'foreach $x (split(/\s+/,$_)) { $sum += sprintf("%.2f",(substr($x,,0,6))) } END { printf "%.2f",$sum }' $outfn`
 }
 

diff --git a/test/lengthy_test_suite.sh b/test/lengthy_test_suite.sh
@@ -18,6 +18,7 @@ testPlinkStandardRelatednessMatrixK() {
 }
 
 testPlinkMultivariateLinearMixedModelMultiplePhenotypes_Issue58() {
+    echo "Long running test!"
     # This test passes, but takes over 30 minutes to run!
     # n=2 is original pheno in fam file
     # n=1 is causal1