我使用 Monte-Carlo 方法使用负二项分布估计来自 GLM 的新观测值的预测区间。我将这种方法用于线性模型并得到了预测区间的可靠估计,但我对统计数据没有足够的信心知道这种方法是否也可以应用于非高斯分布。你能告诉我使用它是否有意义吗?
方法:
- 使用 Wald 方法(在链接上)计算期望值的 se。
- 使用这个 se 从高斯分布中随机简单 1000 个值,期望值为并且作为.
- 对于每个随机值,从负二项分布中随机抽取 1000 个值,随机值的指数为 mu。
- 预测区间计算为给定预期值获得的 1e6 个随机值的 2.5 和 97.5 分位数。
这种方法有意义吗?
问候,马克西姆