gwasurvivr: an R package for genome wide survival analysis

Abstract

SummaryTo address the limited software options for performing survival analyses with millions of SNPs, we developed gwasurvivr, an R/Bioconductor package with a simple interface for conducting genome wide survival analyses using VCF (outputted from Michigan or Sanger imputation servers), IMPUTE2 or PLINK files. To decrease the number of iterations needed for convergence when optimizing the parameter estimates in the Cox model we modified the R package survival; covariates in the model are first fit without the SNP, and those parameter estimates are used as initial points. We benchmarked gwasurvivr with other software capable of conducting genome wide survival analysis (genipe, SurvivalGWAS_SV, and GWASTools). gwasurvivr is significantly faster and shows better scalability as sample size, number of SNPs and number of covariates increases.\n\nAvailability and implementationgwasurvivr, including source code, documentation, and vignette are available at: http://bioconductor.org/packages/gwasurvivr\n\nContactAbbas Rizvi, rizvi.33@osu.edu; Lara E Sucheston-Campbell, suchestoncampbell.1@osu.edu\n\nSupplementary information: Supplementary data are available at https://github.com/suchestoncampbelllab/gwasurvivr_manuscript