0

我正在尝试使用包中的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))

但是,会发生不同的错误: 在此处输入图像描述

这个问题有什么解决办法吗?谢谢你的帮助。

4

1 回答 1

1

试试这个代码,即使它会给你 NA 标准错误:

data_you <- 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)

fit_bur <- fitdist(data_you,"burr", start = list(shape1 = 0.3, shape2 = 1, rate = 1),lower=c(0,0))
于 2018-10-25T12:03:11.877 回答