0

我正在尝试绘制不同的分布,以查看最适合我的数据(“x”)的分布。我正在使用包 fitdistrplus 来做到这一点。除了“tnorm”(截断的正态分布)之外,所有分布都可以工作,它会显示一条错误消息:

fit_tn<-fitdist(x,"tnorm",start=list(mean=mean(x),sd=sd(x)),fix.arg=list(a=0,b=Inf))
$start.arg
$start.arg$mean
[1] 28.61737

$start.arg$sd
[1] 29.45779


$fix.arg
$fix.arg$a
[1] 0

$fix.arg$b
[1] Inf


Error in fitdist(x, "tnorm", start = list(mean = mean(x), sd = sd(x)),  :
  the function mle failed to estimate the parameters,
                with the error code 10

任何想法导致此错误消息以及如何解决它?

非常感谢任何帮助!!!

有效的发行版:

fit_ex<-fitdist(x,"exp")
fit_nm<-fitdist(x,"norm")
fit_ln<-fitdist(x,"lnorm")
fit_gm<-fitdist(x,"gamma")
fit_wb<-fitdist(x,"weibull")

使用的包:

library(fitdistrplus)
library(extraDistr)
library(actuar)

关于我的数据的一些信息:

min(x)
[1] 1
max(x)
[1] 235
length(x)
[1] 426
4

0 回答 0