aboutsummaryrefslogtreecommitdiff
path: root/test
diff options
context:
space:
mode:
Diffstat (limited to 'test')
-rwxr-xr-xtest/dev_test_suite.sh8
1 files changed, 5 insertions, 3 deletions
diff --git a/test/dev_test_suite.sh b/test/dev_test_suite.sh
index 0992635..34a1576 100755
--- a/test/dev_test_suite.sh
+++ b/test/dev_test_suite.sh
@@ -158,15 +158,17 @@ testPlinkUnivariateLinearMixedModelLOCO1() {
testCorrelatedPhenotypesMvLLM() {
+ # https://github.com/genetics-statistics/GEMMA/issues/179
+ outn=corrpheno
$gemma $gemmaopts -p data/correlated_phenotypes/Ysim_reg_gemma.txt \
-g data/correlated_phenotypes/Genotypes_gemma.csv \
-d data/correlated_phenotypes/Kinship_eigenval_gemma.txt \
-u data/correlated_phenotypes/Kinship_eigenvec_gemma.txt \
- -lmm 2 -n 1 9 4 6 10 -o corrpheno
+ -lmm 2 -n 1 9 4 6 10 -o $outn
assertEquals 0 $?
- # outfn=output/$outn.assoc.txt
+ outfn=output/$outn.assoc.txt
# assertEquals "68" `wc -l < $outfn`
- # assertEquals "15465346.22" `perl -nle 'foreach $x (split(/\s+/,$_)) { $sum += sprintf("%.2f",(substr($x,,0,6))) } END { printf "%.2f",$sum }' $outfn`
+ assertEquals "777.32" `perl -nle 'foreach $x (split(/\s+/,$_)) { $sum += sprintf("%.2f",(substr($x,,0,6))) } END { printf "%.2f",$sum }' $outfn`
}