diff --git a/R/bayesfactor_parameters.R b/R/bayesfactor_parameters.R
index 34139f3d6..1251c83d5 100644
--- a/R/bayesfactor_parameters.R
+++ b/R/bayesfactor_parameters.R
@@ -132,8 +132,8 @@
 #'
 #'   # emmGrid objects
 #'   # ---------------
-#'   # group_diff <- pairs(emmeans(stan_model, ~group))
-#'   # bayesfactor_parameters(group_diff, prior = stan_model, verbose = FALSE)
+#'   group_diff <- pairs(emmeans(stan_model, ~group, data = sleep))
+#'   bayesfactor_parameters(group_diff, prior = stan_model, verbose = FALSE)
 #'
 #'   # Or
 #'   # group_diff_prior <- pairs(emmeans(unupdate(stan_model), ~group))
diff --git a/R/bayesfactor_restricted.R b/R/bayesfactor_restricted.R
index db9b3e246..77ae4fc7e 100644
--- a/R/bayesfactor_restricted.R
+++ b/R/bayesfactor_restricted.R
@@ -94,7 +94,7 @@
 #' contrasts(disgust$condition) <- contr.equalprior_pairs # see vignette
 #' fit_model <- rstanarm::stan_glm(score ~ condition, data = disgust, family = gaussian())
 #'
-#' em_condition <- emmeans::emmeans(fit_model, ~condition)
+#' em_condition <- emmeans::emmeans(fit_model, ~condition, data = disgust)
 #' hyps <- c("lemon < control & control < sulfur")
 #'
 #' bayesfactor_restricted(em_condition, prior = fit_model, hypothesis = hyps)
diff --git a/man/bayesfactor_parameters.Rd b/man/bayesfactor_parameters.Rd
index 76cf4cd3b..b5726ecbf 100644
--- a/man/bayesfactor_parameters.Rd
+++ b/man/bayesfactor_parameters.Rd
@@ -278,8 +278,8 @@ if (require("rstanarm") && require("emmeans") && require("logspline")) {
 
   # emmGrid objects
   # ---------------
-  # group_diff <- pairs(emmeans(stan_model, ~group))
-  # bayesfactor_parameters(group_diff, prior = stan_model, verbose = FALSE)
+  group_diff <- pairs(emmeans(stan_model, ~group, data = sleep))
+  bayesfactor_parameters(group_diff, prior = stan_model, verbose = FALSE)
 
   # Or
   # group_diff_prior <- pairs(emmeans(unupdate(stan_model), ~group))
diff --git a/man/bayesfactor_restricted.Rd b/man/bayesfactor_restricted.Rd
index 3e83abd4f..1fe715e68 100644
--- a/man/bayesfactor_restricted.Rd
+++ b/man/bayesfactor_restricted.Rd
@@ -208,7 +208,7 @@ data("disgust")
 contrasts(disgust$condition) <- contr.equalprior_pairs # see vignette
 fit_model <- rstanarm::stan_glm(score ~ condition, data = disgust, family = gaussian())
 
-em_condition <- emmeans::emmeans(fit_model, ~condition)
+em_condition <- emmeans::emmeans(fit_model, ~condition, data = disgust)
 hyps <- c("lemon < control & control < sulfur")
 
 bayesfactor_restricted(em_condition, prior = fit_model, hypothesis = hyps)
diff --git a/vignettes/bayes_factors.Rmd b/vignettes/bayes_factors.Rmd
index b2e791187..479475be2 100644
--- a/vignettes/bayes_factors.Rmd
+++ b/vignettes/bayes_factors.Rmd
@@ -969,7 +969,7 @@ and the difference between them:
 set.seed(1)
 ```
 
-```{r}
+```{r eval=FALSE}
 library(emmeans)
 
 (group_diff <- emmeans(model, pairwise ~ group, data = sleep))
@@ -977,6 +977,14 @@ library(emmeans)
 # pass the original model via prior
 bayesfactor_parameters(group_diff, prior = model)
 ```
+```{r echo=FALSE, message=FALSE, warning=FALSE}
+library(emmeans)
+
+suppressWarnings((group_diff <- emmeans(model, pairwise ~ group, data = sleep)))
+
+# pass the original model via prior
+suppressWarnings(bayesfactor_parameters(group_diff, prior = model))
+```
 
 That is strong evidence for the mean of group 1 being 0, and for group 2 for not
 being 0, but hardly any evidence for the difference between them being not 0.