我想将概率标度上的值转换为 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 概率的存在?