diff options
Diffstat (limited to 'test')
| -rwxr-xr-x | test/dev_test_suite.sh | 8 |
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` } |