QrhoIC {ICSKAT}

QrhoIC.R

Description

Calculate the test statistic, distribution, and p-value for each value of Krho in SKATO.

Usage

QrhoIC(
  rhoVec,
  icskatOut,
  liu = TRUE,
  bootstrapOut = NULL,
  alwaysCentral = FALSE
)

Arguments

rhoVec

Numeric vector of the rho values to use in SKATO.

icskatOut

The output list returned from a call to ICSKAT().

liu

Boolean for whether to use Liu (TRUE) or Davies (FALSE) method in calculating p-values for each Qrho. Default is Liu, following SKAT package. If wanting to use bootstrap moments for Qrho, need to use Liu method.

bootstrapOut

The output (a list) from a call the ICSKATO_bootstrap() function, holding moments for Qrho.

alwaysCentral

A boolean, if TRUE, follow SKAT package practice of always setting delta=0 in chi-square moment matching.

Value

Data frame holding the SKAT pvalue + test statistic for each fixed rho, the matched noncentrality + degree of freedom parameters for each fixed rho (using both bootstrap and analytic calculation), and the mean and variance of each Qrho using both bootstrap and analytic calculation.