aboutsummaryrefslogtreecommitdiff
path: root/test
diff options
context:
space:
mode:
authorPjotr Prins2017-11-09 11:43:53 +0000
committerPjotr Prins2017-11-09 11:43:53 +0000
commitdbce439e7196ab6882634b270f2675b3ebcb33cd (patch)
tree4c593fa5594d1482ccbb57bf4d799e9d0529aebe /test
parent4c158a9bb45e0dec14436a8a19707f485015bf2e (diff)
downloadpangemma-dbce439e7196ab6882634b270f2675b3ebcb33cd.tar.gz
Adding tests for bslmm
Diffstat (limited to 'test')
-rwxr-xr-xtest/test_suite.sh23
1 files changed, 23 insertions, 0 deletions
diff --git a/test/test_suite.sh b/test/test_suite.sh
index 0d56e31..58e5574 100755
--- a/test/test_suite.sh
+++ b/test/test_suite.sh
@@ -16,6 +16,29 @@ testBslmm() {
assertEquals "4043967139.42" `perl -nle 'foreach $x (split(/\s+/,$_)) { $sum += sprintf("%.2f",(substr($x,,0,6))) } END { printf "%.2f",$sum }' $outfn2`
}
+testBslmm2() {
+ $gemma -g ../example/mouse_hs1940.geno.txt.gz \
+ -p ../example/mouse_hs1940.pheno.txt \
+ -n 2 -debug \
+ -a ../example/mouse_hs1940.anno.txt \
+ -gk 1 -o mouse_hs1940_CD8_train
+ assertEquals 0 $?
+ $gemma -g ../example/mouse_hs1940.geno.txt.gz \
+ -p ../example/mouse_hs1940.pheno.txt \
+ -n 2 -debug \
+ -epm ./output/mouse_hs1940_CD8_bslmm.param.txt \
+ -emu ./output/mouse_hs1940_CD8_bslmm.log.txt \
+ -ebv ./output/mouse_hs1940_CD8_bslmm.bv.txt \
+ -k ./output/mouse_hs1940_CD8_train.cXX.txt \
+ -predict \
+ -o mouse_hs1940_CD8_prdt_k
+ assertEquals 0 $?
+ outfn=output/$outn.hyp.txt
+ assertEquals "45181.93" `perl -nle 'foreach $x (split(/\s+/,$_)) { $sum += sprintf("%.2f",(substr($x,,0,6))) } END { printf "%.2f",$sum }' $outfn`
+ ../bin/gemma -g mouse_hs1940.geno.txt.gz -p mouse_hs1940.pheno.txt -n 4 -epm ./output/mouse_hs1940_CD8_bslmm_cc1.param.txt -emu ./output/mouse_hs1940_CD8_bslmm_cc1.log.txt -predict -o mouse_hs1940_CD8_prdt_cc1
+ exit 1
+}
+
testCenteredRelatednessMatrixKFullLOCO1() {
outn=mouse_hs1940_full_LOCO1
$gemma -g ../example/mouse_hs1940.geno.txt.gz \