From 186da1d636f88421d1bbca0033e2d99a3b15e6cd Mon Sep 17 00:00:00 2001
From: Pjotr Prins
Date: Fri, 24 Nov 2023 09:43:17 +0100
Subject: Handle genotypes in gemma module

---
 gn/runner/gemma.scm | 4 ++++
 1 file changed, 4 insertions(+)

(limited to 'gn')

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)
-- 
cgit v1.2.3