aboutsummaryrefslogtreecommitdiff
path: root/gn/runner
diff options
context:
space:
mode:
authorPjotr Prins2023-11-24 09:43:17 +0100
committerPjotr Prins2023-11-24 09:43:17 +0100
commit186da1d636f88421d1bbca0033e2d99a3b15e6cd (patch)
tree9599ae1cc0bed905bb2ff95a33dd9d1cbccca55a /gn/runner
parentfb22ae84c143a27ddc4699327e8aff6c93a0bf6c (diff)
downloadgn-guile-186da1d636f88421d1bbca0033e2d99a3b15e6cd.tar.gz
Handle genotypes in gemma module
Diffstat (limited to 'gn/runner')
-rw-r--r--gn/runner/gemma.scm4
1 files changed, 4 insertions, 0 deletions
diff --git a/gn/runner/gemma.scm b/gn/runner/gemma.scm
index d9e43eb..67dae54 100644
--- a/gn/runner/gemma.scm
+++ b/gn/runner/gemma.scm
@@ -1,16 +1,20 @@
(define-module (gn runner gemma)
#:use-module (json)
+ #:use-module (gn data genotype)
#:use-module (ice-9 match)
#:use-module (ice-9 format)
#:use-module (ice-9 iconv)
#:use-module (ice-9 receive)
#:use-module (ice-9 string-fun)
+ #:use-module (rnrs base)
#:export (
run-gemma
))
(define (run-gemma data-id name trait-name traits)
+ (define bxd-inds (geno-inds-bxd "BXD.json"))
+ (assert (= 235 (length bxd-inds)))
(if name
(display (string-append "WE HAVE OUR BXD DATASET " name " and trait " trait-name " for precompute!\n")))
(display data-id)