diff --git a/R/PLNfit-class.R b/R/PLNfit-class.R
index b2dc83f3..64786a92 100644
--- a/R/PLNfit-class.R
+++ b/R/PLNfit-class.R
@@ -205,6 +205,7 @@ PLNfit <- R6Class(
     },
 
     variance_jackknife = function(Y, X, O, w, config = config_default_nlopt) {
+      cat("Computing jackknife variance estimator", sep = "\n")
       jacks <- future.apply::future_lapply(seq_len(self$n), function(i) {
         data <- list(Y = Y[-i, , drop = FALSE],
                      X = X[-i, , drop = FALSE],