Premium
FastSKAT: Sequence kernel association tests for very large sets of markers
Author(s) -
Lumley Thomas,
Brody Jennifer,
Peloso Gina,
Morrison Alanna,
Rice Kenneth
Publication year - 2018
Publication title -
genetic epidemiology
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.301
H-Index - 98
eISSN - 1098-2272
pISSN - 0741-0395
DOI - 10.1002/gepi.22136
Subject(s) - eigenvalues and eigenvectors , kernel (algebra) , mathematics , covariance matrix , set (abstract data type) , quantile , sequence (biology) , matrix (chemical analysis) , algorithm , statistics , computer science , combinatorics , genetics , biology , physics , materials science , quantum mechanics , composite material , programming language
The sequence kernel association test (SKAT) is widely used to test for associations between a phenotype and a set of genetic variants that are usually rare. Evaluating tail probabilities or quantiles of the null distribution for SKAT requires computing the eigenvalues of a matrix related to the genotype covariance between markers. Extracting the full set of eigenvalues of this matrix (an n × n matrix, for n subjects) has computational complexity proportional to n 3 . As SKAT is often used when n > 10 4 , this step becomes a major bottleneck in its use in practice. We therefore propose fastSKAT, a new computationally inexpensive but accurate approximations to the tail probabilities, in which the k largest eigenvalues of a weighted genotype covariance matrix or the largest singular values of a weighted genotype matrix are extracted, and a single term based on the Satterthwaite approximation is used for the remaining eigenvalues. While the method is not particularly sensitive to the choice of k , we also describe how to choose its value, and show how fastSKAT can automatically alert users to the rare cases where the choice may affect results. As well as providing faster implementation of SKAT, the new method also enables entirely new applications of SKAT that were not possible before; we give examples grouping variants by topologically associating domains, and comparing chromosome‐wide association by class of histone marker.