From 2c853b26030dc5e3ff89d7736718abde56c26ca5 Mon Sep 17 00:00:00 2001
From: Pjotr Prins
Date: Thu, 9 Nov 2017 11:00:47 +0000
Subject: Tests: added bslmm test
---
test/test_suite.sh | 14 ++++++++++++++
1 file changed, 14 insertions(+)
(limited to 'test/test_suite.sh')
diff --git a/test/test_suite.sh b/test/test_suite.sh
index dc6053a..0d56e31 100755
--- a/test/test_suite.sh
+++ b/test/test_suite.sh
@@ -2,6 +2,20 @@
gemma=../bin/gemma
+testBslmm() {
+ outn=mouse_hs1940_CD8_bslmm
+ $gemma -g ../example/mouse_hs1940.geno.txt.gz \
+ -p ../example/mouse_hs1940.pheno.txt \
+ -n 2 -a ../example/mouse_hs1940.anno.txt \
+ -bslmm -debug \
+ -o $outn -w 1000 -s 10000 -seed 1
+ assertEquals 0 $?
+ outfn1=output/$outn.hyp.txt
+ outfn2=output/$outn.param.txt
+ assertEquals "45181.93" `perl -nle 'foreach $x (split(/\s+/,$_)) { $sum += sprintf("%.2f",(substr($x,,0,6))) } END { printf "%.2f",$sum }' $outfn1`
+ assertEquals "4043967139.42" `perl -nle 'foreach $x (split(/\s+/,$_)) { $sum += sprintf("%.2f",(substr($x,,0,6))) } END { printf "%.2f",$sum }' $outfn2`
+}
+
testCenteredRelatednessMatrixKFullLOCO1() {
outn=mouse_hs1940_full_LOCO1
$gemma -g ../example/mouse_hs1940.geno.txt.gz \
--
cgit v1.2.3