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(LBI) > # Extreme value test > > LIbin(0, 1) y n PE LL UL k 0.000000 1.000000 0.000000 0.000000 0.975008 40.000000 > LIbin(0, 2) y n PE LL UL k 0.000000 2.000000 0.000000 0.000000 0.841889 40.000000 > LIbin(0, 10) y n PE LL UL k 0.0000000 10.0000000 0.0000000 0.0000000 0.2015402 9.4962937 > LIbin(0, 1e308) y n PE LL UL k 0e+00 1e+308 0e+00 0e+00 -1e-08 1e+00 > LIbin(1, 1) y n PE LL UL k 1.00000000 1.00000000 1.00000000 0.02499196 1.00000000 40.00000000 > LIbin(10, 10) y n PE LL UL k 10.0000000 10.0000000 1.0000000 0.7984598 1.0000000 9.4962937 > LIbin(1e308, 1e308) y n PE LL UL k 1e+308 1e+308 1e+00 1e+00 1e+00 1e+00 > > LIpois(0) PE LL UL k 0.000000 0.000000 3.688879 40.000000 > LIpois(1) PE LL UL k 1.000000000 0.009309443 6.571648676 40.000000000 > LIpois(2) PE LL UL k 2.0000000 0.1237596 8.6079586 40.0000000 > LIpois(3) PE LL UL k 3.0000000 0.3643815 10.4259307 40.0000000 > > #LInorm > LInorm(c(0, 1, 2)) PE LL UL mean 1.0000000 -0.2804093 2.280409 sd 0.8164966 0.3620906 4.529665 var 0.6666667 0.1311096 20.517862 attr(,"n") [1] 3 attr(,"k") [1] 40 attr(,"log(k)") [1] 3.688879 attr(,"maxLogLik") [1] -3.648618 > > RDLI(0, 1, 0, 1) p1 p2 RD lower upper k 0.000000 0.000000 0.000000 -0.975008 0.975008 40.000000 > RDLI(0, 1, 1, 1) p1 p2 RD lower upper k 0.0000000 1.0000000 -1.0000000 -1.0000000 0.6837781 40.0000000 > RDLI(0, 10, 0, 10) p1 p2 RD lower upper k 0.0000000 0.0000000 0.0000000 -0.1873537 0.1873536 7.9624793 > RDLI(0, 10, 10, 10) p1 p2 RD lower upper k 0.0000000 1.0000000 -1.0000000 -1.0000000 -0.8029295 7.9624793 > RDLI(1, 1, 0, 1) p1 p2 RD lower upper k 1.0000000 0.0000000 1.0000000 -0.6837781 1.0000000 40.0000000 > RDLI(1, 1, 1, 1) p1 p2 RD lower upper k 1.000000 1.000000 0.000000 -0.975008 0.975008 40.000000 > RDLI(10, 10, 0, 10) p1 p2 RD lower upper k 1.0000000 0.0000000 1.0000000 0.8029295 1.0000000 7.9624793 > RDLI(10, 10, 10, 10) p1 p2 RD lower upper k 1.0000000 1.0000000 0.0000000 -0.1873536 0.1873537 7.9624793 > > proc.time() user system elapsed 0.21 0.06 0.21