From 96bc4c74f2e5b7a8ffa2bf25bdb6680197fd53ed Mon Sep 17 00:00:00 2001
From: BonfaceKilz
Date: Thu, 18 Feb 2021 16:45:30 +0300
Subject: Create a new method for generating a valid phenotype file on the fly

---
 gn3/computations/__init__.py |  0
 gn3/computations/gemma.py    | 14 ++++++++++++++
 2 files changed, 14 insertions(+)
 create mode 100644 gn3/computations/__init__.py
 create mode 100644 gn3/computations/gemma.py

(limited to 'gn3')

diff --git a/gn3/computations/__init__.py b/gn3/computations/__init__.py
new file mode 100644
index 0000000..e69de29
diff --git a/gn3/computations/gemma.py b/gn3/computations/gemma.py
new file mode 100644
index 0000000..9a972ba
--- /dev/null
+++ b/gn3/computations/gemma.py
@@ -0,0 +1,14 @@
+"""Procedures related gemma computations"""
+
+
+def generate_pheno_txt_file(trait_filename: str,
+                            values: str,
+                            tmpdir: str = "/tmp") -> str:
+    """Given VALUES, and TMPDIR, generate a valide traits file"""
+    with open(f"{tmpdir}/gn2/{trait_filename}", "w") as _file:
+        for value in values:
+            if value == "x":
+                _file.write("NA\n")
+            else:
+                _file.write(f"{value}\n")
+    return f"{tmpdir}/gn2/{trait_filename}"
-- 
cgit v1.2.3