aboutsummaryrefslogtreecommitdiff
path: root/tests/unit/computations/test_gemma.py
blob: 2d2fbf9a363cf119d36759e84cd2eb16baedbc6f (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
"""Test cases for procedures defined in computations.gemma"""
import unittest

from unittest import mock
from gn3.computations.gemma import generate_hash_of_string
from gn3.computations.gemma import generate_pheno_txt_file
from gn3.computations.gemma import generate_gemma_computation_cmd


class TestGemma(unittest.TestCase):
    """Test cases for computations.gemma module"""
    def test_generate_pheno_txt_file(self):
        """Test that the pheno text file is generated correctly"""
        open_mock = mock.mock_open()
        with mock.patch("gn3.computations.gemma.open", open_mock, create=True):
            _file = generate_pheno_txt_file(tmpdir="/tmp",
                                            trait_filename="phenotype.txt",
                                            values=["x", "x", "BXD07 438.700"])
            self.assertEqual(_file, ("/tmp/gn2/phenotype_"
                                     "P7y6QWnwBPedSZdL0+m/GQ.txt"))
        open_mock.assert_called_with(("/tmp/gn2/phenotype_"
                                      "P7y6QWnwBPedSZdL0+m/GQ.txt"), "w")
        open_mock.return_value.write.assert_has_calls([
            mock.call("NA\n"),
            mock.call("NA\n"),
            mock.call("BXD07 438.700\n")
        ])

    def test_generate_hash_of_string(self):
        """Test that a string is hashed correctly"""
        self.assertEqual(generate_hash_of_string("I^iQP&TlSR^z"),
                         "hMVRw8kbEp49rOmoIkhMjA")

    @mock.patch("gn3.computations.gemma.do_paths_exist")
    def test_compose_k_computation_cmd(self, mock_pathsp):
        """Test that a K computation cmd is constructed properly"""
        mock_pathsp.return_value = True
        self.assertEqual(
            generate_gemma_computation_cmd(
                gemma_cmd="gemma-wrapper",
                gemma_wrapper_kwargs=None,
                gemma_kwargs={
                    "geno_filename": "genofile.txt",
                    "trait_filename": "test.txt",
                    "covar_filename": "genofile_snps.txt"
                },
                output_file="/tmp/gn2/k_output_gUFhGu4rLG7k+CXLPk1OUg.txt",
            ), ("gemma-wrapper --json -- "
                "-g genofile.txt -p "
                "test.txt -a genofile_snps.txt "
                "-gk > /tmp/gn2/"
                "k_output_gUFhGu4rLG7k+CXLPk1OUg.txt"))