Handle genotypes in gemma module

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)