hrStat {PoolBal}

UMP beta p-value pooled statistic

Description

Computes the UMP p-value pooling statistic for a restricted beta family.

Usage

hrStat(p, w = 1)

Arguments

p

numeric vector of p-values between 0 and 1

w

numeric value between 0 and 1

Details

To test the null hypotheses that all p-values are uniform against a restricted beta family 0 < a <= 1 <= b, the most powerful pooled p-value linearly combines upper and lower tail probabilities of the chi-squared distribution with two degrees of freedom with weights w and (1 - w) where w = (1 - a)/(b - a).

This function computes the statistic given by this combination for a collection of p-values, simulation or approximation is required to convert this to a p-value.

Value

A numeric value giving the pooled statistic.

Author(s)

Chris Salahub

Examples

p <- c(0.1, 0.5, 0.9)
hrStat(p, 0.2)
hrStat(p, 0.5)
hrStat(p, 0.9)