distr.test {digitTests}

Test of Digits against a Reference Distribution

Description

This function extracts and performs a test of the distribution of (leading) digits in a vector against a reference distribution. By default, the distribution of leading digits is checked against Benford's law.

Usage

distr.test(x, check = 'first', reference = 'benford')

Arguments

x

a numeric vector.

check

location of the digits to analyze. Can be first, firsttwo, or last.

reference

which character string given the reference distribution for the digits, or a vector of probabilities for each digit. Can be benford for Benford's law, uniform for the uniform distribution. An error is given if any entry of reference is negative. Probabilities that do not sum to one are normalized.

Details

Benford's law is defined as p(d) = log10(1/d). The uniform distribution is defined as p(d) = 1/d.

Value

An object of class dt.distr containing:

observed

the observed counts.

expected

the expected counts under the null hypothesis.

n

the number of observations in x.

statistic

the value the chi-squared test statistic.

parameter

the degrees of freedom of the approximate chi-squared distribution of the test statistic.

p.value

the p-value for the test.

check

checked digits.

digits

vector of digits.

reference

reference distribution

data.name

a character string giving the name(s) of the data.

Author(s)

Koen Derks, [e-mail]

References

Benford, F. (1938). The law of anomalous numbers. In Proceedings of the American Philosophical Society, 551-572.

See Also

distr.btest rv.test

Examples

set.seed(1)
x <- rnorm(100)

# Digit analysis against Benford's law
distr.test(x, check = 'first', reference = 'benford')

# Digit analysis against custom distribution
distr.test(x, check = 'last', reference = rep(1/9, 9))