View source: R/byExtremality.R
byExtremality | R Documentation |
max(SD(X_j)) if X_j ~ Beta(a, b) < max(SD(X_j)) if X_j ~ Bernoulli(a/(a+b))
byExtremality(x, k = 500)
x |
a matrix of beta values, or a GenomicRatioSet |
k |
how many rows to return (500) |
for X having a known mean and SD, hence solvable for a + b by MoM, define
extremality = sd(X_j) / bernoulliSD(mean(X_j))
This function selects the k most extremal rows of x and returns their values.
the most extremal k rows of x (returns the same class as x)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.