aboutsummaryrefslogtreecommitdiff
path: root/wqflask/maintenance/get_group_samplelists.py
diff options
context:
space:
mode:
Diffstat (limited to 'wqflask/maintenance/get_group_samplelists.py')
-rw-r--r--wqflask/maintenance/get_group_samplelists.py47
1 files changed, 0 insertions, 47 deletions
diff --git a/wqflask/maintenance/get_group_samplelists.py b/wqflask/maintenance/get_group_samplelists.py
deleted file mode 100644
index 0a450d3f..00000000
--- a/wqflask/maintenance/get_group_samplelists.py
+++ /dev/null
@@ -1,47 +0,0 @@
-import os
-import glob
-import gzip
-
-from base import webqtlConfig
-
-
-def get_samplelist(file_type, geno_file):
- if file_type == "geno":
- return get_samplelist_from_geno(geno_file)
- elif file_type == "plink":
- return get_samplelist_from_plink(geno_file)
-
-
-def get_samplelist_from_geno(genofilename):
- if os.path.isfile(genofilename + '.gz'):
- genofilename += '.gz'
- genofile = gzip.open(genofilename)
- else:
- genofile = open(genofilename)
-
- for line in genofile:
- line = line.strip()
- if not line:
- continue
- if line.startswith(("#", "@")):
- continue
- break
-
- headers = line.split("\t")
-
- if headers[3] == "Mb":
- samplelist = headers[4:]
- else:
- samplelist = headers[3:]
- return samplelist
-
-
-def get_samplelist_from_plink(genofilename):
- genofile = open(genofilename)
-
- samplelist = []
- for line in genofile:
- line = line.split(" ")
- samplelist.append(line[1])
-
- return samplelist