LSKAT is an R program that performs association testing between a set of SNPs and a longitudinal quantitative trait in population samples. The main reference for this program is
Wang Z., Xu K., Zhang X., Wu X., and Wang Z., (2016) Longitudinal SNP-set association analysis of quantitative phenotypes. Genetic Epidemiology.
To do ...
-
Package SKAT, mvtnorm, snpStats, snowfall.
-
PLINK 1.0 or 1.9 (https://www.cog-genomics.org/plink2)
Please install the required R package before you install LSKAT package. After the installation of SKAT, mvtnorm, snpStats and snowfall package, please install the LSKAT as following steps.
git clone https://github.com/ZWang-Lab/LSKAT.git
R CMD INSTALL LSKAT
-
Please download windows package from (https://github.com/ZWang-Lab/LSKAT/blob/master/LSKAT-win.zip)
-
Install the package in R GUI by selecting the menu "Packages|Install package(s) from local zip files..."
##Usage instructions##
LSKAT is an R package which provides the main functions:
- NULL model estimation to estimate the model parameters and the residuals.
- Gene association test by LSKAT
- LSKAT pipeline for PLINK data set with longitudial phenotype traits.
In general, two ways are common to do data analysis by LSKAT, one is test the genetic association between the longitudinal phenotype traits and SNP-set genotypic matrix, i.e, test LSKAT using the SNP-set matrix given the estimated NULL model which assumes no genetic effects contibute to phenotype traits.
## NULL model estimation
r.model <- longskat_est_model( phe.long.matrix, phe.cov.matrix, phe.time.matrix);
## Gene association test
r.lskat <- longskat_gene_test( r.model, snp.mat);
Other way is to run LSKAT on whole PLINK data set using the pipeline provided by the function longskat_gene_plink
r.lskat <- longskat_gene_plink( file.plink.bed, file.plink.bim, file.plink.fam,
file.phe.long, file.phe.cov, NULL, file.gene.set );
All functions and examples in the LSKAT are available in the manual (https://github.com/ZWang-Lab/LSKAT/blob/master/manual.pdf).