R Under development (unstable) (2024-06-13 r86735 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.3620787  4.529653
var  0.6666667  0.1311010 20.517752
attr(,"n")
[1] 3
attr(,"k")
[1] 40
attr(,"log(k)")
[1] 3.688879
attr(,"2*log(k)")
[1] 7.377759
> 
> 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.20    0.10    0.25