R Under development (unstable) (2024-07-10 r86888 ucrt) -- "Unsuffered Consequences" Copyright (C) 2024 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. > library(fitdistrplus) Loading required package: MASS Loading required package: survival > library(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 14.7947393 NaN shape2 1.9319992 NaN shape3 0.3211336 NaN scale 1.9923966 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 5.07 0.26 5.32