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

knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

library(RVMMAT)

Use a simulated longitudinal binary phenotype and genotype data with random ascertainment for 1000 subjects, each with 5 repeated measures. Each subject has 2 causal SNPs.

p0 = rvmmat_simu(n.sample = 1000, n.time =5)

Estimate null model:

m0 = rvmmat_est(y.long = p0$phe.long, y.cov = p0$phe.cov.long, time = p0$phe.time,phe.model=p0$phe.model)

Perform VAMMT and RVMMAT tests

p_val = rvmmat_test(m0, G = p0$snp.mat)
tail(p_val)

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