Skip to contents

Create an object of class 'SansSouci'

Source: R/SansSouci-class.R
SansSouci.Rd

Create an object of class 'SansSouci'

Usage

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)