kuiper_stat {twosamples}R Documentation

Kuiper Test

Description

Performs a permutation based two sample test using the Kuiper test statistic (kuiper_stat).

Usage

kuiper_stat(a, b, power = 1)

kuiper_test(a, b, nboots = 2000, p = default.p)

Arguments

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

Details

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.

Value

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.

Functions

Examples

vec1 = rnorm(20)
vec2 = rnorm(20,4)
kuiper_test(vec1,vec2)

[Package twosamples version 1.0.0 Index]