Skip to content

Commit 2c609c1

Browse files
committed
feat(pipelines): use BipEx case/control counts from ht globals
1 parent b68a616 commit 2c609c1

File tree

1 file changed

+4
-3
lines changed

1 file changed

+4
-3
lines changed

data_pipeline/data_pipeline/datasets/bipex/bipex_gene_results.py

Lines changed: 4 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -17,6 +17,8 @@ def prepare_gene_results(test_genes, _output_root):
1717
if test_genes:
1818
results = filter_results_table_to_test_gene(results)
1919

20+
n_cases = hl.eval(results.globals["case_total"])
21+
n_controls = hl.eval(results.globals["control_total"])
2022
results = results.select_globals()
2123

2224
results = results.annotate(
@@ -50,9 +52,8 @@ def prepare_gene_results(test_genes, _output_root):
5052

5153
results = results.annotate(
5254
gene_id=gene_model_ht[results["gene_symbol"]].gene_id,
53-
# FIXME: suggest anlyst include this in input file, remove this when they do
54-
n_cases=65_018,
55-
n_controls=169_631,
55+
n_cases=n_cases,
56+
n_controls=n_controls,
5657
)
5758

5859
results = results.group_by("gene_id").aggregate(group_results=hl.agg.collect(results.row.drop("gene_id")))

0 commit comments

Comments
 (0)