我正在尝试使用包中的fitdist
函数fitdistrplus
来使我的数据适合 Burr 分布。
library(fitdistrplus)
library(actuar)
fit_bur <- fitdist(data=c(6,8,6,3,3,6,6,6,7,8,4,5,1,4,4,0,0,0,0,7,5,6,6,5,3,4,8,6,5,6,3,6,5,4,3,9,4,2,4), method = "mle",distr="burr",
start = list(shape1 = 0.3, shape2 = 1, rate = 1))
受此处发布的问题的启发,我注意到错误可能是我的数据包含 0 值,因此我将 0.1 添加到所有 0 值。
fit_bur <- fitdist(data=c(6,8,6,3,3,6,6,6,7,8,4,5,1,4,4,0.1,0.1,0.1,0.1,7,5,6,6,5,3,4,8,6,5,6,3,6,5,4,3,9,4,2,4), method = "mle",distr="burr",
start = list(shape1 = 0.3, shape2 = 1, rate = 1))
这个问题有什么解决办法吗?谢谢你的帮助。