kalis (doi: 10.1186/s12859-024-05688-8) provides a high performance implementation of the Li & Stephens model (doi: 10.1093/genetics/165.4.2213) for local ancestry inference (local referring to a region of the genome). For a set of N phased haplotypes, kalis computes the posterior marginal probability of each haplotype copying every other haplotype by running N hidden Markov models in parallel. This yields an N x N matrix that summarizes the recent local ancestry at each variant of interest. The package provides functionality for specifying a recombination map, site-specific mutation rates, and differing prior copying probabilities for each recipient haplotype. Extensive use is made of low level threading and CPU vector instructions.
Installation
kalis will appear on CRAN in due course. For now, the current development version can be installed using the remotes package as follows:
install.packages("remotes")
remotes::install_github("louisaslett/kalis")Note that kalis uses various low-level optimisations meaning that you should ensure the compiler is targeting your local CPU architecture. The simplest way to do this is to pass configure variables setting the correct CFLAGS at install time:
remotes::install_github("louisaslett/kalis",
configure.vars = c(kalis = "PKG_CFLAGS='-march=native -mtune=native -O3'"))If you do not have the correct flags set, kalis will fall back to an implementation which does not use the special vector instruction set architecture of your CPU and will provide a warning when you load the package in your R session.
Citation
If you make use of this software, please cite:
Aslett L. J. M., Christ R. R. (2024). “kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R.” BMC Bioinformatics, 25(1), 1-18. doi:10.1186/s12859-024-05688-8.
bibTeX:
@Article{,
title = {kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R},
author = {Aslett, L. J. M. and Christ, R. R.},
journal = {BMC Bioinformatics},
year = {2024},
volume = {25},
number = {1},
pages = {1--18},
doi = {10.1186/s12859-024-05688-8}
}