diff --git a/DESCRIPTION b/DESCRIPTION
index a4e12df..ec34bcf 100644
--- a/DESCRIPTION
+++ b/DESCRIPTION
@@ -46,7 +46,7 @@ Imports:
     coda,
     purrr,
     readr,
-    limSolve,
+    MASS,
     utils
 Suggests:
     knitr,
diff --git a/R/mptinr.R b/R/mptinr.R
index 16a7538..562aad6 100644
--- a/R/mptinr.R
+++ b/R/mptinr.R
@@ -801,7 +801,7 @@ mpt_mptinr_complete <- function(dataset,
     variance_covariance <- tryCatch(
       solve(fit_mptinr_tmp$hessian[[1]])
       , error = function(x){
-        limSolve::Solve(fit_mptinr_tmp$hessian[[1]])
+        MASS::ginv(fit_mptinr_tmp$hessian[[1]])
       })
     par_se <- rep(NA_real_, length(rownames(fit_mptinr_tmp$parameters)))
     par_se <- sqrt(diag(variance_covariance))