diff --git a/R/CacheHaplotypes.R b/R/CacheHaplotypes.R
index dcc6b3c..fe5e967 100644
--- a/R/CacheHaplotypes.R
+++ b/R/CacheHaplotypes.R
@@ -48,6 +48,9 @@ assign("L", NA, envir = pkgVars) # must be integer
 #'
 #'
 #'
+#' @references
+#' Aslett, L.J.M. and Christ, R.R. (2024) "kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R", *BMC Bioinformatics*, **25**(1). Available at: \doi{10.1186/s12859-024-05688-8}.
+#'
 #' @param haps can be the name of a file from which the haplotypes are to be read, or can be an R matrix containing only 0/1s.
 #'   See Details section for supported file types.
 #' @param loci.idx an optional vector of indices specifying the variants to load into the cache, indexed from 1.
@@ -219,6 +222,9 @@ CacheHaplotypes.err <- function(err) {
 #' To achieve higher performance, kalis internally represents haplotypes in an efficient raw binary format in memory which cannot be directly viewed or manipulated in R.
 #' This function enables you to copy whole or partial views of haplotypes/variants out of this low-level format and into a standard R matrix of 0's and 1's.
 #'
+#' @references
+#' Aslett, L.J.M. and Christ, R.R. (2024) "kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R", *BMC Bioinformatics*, **25**(1). Available at: \doi{10.1186/s12859-024-05688-8}.
+#'
 #' @param loci.idx which variants to retrieve from the cache, specified as a (vector) index.
 #'   This enables specifying variants by offset in the order they were loaded into the cache (from 1 to the number of variants).
 #' @param hap.idx which haplotypes to retrieve from the cache, specified as a (vector) index.
@@ -295,6 +301,9 @@ QueryCache <- function(loci.idx = NULL, hap.idx = NULL) {
 #' In particular, this cache sits outside R's memory management and will never be garbage collected (unless R is quit or the package is unloaded).
 #' Therefore, this function is provided to enable freeing the memory used by this cache.
 #'
+#' @references
+#' Aslett, L.J.M. and Christ, R.R. (2024) "kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R", *BMC Bioinformatics*, **25**(1). Available at: \doi{10.1186/s12859-024-05688-8}.
+#'
 #' @return Nothing is returned.
 #'
 #' @seealso [CacheHaplotypes()] to create a haplotype cache;
diff --git a/R/CacheSummary.R b/R/CacheSummary.R
index 8c5b77a..c923096 100644
--- a/R/CacheSummary.R
+++ b/R/CacheSummary.R
@@ -1,5 +1,8 @@
 #' Retrieve information about the haplotype cache
 #'
+#' @references
+#' Aslett, L.J.M. and Christ, R.R. (2024) "kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R", *BMC Bioinformatics*, **25**(1). Available at: \doi{10.1186/s12859-024-05688-8}.
+#'
 #' @return
 #'   `CacheSummary()` prints information about the current state of the kalis cache.
 #'     Also invisibly returns a vector giving the dimensions of the cached haplotype data (num variants, num haplotypes), or `NULL` if the cache is empty.
diff --git a/R/IndividualSequenceIO_H5.R b/R/IndividualSequenceIO_H5.R
index ac3f74c..d6b593c 100644
--- a/R/IndividualSequenceIO_H5.R
+++ b/R/IndividualSequenceIO_H5.R
@@ -11,6 +11,9 @@
 #'
 #' Note that if `hdf5.file` exists but does not contain a dataset named `haps`, then `WriteHaplotypes` will simply create a `haps` dataset within the existing file.
 #'
+#' @references
+#' Aslett, L.J.M. and Christ, R.R. (2024) "kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R", *BMC Bioinformatics*, **25**(1). Available at: \doi{10.1186/s12859-024-05688-8}.
+#'
 #' @param hdf5.file the name of the file which the haplotypes are to be written to.
 #' @param haps a vector or a matrix where each column is a haplotype to be stored in the file `hdf5.file`.
 #' @param hap.ids a character vector naming haplotypes when writing, or which haplotypes are to be read.
diff --git a/R/Parameters.R b/R/Parameters.R
index 846ef91..187194e 100644
--- a/R/Parameters.R
+++ b/R/Parameters.R
@@ -6,11 +6,16 @@
 #'
 #' **NOTE:** the corresponding haplotype data *must* have already been inserted into the kalis cache by a call to [CacheHaplotypes()], since this function performs checks to confirm the dimensionality matches.
 #'
-#' TODO: add kalis paper cross ref.
 #' See page 3 in Supplemental Information for the original ChromoPainter paper (Lawson et al., 2012) for motivation behind our parameterisation, which is as follows:
 #'
 #' \deqn{\rho = 1 - \exp(-s \times cM^\gamma)}{\rho = 1 - exp(-s * cM^\gamma)}
 #'
+#' For a complete description, see the main kalis paper, Aslett and Christ (2024).
+#'
+#' @references
+#' Aslett, L.J.M. and Christ, R.R. (2024) "kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R", *BMC Bioinformatics*, **25**(1). Available at: \doi{10.1186/s12859-024-05688-8}.
+#'
+#' Lawson, D.J., Hellenthal, G., Myers, S. and Falush, D. (2012). "Inference of population structure using dense haplotype data", *PLoS genetics*, **8**(1). Available at: \doi{10.1371/journal.pgen.1002453}.
 #'
 #' @param cM a vector specifying the recombination distance between variants in centimorgans.
 #'   Note element i of this vector should be the distance between variants `i` and `i+1` (not `i` and `i-1`), and thus length one less than the number of variants.
@@ -25,10 +30,6 @@
 #'
 #' @seealso [Parameters()] to use the resulting recombination probabilities to construct a `kalisParameters` object.
 #'
-#' @references
-#'   Lawson, D. J., Hellenthal, G., Myers, S., & Falush, D. (2012). Inference of
-#'   population structure using dense haplotype data. *PLoS genetics*, **8**(1).
-#'
 #' @examples
 #' # Load the mini example data and recombination map from the package built-in #' # dataset
 #' data("SmallHaps")
@@ -119,6 +120,13 @@ CalcRho <- function(cM = 0, s = 1, gamma = 1, floor = TRUE) {
 #'
 #' Note that there is a computational cost associated with non-uniform copying probabilities, so it is recommended to leave the default of uniform probabilities when appropriate (**Note:** *do not* specify a uniform matrix when uniform probabilities are intended, since this would end up incurring the computational cost of non-uniform probabilities).
 #'
+#' @references
+#' Aslett, L.J.M. and Christ, R.R. (2024) "kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R", *BMC Bioinformatics*, **25**(1). Available at: \doi{10.1186/s12859-024-05688-8}.
+#'
+#' Lawson, D.J., Hellenthal, G., Myers, S.R. and Falush, D. (2012). "Inference of population structure using dense haplotype data", *PLoS Genetics*, **8**(1). Available at: \doi{10.1371/journal.pgen.1002453}.
+#'
+#' Speidel, L., Forest, M., Shi, S. and Myers, S.R. (2019). "A method for genome-wide genealogy estimation for thousands of samples", *Nature Genetics*, **51**, p. 1321-1329. Available at: \doi{10.1038/s41588-019-0484-x}.
+#'
 #' @param rho recombination probability vector (must be \eqn{L-1} long).
 #'   See [CalcRho()] for assistance constructing this from a recombination
 #'   map/distances.
@@ -140,13 +148,6 @@ CalcRho <- function(cM = 0, s = 1, gamma = 1, floor = TRUE) {
 #' @seealso [MakeForwardTable()] and [MakeBackwardTable()] which construct table objects which internally reference a parameters environment;
 #'   [Forward()] and [Backward()] which propagate those tables according to the Li and Stephens model.
 #'
-#' @references
-#'   Lawson, D. J., Hellenthal, G., Myers, S., & Falush, D. (2012). Inference of
-#'   population structure using dense haplotype data. *PLoS genetics*, **8**(1).
-#'
-#'   Speidel, L., Forest, M., Shi, S., & Myers, S. (2019). A method for
-#'   genome-wide genealogy estimation for thousands of samples. *Nature Genetics*, **51**(1321–1329).
-#'
 #' @examples
 #' # Load the mini example data and recombination map from the package built-in #' # dataset
 #' data("SmallHaps")
diff --git a/R/Probs.R b/R/Probs.R
index a97303e..77f6024 100644
--- a/R/Probs.R
+++ b/R/Probs.R
@@ -43,6 +43,9 @@
 #' Typically, that is simply \eqn{N \times N}{N x N} for \eqn{N} haplotypes.
 #' However, if kalis is being run in a distributed manner, `M` will be a \eqn{N \times R}{N x R} matrix where \eqn{R} is the number of recipient haplotypes on the current machine.
 #'
+#' @references
+#' Aslett, L.J.M. and Christ, R.R. (2024) "kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R", *BMC Bioinformatics*, **25**(1). Available at: \doi{10.1186/s12859-024-05688-8}.
+#'
 #' @param fwd a forward table as returned by [MakeForwardTable()] and propagated to a target variant by [Forward()].
 #'   Must be at the same variant as `bck` (unless `bck` is in "beta-theta space" in which case if must be downstream ... see [Backward()] for details).
 #' @param bck a backward table as returned by [MakeBackwardTable()] and propagated to a target variant by [Backward()].
@@ -164,6 +167,11 @@ PostProbs <- function(fwd, bck, unif.on.underflow = FALSE, M = NULL, beta.theta.
 #' Typically, that is simply \eqn{N \times N}{N x N} for \eqn{N} haplotypes.
 #' However, if kalis is being run in a distributed manner, `M` will be a \eqn{N \times R}{N x R} matrix where \eqn{R} is the number of recipient haplotypes on the current machine.
 #'
+#' @references
+#' Aslett, L.J.M. and Christ, R.R. (2024) "kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R", *BMC Bioinformatics*, **25**(1). Available at: \doi{10.1186/s12859-024-05688-8}.
+#'
+#' Speidel, L., Forest, M., Shi, S. and Myers, S.R. (2019). "A method for genome-wide genealogy estimation for thousands of samples", *Nature Genetics*, **51**, p. 1321-1329. Available at: \doi{10.1038/s41588-019-0484-x}.
+#'
 #' @param fwd a forward table as returned by [MakeForwardTable()] and propagated to a target variant by [Forward()].
 #'   Must be at the same variant as `bck` (unless `bck` is in "beta-theta space" in which case if must be downstream ... see [Backward()] for details).
 #' @param bck a backward table as returned by [MakeBackwardTable()] and propagated to a target variant by [Backward()].
@@ -184,9 +192,6 @@ PostProbs <- function(fwd, bck, unif.on.underflow = FALSE, M = NULL, beta.theta.
 #'
 #'   If you wish to plot this matrix or perform clustering, you may want to symmetrize the matrix first.
 #'
-#' @references
-#'   Speidel, L., Forest, M., Shi, S., & Myers, S. (2019). A method for genome-wide genealogy estimation for thousands of samples. *Nature Genetics*, **51**(1321–1329).
-#'
 #' @seealso
 #'   [PostProbs()] to calculate the posterior marginal probabilities \eqn{p_{ji}}{p_(j,i)};
 #'   [Forward()] to propagate a Forward table to a new variant;
@@ -319,6 +324,9 @@ input_checks_for_probs_and_dist_mat <-  function(fwd,bck,beta.theta.opts = NULL)
 #'
 #' Clusters the given distance matrix and generates a heatmap to display it.
 #'
+#' @references
+#' Aslett, L.J.M. and Christ, R.R. (2024) "kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R", *BMC Bioinformatics*, **25**(1). Available at: \doi{10.1186/s12859-024-05688-8}.
+#'
 #' @param x
 #'        a distance matrix, such as returned by [DistMat()].
 #' @param cluster.method
diff --git a/R/SmallHaps-data.R b/R/SmallHaps-data.R
index 26dfcfa..fe57838 100644
--- a/R/SmallHaps-data.R
+++ b/R/SmallHaps-data.R
@@ -14,7 +14,9 @@
 #' @keywords datasets
 #'
 #' @references
-#'   Kelleher, J., Etheridge, A. M., & McVean, G. (2016). Efficient coalescent simulation and genealogical analysis for large sample sizes. *PLoS computational biology*, **12**(5).
+#' Aslett, L.J.M. and Christ, R.R. (2024) "kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R", *BMC Bioinformatics*, **25**(1). Available at: \doi{10.1186/s12859-024-05688-8}.
+#'
+#' Kelleher, J., Etheridge, A.M. and McVean, G. (2016) "Efficient coalescent simulation and genealogical analysis for large sample sizes", *PLoS Computational Biology*, **12**(5). Available at: \doi{10.1371/journal.pcbi.1004842}.
 #'
 #' @examples
 #' data("SmallHaps")
diff --git a/R/TableMaker.R b/R/TableMaker.R
index 8792c21..b5fc56c 100644
--- a/R/TableMaker.R
+++ b/R/TableMaker.R
@@ -10,6 +10,9 @@
 #'
 #' Since each column corresponds to an independent Li and Stephens hidden Markov model (ie for each recipient), it is possible to create a partial forward table object which corresponds to a subset of recipients using the `from_recipient` and `to_recipient` arguments.
 #'
+#' @references
+#' Aslett, L.J.M. and Christ, R.R. (2024) "kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R", *BMC Bioinformatics*, **25**(1). Available at: \doi{10.1186/s12859-024-05688-8}.
+#'
 #' @param pars a `kalisParameters` object specifying the genetics parameters to be associated with this forward table.
 #'   These parameters can be set up by using the [Parameters()] function.
 #' @param from_recipient first recipient haplotype included if creating a partial forward table.
@@ -130,6 +133,9 @@ print.kalisForwardTable <- function(x, ...) {
 #'
 #' Since each column corresponds to an independent Li and Stephens hidden Markov model (ie for each recipient), it is possible to create a partial backward table object which corresponds to a subset of recipients using the `from_recipient` and `to_recipient` arguments.
 #'
+#' @references
+#' Aslett, L.J.M. and Christ, R.R. (2024) "kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R", *BMC Bioinformatics*, **25**(1). Available at: \doi{10.1186/s12859-024-05688-8}.
+#'
 #' @param pars a `kalisParameters` object specifying the genetics parameters to be associated with this backward table.
 #'   These parameters can be set up by using the [Parameters()] function.
 #' @param from_recipient first recipient haplotype included if creating a partial backward table.
@@ -249,6 +255,9 @@ print.kalisBackwardTable <- function(x, ...) {
 #'
 #' This function is therefore designed to enable explicit copying of tables.
 #'
+#' @references
+#' Aslett, L.J.M. and Christ, R.R. (2024) "kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R", *BMC Bioinformatics*, **25**(1). Available at: \doi{10.1186/s12859-024-05688-8}.
+#'
 #' @param to a `kalisForwardTable` or `kalisBackwardTable` object which is to be copied into.
 #' @param from a `kalisForwardTable` or `kalisBackwardTable` object which is to be copied from.
 #'
@@ -340,6 +349,9 @@ CopyTable <- function(to, from) {
 #' It is *much* faster to reset a forward/backward table rather than remove and make a new one.
 #' This function marks a table as reset so that it will be propagated as if freshly allocated.
 #'
+#' @references
+#' Aslett, L.J.M. and Christ, R.R. (2024) "kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R", *BMC Bioinformatics*, **25**(1). Available at: \doi{10.1186/s12859-024-05688-8}.
+#'
 #' @param tbl a `kalisForwardTable` or `kalisBackwardTable` object
 #'   which is to be reset.
 #'
diff --git a/man/CacheHaplotypes.Rd b/man/CacheHaplotypes.Rd
index b22bed1..98f5f69 100644
--- a/man/CacheHaplotypes.Rd
+++ b/man/CacheHaplotypes.Rd
@@ -143,6 +143,9 @@ CacheHaplotypes(haps)
 data("SmallHaps")
 CacheHaplotypes(SmallHaps)
 
+}
+\references{
+Aslett, L.J.M. and Christ, R.R. (2024) "kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R", \emph{BMC Bioinformatics}, \strong{25}(1). Available at: \doi{10.1186/s12859-024-05688-8}.
 }
 \seealso{
 \code{\link[=CacheSummary]{CacheSummary()}} for a list detailing the current cache status;
diff --git a/man/CacheSummary.Rd b/man/CacheSummary.Rd
index d53cddd..da7a847 100644
--- a/man/CacheSummary.Rd
+++ b/man/CacheSummary.Rd
@@ -40,3 +40,6 @@ N()
 L()
 
 }
+\references{
+Aslett, L.J.M. and Christ, R.R. (2024) "kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R", \emph{BMC Bioinformatics}, \strong{25}(1). Available at: \doi{10.1186/s12859-024-05688-8}.
+}
diff --git a/man/CalcRho.Rd b/man/CalcRho.Rd
index 3f87fbd..4714f4b 100644
--- a/man/CalcRho.Rd
+++ b/man/CalcRho.Rd
@@ -30,10 +30,11 @@ This is a utility function to calculate the recombination probabilities rho (the
 
 \strong{NOTE:} the corresponding haplotype data \emph{must} have already been inserted into the kalis cache by a call to \code{\link[=CacheHaplotypes]{CacheHaplotypes()}}, since this function performs checks to confirm the dimensionality matches.
 
-TODO: add kalis paper cross ref.
 See page 3 in Supplemental Information for the original ChromoPainter paper (Lawson et al., 2012) for motivation behind our parameterisation, which is as follows:
 
 \deqn{\rho = 1 - \exp(-s \times cM^\gamma)}{\rho = 1 - exp(-s * cM^\gamma)}
+
+For a complete description, see the main kalis paper, Aslett and Christ (2024).
 }
 \examples{
 # Load the mini example data and recombination map from the package built-in #' # dataset
@@ -49,8 +50,9 @@ rho
 
 }
 \references{
-Lawson, D. J., Hellenthal, G., Myers, S., & Falush, D. (2012). Inference of
-population structure using dense haplotype data. \emph{PLoS genetics}, \strong{8}(1).
+Aslett, L.J.M. and Christ, R.R. (2024) "kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R", \emph{BMC Bioinformatics}, \strong{25}(1). Available at: \doi{10.1186/s12859-024-05688-8}.
+
+Lawson, D.J., Hellenthal, G., Myers, S. and Falush, D. (2012). "Inference of population structure using dense haplotype data", \emph{PLoS genetics}, \strong{8}(1). Available at: \doi{10.1371/journal.pgen.1002453}.
 }
 \seealso{
 \code{\link[=Parameters]{Parameters()}} to use the resulting recombination probabilities to construct a \code{kalisParameters} object.
diff --git a/man/ClearHaplotypeCache.Rd b/man/ClearHaplotypeCache.Rd
index 7afce6c..9e832b5 100644
--- a/man/ClearHaplotypeCache.Rd
+++ b/man/ClearHaplotypeCache.Rd
@@ -31,6 +31,9 @@ ClearHaplotypeCache()
 # Verify it is gone
 CacheSummary()
 
+}
+\references{
+Aslett, L.J.M. and Christ, R.R. (2024) "kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R", \emph{BMC Bioinformatics}, \strong{25}(1). Available at: \doi{10.1186/s12859-024-05688-8}.
 }
 \seealso{
 \code{\link[=CacheHaplotypes]{CacheHaplotypes()}} to create a haplotype cache;
diff --git a/man/CopyTable.Rd b/man/CopyTable.Rd
index 1020f2f..a8d044d 100644
--- a/man/CopyTable.Rd
+++ b/man/CopyTable.Rd
@@ -61,6 +61,9 @@ Forward(fwd, pars, 30)
 fwd
 fwd2
 
+}
+\references{
+Aslett, L.J.M. and Christ, R.R. (2024) "kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R", \emph{BMC Bioinformatics}, \strong{25}(1). Available at: \doi{10.1186/s12859-024-05688-8}.
 }
 \seealso{
 \code{\link[=MakeForwardTable]{MakeForwardTable()}}, \code{\link[=MakeBackwardTable]{MakeBackwardTable()}} to create tables which can be copied into.
diff --git a/man/DistMat.Rd b/man/DistMat.Rd
index ffe3e6c..cf38ad3 100644
--- a/man/DistMat.Rd
+++ b/man/DistMat.Rd
@@ -129,7 +129,9 @@ plot(d)
 
 }
 \references{
-Speidel, L., Forest, M., Shi, S., & Myers, S. (2019). A method for genome-wide genealogy estimation for thousands of samples. \emph{Nature Genetics}, \strong{51}(1321–1329).
+Aslett, L.J.M. and Christ, R.R. (2024) "kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R", \emph{BMC Bioinformatics}, \strong{25}(1). Available at: \doi{10.1186/s12859-024-05688-8}.
+
+Speidel, L., Forest, M., Shi, S. and Myers, S.R. (2019). "A method for genome-wide genealogy estimation for thousands of samples", \emph{Nature Genetics}, \strong{51}, p. 1321-1329. Available at: \doi{10.1038/s41588-019-0484-x}.
 }
 \seealso{
 \code{\link[=PostProbs]{PostProbs()}} to calculate the posterior marginal probabilities \eqn{p_{ji}}{p_(j,i)};
diff --git a/man/MakeBackwardTable.Rd b/man/MakeBackwardTable.Rd
index 5b78dcb..4718432 100644
--- a/man/MakeBackwardTable.Rd
+++ b/man/MakeBackwardTable.Rd
@@ -62,6 +62,9 @@ bck
 
 # Now ready to run the HMM backward recursions with Backward() ...
 
+}
+\references{
+Aslett, L.J.M. and Christ, R.R. (2024) "kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R", \emph{BMC Bioinformatics}, \strong{25}(1). Available at: \doi{10.1186/s12859-024-05688-8}.
 }
 \seealso{
 \code{\link[=Backward]{Backward()}} to propagate the newly created \code{kalisBackwardTable};
diff --git a/man/MakeForwardTable.Rd b/man/MakeForwardTable.Rd
index a93c354..9e60ee9 100644
--- a/man/MakeForwardTable.Rd
+++ b/man/MakeForwardTable.Rd
@@ -61,6 +61,9 @@ fwd
 
 # Now ready to run the HMM forward recursions with Forward() ...
 
+}
+\references{
+Aslett, L.J.M. and Christ, R.R. (2024) "kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R", \emph{BMC Bioinformatics}, \strong{25}(1). Available at: \doi{10.1186/s12859-024-05688-8}.
 }
 \seealso{
 \code{\link[=Forward]{Forward()}} to propagate the newly created \code{kalisForwardTable};
diff --git a/man/Parameters.Rd b/man/Parameters.Rd
index 0b14657..f34bdf7 100644
--- a/man/Parameters.Rd
+++ b/man/Parameters.Rd
@@ -98,11 +98,11 @@ fwd
 
 }
 \references{
-Lawson, D. J., Hellenthal, G., Myers, S., & Falush, D. (2012). Inference of
-population structure using dense haplotype data. \emph{PLoS genetics}, \strong{8}(1).
+Aslett, L.J.M. and Christ, R.R. (2024) "kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R", \emph{BMC Bioinformatics}, \strong{25}(1). Available at: \doi{10.1186/s12859-024-05688-8}.
 
-Speidel, L., Forest, M., Shi, S., & Myers, S. (2019). A method for
-genome-wide genealogy estimation for thousands of samples. \emph{Nature Genetics}, \strong{51}(1321–1329).
+Lawson, D.J., Hellenthal, G., Myers, S.R. and Falush, D. (2012). "Inference of population structure using dense haplotype data", \emph{PLoS Genetics}, \strong{8}(1). Available at: \doi{10.1371/journal.pgen.1002453}.
+
+Speidel, L., Forest, M., Shi, S. and Myers, S.R. (2019). "A method for genome-wide genealogy estimation for thousands of samples", \emph{Nature Genetics}, \strong{51}, p. 1321-1329. Available at: \doi{10.1038/s41588-019-0484-x}.
 }
 \seealso{
 \code{\link[=MakeForwardTable]{MakeForwardTable()}} and \code{\link[=MakeBackwardTable]{MakeBackwardTable()}} which construct table objects which internally reference a parameters environment;
diff --git a/man/PostProbs.Rd b/man/PostProbs.Rd
index 389927a..8b15103 100644
--- a/man/PostProbs.Rd
+++ b/man/PostProbs.Rd
@@ -104,6 +104,9 @@ Backward(bck, pars, 100)
 p <- PostProbs(fwd, bck)
 d <- DistMat(fwd, bck)
 
+}
+\references{
+Aslett, L.J.M. and Christ, R.R. (2024) "kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R", \emph{BMC Bioinformatics}, \strong{25}(1). Available at: \doi{10.1186/s12859-024-05688-8}.
 }
 \seealso{
 \code{\link[=DistMat]{DistMat()}} to generate calculate \eqn{d_{ji}}{d_(j,i)} distances directly;
diff --git a/man/QueryCache.Rd b/man/QueryCache.Rd
index dc5de1d..67621a1 100644
--- a/man/QueryCache.Rd
+++ b/man/QueryCache.Rd
@@ -36,6 +36,9 @@ CacheHaplotypes(haps)
 res <- QueryCache(50:150, 10:11)
 all(res == haps[50:150, 10:11])
 
+}
+\references{
+Aslett, L.J.M. and Christ, R.R. (2024) "kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R", \emph{BMC Bioinformatics}, \strong{25}(1). Available at: \doi{10.1186/s12859-024-05688-8}.
 }
 \seealso{
 \code{\link[=CacheHaplotypes]{CacheHaplotypes()}} to fill the memory cache with haplotype data.
diff --git a/man/ResetTable.Rd b/man/ResetTable.Rd
index 9797a48..024ef5c 100644
--- a/man/ResetTable.Rd
+++ b/man/ResetTable.Rd
@@ -42,6 +42,9 @@ ResetTable(fwd)
 Forward(fwd, pars, 5)
 fwd
 
+}
+\references{
+Aslett, L.J.M. and Christ, R.R. (2024) "kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R", \emph{BMC Bioinformatics}, \strong{25}(1). Available at: \doi{10.1186/s12859-024-05688-8}.
 }
 \seealso{
 \code{\link[=MakeForwardTable]{MakeForwardTable()}}, \code{\link[=MakeBackwardTable]{MakeBackwardTable()}} to create fresh tables.
diff --git a/man/SmallHaps.Rd b/man/SmallHaps.Rd
index 6ff1916..5b540b1 100644
--- a/man/SmallHaps.Rd
+++ b/man/SmallHaps.Rd
@@ -40,6 +40,8 @@ pars
 
 }
 \references{
-Kelleher, J., Etheridge, A. M., & McVean, G. (2016). Efficient coalescent simulation and genealogical analysis for large sample sizes. \emph{PLoS computational biology}, \strong{12}(5).
+Aslett, L.J.M. and Christ, R.R. (2024) "kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R", \emph{BMC Bioinformatics}, \strong{25}(1). Available at: \doi{10.1186/s12859-024-05688-8}.
+
+Kelleher, J., Etheridge, A.M. and McVean, G. (2016) "Efficient coalescent simulation and genealogical analysis for large sample sizes", \emph{PLoS Computational Biology}, \strong{12}(5). Available at: \doi{10.1371/journal.pcbi.1004842}.
 }
 \keyword{datasets}
diff --git a/man/WriteHaplotypes.Rd b/man/WriteHaplotypes.Rd
index 5d2d0fe..a27eed8 100644
--- a/man/WriteHaplotypes.Rd
+++ b/man/WriteHaplotypes.Rd
@@ -94,6 +94,9 @@ CacheHaplotypes("~/myhaps.h5")
 all(haps == QueryCache())
 }
 
+}
+\references{
+Aslett, L.J.M. and Christ, R.R. (2024) "kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R", \emph{BMC Bioinformatics}, \strong{25}(1). Available at: \doi{10.1186/s12859-024-05688-8}.
 }
 \seealso{
 \code{\link[=CacheHaplotypes]{CacheHaplotypes()}} to fill the kalis cache with haplotypes.
diff --git a/man/plot.kalisDistanceMatrix.Rd b/man/plot.kalisDistanceMatrix.Rd
index 5246eb4..9ca8ec3 100644
--- a/man/plot.kalisDistanceMatrix.Rd
+++ b/man/plot.kalisDistanceMatrix.Rd
@@ -24,3 +24,6 @@ Clusters the given distance matrix and generates a heatmap to display it.
 # TODO
 
 }
+\references{
+Aslett, L.J.M. and Christ, R.R. (2024) "kalis: a modern implementation of the Li & Stephens model for local ancestry inference in R", \emph{BMC Bioinformatics}, \strong{25}(1). Available at: \doi{10.1186/s12859-024-05688-8}.
+}