The sequence kernel association test for the proportional odds model

Motivation The Sequence Kernel Association Test (SKAT) and its extensions are the most popular methods for studying the association between phenotypes and a set of single nucleotide polymorphisms. Their practical application is very wide, but most of these methods are designed for continuous and binary phenotypes. Ordered categorical phenotypes are also very common in practice, so there is an urgent need to develop SKAT-type tests for proportional odds model. Results To accommodate ordered categorical phenotypes, we propose a test named the Sequence Kernel Association Test for the Proportional Odds Model (POM-SKAT). It constructs a score test for the variance of the coefficients of interest using a quasi-likelihood and the P-value is evaluated by approximating the asymptotic distribution of the test statistic with the Pearson Type III distribution. Simulation studies demonstrate that our method performs well and achieves high power in detecting gene-phenotype associations. We apply POM-SKAT to rheumatoid arthritis data provided by Genetic Analysis Workshop 16, identifying multiple relevant gene variants. Availability and implementation Code is available at GitHub (http://github.com/amss-stat/POM-SKAT).