The hardware and bandwidth for this mirror is donated by METANET, the Webhosting and Full Service-Cloud Provider.
If you wish to report a bug, or if you are interested in having us mirror your free-software or open-source project, please feel free to contact us at mirror[@]metanet.ch.
elgbd performs hypothesis testing for general block designs with empirical likelihood. The core computational routines are implemented with the ‘Eigen’ ‘C++’ library and ‘RcppEigen’ interface, with ‘OpenMP’ for parallel computation. Details of the testing procedures are given in Kim, MacEachern, and Peruggia (2023). This work was supported by the U.S. National Science Foundation under Grants No. SES-1921523 and DMS-2015552.
You can install the latest stable release of elgbd from CRAN.
You can install the development version of elgbd from GitHub.
library(elgbd)
set.seed(58234)
# Analysis of variance
el_aov(formula = Sepal.Length ~ Species, data = iris)
#> Call:
#> el_aov(formula = Sepal.Length ~ Species, data = iris)
#>
#> Minimizer:
#> 5.9595 5.9595 5.9595
#>
#> Statistic:
#> 94.9054
# All pairwise comparisons
data("clothianidin")
el_pairwise(clo ~ trt | blk, data = clothianidin, B = 1000)
#>
#> Empirical Likelihood Multiple Tests
#>
#> All pairwise comparisons
#>
#> Estimate Chisq Lwr.ci Upr.ci p.adj
#> Naked - Fungicide -1.0525 14.098 -2.1150 -0.2794 0.002 **
#> Naked - Low -1.6794 15.935 -2.6705 -0.6435 0.001 ***
#> Naked - High -3.1726 42.803 -4.5731 -1.6940 <0.001 ***
#> Fungicide - Low -0.6269 1.834 -1.6303 0.6448 0.523
#> Fungicide - High -2.1201 14.085 -3.2313 -0.7247 0.002 **
#> Low - High -1.4932 11.715 -2.7228 -0.3954 0.003 **
#> ---
#> Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
#>
#> k: 1, level: 0.05, method: AMC, cutoff: 6.6737
# Comparisons with control
el_pairwise(clo ~ trt | blk,
data = clothianidin, control = "Naked", method = "NB", B = 1000
)
#>
#> Empirical Likelihood Multiple Tests
#>
#> Comparisons with control
#>
#> Estimate Chisq Lwr.ci Upr.ci p.adj
#> Fungicide - Naked 1.0525 14.10 0.3506 1.9900 0.005 **
#> Low - Naked 1.6794 15.94 0.7586 2.5646 0.004 **
#> High - Naked 3.1726 42.80 1.8489 4.4282 0.003 **
#> ---
#> Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
#>
#> k: 1, level: 0.05, method: NB, cutoff: 6.3291
# Test for equal means
el_test(clo ~ trt | blk, clothianidin,
lhs = matrix(c(
1, -1, 0, 0,
0, 1, -1, 0,
0, 0, 1, -1
), byrow = TRUE, nrow = 3L)
)
#>
#> Empirical Likelihood Test
#>
#> General block designs
#>
#> Maximum EL estimates:
#> [1] -4.479086 -3.426618 -2.799700 -1.306530
#>
#> Statistic: 28.055These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.