# test for wqflask/marker_regression/plink_mapping.py import unittest from unittest import mock from gn2.wqflask.marker_regression.plink_mapping import build_line_list from gn2.wqflask.marker_regression.plink_mapping import get_samples_from_ped_file from gn2.wqflask.marker_regression.plink_mapping import flat_files from gn2.wqflask.marker_regression.plink_mapping import gen_pheno_txt_file_plink from gn2.wqflask.marker_regression.plink_mapping import parse_plink_output class AttributeSetter: def __init__(self, obj): for key, val in obj.items(): setattr(self, key, val) class TestPlinkMapping(unittest.TestCase): def test_build_line_list(self): """test for building line list""" line_1 = "this is line one test" irregular_line = " this is an, irregular line " exp_line1 = ["this", "is", "line", "one", "test"] results = build_line_list(irregular_line) self.assertEqual(exp_line1, build_line_list(line_1)) self.assertEqual([], build_line_list()) self.assertEqual(["this", "is", "an,", "irregular", "line"], results) @mock.patch("gn2.wqflask.marker_regression.plink_mapping.flat_files") def test_get_samples_from_ped_file(self, mock_flat_files): """test for getting samples from ped file""" dataset = AttributeSetter({"group": AttributeSetter({"name": "n_1"})}) file_sample = """Expected_1\tline test Expected_2\there Expected_3\tthree""" mock_flat_files.return_value = "/home/user/" with mock.patch("builtins.open", mock.mock_open(read_data=file_sample)) as mock_open: results = get_samples_from_ped_file(dataset) mock_flat_files.assert_called_once_with("mapping") mock_open.assert_called_once_with("/home/user/n_1.ped", "r") self.assertEqual( ["Expected_1", "Expected_2", "Expected_3"], results) @mock.patch("gn2.wqflask.marker_regression.plink_mapping.TMPDIR", "/home/user/data/") @mock.patch("gn2.wqflask.marker_regression.plink_mapping.get_samples_from_ped_file") def test_gen_pheno_txt_file_plink(self, mock_samples): """test for getting gen_pheno txt file""" mock_samples.return_value = ["Expected_1", "Expected_2", "Expected_3"] trait = AttributeSetter({"name": "TX"}) dataset = AttributeSetter({"group": AttributeSetter({"name": "n_1"})}) vals = ["value=K1", "value=K2", "value=K3"] with mock.patch("builtins.open", mock.mock_open()) as mock_open: results = gen_pheno_txt_file_plink(this_trait=trait, dataset=dataset, vals=vals, pheno_filename="ph_file") mock_open.assert_called_once_with( "/home/user/data/ph_file.txt", "wb") filehandler = mock_open() calls_expected = [mock.call('FID\tIID\tTX\n'), mock.call('Expected_1\tExpected_1\tK1\nExpected_2\tExpected_2\tK2\nExpected_3\tExpected_3\tK3\n')] filehandler.write.assert_has_calls(calls_expected) filehandler.close.assert_called_once() @mock.patch("gn2.wqflask.marker_regression.plink_mapping.TMPDIR", "/home/user/data/") @mock.patch("gn2.wqflask.marker_regression.plink_mapping.build_line_list") def test_parse_plink_output(self, mock_line_list): """test for parsing plink output""" chromosomes = [0, 34, 110, 89, 123, 23, 2] species = AttributeSetter( {"name": "S1", "chromosomes": AttributeSetter({"chromosomes": chromosomes})}) fake_file = """0 AACCAT T98.6 0.89\n2 AATA B45 0.3\n121 ACG B56.4 NA""" mock_line_list.side_effect = [["0", "AACCAT", "T98.6", "0.89"], [ "2", "AATA", "B45", "0.3"], ["121", "ACG", "B56.4", "NA"]] with mock.patch("builtins.open", mock.mock_open(read_data=fake_file)) as mock_open: parse_results = parse_plink_output( output_filename="P1_file", species=species) mock_open.assert_called_once_with( "/home/user/data/P1_file.qassoc", "rb") expected = (2, {'AACCAT': 0.89, 'AATA': 0.3}) self.assertEqual(parse_results, expected)