diff --git a/R/bcf.R b/R/bcf.R
index 2d02cb6b..2d1564fa 100644
--- a/R/bcf.R
+++ b/R/bcf.R
@@ -2707,6 +2707,7 @@ predict.bcfmodel <- function(
     }
 
     # Add propensities to covariate set if necessary
+    X_combined <- X
     if (object$model_params$propensity_covariate != "none") {
         X_combined <- cbind(X, propensity)
     }