MiHC.stat: uHC(h)'s and wHC(h)'s

In hk1785/MiHC: Microbiome Higher Criticism Analysis

uHC(h)'s and wHC(h)'s

Description

This function produces the values of the test statistics of uHC(h)'s and wHC(h)'s.

Usage

MiHC.stat(Zs, hs=hs, Ws=Ws)

Arguments

Zs	The Z-score values for individual OTUs.
hs	A vector of the candidate modulation schema for lower sparsity levels. Default is hc=c(1,3,5,7,9).
Ws	A vector of the weights for invidual OTUs.

Value

The values of the test statistics of uHC(h)'s and wHC(h)'s.

Author(s)

Hyunwook Koh

References

Koh and Zhao. A powerful microbial group association test based on the higher criticism analysis for sparse microbial association signals. (Under revision).

Examples

# Import requisite R packages
require(cluster)
require(permute)
require(phyloseq)

Zs <- rnorm(50)
hs <- c(1,3,5,7,9)
Ws <- runif(50)

MiHC.stat(Zs, hs=hs, Ws=Ws)

hk1785/MiHC documentation built on Dec. 21, 2024, 4:29 a.m.