diff --git a/src/bayesvalidrox/bayes_inference/bayes_inference.py b/src/bayesvalidrox/bayes_inference/bayes_inference.py
index 2e44a737c20cce12b167887d046ae3fab5e8e1ae..061283d3cc2f51ffe6d3a56067d25cc0b6b2e4e0 100644
--- a/src/bayesvalidrox/bayes_inference/bayes_inference.py
+++ b/src/bayesvalidrox/bayes_inference/bayes_inference.py
@@ -744,7 +744,7 @@ class BayesInference:
 
                 # Bayes factor Eq. 11
                 ln_BF[i] = -(n_samples**2/(2*n_samples+2))
-                ln_BF[i] *= np.dot(np.dot(y_bar, np.linalg.inv(covMatrix)),
+                ln_BF[i] *= np.dot(np.dot(y_bar, np.linalg.pinv(covMatrix)),
                                    np.transpose(y_bar))
                 ln_BF[i] += np.log(n_samples+1)/2