aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorPjotr Prins2023-11-24 09:43:17 +0100
committerPjotr Prins2023-11-24 09:43:17 +0100
commit186da1d636f88421d1bbca0033e2d99a3b15e6cd (patch)
tree9599ae1cc0bed905bb2ff95a33dd9d1cbccca55a
parentfb22ae84c143a27ddc4699327e8aff6c93a0bf6c (diff)
downloadgn-guile-186da1d636f88421d1bbca0033e2d99a3b15e6cd.tar.gz
Handle genotypes in gemma module
-rw-r--r--gn/runner/gemma.scm4
-rwxr-xr-xscripts/precompute/precompute-hits.scm4
2 files changed, 4 insertions, 4 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)
diff --git a/scripts/precompute/precompute-hits.scm b/scripts/precompute/precompute-hits.scm
index e8ba0a2..b9dd746 100755
--- a/scripts/precompute/precompute-hits.scm
+++ b/scripts/precompute/precompute-hits.scm
@@ -7,7 +7,6 @@
(use-modules (dbi dbi)
(gn db mysql)
(gn data dataset)
- (gn data genotype)
(gn data hits)
(gn data strains)
(gn util convert)
@@ -18,9 +17,6 @@
)
-(define inds (geno-inds-bxd "BXD.json"))
-(display inds)
-(assert (= 235 (length inds)))
;; potentially you want to test connection with mysql client:
;;