diff options
| author | Arun Isaac | 2025-11-13 02:20:27 +0000 |
|---|---|---|
| committer | Arun Isaac | 2025-11-13 02:20:27 +0000 |
| commit | 530652ab438a977c500df896d30b163c511dfc4a (patch) | |
| tree | 0234f74c1508409eb37db311081021b7a9beb6e2 /e2e-tests/hsmice/check-qtl.py | |
| parent | fc0e13275d5d96cc6dd0974582d888156b91734f (diff) | |
| download | pyhegp-530652ab438a977c500df896d30b163c511dfc4a.tar.gz pyhegp-530652ab438a977c500df896d30b163c511dfc4a.tar.lz pyhegp-530652ab438a977c500df896d30b163c511dfc4a.zip | |
Parameterize query when checking for the QTL.
This will allow generalization to JWAS later.
Diffstat (limited to 'e2e-tests/hsmice/check-qtl.py')
| -rw-r--r-- | e2e-tests/hsmice/check-qtl.py | 5 |
1 files changed, 4 insertions, 1 deletions
diff --git a/e2e-tests/hsmice/check-qtl.py b/e2e-tests/hsmice/check-qtl.py index 6e342b1..f7dd739 100644 --- a/e2e-tests/hsmice/check-qtl.py +++ b/e2e-tests/hsmice/check-qtl.py @@ -20,9 +20,12 @@ import sys import pandas as pd +data_file = sys.argv[1] +query_expression = sys.argv[2] + if __name__ == "__main__": df = pd.read_csv(sys.argv[1], sep="\t") - qtl = df.query("p < 1e-10") + qtl = df.query(query_expression) # Assert that the QTL is on chromosome 4. assert (qtl.chromosome == 4).all() # Assert that the QTL is within 2 Mb of the expected position. |