使反伽马分布适合R

时间:2018-12-04 11:56:04

标签: r distribution curve-fitting gamma-distribution

假设我有一个样本,该样本可以遵循反伽马分布(请参见Empirical PDF)。

我想用fitdistr之类的方法估算形状参数alpha和比例参数beta。有可能吗?

我尝试了以下解决方案(在https://stats.stackexchange.com/questions/31934/maximum-likelihood-estimation-of-inverse-gamma-distribution-in-r-or-rpy之后):

f <- function(x, a, b){
  ((b^a)/gamma(a))*((1/x)^(a-1))*exp(-b/x) #PDF Inv. Gamma
}
fitdistr(x, f, list(a=.01, b=.01))

但这对我不起作用。它说:非有限的有限差分值[2]。

可以在https://www.dropbox.com/s/j4n09w1sszcv0ud/data.txt?dl=0上找到数据。

1 个答案:

答案 0 :(得分:0)

我有答案。就像

一样简单
fit = MASS::fitdistr(1/x1, "gamma")