从概率尺度转换为 R 中的 logit 尺度？

Question

我想将概率标度上的值转换为 logit 标度。例如，如果我创建一些概率，我可以像这样转换它们：

# Create some probabilitws between 0 and 1
probs <- round(runif(5, 0.0, 1.0), 2)

# convert to logit
prob_logit <- function(x) {
  out <- log(x / (1 - x))
  return(out)
}

prob_logit(probs)

问题是，如果概率是 0 或 1，那么我的函数将返回 -Inf 或 Inf，因为 logit 将接近 ± 无穷大。所以，我的问题是，如何从概率尺度转换为 logit 尺度，同时还要考虑 0 或 1 概率的存在？