From fb22ae84c143a27ddc4699327e8aff6c93a0bf6c Mon Sep 17 00:00:00 2001
From: Pjotr Prins
Date: Fri, 24 Nov 2023 09:07:55 +0100
Subject: Fetch genotype inds from json file

---
 gn/runner/gemma.scm | 9 ++++++---
 1 file changed, 6 insertions(+), 3 deletions(-)

(limited to 'gn/runner')

diff --git a/gn/runner/gemma.scm b/gn/runner/gemma.scm
index 31a14f7..d9e43eb 100644
--- a/gn/runner/gemma.scm
+++ b/gn/runner/gemma.scm
@@ -17,13 +17,16 @@
   (display traits)
   (newline)
   ;; ---- write phenotype file
-  (call-with-output-file "test.tmp"
+  (call-with-output-file "pheno.txt"
     (lambda (port)
       (write 12 port)
       (newline port)
       (write "HELLO" port)
       ))
+  ;; set up with ./.guix-shell -- guile -L . -s ./scripts/precompute/precompute-hits.scm
 
-  ;; ---- start GEMMA precompute
-  (system "gemma")
+  ;; ---- to start GEMMA precompute inside container
+  ;; env LD_LIBRARY_PATH=$GUIX_ENVIRONMENT/lib/ guile -L . -s ./scripts/precompute/precompute-hits.scm
+  ;; --- First we compute K
+  (system (string-append "env GEMMA_COMMAND=/gemma/bin/gemma /gemma-wrapper/bin/gemma-wrapper --debug -- -gk -g BXD.8_geno.txt.gz -p pheno.txt -a BXD.8_snps.txt"  ))
   )
-- 
cgit v1.2.3