diff --git a/R/stray.R b/R/stray.R
index 1fdec05..25868dd 100644
--- a/R/stray.R
+++ b/R/stray.R
@@ -8,7 +8,7 @@
 #' # Univariate data
 #' y <- c(6, rnorm(49))
 #' scores <- stray_scores(y)
-#' threshold <- stray::find_threshold(scores)
+#' threshold <- stray::find_threshold(scores, outtail = "max")
 #' which(scores > threshold)
 #' stray_anomalies(y)
 #' # Bivariate data
diff --git a/man/stray_scores.Rd b/man/stray_scores.Rd
index ea3b6b8..f0d45c0 100644
--- a/man/stray_scores.Rd
+++ b/man/stray_scores.Rd
@@ -27,7 +27,7 @@ Compute stray scores indicating how anomalous each observation is.
 # Univariate data
 y <- c(6, rnorm(49))
 scores <- stray_scores(y)
-threshold <- stray::find_threshold(scores)
+threshold <- stray::find_threshold(scores, outtail = "max")
 which(scores > threshold)
 stray_anomalies(y)
 # Bivariate data