test.hc: Multiple comparison test using Higher Criticism (HC)...

View source: R/test.hc.R

test.hcR Documentation

Multiple comparison test using Higher Criticism (HC) statitics.

Description

Multiple comparison test using Higher Criticism (HC) statitics.

Usage

test.hc(prob, M, k0, k1, onesided = FALSE, method = "ecc", ei = NULL)

Arguments

prob

- vector of input p-values.

M

- correlation matrix of input statistics (of the input p-values).

k0

- search range starts from the k0th smallest p-value.

k1

- search range ends at the k1th smallest p-value.

onesided

- TRUE if the input p-values are one-sided.

method

- default = "ecc": the effective correlation coefficient method in reference 2. "ave": the average method in reference 3, which is an earlier version of reference 2. The "ecc" method is more accurate and numerically stable than "ave" method.

ei

- the eigenvalues of M if available.

Value

pvalue - The p-value of the HC test.

hcstat - HC statistic.

location - the order of the input p-values to obtain HC statistic.

References

1. Hong Zhang, Jiashun Jin and Zheyang Wu. "Distributions and power of optimal signal-detection statistics in finite case", IEEE Transactions on Signal Processing (2020) 68, 1021-1033 2. Hong Zhang and Zheyang Wu. "The general goodness-of-fit tests for correlated data", Computational Statistics & Data Analysis (2022) 167, 107379 3. Hong Zhang and Zheyang Wu. "Generalized Goodness-Of-Fit Tests for Correlated Data", arXiv:1806.03668. 4. Donoho, David; Jin, Jiashun. "Higher criticism for detecting sparse heterogeneous mixtures". Annals of Statistics 32 (2004).

See Also

stat.hc for the definition of the statistic.

Examples

pval.test = runif(10)
test.hc(pval.test, M=diag(10), k0=1, k1=10)
#When the input are statistics#
stat.test = rnorm(20)
p.test = 2*(1 - pnorm(abs(stat.test)))
test.hc(p.test, M=diag(20), k0=1, k1=10)

SetTest documentation built on Sept. 12, 2024, 7:41 a.m.

Related to test.hc in SetTest...