diff options
author | Pjotr Prins | 2017-11-09 11:43:53 +0000 |
---|---|---|
committer | Pjotr Prins | 2017-11-09 11:43:53 +0000 |
commit | dbce439e7196ab6882634b270f2675b3ebcb33cd (patch) | |
tree | 4c593fa5594d1482ccbb57bf4d799e9d0529aebe /test | |
parent | 4c158a9bb45e0dec14436a8a19707f485015bf2e (diff) | |
download | pangemma-dbce439e7196ab6882634b270f2675b3ebcb33cd.tar.gz |
Adding tests for bslmm
Diffstat (limited to 'test')
-rwxr-xr-x | test/test_suite.sh | 23 |
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 \ |