R Under development (unstable) (2024-12-17 r87446 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.13
Date:     
BugReport: https://github.com/spedygiorgio/mbbefd/issues


> library(fitdistrplus)
> 
> 
> #oiunif
> n <- 1e3
> nboot <- 1000
> nboot <- 10
> x <- roiunif(n, 1/6)
> f1 <- fitDR(x, "oiunif", method="mle")
Warning message:
In fitdist(x, distr = dist, method = method, start = start, calcvcov = TRUE,  :
  The doiunif function should return a vector of with NaN values when input has inconsistent values and not raise an error
> summary(f1)
Fitting of the distribution ' oiunif ' by maximum likelihood 
Parameters : 
   estimate Std. Error
p1    0.155 0.01144403
Loglikelihood:  -431.2849   AIC:  864.5699   BIC:  869.4776 
> summary(fitdist(x, "oiunif", method="mle", start=list(p1=1/2))) #check
Fitting of the distribution ' oiunif ' by maximum likelihood 
Parameters : 
    estimate Std. Error
p1 0.1550017 0.01144413
Loglikelihood:  -431.2849   AIC:  864.5699   BIC:  869.4776 
Warning message:
In fitdist(x, "oiunif", method = "mle", start = list(p1 = 1/2)) :
  The doiunif function should return a vector of with NaN values when input has inconsistent values and not raise an error
> 
> b1 <- bootDR(f1, niter=nboot, silent=FALSE)
Warning message:
In fitdist(x, distr = dist, method = method, start = start, calcvcov = TRUE,  :
  The doiunif function should return a vector of with NaN values when input has inconsistent values and not raise an error
> summary(b1)
Parametric bootstrap medians and 95% percentile CI 
 Median    2.5%   97.5% 
0.16100 0.13645 0.17055 
> 
> plot(b1)
> abline(v=1/6, col="red")
> 
> hist(b1$estim[,1])
> abline(v=1/6, col="red")
> 
> 
> f2 <- fitDR(x, "oiunif", method="tlmme")
Warning message:
In fitdist(x, distr = dist, method = method, start = start, calcvcov = TRUE,  :
  The doiunif function should return a vector of with NaN values when input has inconsistent values and not raise an error
> 
> 
> 
> proc.time()
   user  system elapsed 
   1.85    0.31    2.18