Title: | Two and Three Dimensional Kolmogorov-Smirnov Two-Sample Tests |
---|---|
Description: | The original definition of the two and three dimensional Kolmogorov-Smirnov two-sample test statistics given by Peacock (1983) is implemented. Two R-functions: peacock2 and peacock3, are provided to compute the test statistics in two and three dimensional spaces, respectively. Note the Peacock test is different from the Fasano and Franceschini test (1987). The latter is a variant of the Peacock test. |
Authors: | Yuanhui Xiao |
Maintainer: | Yuanhui Xiao <[email protected]> |
License: | GPL-2 |
Version: | 1.0 |
Built: | 2024-12-08 06:58:55 UTC |
Source: | CRAN |
Two R-functions: peacock2 and peacock3 are provided to compute the two dimensional and three dimensional KS two-sample tests, respectively. The famous KS two sample test was generalized to multidimensional spaces by Peacock (1983). Hence, it is also called the Peacock test. The Peacock is different from the widely used Fasano and Franceschini test (1987). The latter is a variant of the KS test, invented to alleviate the computational intensity of the former. The two R-functions implement the original definition of the KS test given by Peacock (1983).
Package: | Peacock.test |
Type: | Package |
Version: | 1.0 |
Date: | 2016-07-13 |
License: | GPL-2 |
The two functions: peacock2 and peacock3, provided in this package are self-explanatory and their usage is straightforward.
Yuanhui Xiao
Maintainer: Yuanhui Xiao <xiao\[email protected]> ~~ The author and/or maintainer of the package ~~
Fasano, G., Franceschini, A. (1987)<DOI:10.1093/mnras/225.1.155>. A multidimensional version of the Kolmogorov-Smirnov test. Monthly Notices of the Royal Astronomical Society 225:155-170.
Peacock J.A. (1983) <DOI:10.1093/mnras/202.3.615>. Two-dimensional goodness-of-fit testing in astronomy. Monthly Notices of the Royal Astronomical Society 202:615-627.
Xiao Y(2016). A fast algorithm for two-dimensional Kolmogorov-Smirnov two-sample tests, Journal Computational Statistics and Data Analysis, under revision.
# two-dimensional case x2 <- matrix(rnorm(12, 0, 1), ncol=2) y2 <- matrix(rnorm(16, 0, 1), ncol=2) ks2 <- peacock2(x2, y2) ks2 # three-dimensional case x3 <- matrix(rnorm(12, 0, 1), ncol=3) y3 <- matrix(rnorm(18, 0, 1), ncol=3) ks3 <- peacock3(x3, y3) ks3
# two-dimensional case x2 <- matrix(rnorm(12, 0, 1), ncol=2) y2 <- matrix(rnorm(16, 0, 1), ncol=2) ks2 <- peacock2(x2, y2) ks2 # three-dimensional case x3 <- matrix(rnorm(12, 0, 1), ncol=3) y3 <- matrix(rnorm(18, 0, 1), ncol=3) ks3 <- peacock3(x3, y3) ks3
This function returns the greatest common divisor of two integers
gcd(x, y)
gcd(x, y)
x |
a nonnegative integer |
y |
a nonnegative integer |
The function utilizes the Euclidean algorithm to compute the greatest common divisor of two integers
An integer, which is the greatest common divisor. If both arguments are zero, then the returned value is 1.
Yuanhui Xiao
x <- 4 y <- 6 d <- gcd(x, y) d
x <- 4 y <- 6 d <- gcd(x, y) d
This function implements the original definition of the two-dimensional Kolmogorov-Smirnov test by Peacock (1983). This test is not the widely used Fasano-Franceschini test (1987). The latter is a variant the Peacock test.
peacock2(x, y)
peacock2(x, y)
x |
x is the object representing the first sample. x should be able to be converted a matrix, where each row represents a sample point. If the object cannot be converted to a matrix, the function will stop and throw an error message. Pleaste note that only the first two columns of the matrix will be used, and the rest columns are just ingored. |
y |
Similar to x, y is the object representing the second sample. |
the value of the test statistic
Yuanhui Xiao
x <- matrix(rnorm(12, 0, 1), ncol=2) y <- matrix(rnorm(16, 0, 1), ncol=2) ks <- peacock2(x, y) ks
x <- matrix(rnorm(12, 0, 1), ncol=2) y <- matrix(rnorm(16, 0, 1), ncol=2) ks <- peacock2(x, y) ks
This function implements the original definition of the three-dimensional Kolmogorov-Smirnov test by Peacock (1983). This test is not the widely used Fasano-Franceschini test (1987). The latter is a variant the Peacock test.
peacock3(x, y)
peacock3(x, y)
x |
x is the object representing the first sample. x should be able to be converted a matrix, where each row represents a sample point. If the object cannot be converted to a matrix, the function will stop and throw an error message. Pleaste note that only the first three columns of the matrix will be used, and the rest columns are just ingored. |
y |
Similar to x, y is the object representing the second sample. |
the value of the test statistic
Yuanhui Xiao
x <- matrix(rnorm(12, 0, 1), ncol=3) y <- matrix(rnorm(18, 0, 1), ncol=3) ks <- peacock3(x, y) ks
x <- matrix(rnorm(12, 0, 1), ncol=3) y <- matrix(rnorm(18, 0, 1), ncol=3) ks <- peacock3(x, y) ks