aboutsummaryrefslogtreecommitdiff
path: root/r_qtl/r_qtl2.py
blob: 94fa842d4f83174c742728c1de1bb601c2a8c347 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
"""The R/qtl2 parsing and processing code."""
import io
import csv
import json
import yaml
from pathlib import Path
from functools import reduce
from typing import Any, List, Union, Iterator
from zipfile import ZipFile, ZipInfo, is_zipfile

from r_qtl.errors import InvalidFormat

def thread_op(value, *functions):
    """Thread the `value` through the sequence of `functions`."""
    return reduce(lambda result, func: func(result), functions, value)

def control_data(zfile: ZipFile) -> dict:
    """Retrieve the control file from the zip file info."""
    files = tuple(filename
                  for filename in zfile.namelist()
                  if (filename.endswith(".yaml") or filename.endswith(".json")))
    num_files = len(files)
    if num_files == 0:
        raise InvalidFormat("Expected a json or yaml control file.")

    if num_files > 1:
        raise InvalidFormat("Found more than one possible control file.")

    return (json.loads(zfile.read(files[0]))
            if files[0].endswith(".json")
            else yaml.safe_load(zfile.read(files[0])))

def genotype_metadata(zfile: ZipFile, cdata: dict) -> dict:
    """Read Individual ID key and the marker names."""
    # TODO: Handle transposed files
    line_num = 0
    with zfile.open(cdata["geno"]) as genofile:
        for line in filter(lambda line: not line.startswith("#"),
                           io.TextIOWrapper(genofile)):
            line_parts = line.strip().split(cdata.get("sep", ","))
            return {
                "individual_id_key": line_parts[0].strip(),
                "markers": tuple(marker.strip() for marker in line_parts[1:])
            }

def genotype_data(zfile: ZipFile, cdata: dict) -> Iterator[dict]:
    """Load the genotype file, making use of the control data."""
    # TODO: Handle transposed files
    with zfile.open(cdata["geno"]) as genofile:
        reader = csv.DictReader(filter(lambda line: not line.startswith("#"),
                                       io.TextIOWrapper(genofile)),
                                delimiter=cdata.get("sep", ","))
        for row in reader:
            yield {
                key: thread_op(
                    value,
                    # replace genotype codes
                    lambda val: cdata["genotypes"].get(val, val),
                    # replace N/A strings
                    lambda val: (None if val in cdata["na.strings"] else val))
                for key,value
                in row.items()
            }

def genetic_map_data(zfile: ZipFile, cdata: dict) -> dict:
    """Read gmap files to get the genome mapping data"""
    if not cdata.get("geno_transposed", False):
        with zfile.open(cdata["gmap"]) as gmapfile:
            reader = csv.DictReader(
                filter(lambda line: not line.startswith("#"),
                       io.TextIOWrapper(gmapfile)),
                delimiter=cdata.get("sep", ","))
            return {
                line["marker"]: {
                    key: value for key,value in line.items() if key != "marker"
                } for line in reader
            }

    with zfile.open(cdata["gmap"]) as gmapfile:
        lines = [[field.strip() for field in
                  line.strip().split(cdata.get("sep", ","))]
                 for line in
                 filter(lambda line: not line.startswith("#"),
                        io.TextIOWrapper(gmapfile))]

    headers = tuple(line[0] for line in lines)
    return reduce(
        lambda gmap, row: {
            **gmap,
            row[0]: dict(zip(headers[1:], row[1:]))},
        zip(*(line[1:] for line in lines)),
        {})

def read_r_qtl2_files(filepath: Path):
    """Read R/qtl2 format zip files."""
    with ZipFile(filepath, "r") as zfile:
        cf = control_data(zfile)
        raise NotImplementedError("Implementation is incomplete.")