diff --git a/Scripts/filterPausing.py b/Scripts/filterPausing.py index aca0705..04dcb60 100644 --- a/Scripts/filterPausing.py +++ b/Scripts/filterPausing.py @@ -2,6 +2,7 @@ # -*- coding: utf-8 -*- import os import pandas as pd +import statsmodels.stats.multitest as mtest dfs=[] for fn in snakemake.input: @@ -10,5 +11,7 @@ df=pd.read_csv(fn,delimiter='\t',header=None) dfs.append(df) df=pd.concat(dfs,sort=[0,1],ignore_index=True) +df[4]=1.-df[4]/100 +df[4]=mtest.fdrcorrection(df[4])[1] df.to_csv(snakemake.output[1],sep='\t',index=False,header=False) -df[df[4]>=snakemake.params[0]].to_csv(snakemake.output[0],sep='\t',index=False,header=False) +df[df[4]