diff --git a/bin/merge.R b/bin/merge.R index 0d349b5..129e41d 100644 --- a/bin/merge.R +++ b/bin/merge.R @@ -32,5 +32,5 @@ difference = formatC(sumb/suma, digits=4) loss = formatC(1 - suma/sumb, digits=2) lengthb = formatC(mean(before$V3-before$V2), digits=4) lengtha = formatC(mean(splitted$length), digits=4) -stats=data.frame(sum_nt_input=sumb, sum_nt_filtered=suma, factor=difference, loss=loss, mean_length_input=lengthb, mean_length_filtered=lengtha) +stats=data.frame(sum_nt_input=sumb, sum_nt_filtered=suma, factor=difference, loss=loss, mean_length_input=lengthb, mean_length_filtered=lengtha, flag_1_ratio=length(which(splitted$containsMaxpos == 1))/dim(splitted)[1]) write.table(stats, "../FilterMotifs.stats", row.names=FALSE, quote=FALSE, sep='\t')