在R中使用fitdist时出错-必须定义dllogis函数

Question

通过运行以下代码：

set.seed(12)
eta <- 10
nu <- 2
n <- 100
t <- eta*(((-1)*log(1-runif(n)))^(1/nu))# generates T_i
t
library("fitdistrplus")
fitll <- fitdist(t, "llogis")

在R中，我收到了这个错误：

Error in fitdist(t, "llogis") : The  dllogis  function must be defined

我之前安装过包"fitdistrplus“。对如何修复这个错误有什么想法吗？对于其他发行版，例如"gamma“，此代码可以很好地工作。

Answer 1

尝尝这个。看起来需要调用包含所需发行版的包(flexsurv)：

set.seed(12)
eta <- 10
nu <- 2
n <- 100
t <- eta*(((-1)*log(1-runif(n)))^(1/nu))# generates T_i
t
library("fitdistrplus")
library(flexsurv)
fitll <- fitdist(t, "llogis")
fitll

输出：

fitll
Fitting of the distribution ' llogis ' by maximum likelihood 
Parameters:
      estimate Std. Error
shape 3.141862  0.2609741
scale 7.659338  0.4246839