“fitdistrplus”包有问题，t-distribution

Question

我正在尝试将 t 分布拟合到我的数据中，但无法这样做。我的第一次尝试是

fitdistr(myData, "t")

有 41 个警告，都说产生了 NaN。我不知道如何，似乎涉及对数。所以我稍微调整了我的数据，使所有数据都> 0，但我仍然有同样的问题（虽然警告减少了 9 个......）。sstdFit() 也有同样的问题，会产生 NaN。

因此，我尝试使用我在 stackoverflow 和 CrossValidated 上看到的 fitdist：

fitdist(myData, "t")

然后我得到

Error in mledist(data, distname, start, fix.arg, ...) :
'start' must be defined as a named list for this distribution

这是什么意思？我试图查看文档，但没有告诉我什么。我只想拟合一个 t 分布，这太令人沮丧了：P

谢谢！

Answer 1

Start 是对分布参数的初始猜测。涉及对数，因为它使用最大似然，因此使用对数似然。

library(fitdistrplus)
dat <- rt(100, df=10)
fit <- fitdist(dat, "t", start=list(df=2))