Simulation for LBRAT test

Source: R/lbrat_simu.R

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

lbrat_simu(n.sample = 1000, n.time = 5, par = list(),
  time_cov = TRUE, snp.count = 1000, 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 covaraites and individual specific time dependent random effects
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 LBRAT test