diff --git a/scripts/filter_and_BH.R b/scripts/filter_and_BH.R index ba35bae..f1a24b3 100644 --- a/scripts/filter_and_BH.R +++ b/scripts/filter_and_BH.R @@ -33,7 +33,7 @@ regtools_data_filtered_sorted = regtools_data_filtered[order(adjusted_p)] write_file = paste("compare_junctions/hist/", "junction_pvalues_filtered_BH", tag, ".tsv", sep="") write.table(regtools_data_filtered_sorted, file=write_file, quote=FALSE, sep='\t', row.names = FALSE) -threshold = 0.025 +threshold = 0.05 is_significant = regtools_data_filtered_sorted$adjusted_p < threshold regtools_data_significant_filtered_sorted = regtools_data_filtered_sorted[is_significant]