R-使用fitdistrplus软件包拟合广义伽玛分布

编程问答 2022-08-15

问题描述

我正在尝试将通用的伽马分布（威布尔，正态，指数和对数正态分布）适合我国的Covid-19数据的日常情况（五个月）。我以前运行的代码如下：

#required library
library("MASS")
library("survival")
library("fitdistrplus")

但是，当我尝试拟合每个发行版时，出现以下错误。

fn_w <- fitdist(daily_cases,"weibull") #weibull distribution
Error in checkparamlist(arg_startfix$start.arg,arg_startfix$fix.arg,: 
  'start' should not have NA or NaN values.

fn_g <- fitdist(daily_cases,"gamma") #for gamma distribution
<simpleError in optim(par = vstart,fn = fnobj,fix.arg = fix.arg,obs = data,gr = gradient,ddistnam = ddistname,hessian = TRUE,method = meth,lower = lower,upper = upper,...): function cannot be evaluated at initial parameters>
Error in fitdist(daily_cases,"gamma") : 
  the function mle Failed to estimate the parameters,with the error code 100

> fn_ln <- fitdist(daily_cases,"lnorm") #for log normal distribution
Error in computing default starting values.
Error in manageparam(start.arg = start,: 
  Error in start.arg.default(obs,distname) : 
  values must be positive to fit a lognormal distribution

只有指数分布才能有效运行。

> fn_exp <- fitdist(daily_cases,"exp")
> summary(fn_exp)
Fitting of the distribution ' exp ' by maximum likelihood 
Parameters : 
        estimate   Std. Error
rate 0.002949398 0.0002060355
Loglikelihood:  -1044.402   AIC:  2090.803   BIC:  2093.834

我该如何纠正？ 预先感谢

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）

distributiondistributionrr