时间戳

首页 > 解决方案 > 将逆伽马分布拟合到 R 中的数据

r - 将逆伽马分布拟合到 R 中的数据

问题描述

假设我有一个可以遵循逆伽马分布的样本（请参阅Empirical PDF）。

我想用类似的东西估计形状参数 alpha 和比例参数 beta fitdistr。可能吗？

我已经尝试过这个解决方案（遵循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找到。

标签： rdistributioncurve-fittinggamma-distribution

解决方案

我有答案。这很简单

fit = MASS::fitdistr(1/x1, "gamma")

推荐阅读