增加伽马分布密度的边界

时间:2015-10-22 02:21:18

标签: r

dgamma除了这些之外还有其他元素吗?

dgamma(x, shape, rate = 1, scale = 1/rate, log = FALSE)

例如:下限/上限等

2 个答案:

答案 0 :(得分:1)

我不知道为dgamma设置下限和上限的任何选项,但一般来说,你可以编写一个函数来扩展你的发行版,如下所示:

scale <- function(dist){(dist-min(dist))/(max(dist)-min(dist))*(upper-lower)+lower}

如果dist是你的分布(dist&lt; -dgamma(...)),那么上面和下面是分布的理想上下边界。请注意,尽管分布的一般形状是守恒的,但统计特性可能会发生变化。

答案 1 :(得分:1)

如果你想要一个截断的伽玛分布,那么只需将dgamma(在指定的范围内)除以截断限制末端的pgamma值之间的差值。

 dgamma.trunc <- function(x, ..., lower=0, upper=Inf){if( x < lower|| x > upper){0} else {
                   dgamma(x, ...)/( pgamma(upper, ...)-pgamma(lower, ...) )}   }