plot_mhg: Plot the results of the minimum hypergeometric test.

In slowkow/mhg: Non-parametric rank enrichment test for binary data.

Description

Plot the results of the minimum hypergeometric test.

Usage

plot_mhg(values, x, res, n = 500, main = "", value = "Value", cex = 1.5,
  cex.lab = 1.5, cex.axis = 1.5, cex.main = 2)

Arguments

values	The values corresponding to the ranked vector.
x	A binary vector of successes and failures.
res	A list returned by mhg_test.
n	Plot the first n items in the list.
main	The title of the plot.
value	The name of the value used to rank the items.
cex	The relative size of the legend text.
cex.lab	The relative size of the axis labels.
cex.axis	The relative size of the axis tickmarks.
cex.main	The relative size of the title.

See Also

mhg_test

Examples

# Size of the population.
N <- 5000L
# Successes in the population.
K <- 100L
# Only consider enrichments in the first L observations.
L <- N / 4L
# Require at least X successes in the first L observations.
X <- 5L

set.seed(42)

# Binary vector of successes and failures.
x <- rep(0, N)
x[sample(100, 5)] <- 1
x[sample(200, 10)] <- 1

res <- mhg_test(x, N, K, L, X)

abs(res$pvalue - 1.810658e-05) < 1e-6 # TRUE

# Plot the result.
plot_mhg(sort(rnorm(N)), x, res, L)

slowkow/mhg documentation built on May 30, 2019, 3:06 a.m.