R语言ARIMA带0-100%区间限制的预测方法及对数变换问题求助

问题2：ARIMA模型中实现带0-100%区间限制的预测 #

思路和ETS模型一致：先对百分比数据做logit变换（专门针对0-1比例的对数变换，百分比转成比例后用更方便），拟合ARIMA模型后再逆变换回百分比，这样置信区间会自动被限制在0-100%之间。

完整实现代码：

# 1. 将百分比转成0-1的比例
room$prop <- room$rate / 100

# 2. 预处理比例数据：避免0或1（logit变换不允许边界值）
room$prop <- ifelse(room$prop <= 0, 0.001, room$prop)
room$prop <- ifelse(room$prop >= 1, 0.999, room$prop)

# 3. logit变换（等价于log(p/(1-p))，用R内置函数qlogis更简洁）
transformed_prop <- qlogis(room$prop)

# 4. 拟合ARIMA模型（用auto.arima自动选择最优模型参数）
fit_arima <- auto.arima(transformed_prop)
fc_arima <- forecast(fit_arima, h=50)

# 5. 逆变换：转回百分比尺度（用plogis做logit逆变换）
fc_arima$mean <- 100 * plogis(fc_arima$mean)
fc_arima$lower <- 100 * plogis(fc_arima$lower)
fc_arima$upper <- 100 * plogis(fc_arima$upper)

# 可视化预测结果
plot(fc_arima)

额外提示：

• 处理边界值时，尽量用接近边界的小值（比如0.001、0.999），避免引入过大的偏差；
• 如果你的数据中边界值很少，也可以考虑直接删除这些异常行，但不建议（会丢失数据）；
• qlogis()和plogis()是R中专门用于logit变换的函数，比手动写log(p/(1-p))更简洁且不易出错。

内容的提问来源于stack exchange，提问作者WWYY