diff --git a/R/mc_main_function.R b/R/mc_main_function.R
index a3e4232dc7a5723d629f8ea8739cc8e4e577b3fb..c884b3afda0a31ec8436224032b6cc01a32c0c3d 100644
--- a/R/mc_main_function.R
+++ b/R/mc_main_function.R
@@ -114,6 +114,7 @@ mcglm <- function(linear_pred, matrix_pred, link, variance, covariance, offset,
     model_fit$offset <- offset
     model_fit$power_fixed
     model_fit$sparse <- sparse
+    model_fit$data <- data
     class(model_fit) <- "mcglm"
   }
   return(model_fit)