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: In sqrt(diag(varcovar)): NaN が生成されました 2: In sqrt (1/diag(V)) : 生成された NaN 3: cov2cor(varcovar) : diag(.) には 0 または NA エントリがありました。非有限の結果は疑わしい

また、出力には 62 個の値しか含まれていませんが、336 セットの分布があるため、336 個の値を期待しています。