R Under development (unstable) (2019-12-01 r77492) -- "Unsuffered Consequences" Copyright (C) 2019 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 (64-bit) 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. > ### Regression tests for naming of results > > set.seed(290875) > library("coin") Loading required package: survival > > y1 <- sample(1:20) > y2 <- rnorm(20) > x <- gl(2, 10, labels = c("A", "B")) > > > ### > ### Asymptotic > ### > > ### Scalar > asy_scl <- independence_test(y1 ~ x, distr = "asymptotic", teststat = "scalar") > s <- statistic(asy_scl) > > stopifnot(identical(dimnames( expectation(asy_scl) ), list("A", "") )) # non-matrix in < 1.4-0 > stopifnot(identical(dimnames( variance(asy_scl) ), list("A", "") )) # non-matrix in < 1.4-0 > stopifnot(identical(dimnames( covariance(asy_scl) ), list("A", "A") )) > stopifnot( is.null( names( statistic(asy_scl, type = "test") ) )) # named in < 1.3-0 > stopifnot(identical(dimnames( statistic(asy_scl, type = "linear") ), list("A", "") )) > stopifnot(identical(dimnames( statistic(asy_scl, type = "centered") ), list("A", "") )) > stopifnot(identical(dimnames( statistic(asy_scl, type = "standardized") ), list("A", "") )) > ## stopifnot( is.null( names( support(asy_scl) ) )) # NA > stopifnot( is.null( names( dperm(asy_scl, x = s) ) )) > stopifnot(identical( names( dperm(asy_scl, x = c(s = s)) ), "s" )) > stopifnot( is.null(dimnames( dperm(asy_scl, x = matrix(s, nrow = 1)) ) )) # non-matrix in < 1.3-1 > stopifnot(identical(dimnames( dperm(asy_scl, x = matrix(s, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # non-matrix in < 1.3-1 > stopifnot( is.null( names( pperm(asy_scl, q = s) ) )) > stopifnot(identical( names( pperm(asy_scl, q = c(s = s)) ), "s" )) > stopifnot( is.null(dimnames( pperm(asy_scl, q = matrix(s, nrow = 1)) ) )) # non-matrix in < 1.3-1 > stopifnot(identical(dimnames( pperm(asy_scl, q = matrix(s, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # non-matrix in < 1.3-1 > stopifnot( is.null( names( qperm(asy_scl, p = 0.75) ) )) > stopifnot(identical( names( qperm(asy_scl, p = c(s = 0.75)) ), "s" )) > stopifnot( is.null(dimnames( qperm(asy_scl, p = matrix(0.75, nrow = 1)) ) )) # non-matrix in < 1.3-1 > stopifnot(identical(dimnames( qperm(asy_scl, p = matrix(0.75, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # non-matrix in < 1.3-1 > stopifnot( is.null( names( rperm(asy_scl, n = 5) ) )) > stopifnot( is.null( names( pvalue(asy_scl, method = "global") ) )) > stopifnot( is.null( names( pvalue(asy_scl, method = "single-step") ) )) > stopifnot( is.null( names( pvalue(asy_scl, method = "step-down") ) )) > stopifnot( is.null( names( pvalue(asy_scl, method = "unadjusted") ) )) > ## stopifnot( is.null( names( midpvalue(asy_scl) ) )) # NA > ## stopifnot(identical( names(pvalue_interval(asy_scl) ), c("p_0", "p_1") )) # NA > ## stopifnot( is.null( names( size(asy_scl, alpha = 0.05) ) )) # NA > rm(asy_scl, s) > > ### Quadratic univariate > asy_qdr_u <- independence_test(y1 ~ x, distr = "asymptotic", teststat = "quadratic") > s <- statistic(asy_qdr_u) > > stopifnot(identical(dimnames( expectation(asy_qdr_u) ), list("A" , "") )) # non-matrix in < 1.4-0 > stopifnot(identical(dimnames( variance(asy_qdr_u) ), list("A" , "") )) # non-matrix in < 1.4-0 > stopifnot(identical(dimnames( covariance(asy_qdr_u) ), list("A" , "A") )) > stopifnot( is.null( names( statistic(asy_qdr_u, type = "test") ) )) > stopifnot(identical(dimnames( statistic(asy_qdr_u, type = "linear") ), list("A" , "") )) > stopifnot(identical(dimnames( statistic(asy_qdr_u, type = "centered") ), list("A" , "") )) > stopifnot(identical(dimnames( statistic(asy_qdr_u, type = "standardized") ), list("A" , "") )) > ## stopifnot( is.null( names( support(asy_qdr_u) ) )) # NA > stopifnot( is.null( names( dperm(asy_qdr_u, x = s) ) )) > stopifnot(identical( names( dperm(asy_qdr_u, x = c(s = s)) ), "s" )) > stopifnot( is.null(dimnames( dperm(asy_qdr_u, x = matrix(s, nrow = 1)) ) )) # non-matrix in < 1.3-1 > stopifnot(identical(dimnames( dperm(asy_qdr_u, x = matrix(s, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # non-matrix in < 1.3-1 > stopifnot( is.null( names( pperm(asy_qdr_u, q = s) ) )) > stopifnot(identical( names( pperm(asy_qdr_u, q = c(s = s)) ), "s" )) > stopifnot( is.null(dimnames( pperm(asy_qdr_u, q = matrix(s, nrow = 1)) ) )) # non-matrix in < 1.3-1 > stopifnot(identical(dimnames( pperm(asy_qdr_u, q = matrix(s, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # non-matrix in < 1.3-1 > stopifnot( is.null( names( qperm(asy_qdr_u, p = 0.75) ) )) > stopifnot(identical( names( qperm(asy_qdr_u, p = c(s = 0.75)) ), "s" )) > stopifnot( is.null(dimnames( qperm(asy_qdr_u, p = matrix(0.75, nrow = 1)) ) )) # non-matrix in < 1.3-1 > stopifnot(identical(dimnames( qperm(asy_qdr_u, p = matrix(0.75, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # non-matrix in < 1.3-1 > stopifnot( is.null( names( rperm(asy_qdr_u, n = 5) ) )) > stopifnot( is.null( names( pvalue(asy_qdr_u, method = "global") ) )) > stopifnot( is.null( names( pvalue(asy_qdr_u, method = "single-step") ) )) > stopifnot( is.null( names( pvalue(asy_qdr_u, method = "step-down") ) )) > stopifnot( is.null( names( pvalue(asy_qdr_u, method = "unadjusted") ) )) > ## stopifnot( is.null( names( midpvalue(asy_qdr_u) ) )) # NA > ## stopifnot(identical( names(pvalue_interval(asy_qdr_u) ), c("p_0", "p_1") )) # NA > ## stopifnot( is.null( names( size(asy_qdr_u, alpha = 0.05) ) )) # NA > rm(asy_qdr_u, s) > > ### Maximum > asy_mxm <- independence_test(y1 + y2 ~ x, distribution = "asymptotic", teststat = "maximum") > s <- statistic(asy_mxm) > > stopifnot(identical(dimnames( expectation(asy_mxm) ), list("A", c("y1", "y2")) )) # non-matrix in < 1.4-0 > stopifnot(identical(dimnames( variance(asy_mxm) ), list("A", c("y1", "y2")) )) # non-matrix in < 1.4-0 > stopifnot(identical(dimnames( covariance(asy_mxm) ), list(c("A:y1", "A:y2"), c("A:y1", "A:y2")) )) > stopifnot( is.null( names( statistic(asy_mxm, type = "test") ) )) > stopifnot(identical(dimnames( statistic(asy_mxm, type = "linear") ), list("A", c("y1", "y2")) )) > stopifnot(identical(dimnames( statistic(asy_mxm, type = "centered") ), list("A", c("y1", "y2")) )) > stopifnot(identical(dimnames( statistic(asy_mxm, type = "standardized") ), list("A", c("y1", "y2")) )) > ## stopifnot( is.null( names( support(asy_mxm) ) )) # NA > ## stopifnot( is.null( names( dperm(asy_mxm, x = s) ) )) # NA > ## stopifnot(identical( names( dperm(asy_mxm, x = c(s = s)) ), "s" )) # NA > ## stopifnot( is.null(dimnames( dperm(asy_mxm, x = matrix(s, nrow = 1)) ) )) # NA > ## stopifnot(identical(dimnames( dperm(asy_mxm, x = matrix(s, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # NA > stopifnot( is.null( names( pperm(asy_mxm, q = s) ) )) > stopifnot(identical( names( pperm(asy_mxm, q = c(s = s)) ), "s" )) > stopifnot( is.null(dimnames( pperm(asy_mxm, q = matrix(s, nrow = 1)) ) )) # non-matrix in < 1.3-1 > stopifnot(identical(dimnames( pperm(asy_mxm, q = matrix(s, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # non-matrix in < 1.3-1 > stopifnot( is.null( names( qperm(asy_mxm, p = 0.75) ) )) > stopifnot(identical( names( qperm(asy_mxm, p = c(s = 0.75)) ), "s" )) > stopifnot( is.null(dimnames( qperm(asy_mxm, p = matrix(0.75, nrow = 1)) ) )) # non-matrix in < 1.3-1 > stopifnot(identical(dimnames( qperm(asy_mxm, p = matrix(0.75, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # non-matrix in < 1.3-1 > ## stopifnot( is.null( names( rperm(asy_mxm, n = 5) ) )) # NA > stopifnot( is.null( names( pvalue(asy_mxm, method = "global") ) )) > stopifnot(identical(dimnames( pvalue(asy_mxm, method = "single-step") ), list("A", c("y1", "y2")) )) > stopifnot(identical(dimnames( pvalue(asy_mxm, method = "step-down") ), list("A", c("y1", "y2")) )) > stopifnot(identical(dimnames( pvalue(asy_mxm, method = "unadjusted") ), list("A", c("y1", "y2")) )) > ## stopifnot( is.null( names( midpvalue(asy_mxm) ) )) # NA > ## stopifnot(identical( names(pvalue_interval(asy_mxm) ), c("p_0", "p_1") )) # NA > ## stopifnot( is.null( names( size(asy_mxm, alpha = 0.05) ) )) # NA > rm(asy_mxm, s) > > ### Quadratic multivariate > asy_qdr_m <- independence_test(y1 + y2 ~ x, distr = "asymptotic", teststat = "quadratic") > s <- statistic(asy_qdr_m) > > stopifnot(identical(dimnames( expectation(asy_qdr_m) ), list("A", c("y1", "y2")) )) # non-matrix in < 1.4-0 > stopifnot(identical(dimnames( variance(asy_qdr_m) ), list("A", c("y1", "y2")) )) # non-matrix in < 1.4-0 > stopifnot(identical(dimnames( covariance(asy_qdr_m) ), list(c("A:y1", "A:y2"), c("A:y1", "A:y2")) )) > stopifnot( is.null( names( statistic(asy_qdr_m, type = "test") ) )) > stopifnot(identical(dimnames( statistic(asy_qdr_m, type = "linear") ), list("A", c("y1", "y2")) )) > stopifnot(identical(dimnames( statistic(asy_qdr_m, type = "centered") ), list("A", c("y1", "y2")) )) > stopifnot(identical(dimnames( statistic(asy_qdr_m, type = "standardized") ), list("A", c("y1", "y2")) )) > ## stopifnot( is.null( names( support(asy_qdr_m) ) )) # NA > stopifnot( is.null( names( dperm(asy_qdr_m, x = s) ) )) > stopifnot(identical( names( dperm(asy_qdr_m, x = c(s = s)) ), "s" )) > stopifnot( is.null(dimnames( dperm(asy_qdr_m, x = matrix(s, nrow = 1)) ) )) # non-matrix in < 1.3-1 > stopifnot(identical(dimnames( dperm(asy_qdr_m, x = matrix(s, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # non-matrix in < 1.3-1 > stopifnot( is.null( names( pperm(asy_qdr_m, q = s) ) )) > stopifnot(identical( names( pperm(asy_qdr_m, q = c(s = s)) ), "s" )) > stopifnot( is.null(dimnames( pperm(asy_qdr_m, q = matrix(s, nrow = 1)) ) )) # non-matrix in < 1.3-1 > stopifnot(identical(dimnames( pperm(asy_qdr_m, q = matrix(s, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # non-matrix in < 1.3-1 > stopifnot( is.null( names( qperm(asy_qdr_m, p = 0.75) ) )) > stopifnot(identical( names( qperm(asy_qdr_m, p = c(s = 0.75)) ), "s" )) > stopifnot( is.null(dimnames( qperm(asy_qdr_m, p = matrix(0.75, nrow = 1)) ) )) # non-matrix in < 1.3-1 > stopifnot(identical(dimnames( qperm(asy_qdr_m, p = matrix(0.75, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # non-matrix in < 1.3-1 > stopifnot( is.null( names( rperm(asy_qdr_m, n = 5) ) )) > stopifnot( is.null( names( pvalue(asy_qdr_m, method = "global") ) )) > stopifnot( is.null( names( pvalue(asy_qdr_m, method = "single-step") ) )) > stopifnot( is.null( names( pvalue(asy_qdr_m, method = "step-down") ) )) > stopifnot( is.null( names( pvalue(asy_qdr_m, method = "unadjusted") ) )) > ## stopifnot( is.null( names( midpvalue(asy_qdr_m) ) )) # NA > ## stopifnot(identical( names(pvalue_interval(asy_qdr_m) ), c("p_0", "p_1") )) # NA > ## stopifnot( is.null( names( size(asy_qdr_m, alpha = 0.05) ) )) # NA > rm(asy_qdr_m, s) > > > ### > ### Approximate > ### > > ### Scalar > app_scl <- independence_test(y1 ~ x, distr = "approximate") > s <- statistic(app_scl) > > stopifnot(identical(dimnames( expectation(app_scl) ), list("A", "") )) # non-matrix in < 1.4-0 > stopifnot(identical(dimnames( variance(app_scl) ), list("A", "") )) # non-matrix in < 1.4-0 > stopifnot(identical(dimnames( covariance(app_scl) ), list("A", "A") )) > stopifnot( is.null( names( statistic(app_scl, type = "test") ) )) # named in < 1.3-0 > stopifnot(identical(dimnames( statistic(app_scl, type = "linear") ), list("A", "") )) > stopifnot(identical(dimnames( statistic(app_scl, type = "centered") ), list("A", "") )) > stopifnot(identical(dimnames( statistic(app_scl, type = "standardized") ), list("A", "") )) > stopifnot( is.null( names( support(app_scl) ) )) > stopifnot( is.null( names( dperm(app_scl, x = s) ) )) > stopifnot(identical( names( dperm(app_scl, x = c(s = s)) ), "s" )) > stopifnot( is.null(dimnames( dperm(app_scl, x = matrix(s, nrow = 1)) ) )) # non-matrix in < 1.3-1 > stopifnot(identical(dimnames( dperm(app_scl, x = matrix(s, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # non-matrix in < 1.3-1 > stopifnot( is.null( names( pperm(app_scl, q = s) ) )) > stopifnot(identical( names( pperm(app_scl, q = c(s = s)) ), "s" )) > stopifnot( is.null(dimnames( pperm(app_scl, q = matrix(s, nrow = 1)) ) )) # non-matrix in < 1.3-1 > stopifnot(identical(dimnames( pperm(app_scl, q = matrix(s, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # non-matrix in < 1.3-1 > stopifnot( is.null( names( qperm(app_scl, p = 0.75) ) )) > stopifnot(identical( names( qperm(app_scl, p = c(s = 0.75)) ), "s" )) # unnamed in < 1.3-0 > stopifnot( is.null(dimnames( qperm(app_scl, p = matrix(0.75, nrow = 1)) ) )) # non-matrix in < 1.3-1 > stopifnot(identical(dimnames( qperm(app_scl, p = matrix(0.75, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # non-matrix in < 1.3-1 > stopifnot( is.null( names( rperm(app_scl, n = 5) ) )) > stopifnot( is.null( names( pvalue(app_scl, method = "global") ) )) > stopifnot( is.null( names( pvalue(app_scl, method = "single-step") ) )) > stopifnot( is.null( names( pvalue(app_scl, method = "step-down") ) )) > stopifnot( is.null( names( pvalue(app_scl, method = "unadjusted") ) )) > stopifnot( is.null( names( midpvalue(app_scl) ) )) > stopifnot(identical( names(pvalue_interval(app_scl) ), c("p_0", "p_1") )) > stopifnot( is.null( names( size(app_scl, alpha = 0.05) ) )) > rm(app_scl, s) > > ### Quadratic univariate > app_qdr_u <- independence_test(y1 ~ x, distr = "approximate", teststat = "quadratic") > s <- statistic(app_qdr_u) > > stopifnot(identical(dimnames( expectation(app_qdr_u) ), list("A", "") )) # non-matrix in < 1.4-0 > stopifnot(identical(dimnames( variance(app_qdr_u) ), list("A", "") )) # non-matrix in < 1.4-0 > stopifnot(identical(dimnames( covariance(app_qdr_u) ), list("A", "A") )) > stopifnot( is.null( names( statistic(app_qdr_u, type = "test") ) )) > stopifnot(identical(dimnames( statistic(app_qdr_u, type = "linear") ), list("A", "") )) > stopifnot(identical(dimnames( statistic(app_qdr_u, type = "centered") ), list("A", "") )) > stopifnot(identical(dimnames( statistic(app_qdr_u, type = "standardized") ), list("A", "") )) > stopifnot( is.null( names( support(app_qdr_u) ) )) > stopifnot( is.null( names( dperm(app_qdr_u, x = s) ) )) > stopifnot(identical( names( dperm(app_qdr_u, x = c(s = s)) ), "s" )) > stopifnot( is.null(dimnames( dperm(app_qdr_u, x = matrix(s, nrow = 1)) ) )) # non-matrix in < 1.3-1 > stopifnot(identical(dimnames( dperm(app_qdr_u, x = matrix(s, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # non-matrix in < 1.3-1 > stopifnot( is.null( names( pperm(app_qdr_u, q = s) ) )) > stopifnot(identical( names( pperm(app_qdr_u, q = c(s = s)) ), "s" )) > stopifnot( is.null(dimnames( pperm(app_qdr_u, q = matrix(s, nrow = 1)) ) )) # non-matrix in < 1.3-1 > stopifnot(identical(dimnames( pperm(app_qdr_u, q = matrix(s, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # non-matrix in < 1.3-1 > stopifnot( is.null( names( qperm(app_qdr_u, p = 0.75) ) )) > stopifnot(identical( names( qperm(app_qdr_u, p = c(s = 0.75)) ), "s" )) # unnamed in < 1.3-0 > stopifnot( is.null(dimnames( qperm(app_qdr_u, p = matrix(0.75, nrow = 1)) ) )) # non-matrix in < 1.3-1 > stopifnot(identical(dimnames( qperm(app_qdr_u, p = matrix(0.75, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # non-matrix in < 1.3-1 > stopifnot( is.null( names( rperm(app_qdr_u, n = 5) ) )) > stopifnot( is.null( names( pvalue(app_qdr_u, method = "global") ) )) > stopifnot( is.null( names( pvalue(app_qdr_u, method = "single-step") ) )) > stopifnot( is.null( names( pvalue(app_qdr_u, method = "step-down") ) )) > stopifnot( is.null( names( pvalue(app_qdr_u, method = "unadjusted") ) )) > stopifnot( is.null( names( midpvalue(app_qdr_u) ) )) > stopifnot(identical( names(pvalue_interval(app_qdr_u) ), c("p_0", "p_1") )) > stopifnot( is.null( names( size(app_qdr_u, alpha = 0.05) ) )) > rm(app_qdr_u, s) > > ### Maximum > app_mxm <- independence_test(y1 + y2 ~ x, distr = "approximate", teststat = "maximum") > s <- statistic(app_mxm) > > stopifnot(identical(dimnames( expectation(app_mxm) ), list("A", c("y1", "y2")) )) # non-matrix in < 1.4-0 > stopifnot(identical(dimnames( variance(app_mxm) ), list("A", c("y1", "y2")) )) # non-matrix in < 1.4-0 > stopifnot(identical(dimnames( covariance(app_mxm) ), list(c("A:y1", "A:y2"), c("A:y1", "A:y2")) )) > stopifnot( is.null( names( statistic(app_mxm, type = "test") ) )) > stopifnot(identical(dimnames( statistic(app_mxm, type = "linear") ), list("A", c("y1", "y2")) )) > stopifnot(identical(dimnames( statistic(app_mxm, type = "centered") ), list("A", c("y1", "y2")) )) > stopifnot(identical(dimnames( statistic(app_mxm, type = "standardized") ), list("A", c("y1", "y2")) )) > stopifnot( is.null( names( support(app_mxm) ) )) > stopifnot( is.null( names( dperm(app_mxm, x = s) ) )) > stopifnot(identical( names( dperm(app_mxm, x = c(s = s)) ), "s" )) > stopifnot( is.null(dimnames( dperm(app_mxm, x = matrix(s, nrow = 1)) ) )) # non-matrix in < 1.3-1 > stopifnot(identical(dimnames( dperm(app_mxm, x = matrix(s, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # non-matrix in < 1.3-1 > stopifnot( is.null( names( pperm(app_mxm, q = s) ) )) > stopifnot(identical( names( pperm(app_mxm, q = c(s = s)) ), "s" )) > stopifnot( is.null(dimnames( pperm(app_mxm, q = matrix(s, nrow = 1)) ) )) # non-matrix in < 1.3-1 > stopifnot(identical(dimnames( pperm(app_mxm, q = matrix(s, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # non-matrix in < 1.3-1 > stopifnot( is.null( names( qperm(app_mxm, p = 0.75) ) )) > stopifnot(identical( names( qperm(app_mxm, p = c(s = 0.75)) ), "s" )) # unnamed in < 1.3-0 > stopifnot( is.null(dimnames( qperm(app_mxm, p = matrix(0.75, nrow = 1)) ) )) # non-matrix in < 1.3-1 > stopifnot(identical(dimnames( qperm(app_mxm, p = matrix(0.75, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # non-matrix in < 1.3-1 > stopifnot( is.null( names( rperm(app_mxm, n = 5) ) )) > stopifnot( is.null( names( pvalue(app_mxm, method = "global") ) )) > stopifnot(identical(dimnames( pvalue(app_mxm, method = "single-step") ), list("A", c("y1", "y2")) )) > stopifnot(identical(dimnames( pvalue(app_mxm, method = "step-down") ), list("A", c("y1", "y2")) )) > stopifnot(identical(dimnames( pvalue(app_mxm, method = "unadjusted") ), list("A", c("y1", "y2")) )) > stopifnot( is.null( names( midpvalue(app_mxm) ) )) > stopifnot(identical( names(pvalue_interval(app_mxm) ), c("p_0", "p_1") )) > stopifnot( is.null( names( size(app_mxm, alpha = 0.05) ) )) > rm(app_mxm, s) > > ### Quadratic multivariate > app_qdr_m <- independence_test(y1 + y2 ~ x, distr = "approximate", teststat = "quadratic") > s <- statistic(app_qdr_m) > > stopifnot(identical(dimnames( expectation(app_qdr_m) ), list("A", c("y1", "y2")) )) # non-matrix in < 1.4-0 > stopifnot(identical(dimnames( variance(app_qdr_m) ), list("A", c("y1", "y2")) )) # non-matrix in < 1.4-0 > stopifnot(identical(dimnames( covariance(app_qdr_m) ), list(c("A:y1", "A:y2"), c("A:y1", "A:y2")) )) > stopifnot( is.null( names( statistic(app_qdr_m, type = "test") ) )) > stopifnot(identical(dimnames( statistic(app_qdr_m, type = "linear") ), list("A", c("y1", "y2")) )) > stopifnot(identical(dimnames( statistic(app_qdr_m, type = "centered") ), list("A", c("y1", "y2")) )) > stopifnot(identical(dimnames( statistic(app_qdr_m, type = "standardized") ), list("A", c("y1", "y2")) )) > stopifnot( is.null( names( support(app_qdr_m) ) )) > stopifnot( is.null( names( dperm(app_qdr_m, x = s) ) )) > stopifnot(identical( names( dperm(app_qdr_m, x = c(s = s)) ), "s" )) > stopifnot( is.null(dimnames( dperm(app_qdr_m, x = matrix(s, nrow = 1)) ) )) # non-matrix in < 1.3-1 > stopifnot(identical(dimnames( dperm(app_qdr_m, x = matrix(s, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # non-matrix in < 1.3-1 > stopifnot( is.null( names( pperm(app_qdr_m, q = s) ) )) > stopifnot(identical( names( pperm(app_qdr_m, q = c(s = s)) ), "s" )) > stopifnot( is.null(dimnames( pperm(app_qdr_m, q = matrix(s, nrow = 1)) ) )) # non-matrix in < 1.3-1 > stopifnot(identical(dimnames( pperm(app_qdr_m, q = matrix(s, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # non-matrix in < 1.3-1 > stopifnot( is.null( names( qperm(app_qdr_m, p = 0.75) ) )) > stopifnot(identical( names( qperm(app_qdr_m, p = c(s = 0.75)) ), "s" )) # unnamed in < 1.3-0 > stopifnot( is.null(dimnames( qperm(app_qdr_m, p = matrix(0.75, nrow = 1)) ) )) # non-matrix in < 1.3-1 > stopifnot(identical(dimnames( qperm(app_qdr_m, p = matrix(0.75, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # non-matrix in < 1.3-1 > stopifnot( is.null( names( rperm(app_qdr_m, n = 5) ) )) > stopifnot( is.null( names( pvalue(app_qdr_m, method = "global") ) )) > stopifnot( is.null( names( pvalue(app_qdr_m, method = "single-step") ) )) > stopifnot( is.null( names( pvalue(app_qdr_m, method = "step-down") ) )) > stopifnot( is.null( names( pvalue(app_qdr_m, method = "unadjusted") ) )) > stopifnot( is.null( names( midpvalue(app_qdr_m) ) )) > stopifnot(identical( names(pvalue_interval(app_qdr_m) ), c("p_0", "p_1") )) > stopifnot( is.null( names( size(app_qdr_m, alpha = 0.05) ) )) > rm(app_qdr_m, s) > > > ### > ### Exact using shift algorithm > ### > > ### Scalar > shf_scl <- independence_test(y1 ~ x, distr = exact(algo = "shift")) > s <- statistic(shf_scl) > > stopifnot(identical(dimnames( expectation(shf_scl) ), list("A", "") )) # non-matrix in < 1.4-0 > stopifnot(identical(dimnames( variance(shf_scl) ), list("A", "") )) # non-matrix in < 1.4-0 > stopifnot(identical(dimnames( covariance(shf_scl) ), list("A", "A") )) > stopifnot( is.null( names( statistic(shf_scl, type = "test") ) )) # named in < 1.3-0 > stopifnot(identical(dimnames( statistic(shf_scl, type = "linear") ), list("A", "") )) > stopifnot(identical(dimnames( statistic(shf_scl, type = "centered") ), list("A", "") )) > stopifnot(identical(dimnames( statistic(shf_scl, type = "standardized") ), list("A", "") )) > stopifnot( is.null( names( support(shf_scl) ) )) > stopifnot( is.null( names( dperm(shf_scl, x = s) ) )) > stopifnot(identical( names( dperm(shf_scl, x = c(s = s)) ), "s" )) > stopifnot( is.null(dimnames( dperm(shf_scl, x = matrix(s, nrow = 1)) ) )) # non-matrix in < 1.3-1 > stopifnot(identical(dimnames( dperm(shf_scl, x = matrix(s, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # non-matrix in < 1.3-1 > stopifnot( is.null( names( pperm(shf_scl, q = s) ) )) > stopifnot(identical( names( pperm(shf_scl, q = c(s = s)) ), "s" )) > stopifnot( is.null(dimnames( pperm(shf_scl, q = matrix(s, nrow = 1)) ) )) # non-matrix in < 1.3-1 > stopifnot(identical(dimnames( pperm(shf_scl, q = matrix(s, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # non-matrix in < 1.3-1 > stopifnot( is.null( names( qperm(shf_scl, p = 0.75) ) )) > stopifnot(identical( names( qperm(shf_scl, p = c(s = 0.75)) ), "s" )) > stopifnot( is.null(dimnames( qperm(shf_scl, p = matrix(0.75, nrow = 1)) ) )) # non-matrix in < 1.3-1 > stopifnot(identical(dimnames( qperm(shf_scl, p = matrix(0.75, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # non-matrix in < 1.3-1 > stopifnot( is.null( names( rperm(shf_scl, n = 5) ) )) > stopifnot( is.null( names( pvalue(shf_scl, method = "global") ) )) > stopifnot( is.null( names( pvalue(shf_scl, method = "single-step") ) )) > stopifnot( is.null( names( pvalue(shf_scl, method = "step-down") ) )) > stopifnot( is.null( names( pvalue(shf_scl, method = "unadjusted") ) )) > stopifnot( is.null( names( midpvalue(shf_scl) ) )) > stopifnot(identical( names(pvalue_interval(shf_scl) ), c("p_0", "p_1") )) > stopifnot( is.null( names( size(shf_scl, alpha = 0.05) ) )) > rm(shf_scl, s) > > ### Quadratic univarite > shf_qdr_u <- independence_test(y1 ~ x, distr = exact(algo = "shift"), teststat = "quadratic") > s <- statistic(shf_qdr_u) > > stopifnot(identical(dimnames( expectation(shf_qdr_u) ), list("A", "") )) # non-matrix in < 1.4-0 > stopifnot(identical(dimnames( variance(shf_qdr_u) ), list("A", "") )) # non-matrix in < 1.4-0 > stopifnot(identical(dimnames( covariance(shf_qdr_u) ), list("A", "A") )) > stopifnot( is.null( names( statistic(shf_qdr_u, type = "test") ) )) > stopifnot(identical(dimnames( statistic(shf_qdr_u, type = "linear") ), list("A", "") )) > stopifnot(identical(dimnames( statistic(shf_qdr_u, type = "centered") ), list("A", "") )) > stopifnot(identical(dimnames( statistic(shf_qdr_u, type = "standardized") ), list("A", "") )) > stopifnot( is.null( names( support(shf_qdr_u) ) )) > stopifnot( is.null( names( dperm(shf_qdr_u, x = s) ) )) > stopifnot(identical( names( dperm(shf_qdr_u, x = c(s = s)) ), "s" )) > stopifnot( is.null(dimnames( dperm(shf_qdr_u, x = matrix(s, nrow = 1)) ) )) # non-matrix in < 1.3-1 > stopifnot(identical(dimnames( dperm(shf_qdr_u, x = matrix(s, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # non-matrix in < 1.3-1 > stopifnot( is.null( names( pperm(shf_qdr_u, q = s) ) )) > stopifnot(identical( names( pperm(shf_qdr_u, q = c(s = s)) ), "s" )) > stopifnot( is.null(dimnames( pperm(shf_qdr_u, q = matrix(s, nrow = 1)) ) )) # non-matrix in < 1.3-1 > stopifnot(identical(dimnames( pperm(shf_qdr_u, q = matrix(s, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # non-matrix in < 1.3-1 > stopifnot( is.null( names( qperm(shf_qdr_u, p = 0.75) ) )) > stopifnot(identical( names( qperm(shf_qdr_u, p = c(s = 0.75)) ), "s" )) > stopifnot( is.null(dimnames( qperm(shf_qdr_u, p = matrix(0.75, nrow = 1)) ) )) # non-matrix in < 1.3-1 > stopifnot(identical(dimnames( qperm(shf_qdr_u, p = matrix(0.75, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # non-matrix in < 1.3-1 > stopifnot( is.null( names( rperm(shf_qdr_u, n = 5) ) )) > stopifnot( is.null( names( pvalue(shf_qdr_u, method = "global") ) )) > stopifnot( is.null( names( pvalue(shf_qdr_u, method = "single-step") ) )) > stopifnot( is.null( names( pvalue(shf_qdr_u, method = "step-down") ) )) > stopifnot( is.null( names( pvalue(shf_qdr_u, method = "unadjusted") ) )) > stopifnot( is.null( names( midpvalue(shf_qdr_u) ) )) > stopifnot(identical( names(pvalue_interval(shf_qdr_u) ), c("p_0", "p_1") )) > stopifnot( is.null( names( size(shf_qdr_u, alpha = 0.05) ) )) > rm(shf_qdr_u, s) > > > ### > ### Exact using split-up algorithm > ### > > ### Scalar > spl_scl <- independence_test(y1 ~ x, distr = exact(algo = "split")) > s <- statistic(spl_scl) > > stopifnot(identical(dimnames( expectation(spl_scl) ), list("A", "") )) # non-matrix in < 1.4-0 > stopifnot(identical(dimnames( variance(spl_scl) ), list("A", "") )) # non-matrix in < 1.4-0 > stopifnot(identical(dimnames( covariance(spl_scl) ), list("A", "A") )) > stopifnot( is.null( names( statistic(spl_scl, type = "test") ) )) # named in < 1.3-0 > stopifnot(identical(dimnames( statistic(spl_scl, type = "linear") ), list("A", "") )) > stopifnot(identical(dimnames( statistic(spl_scl, type = "centered") ), list("A", "") )) > stopifnot(identical(dimnames( statistic(spl_scl, type = "standardized") ), list("A", "") )) > ## stopifnot( is.null( names( support(spl_scl) ) )) # NA > ## stopifnot( is.null( names( dperm(spl_scl, x = s) ) )) # NA > ## stopifnot(identical( names( dperm(spl_scl, x = c(s = s)) ), "s" )) # NA > ## stopifnot( is.null(dimnames( dperm(spl_scl, x = matrix(s, nrow = 1)) ) )) # NA > ## stopifnot(identical(dimnames( dperm(spl_scl, x = matrix(s, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # NA > stopifnot( is.null( names( pperm(spl_scl, q = s) ) )) > stopifnot(identical( names( pperm(spl_scl, q = c(s = s)) ), "s" )) > stopifnot( is.null(dimnames( pperm(spl_scl, q = matrix(s, nrow = 1)) ) )) # non-matrix in < 1.3-1 > stopifnot(identical(dimnames( pperm(spl_scl, q = matrix(s, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # non-matrix in < 1.3-1 > stopifnot( is.null( names( qperm(spl_scl, p = 0.75) ) )) > stopifnot(identical( names( qperm(spl_scl, p = c(s = 0.75)) ), "s" )) > stopifnot( is.null(dimnames( qperm(spl_scl, p = matrix(0.75, nrow = 1)) ) )) # non-matrix in < 1.3-1 > stopifnot(identical(dimnames( qperm(spl_scl, p = matrix(0.75, nrow = 1, dimnames = list("s", NULL))) ), list("s", NULL) )) # non-matrix in < 1.3-1 > stopifnot( is.null( names( rperm(spl_scl, n = 5) ) )) > stopifnot( is.null( names( pvalue(spl_scl, method = "global") ) )) > stopifnot( is.null( names( pvalue(spl_scl, method = "single-step") ) )) > stopifnot( is.null( names( pvalue(spl_scl, method = "step-down") ) )) > stopifnot( is.null( names( pvalue(spl_scl, method = "unadjusted") ) )) > ## stopifnot( is.null( names( midpvalue(spl_scl) ) )) # NA > ## stopifnot(identical( names(pvalue_interval(spl_scl) ), c("p_0", "p_1") )) # NA > ## stopifnot( is.null( names( size(spl_scl, alpha = 0.05) ) )) # NA > rm(spl_scl, s) > > proc.time() user system elapsed 2.15 0.15 2.29
Generated by dwww version 1.15 on Sun Jun 16 11:54:57 CEST 2024.