kuiper_stat {twosamples} | R Documentation |
Performs a permutation based two sample test using the Kuiper test statistic (kuiper_stat).
kuiper_stat(a, b, power = 1) kuiper_test(a, b, nboots = 2000, p = default.p)
a |
a vector of numbers |
b |
a vector of numbers |
power |
power to raise test stat to |
nboots |
Number of bootstrap iterations |
p |
power to raise test stat to |
The Kuiper test compares two ECDFs by looking at the maximum positive and negative difference between them. Formally – if E is the ECDF of sample 1 and F is the ECDF of sample 2, then KUIPER = max_x E(x)-F(x) + max_x F(x)-E(x). The test p-value is calculated by randomly resampling two samples of the same size using the combined sample.
Output is a length 2 Vector with test stat and p-value in that order. That vector has 3 attributes – the sample sizes of each sample, and the number of bootstraps performed for the pvalue.
kuiper_stat
: Kuiper Test statistic
kuiper_test
: Permutation based two sample Kuiper test
vec1 = rnorm(20) vec2 = rnorm(20,4) kuiper_test(vec1,vec2)