Closed mdozmorov closed 9 years ago
Quantiles are better than % threshold.
mtx <- read.table("data/E001-H3K4me1.broadPeak.gz") hist(mtx$V7, n=100)
summary(mtx$V7) quantile(mtx$V7)
print(c(max(mtx$V7) * 0.25, max(mtx$V7) * 0.5, max(mtx$V7) * 0.75))
Quantiles are better than % threshold.
mtx <- read.table("data/E001-H3K4me1.broadPeak.gz") hist(mtx$V7, n=100)
summary(mtx$V7) quantile(mtx$V7)
print(c(max(mtx$V7) * 0.25, max(mtx$V7) * 0.5, max(mtx$V7) * 0.75))