r - 在 for 循环中使用 fitdist 的 fitdistr 包中的错误

Question

我正在运行以下代码，即将负二项式分布的 AIC 值保存在向量中。

max_nbinom_aic <-data.frame(NBINOM_AIC=numeric(0))
 for(i in 1:336) {
     a <- max_samples[,i]
     d <- fitdist(a,distr ="nbinom") 
 max_nbinom_aic <- rbind(max_nbinom_aic,d$aic)
}

我收到以下错误消息：

fitdist(a, distr = "nbinom") 中的错误：函数 mle 无法估计参数，错误代码为 100 另外：警告消息：1：在 sqrt(diag(varcovar)) 中：产生了 NaN 2：在 sqrt (1/diag(V)) : NaNs 产生 3: 在 cov2cor(varcovar) : diag(.) 有 0 或 NA 条目；非限定结果值得怀疑

此外，输出仅包含 62 个值，而我期望 336 个值，因为我有 336 组分布。