Bug with Pareto distribution and eval quantile

[vincenzocoia]

library(tidyverse)
library(distplyr)
dst_pareto <- local({
  ppareto <- function(q, xi, lower.tail = TRUE) {
    res <- pmin(1, q ^ (-1 / xi))
    if (lower.tail) res <- 1 - res
    res
  }
  dpareto <- function(q, xi) sapply(q, function(q_) if (q_ >= 1) q_ ^ (-1 / xi - 1) / xi else 0)
  qpareto <- function(p, xi) (1 - p) ^ (-xi)
  function(xi) dst_parametric("pareto", xi = xi, .variable = "continuous")
})
eval_quantile(maximise(dst_pareto(1), dst_pareto(2) + 5), at = 0.5)

Error in while (cdf_gt(cdf_left, p_min)) { :
missing value where TRUE/FALSE needed