R Under development (unstable) (2025-07-02 r88374 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(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 > require(fitdistrplus) Loading required package: fitdistrplus Loading required package: MASS Loading required package: survival > n <- 1e4 > > > #### Inverse exponential #### > > x <- rinvexp(n, 3) > > if(FALSE) + { + fitdistrplus:::startarg_invtransgamma_family(x, "invexp") + fitdistrplus:::startarg_transgamma_family(1/x, "exp") + } > > cdfcomp(fitdist(x, "invexp"), xlogscale = TRUE, do.points = FALSE) > > > #### Inverse transformed gamma #### > > x <- rinvtrgamma(n, 3, 3, 10) > > if(FALSE) + { + fitdistrplus:::startarg_invtransgamma_family(x, "invtrgamma") + fitdistrplus:::startarg_transgamma_family(1/x, "trgamma") + + cutshapeparam(list("shape0"=1e3, "shape33"=1e-10, "theta"=22)) + cutshapeparam(list("shape0"=1e3, "shape33"=1e3, "theta"=22)) + + cutshapeparam(fitdistrplus:::startarg_invtransgamma_family(x, "invtrgamma")) + } > > cdfcomp(fitdist(x, "invtrgamma", lower=0), xlogscale = TRUE, do.points = FALSE) 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 > cdfcomp(fitdist(x, "invtrgamma", lower=0, start=list("shape1"=100, "shape2"=1, "scale"=1/2)), + xlogscale = TRUE, do.points = FALSE) 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 > curve(pinvtrgamma(x, 10, 1, 1)) > curve(pinvtrgamma(x, 25, 1, 1/2)) > curve(pinvtrgamma(x, 100, 1, 1/2)) > > x <- rinvparalogis(n, 2, 10) > > cdfcomp(fitdist(x, "invparalogis"), xlogscale = TRUE, do.points = FALSE) > cdfcomp(fitdist(x, "invparalogis", lower=0, start=list("shape"=1e-5, "scale"=1/2), + control=list(trace=1, REPORT=1)), + xlogscale = TRUE, do.points = FALSE) Nelder-Mead direct search function minimizer function value for initial parameters = 21.923456 Scaled convergence tolerance is 3.26685e-07 Stepsize computed as 0.050000 BUILD 3 21.923456 4.897738 EXTENSION 5 21.923456 4.088382 EXTENSION 7 4.897738 2.299196 EXTENSION 9 4.088382 1.445241 EXTENSION 11 2.299196 0.389461 EXTENSION 13 1.445241 0.169056 LO-REDUCTION 15 0.389461 0.137860 HI-REDUCTION 17 0.169056 0.132500 REFLECTION 19 0.137860 0.067586 HI-REDUCTION 21 0.132500 0.067586 EXTENSION 23 0.086832 0.012889 EXTENSION 25 0.067586 -0.158650 EXTENSION 27 0.012889 -0.522993 HI-REDUCTION 29 -0.125411 -0.522993 LO-REDUCTION 31 -0.158650 -0.522993 HI-REDUCTION 33 -0.335275 -0.522993 HI-REDUCTION 35 -0.433995 -0.522993 LO-REDUCTION 37 -0.455018 -0.522993 REFLECTION 39 -0.513602 -0.553659 HI-REDUCTION 41 -0.522993 -0.553659 EXTENSION 43 -0.548770 -0.607820 LO-REDUCTION 45 -0.553659 -0.607820 REFLECTION 47 -0.607487 -0.634187 HI-REDUCTION 49 -0.607820 -0.634187 EXTENSION 51 -0.627939 -0.667283 HI-REDUCTION 53 -0.634187 -0.667283 LO-REDUCTION 55 -0.648504 -0.667283 REFLECTION 57 -0.664409 -0.667775 HI-REDUCTION 59 -0.667283 -0.668500 HI-REDUCTION 61 -0.667775 -0.669231 HI-REDUCTION 63 -0.668500 -0.669438 LO-REDUCTION 65 -0.669231 -0.669438 HI-REDUCTION 67 -0.669299 -0.669507 LO-REDUCTION 69 -0.669438 -0.669541 HI-REDUCTION 71 -0.669507 -0.669541 HI-REDUCTION 73 -0.669535 -0.669554 LO-REDUCTION 75 -0.669541 -0.669554 HI-REDUCTION 77 -0.669553 -0.669555 LO-REDUCTION 79 -0.669554 -0.669557 HI-REDUCTION 81 -0.669555 -0.669557 HI-REDUCTION 83 -0.669557 -0.669557 LO-REDUCTION 85 -0.669557 -0.669558 Exiting from Nelder Mead minimizer 87 function evaluations used Nelder-Mead direct search function minimizer function value for initial parameters = -0.669784 Scaled convergence tolerance is 9.98057e-09 Stepsize computed as 0.197906 BUILD 3 0.535755 -0.669784 HI-REDUCTION 5 -0.148439 -0.669784 HI-REDUCTION 7 -0.480218 -0.669784 HI-REDUCTION 9 -0.606293 -0.669784 HI-REDUCTION 11 -0.647645 -0.669784 LO-REDUCTION 13 -0.657091 -0.669784 LO-REDUCTION 15 -0.664877 -0.669784 LO-REDUCTION 17 -0.667062 -0.669784 LO-REDUCTION 19 -0.668068 -0.669784 LO-REDUCTION 21 -0.669202 -0.669784 HI-REDUCTION 23 -0.669476 -0.669784 LO-REDUCTION 25 -0.669585 -0.669784 LO-REDUCTION 27 -0.669704 -0.669784 LO-REDUCTION 29 -0.669744 -0.669784 LO-REDUCTION 31 -0.669758 -0.669784 LO-REDUCTION 33 -0.669775 -0.669784 LO-REDUCTION 35 -0.669781 -0.669784 HI-REDUCTION 37 -0.669782 -0.669784 LO-REDUCTION 39 -0.669783 -0.669784 HI-REDUCTION 41 -0.669784 -0.669784 HI-REDUCTION 43 -0.669784 -0.669784 LO-REDUCTION 45 -0.669784 -0.669784 HI-REDUCTION 47 -0.669784 -0.669784 HI-REDUCTION 49 -0.669784 -0.669784 Exiting from Nelder Mead minimizer 51 function evaluations used Nelder-Mead direct search function minimizer function value for initial parameters = -0.669784 Scaled convergence tolerance is 9.98057e-09 Stepsize computed as 0.197905 BUILD 3 0.535833 -0.669784 HI-REDUCTION 5 -0.148352 -0.669784 HI-REDUCTION 7 -0.480146 -0.669784 HI-REDUCTION 9 -0.606244 -0.669784 HI-REDUCTION 11 -0.647615 -0.669784 LO-REDUCTION 13 -0.657081 -0.669784 LO-REDUCTION 15 -0.664888 -0.669784 LO-REDUCTION 17 -0.667074 -0.669784 LO-REDUCTION 19 -0.668072 -0.669784 LO-REDUCTION 21 -0.669200 -0.669784 HI-REDUCTION 23 -0.669472 -0.669784 LO-REDUCTION 25 -0.669583 -0.669784 LO-REDUCTION 27 -0.669705 -0.669784 LO-REDUCTION 29 -0.669745 -0.669784 LO-REDUCTION 31 -0.669759 -0.669784 LO-REDUCTION 33 -0.669775 -0.669784 LO-REDUCTION 35 -0.669780 -0.669784 LO-REDUCTION 37 -0.669781 -0.669784 LO-REDUCTION 39 -0.669783 -0.669784 LO-REDUCTION 41 -0.669784 -0.669784 LO-REDUCTION 43 -0.669784 -0.669784 LO-REDUCTION 45 -0.669784 -0.669784 HI-REDUCTION 47 -0.669784 -0.669784 HI-REDUCTION 49 -0.669784 -0.669784 Exiting from Nelder Mead minimizer 51 function evaluations used > > > proc.time() user system elapsed 27.28 0.76 28.04