Create an object of class 'SansSouci'
SansSouci(Y, groups, truth = NULL)Arguments
- Y
A matrix of \(m\) variables (hypotheses) by \(n\) observations
- groups
A numeric vector of \(n\) values in \(0, 1\), the groups of observations on which to perform two-sample tests
- truth
An optional numeric vector of $m$ values in $0,1$, the status of each null hypothesis (0 means H0 is true, 1 means H1 is true). Typically used in simulations.
Value
An object of class "SansSouci"
Examples
data(expr_ALL, package = "sanssouci.data")
groups <- ifelse(colnames(expr_ALL) == "NEG", 0, 1)
table(groups)
#> groups
#> 0 1
#> 42 37
a <- SansSouci(Y = expr_ALL, groups = groups)
res <- fit(a, B = 100, alpha = 0.1)
plot(res, xmax = 500)
res_beta <- fit(res, B = 100, alpha = 0.1, family = "Beta", K = 10)
plot(res_beta, xmax = 500)
