enrichmentQuality: Best and worst enrichment levels

In benja0x40/MRA.TA: Multi Resolution Analysis of Tiling Arrays

Best and worst enrichment levels

Description

Identify probes most likely associated with best and worst enrichment levels.

Usage

enrichmentQuality(x, y, q.best = 0.015, q.worst = 0.005, plots = F)

Arguments

x	raw Ctr signal intensity
y	raw 4C signal intensity
q.best
q.worst
plots

Author(s)

Benjamin Leblanc

See Also

plotProbeDistanceControls

Examples

x <- 2^pmin(rlnorm(100000, 1, 0.5), 10)
y <- 2^pmin(rlnorm(100000, 0.5, 1), 10)

chk <- enrichmentQuality(x, y, plots=T)

# Visualize thresholds on enrichment quality scores
clr <- with(
  chk, rgb(is.worst, is.best, 0, ifelse(is.worst | is.best, 0.5, 0.1))
)
with(chk, plot(w.score, b.score, pch='.', col=clr))

benja0x40/MRA.TA documentation built on March 13, 2023, 5:15 a.m.