diff --git a/R/listSims.R b/R/listSims.R
index 2fc06d168bdd2552b2e0dbead7e40e633be7187e..4dadb145a5c1b40d69849adaf4544322f50e9abb 100644
--- a/R/listSims.R
+++ b/R/listSims.R
@@ -40,7 +40,7 @@ listSims <- function(print = TRUE) {
                    "Gamma distributed means and negative binomial counts. Cells
                    are given a size factor and differential expression can be
                    simulated with fixed fold changes."),
-                 c("Lun 2", "lun2", "10.1101/073973",
+                 c("Lun 2", "lun2", "10.1093/biostatistics/kxw055",
                    "MarioniLab/PlateEffects2016",
                    "Negative binomial counts where the means and dispersions
                    have been sampled from a real dataset. The core feature of
diff --git a/R/lun2-simulate.R b/R/lun2-simulate.R
index 35329de75e4dd9f9fff1fb643c349e475c8997b4..c839eed6e1b6f53ff83d6e828d87c67ed5f0b822 100644
--- a/R/lun2-simulate.R
+++ b/R/lun2-simulate.R
@@ -30,9 +30,9 @@
 #'
 #' @references
 #' Lun ATL, Marioni JC. Overcoming confounding plate effects in differential
-#' expression analyses of single-cell RNA-seq data. bioRxiv (2016).
+#' expression analyses of single-cell RNA-seq data. Biostatistics (2017).
 #'
-#' Paper: \url{dx.doi.org/10.1101/073973}
+#' Paper: \url{dx.doi.org/10.1093/biostatistics/kxw055}
 #'
 #' Code: \url{https://github.com/MarioniLab/PlateEffects2016}
 #'
diff --git a/man/lun2Simulate.Rd b/man/lun2Simulate.Rd
index e091ce9d8c739de001334b4f9b8a0f6b2d3975c4..4b6c66dfc317aed96b34d707b7b593f93a7e317c 100644
--- a/man/lun2Simulate.Rd
+++ b/man/lun2Simulate.Rd
@@ -45,9 +45,9 @@ sim <- lun2Simulate()
 }
 \references{
 Lun ATL, Marioni JC. Overcoming confounding plate effects in differential
-expression analyses of single-cell RNA-seq data. bioRxiv (2016).
+expression analyses of single-cell RNA-seq data. Biostatistics (2017).
 
-Paper: \url{dx.doi.org/10.1101/073973}
+Paper: \url{dx.doi.org/10.1093/biostatistics/kxw055}
 
 Code: \url{https://github.com/MarioniLab/PlateEffects2016}
 }