exceedance_profile: Profile the data

In Jiefei-Wang/exceedance: Multiple Hypothesis Testing While Controlling the Exceedance Probability of the False Discovery Proportion

Profile the data

Description

Preprocessing the data and collecting information before performing the exceedance control. This step is developed for reducing the computational burder only. It is not related to statistics.

Usage

exceedance_profile(x, param)

Arguments

x	numeric, p-values from some tests (E.g T-test for high-thoughput gene data.), with each element representing a hypothesis.
param	exceedance_parameters object

Value

an exceedance_profile object

Examples

## The 3rd pvalue statistic
param <- param_fast_GW(statistic = "kth_p", param1 = 3)

## generate p-values
x <- rbeta(10, 1, 10)

## profile the data
profile <- profile_pvalue(x,param)


## compute the 95% confidence envolop
alpha <- 0.05

## reject the first three hypotheses
exceedance_confidence(profile, alpha, ri = 3)

## reject the hypothese which pvalues are equal to
## the first three samples.
## In other word, this is equivalent to reject the first three hypotheses
exceedance_confidence(profile, alpha, rx = x[1:3])

## reject the hypotheses which have the lowest 3 p-values
exceedance_confidence(profile, alpha, sri = 3)


## Determine which hypotheses can be rejected while controlling the
## exceedance rate: P(FDP > bound) < alpha
alpha <- 0.05
bound <- 0.2
exceedance_inference(profile, alpha, bound)

Jiefei-Wang/exceedance documentation built on May 11, 2022, 1:43 a.m.