about summary refs log tree commit diff
diff options
context:
space:
mode:
authorPjotr Prins2017-11-09 11:07:03 +0000
committerPjotr Prins2017-11-09 11:07:03 +0000
commit4c158a9bb45e0dec14436a8a19707f485015bf2e (patch)
tree6cfb2f11f49d78af4f3dc46a3abb6deedf2a834b
parent2c853b26030dc5e3ff89d7736718abde56c26ca5 (diff)
downloadpangemma-4c158a9bb45e0dec14436a8a19707f485015bf2e.tar.gz
Tests: add -lm test
-rwxr-xr-xtest/dev_test_suite.sh17
1 files changed, 17 insertions, 0 deletions
diff --git a/test/dev_test_suite.sh b/test/dev_test_suite.sh
index 284c9aa..0966cba 100755
--- a/test/dev_test_suite.sh
+++ b/test/dev_test_suite.sh
@@ -32,6 +32,23 @@ testBXDStandardRelatednessMatrixK() {
     assertEquals "-116.11" `perl -nle 'foreach $x (split(/\s+/,$_)) { $sum += sprintf("%.2f",(substr($x,,0,6))) } END { printf "%.2f",$sum }' $outfn`
 }
 
+testBXDLMLikelihoodRatio() {
+    outn=BXD_LM_LR
+    $gemma -g ../example/BXD_geno.txt.gz \
+           -p ../example/BXD_pheno.txt \
+           -c ../example/BXD_covariates2.txt \
+           -a ../example/BXD_snps.txt \
+           -k ./output/BXD.cXX.txt \
+           -lm 4 -maf 0.1 \
+           -debug \
+           -o $outn
+    assertEquals 0 $?
+
+    outfn=output/$outn.assoc.txt
+    assertEquals "95134" `wc -w < $outfn`
+    assertEquals "3089042886.28" `perl -nle 'foreach $x (split(/\s+/,$_)) { $sum += sprintf("%.2f",(substr($x,,0,6))) } END { printf "%.2f",$sum }' $outfn`
+}
+
 testBXDLMMLikelihoodRatio() {
     outn=BXD_LMM_LR
     $gemma -g ../example/BXD_geno.txt.gz \