R中Levy分布的参数

Question

有人能帮我弄清楚如何使用R得到征费分布的参数估计吗？与正态分布和学生T分布不同的是，它们分别具有来自fBasics包的nFit(x)函数和来自fGarch包的stdFit(x)函数，而Levy分布似乎没有相同的能力。我正在尝试将我的数据集拟合到levy分布，所以我需要R中的一个函数，它可以给我对增量和伽马参数的良好估计。

我能够使用fBasics和fGarch包来计算正态分布的参数(µ和sigma)和对称学生T分布的参数(µ、sd和nu)。现在，我希望为Levy分布计算以下参数: delta和gamma。增量是位置参数，而gamma是比例参数。

通过执行以下操作，我尝试使用MASS中的fitdistr函数来拟合它：

> fitdistr(x,dlevy,list(m = 0, s = 1))
Error in densfun(x, parm[1], parm[2], ...) : some y <= m

然而，现在我被卡住了。我一直收到上面的“错误”信息。假设Levy发行版没有内置函数，那么如何处理呢？

提前感谢！

Answer 1

看起来VGAM包中的levy函数就是您要查找的(click)。