zwl_sim {highDmean}

Apply the test by Zhang and Wang (2020) to multiple simulated two-sample datasets

Description

Apply the two-sample high-dimensional test by Zhang and Wang (2020) to multiple simulated two-sample high dimensional datasets. This function is useful for Monte Carlo experiments.

Usage

zwl_sim(DATA, order = 0)

Arguments

DATA

The list of dataset lists generated by buildData.

order

The order of the center correction. Possible choices are 0, 2. To use the ZWLm test, set order = 0; to use the ZWL test, set order = 2. THE ZWLm test is recommended for moderate sample sizes.

Value

A dataframe with each row consisting the values of the test statistics, p-values, Tn, and the estimate of Var(Tn).

References

Zhang, H. and Wang, H. (2020). Result consistency of high dimensional two-sample tests applied to gene ontology terms with gene sets. Manuscript in review.

Examples

# Generate 3 simulated two-sample datasets and apply the ZWL test
data <- buildData(n = 45, m =60, p = 300,
          muX = rep(0,300), muY = rep(0,300),
          dep = 'IND', S = 3, innov = rnorm)
zwl_sim(data, order = 2)