R Under development (unstable) (2025-01-06 r87534 ucrt) -- "Unsuffered Consequences"
Copyright (C) 2025 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> require("fitdistrplus")
Loading required package: fitdistrplus
Loading required package: MASS
Loading required package: survival
> require("actuar")
Loading required package: actuar

Attaching package: 'actuar'

The following objects are masked from 'package:stats':

    sd, var

The following object is masked from 'package:grDevices':

    cm

> 
> set.seed(1234)
> 
> n <- 1e3
> 
> x <- rlgamma(n, 1, 1)
> fitdistrplus:::startarg_othercontinuous_actuar_family(x, "lgamma")
$shapelog
[1] 0.9538171

$ratelog
[1] 0.9603243

> fitdist(x, "lgamma")
Fitting of the distribution ' lgamma ' by maximum likelihood 
Parameters:
         estimate Std. Error
shapelog 1.000753 0.03940913
ratelog  1.007445 0.05087335
> 
> 
> x <- rgumbel(n, 1, 1)
> fitdistrplus:::startarg_othercontinuous_actuar_family(x, "gumbel")
$scale
[1] 1.009752

$alpha
[1] 0.997063

> fitdist(x, "gumbel")
Fitting of the distribution ' gumbel ' by maximum likelihood 
Parameters:
       estimate Std. Error
scale 1.0197296  0.0251320
alpha 0.9853615  0.0339739
> 
> 
> x <- rinvgauss(n, pi, 2*pi)
> fitdistrplus:::startarg_othercontinuous_actuar_family(x, "invgauss")
$mean
[1] 3.199838

$dispersion
[1] 0.1439096

> fitdist(x, "invgauss")
Fitting of the distribution ' invgauss ' by maximum likelihood 
Parameters:
            estimate  Std. Error
mean       3.1995736 0.069266474
dispersion 0.1465136 0.006552722
> 
> 
> x <- rgenbeta(n, 2, 2, 2, scale=2)
> fitdistrplus:::startarg_othercontinuous_actuar_family(x, "genbeta")
$shape1
[1] 2.355739

$shape2
[1] 16.23328

$shape3
[1] 3

$scale
[1] 2.914548

> fitdist(x, "genbeta", lower=0)
Fitting of the distribution ' genbeta ' by maximum likelihood 
Parameters:
        estimate Std. Error
shape1 0.5749735        NaN
shape2 0.7786942        NaN
shape3 6.3005086        NaN
scale  1.1554106        NaN
Warning messages:
1: In cov2cor(varcovar) :
  diag(V) had non-positive or NA entries; the non-finite result may be dubious
2: In sqrt(diag(varcovar)) : NaNs produced
> 
> proc.time()
   user  system elapsed 
   1.92    0.25    2.15