【发布时间】：2018-12-04 11:56:03
【问题描述】：
假设我有一个可以遵循逆伽马分布的样本（请参阅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找到。
【问题讨论】：
-
分享您的样本会很有帮助。你试过不同的起始参数值吗？
-
感谢您的重播，朱利叶斯。我已经包含了我正在使用的功能。我也改了参数，但是又出现了一个错误。
-
我的意思是分享你的数据。由于您使用的功能与其他人相同，因此它与您的数据和/或参数有关。因此，为了帮助您，我们需要查看您的数据。
-
好的。我已经分享了数据。
-
你的数据是什么？看起来像一些单列递减值？这些只是要分箱的样本吗？
标签： r distribution curve-fitting gamma-distribution