R Under development (unstable) (2016-11-01 r71616) -- "Unsuffered Consequences"
Copyright (C) 2016 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 the correlation problem, i.e.,
> ### testing the independence of two numeric variables
> ### 'x' and 'y' (possibly blocked)
>
> set.seed(290875)
> library("coin")
Loading required package: survival
> isequal <- coin:::isequal
> options(useFancyQuotes = FALSE)
>
> ### generate data
> dat <- data.frame(x = rnorm(100), y = rnorm(100), block = gl(10, 10))
>
> ### not really the same, T = (rank(x) - rank(y))^2 is used here
> cor.test(~ x + y, data = dat, method = "spearman")$p.value
[1] 0.8474266
> cor.test(~ x + y, data = dat, alternative = "less", method = "spearman")$p.value
[1] 0.5763332
> cor.test(~ x + y, data = dat, alternative = "greater", method = "spearman")$p.value
[1] 0.4237133
>
> ### without blocks
> pvalue(spearman_test(y ~ x, data = dat))
[1] 0.846517
> pvalue(spearman_test(x ~ y, data = dat))
[1] 0.846517
> pvalue(spearman_test( ~ y + x, data = dat))
[1] 0.846517
> pvalue(spearman_test( ~ x + y, data = dat))
[1] 0.846517
>
> pvalue(fisyat_test(y ~ x, data = dat))
[1] 0.9155477
> pvalue(fisyat_test(x ~ y, data = dat))
[1] 0.9155477
> pvalue(fisyat_test( ~ y + x, data = dat))
[1] 0.9155477
> pvalue(fisyat_test( ~ x + y, data = dat))
[1] 0.9155477
>
> pvalue(quadrant_test(y ~ x, data = dat))
[1] 0.6906339
> pvalue(quadrant_test(x ~ y, data = dat))
[1] 0.6906339
> pvalue(quadrant_test( ~ y + x, data = dat))
[1] 0.6906339
> pvalue(quadrant_test( ~ x + y, data = dat))
[1] 0.6906339
>
> pvalue(koziol_test(y ~ x, data = dat))
[1] 0.7825737
> pvalue(koziol_test(x ~ y, data = dat))
[1] 0.7825737
> pvalue(koziol_test( ~ y + x, data = dat))
[1] 0.7825737
> pvalue(koziol_test( ~ x + y, data = dat))
[1] 0.7825737
>
> ### with blocks
> pvalue(spearman_test(y ~ x | block, data = dat))
[1] 0.8901013
> pvalue(spearman_test(x ~ y | block, data = dat))
[1] 0.8901013
> pvalue(spearman_test( ~ y + x | block, data = dat))
[1] 0.8901013
> pvalue(spearman_test( ~ x + y | block, data = dat))
[1] 0.8901013
>
> pvalue(fisyat_test(y ~ x | block, data = dat))
[1] 0.9061762
> pvalue(fisyat_test(x ~ y | block, data = dat))
[1] 0.9061762
> pvalue(fisyat_test( ~ y + x | block, data = dat))
[1] 0.9061762
> pvalue(fisyat_test( ~ x + y | block, data = dat))
[1] 0.9061762
>
> pvalue(quadrant_test(y ~ x | block, data = dat))
[1] 0.7489036
> pvalue(quadrant_test(x ~ y | block, data = dat))
[1] 0.7489036
> pvalue(quadrant_test( ~ y + x | block, data = dat))
[1] 0.7489036
> pvalue(quadrant_test( ~ x + y | block, data = dat))
[1] 0.7489036
>
> pvalue(koziol_test(y ~ x | block, data = dat))
[1] 0.8408785
> pvalue(koziol_test(x ~ y | block, data = dat))
[1] 0.8408785
> pvalue(koziol_test( ~ y + x | block, data = dat))
[1] 0.8408785
> pvalue(koziol_test( ~ x + y | block, data = dat))
[1] 0.8408785
>
> ### sanity checks, those should be errors
> dat <- data.frame(x = gl(2, 50), y = rnorm(100), block = rnorm(100))
>
> try(spearman_test(y ~ x, data = dat))
Error in check(object) :
'object' does not represent a univariate correlation problem
> try(spearman_test(y ~ x | block, data = dat))
Error in .local(.Object, ...) : 'block' is not a factor
>
> try(fisyat_test(y ~ x, data = dat))
Error in check(object) :
'object' does not represent a univariate correlation problem
> try(fisyat_test(y ~ x | block, data = dat))
Error in .local(.Object, ...) : 'block' is not a factor
>
> try(quadrant_test(y ~ x, data = dat))
Error in check(object) :
'object' does not represent a univariate correlation problem
> try(quadrant_test(y ~ x | block, data = dat))
Error in .local(.Object, ...) : 'block' is not a factor
>
> try(koziol_test(y ~ x, data = dat))
Error in check(object) :
'object' does not represent a univariate correlation problem
> try(koziol_test(y ~ x | block, data = dat))
Error in .local(.Object, ...) : 'block' is not a factor
>
> proc.time()
user system elapsed
1.07 0.07 1.14
Generated by dwww version 1.15 on Sun Jun 16 13:10:27 CEST 2024.