R Under development (unstable) (2024-10-16 r87241 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(mbbefd)
Loading required package: fitdistrplus
Loading required package: MASS
Loading required package: survival
Loading required package: alabama
Loading required package: numDeriv
Loading required package: Rcpp
Package:  mbbefd
Version:  0.8.12
Date:     
BugReport: https://github.com/spedygiorgio/mbbefd/issues


> 
> data(beaonre)
> x <- beaonre$ClaimCost/beaonre$SumInsured
> 
> 
> # #Nelder Mead
> # mledist(x[x!=1], "gbeta", start=list(shape0=3.373523e-03, 
> #     shape1=2.913619e+02, shape2=7.198632e+00), control=list(trace=1,REPORT=1))
> # #L-BFGS-B
> # mledist(x[x!=1], "gbeta", lower=0, start=list(shape0=3.373523e-03, 
> #     shape1=2.913619e+02, shape2=7.198632e+00), control=list(trace=1, REPORT=1, fnscale=1e-6))
> # 
> # fitDR(x, "oigbeta", method="mle", control=list(trace=1, REPORT=1))
> # fitDR(x, "MBBEFD", method="mle", control=list(trace=1, REPORT=1))
> 
> 
> dlist <- c("oistpareto", "oibeta", "oigbeta", "mbbefd", "MBBEFD")
> dlist <- c("oistpareto", "oibeta", "mbbefd", "MBBEFD")
> flist <- lapply(dlist, function(d) {print(d);
+   fitDR(x, d, method="mle")})
[1] "oistpareto"
[1] "oibeta"
[1] "mbbefd"
<simpleError in if (obj * sign(mu) > obj.old * sign(mu)) break: missing value where TRUE/FALSE needed>
[1] "MBBEFD"
<simpleError in if (obj * sign(mu) > obj.old * sign(mu)) break: missing value where TRUE/FALSE needed>
Warning message:
In cov2cor(f1$vcov) :
  diag(V) had non-positive or NA entries; the non-finite result may be dubious
> names(flist) <- dlist
> 
> 
> cdfcomp(flist, do.points=FALSE, leg=dlist)
> ppcomp(flist, leg=dlist, fitpch=".", addlegend = FALSE)
> legend("bottomright", fill=c("red", "green", "blue", "cyan"), leg=dlist)
> 
> qqcomp(flist, leg=dlist, use.ppoints=TRUE)
> 
> proc.time()
   user  system elapsed 
   2.68    0.23    2.89