diff options
Diffstat (limited to 'test/lmm-integration-tests.scm')
| -rw-r--r-- | test/lmm-integration-tests.scm | 52 |
1 files changed, 0 insertions, 52 deletions
diff --git a/test/lmm-integration-tests.scm b/test/lmm-integration-tests.scm deleted file mode 100644 index 90858ce..0000000 --- a/test/lmm-integration-tests.scm +++ /dev/null @@ -1,52 +0,0 @@ -#!/bin/sh -# -*- mode: scheme; -*- -exec guile --debug -s "$0" "$@" -!# - -(define-module (test-runner) - #:use-module (ice-9 match) - #:use-module (srfi srfi-1) ; for last - #:use-module (srfi srfi-13) - #:use-module (srfi srfi-64) ; for tests - #:use-modules (ice-9 ftw) ; delete file - #:use-module (ice-9 rdelim) - ) - -(define kinship-fn "./output/mouse_hs1940.cXX.txt") -(define gwa-fn "./output/mouse_hs1940.assoc.txt") - -(test-begin "uvlmm-kinship-run") - -(when (file-exists? kinship-fn) - (delete-file kinship-fn)) -(let [(err (system "./build/bin/Debug/gemma -g ./example/mouse_hs1940.geno.txt.gz -gk -o mouse_hs1940 -debug"))] - (test-eqv 0 err)) - -(test-end "uvlmm-kinship-run") - -(test-begin "uvlmm-gwa-run") - -(when (file-exists? gwa-fn) - (delete-file gwa-fn)) -;; The following integration test runs gemma uvlmm and adds up the output column as a check. -;; It uses the kinship-run matrix from the earlier test -(let [(err (system (string-append "./build/bin/Debug/gemma -g ./example/mouse_hs1940.geno.txt.gz -p ./example/mouse_hs1940.pheno.txt -n 1 -a ./example/mouse_hs1940.anno.txt -k " kinship-fn " -o mouse_hs1940 -lmm 9 -debug")))] - (test-eqv 0 err)) -(call-with-input-file gwa-fn - (lambda (port) - (read-line port) ; skip first line - (let* ((fields (string-split (read-line port) #\tab)) - (last-field (last fields))) - (test-eqv 208.0 (truncate (* 1000 (string->number last-field))))) - (test-eqv 5720672.0 - (let loop ((line (read-line port)) - (sum 208.0)) - (if (eof-object? line) - sum - (let* ((fields (string-split line #\tab)) - (last-field (last fields)) - (value (string->number last-field))) - (loop (read-line port) - (+ sum (truncate (* 1000 value)))))))))) - -(test-end "uvlmm-gwa-run") |
