rvmmat_simu: Simulation for RVMMAT test

In ZWang-Lab/RVMMAT: Retrospective Varying Coefficient Mixed Model Association Test

Simulation for RVMMAT test

Description

This function use pre-defined parameters to make the simulation data for the RVMMAT test (including type I and power test)

Usage

rvmmat_simu(
  n.sample = 1000,
  n.time = 5,
  par = list(),
  genetic.effect = 0.6,
  time_cov = TRUE,
  snp.count = 100,
  intercept = TRUE,
  disease.para = list(),
  onlypower = FALSE,
  phe.model = "logistic",
  oversampling = "random"
)

Arguments

n.sample	Numeric, sample size, number of individuals
n.time	Numeric, number of measurements for each individual
par	List, the parameters for the phenotype traits, including covariates and individual specific time dependent random effects
genetic.effect	Numeric, genetic effect size (gamma) in genetic effect function
time_cov	Logical variable, indicating whether time effect is included in phenotypic traits
snp.count	Numeric, number of SNPs
intercept	Logical variable, indicating whether intercept is used in phenotypic traits
disease.para	List, the parameters for disease allele and its effect size for power simulation
onlypower	Logical variable, indicating whether include disease SNPs in the generated SNPs
phe.model	String, the phenotype model, two optional values: 'logistic', 'liability'
oversampling	String, the ascertainment scheme, three optional value: 'random', 'baseline', 'sum'

Value

A list object is returned to be used as object for RVMMAT test

ZWang-Lab/RVMMAT documentation built on Sept. 27, 2024, 5:19 p.m.