我有一个支持 [4,6] 的 Gamma(shape=50, scale=0.1)。我能够通过将完整的伽马分布除以 F(6) - F(4) 来找到它的分布。
p1 = seq(1,10,length=100)
d1 = dgamma(p1, shape=50, scale=0.1)
p2 = seq(4,6,length=100)
d2.full = dgamma(p2, shape=50, scale=0.1)
d2 = d2.full / (pgamma(6, shape=50, scale=0.1) - pgamma(4, shape=50, scale=0.1))
我如何找到这个截断分布的中心 95 可信区间(即 d2)?
编辑:请注意,我的截断伽玛与标准伽玛没有相同的 pdf。原因是截断的 gamma 必须重新归一化,以便它在支持上积分为 1 [4,6]。这就是为什么 d2 = d2.full / (F(6) - F(4))