R Under development (unstable) (2026-01-19 r89308 ucrt) -- "Unsuffered Consequences"
Copyright (C) 2026 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.14
Date:     
BugReport: https://github.com/spedygiorgio/mbbefd/issues


> 
> testfunc <- function(x)
+   c(summary(x), sd=sd(x), tl=etl(x))
> extensive <- TRUE
> extensive <- FALSE
> 
> 
> # test invalid param
> n <- 5
> a <- 0 
> b <- -1/2
> mbbefd:::rmbbefdCpp(n, a, b)
[1] 1 1 1 1 1
> mbbefd:::rmbbefdR(n, a, b)
[1] NaN NaN NaN NaN NaN
> g <- 1/2 
> b <- 3
> 
> mbbefd:::rMBBEFDCpp(n, g, b)
[1] NaN NaN NaN NaN NaN
> mbbefd:::rMBBEFDR(n, g, b)
[1] NaN NaN NaN NaN NaN
> 
> 
> #test of MBBEFD(a,b) distribution
> 
> n <- 10
> a <- 0 
> b <- 1/2
> 
> mbbefd:::rmbbefdCpp(n, a, b)
 [1] 1 1 1 1 1 1 1 1 1 1
> mbbefd:::rmbbefdR(n, a, b)
 [1] 1 1 1 1 1 1 1 1 1 1
> if(extensive)
+ {
+   n <- 1e6
+   print(testfunc(mbbefd:::rmbbefdCpp(n, a, b)))
+   print(testfunc(mbbefd:::rmbbefdR(n, a, b)))
+ }
> 
> a <- 1/2 
> b <- 1
> n <- 10
> 
> mbbefd:::rmbbefdCpp(n, a, b)
 [1] 1 1 1 1 1 1 1 1 1 1
> mbbefd:::rmbbefdR(n, a, b)
 [1] 1 1 1 1 1 1 1 1 1 1
> if(extensive)
+ {
+   n <- 1e6
+   print(testfunc(mbbefd:::rmbbefdCpp(n, a, b)))
+   print(testfunc(mbbefd:::rmbbefdR(n, a, b)))
+ }
> 
> a <- -1/2 
> b <- 3
> n <- 10
> 
> mbbefd:::rmbbefdCpp(n, a, b)
 [1] 1.0000000 1.0000000 1.0000000 0.1876183 1.0000000 1.0000000 1.0000000
 [8] 1.0000000 0.2660039 1.0000000
> mbbefd:::rmbbefdR(n, a, b)
 [1] 0.90613879 1.00000000 1.00000000 0.05606251 1.00000000 1.00000000
 [7] 1.00000000 0.18630839 1.00000000 0.27373165
> if(extensive)
+ {
+   n <- 1e6
+   print(testfunc(mbbefd:::rmbbefdCpp(n, a, b)))
+   print(testfunc(mbbefd:::rmbbefdR(n, a, b)))
+ }
> 
> a <- Inf
> b <- 1/3
> n <- 10
> 
> mbbefd:::rmbbefdCpp(n, a, b)
 [1] 0.5586340 1.0000000 0.7054731 0.8494513 0.6271367 1.0000000 1.0000000
 [8] 0.1575667 0.2638093 1.0000000
> mbbefd:::rmbbefdR(n, a, b)
 [1] 1.0000000 0.2826278 0.8389927 0.6060484 1.0000000 0.7747242 1.0000000
 [8] 0.1218277 0.4895388 0.1672389
> if(extensive)
+ {
+   n <- 1e6
+   print(testfunc(mbbefd:::rmbbefdCpp(n, a, b)))
+   print(testfunc(mbbefd:::rmbbefdR(n, a, b)))
+ }
> 
> 
> 
> 
> #test of MBBEFD(g,b) distribution
> 
> n <- 10
> g <- 1 
> b <- 1/2
> 
> mbbefd:::rMBBEFDCpp(n, g, b)
 [1] 1 1 1 1 1 1 1 1 1 1
> mbbefd:::rMBBEFDR(n, g, b)
 [1] 1 1 1 1 1 1 1 1 1 1
> if(extensive)
+ {
+   n <- 1e6
+   print(testfunc(mbbefd:::rMBBEFDCpp(n, g, b)))
+   print(testfunc(mbbefd:::rMBBEFDR(n, g, b)))
+ }
> 
> n <- 10
> g <- 2 
> b <- 0
> 
> mbbefd:::rMBBEFDCpp(n, g, b)
 [1] 1 1 1 1 1 1 1 1 1 1
> mbbefd:::rMBBEFDR(n, g, b)
 [1] 1 1 1 1 1 1 1 1 1 1
> if(extensive)
+ {
+   n <- 1e6
+   print(testfunc(mbbefd:::rMBBEFDCpp(n, g, b)))
+   print(testfunc(mbbefd:::rMBBEFDR(n, g, b)))
+ }
> 
> 
> n <- 10
> g <- 2 
> b <- 1/2
> 
> mbbefd:::rMBBEFDCpp(n, g, b)
 [1] 0.04370886 0.88747975 1.00000000 0.04954522 0.42454183 0.02079493
 [7] 0.06945197 0.43622669 1.00000000 0.26760277
> mbbefd:::rMBBEFDR(n, g, b)
 [1] 0.44701923 0.57324119 1.00000000 0.27876252 0.20205631 0.13865153
 [7] 0.02998699 1.00000000 1.00000000 1.00000000
> if(extensive)
+ {
+   n <- 1e6
+   print(testfunc(mbbefd:::rMBBEFDCpp(n, g, b)))
+   print(testfunc(mbbefd:::rMBBEFDR(n, g, b)))
+ }
> 
> n <- 10
> g <- 2
> b <- 1
> 
> mbbefd:::rMBBEFDCpp(n, g, b)
 [1] 1.0000000 0.1916770 1.0000000 1.0000000 1.0000000 0.1916573 0.5437108
 [8] 0.8144119 0.1021862 1.0000000
> mbbefd:::rMBBEFDR(n, g, b)
 [1] 0.82924067 1.00000000 0.40315268 0.66846572 1.00000000 1.00000000
 [7] 0.02716495 0.20187392 1.00000000 0.17658812
> if(extensive)
+ {
+   n <- 1e6
+   print(testfunc(mbbefd:::rMBBEFDCpp(n, g, b)))
+   print(testfunc(mbbefd:::rMBBEFDR(n, g, b)))
+ }
> 
> n <- 10
> g <- 2 
> b <- 3
> 
> mbbefd:::rMBBEFDCpp(n, g, b)
 [1] 1.00000000 1.00000000 1.00000000 0.01923398 0.01033167 1.00000000
 [7] 0.35164123 0.05267305 1.00000000 1.00000000
> mbbefd:::rMBBEFDR(n, g, b)
 [1] 0.575824305 0.034622927 0.001365621 1.000000000 0.544291198 0.133321135
 [7] 1.000000000 1.000000000 0.167691499 0.603285949
> if(extensive)
+ {
+   n <- 1e6
+   print(testfunc(mbbefd:::rMBBEFDCpp(n, g, b)))
+   print(testfunc(mbbefd:::rMBBEFDR(n, g, b)))
+ }
> 
> proc.time()
   user  system elapsed 
   1.54    0.26    1.79