Simu_null_proxy: Simulation of the null statistics

In william-denault/WaveletScreaming: A package to perform Wavelet Screening on GWAS data

Description

Simulation of the null statistics

Usage

Simu_null_proxy(
  Y,
  confounder,
  lev_res,
  size = 10000,
  sigma_b = NA,
  base_shrink,
  coeftype = "c",
  verbose = TRUE
)

Arguments

Y	a vector of numeric values used in the wavelet screening function for association (a phenotype or simulated phenotype with the same distribution).
confounder	the confounding matrix with the same sample order as Y. The intercept should not be included if missing will generate an intercept matrix.
lev_res	the level of resolution in the wavelet transform
size	number of simulation to be performed. If not specified set at 10000
sigma_b	the parameter of the NIG prior used for the Betas computation.
base_shrink	numeric, value used in the thresholding of the proportion of assocation, if non specificed set up as 1/sqrt(2*log(sample_size)
coeftype	type of wavelet coefficient used for the screening (choice "c" or "d"). If missing set as "c"
verbose	logical parameter, set as TRUE by default. ID

Value

The simulation under the null of the two test statistics used to build the final test (i.e., L_h and min(ph,pv))

Examples

## Not run: 
Y <- rnorm(4000)
Sim <- Simu_null_proxy(Y=Y,lev_res = 6,sigma_b = 0.2,size=10000)

## End(Not run)

william-denault/WaveletScreaming documentation built on Jan. 23, 2021, 12:34 p.m.