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 <- rfpareto(n, 1, 2, 2, 2, scale=2)
> 
> fitdistrplus:::startarg_fellerpareto_family(x, "pareto4")
$min
[1] 1.281281

$shape1
[1] 1.970823

$shape2
[1] 2.278414

$scale
[1] 2.62346

> 
> fitdistrplus:::startarg_fellerpareto_family(x, "pareto3")
$min
[1] 1.281281

$shape
[1] 2.682464

$scale
[1] 1.741837

> 
> fitdistrplus:::startarg_fellerpareto_family(x, "pareto2")
$min
[1] 1.281281

$shape
[1] 1.362517

$scale
[1] 1.735447

> 
> fitdistrplus:::startarg_fellerpareto_family(x, "pareto1")
$shape
[1] 1.117406

$min
[1] 1.281281

> 
> fitdistrplus:::startarg_fellerpareto_family(x, "pareto")
$shape
[1] 0.9551719

$scale
[1] 1.735447

> 
> fitdistrplus:::startarg_fellerpareto_family(x, "llogis")
$shape
[1] 2.942622

$scale
[1] 4.759425

> 
> fitdist(x, "pareto4")
Fitting of the distribution ' pareto4 ' by maximum likelihood 
Parameters:
       estimate Std. Error
min    1.294655 0.04590328
shape1 1.768239 0.31737902
shape2 2.276273 0.16465235
scale  2.423419 0.25945779
> fitdist(x, "pareto3")
Fitting of the distribution ' pareto3 ' by maximum likelihood 
Parameters:
      estimate Std. Error
min   1.089242 0.09150735
shape 3.057758 0.18642591
scale 1.945448 0.10445751
> fitdist(x, "pareto2", lower=0)
Fitting of the distribution ' pareto2 ' by maximum likelihood 
Parameters:
         estimate Std. Error
min      1.348716        NaN
shape 1228.519336 0.07186649
scale 2485.640895 0.10290494
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
> fitdist(x, "pareto1", lower=0)
Fitting of the distribution ' pareto1 ' by maximum likelihood 
Parameters:
      estimate Std. Error
shape 1.185491        NaN
min   1.348716        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
> fitdist(x, "pareto", lower=0)
Fitting of the distribution ' pareto ' by maximum likelihood 
Parameters:
       estimate Std. Error
shape  9793.147 0.04461554
scale 33029.634 0.09771598
> 
> fitdist(x, "llogis", lower=0)
Fitting of the distribution ' llogis ' by maximum likelihood 
Parameters:
      estimate Std. Error
shape 4.870201        NaN
scale 3.084839        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
> 
> 
> fitdistrplus:::startarg_fellerpareto_family(x, "fpareto")
$min
[1] 1.281281

$shape1
[1] 4.817928

$shape2
[1] 2.064787

$shape3
[1] 5.237083

$scale
[1] 1.557258

> 
> fitdist(x, "fpareto", lower=0)
Fitting of the distribution ' fpareto ' by maximum likelihood 
Parameters:
        estimate Std. Error
min    1.3164307 0.01408560
shape1 1.3425857        NaN
shape2 2.7061466 0.01712428
shape3 0.7581854        NaN
scale  2.3326050        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
> 
> fitdistrplus:::startarg_fellerpareto_family(x-1, "trbeta")
$shape1
[1] 5.704643

$shape2
[1] 2.001012

$shape3
[1] 9.109486

$scale
[1] 1.233544

> 
> fitdist(x-1, "trbeta", lower=0)
Fitting of the distribution ' trbeta ' by maximum likelihood 
Parameters:
       estimate Std. Error
shape1 2.694614 0.08895979
shape2 1.750024 0.04261088
shape3 2.524013 0.32967326
scale  2.129175 0.05071407
> 
> 
> fitdistrplus:::startarg_fellerpareto_family(x-1, "genpareto")
$shape1
[1] 4.533974

$shape2
[1] 13.03459

$scale
[1] 0.6603603

> 
> fitdist(x-1, "genpareto", lower=0)
Fitting of the distribution ' genpareto ' by maximum likelihood 
Parameters:
       estimate Std. Error
shape1 7.189862  0.4644227
shape2 6.965237  3.6829652
scale  2.105227  0.2600406
> 
> fitdistrplus:::startarg_fellerpareto_family(x-1, "paralogis")
$shape
[1] 0.6901228

$scale
[1] 2.662399

> 
> fitdist(x-1, "paralogis", lower=0)
Fitting of the distribution ' paralogis ' by maximum likelihood 
Parameters:
      estimate Std. Error
shape 2.458679  0.0154506
scale 3.308168  0.4577320
> 
> 
> x <- rfpareto(n, 0, 1, 2, 2, scale=2)
> 
> fitdistrplus:::startarg_fellerpareto_family(x, "invburr")
$shape1
[1] 11.76437

$shape2
[1] 1.791445

$scale
[1] 1.529135

> fitdistrplus:::startarg_fellerpareto_family(x, "invpareto")
$shape
[1] 1.440887

$scale
[1] 0.3434468

> fitdistrplus:::startarg_fellerpareto_family(x, "invparalogis")
$shape
[1] 1.203247

$scale
[1] 0.2772426

> 
> 
> proc.time()
   user  system elapsed 
   6.81    0.34    7.14