fix header im variant correlation

kircherlab · Nov 25, 2024 · 0f736e3 · 0f736e3
1 parent 271f16e
commit 0f736e3
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/workflow/scripts/variants/correlateVariantTables.py b/workflow/scripts/variants/correlateVariantTables.py
@@ -57,7 +57,7 @@ def filterOnThreshold(variants, threshold):
                 variants_2 = filterOnThreshold(variants_2, bc_threshold)
 
             click.echo("Join variants file...")
-            variants_join = variants_1.join(variants_2, how="inner", lsuffix='_A', rsuffix='_B')[["log2_expression_A", "log2_expression_B"]]
+            variants_join = variants_1.join(variants_2, how="inner", lsuffix='_A', rsuffix='_B')[["log2FoldChange_expression_A", "log2FoldChange_expression_B"]]
 
             output = pd.concat([output, pd.DataFrame([[condition, rep_1, rep_2, variants_join.shape[0], bc_threshold, variants_join.corr(method="pearson").iloc[0,1],variants_join.corr(method="spearman").iloc[0,1]]])], ignore_index=True)