我想了解绘制 fitdistrplus 对象的问题:
ftnb = tryCatch(expr = fitdist(unlist(x[,ftr]),distr="nbinom",discrete = TRUE),
error = function(err) {})
if (class(ftnb) %in% c("fitdist")) {
print(summary(ftnb))
plot(ftnb))
}
我的分布拟合总结如下:
Fitting of the distribution ' nbinom ' by maximum likelihood Parameters :
estimate Std. Error
size 5.915112e-01 NA
mu 1.330646e+11 NA
Loglikelihood: -2543.973 AIC: 5091.946 BIC: 5097.074
Correlation matrix:
[1] NA
我显然看不出合适,因为我在 2 个参数估计的标准误差中看到了“NA”。plot 命令出现以下错误,我想了解一下:
Error in seq.default(min(xval), max(xval), by = 1) : 'by' argument is much too small
请帮忙。