moezipf: The Marshal-Olkin Extended Zipf Distribution (MOEZipf).

In zipfextR: Zipf Extended Distributions

Description

Probability mass function, cumulative distribution function, quantile function and random number generation for the MOEZipf distribution with parameters α and β. The support of the MOEZipf distribution are the strictly positive integer numbers large or equal than one.

Usage

dmoezipf(x, alpha, beta, log = FALSE)

pmoezipf(q, alpha, beta, log.p = FALSE, lower.tail = TRUE)

qmoezipf(p, alpha, beta, log.p = FALSE, lower.tail = TRUE)

rmoezipf(n, alpha, beta)

Arguments

x, q	Vector of positive integer values.
alpha	Value of the α parameter (α > 1 ).
beta	Value of the β parameter (β > 0 ).
log, log.p	Logical; if TRUE, probabilities p are given as log(p).
lower.tail	Logical; if TRUE (default), probabilities are P[X ≤q x], otherwise, P[X > x].
p	Vector of probabilities.
n	Number of random values to return.

Details

The probability mass function at a positive integer value x of the MOEZipf distribution with parameters α and β is computed as follows:

p(x | α, β) = \frac{x^{-α} β ζ(α) }{[ζ(α) - \bar{β} ζ (α, x)] [ζ (α) - \bar{β} ζ (α, x + 1)]},\, x = 1,2,...,\, α > 1, β > 0,

where ζ(α) is the Riemann-zeta function at α, ζ(α, x) is the Hurtwitz zeta function with arguments α and x, and \bar{β} = 1 - β.

The cumulative distribution function, at a given positive integer value x, is computed as F(x) = 1 - S(x), where the survival function S(x) is equal to:

S(x) = \frac{β\, ζ(α, x + 1)}{ζ(α) - \bar{β}\,ζ(α, x + 1)},\, x = 1, 2, ..

The quantile of the MOEZipf(α, β) distribution of a given probability value p is equal to the quantile of the Zipf(α) distribution at the value:

p\prime = \frac{p\,β}{1 + p\,(β - 1)}

The quantiles of the Zipf(α) distribution are computed by means of the tolerance package.

To generate random data from a MOEZipf one applies the quantile function over n values randomly generated from an Uniform distribution in the interval (0, 1).

Value

dmoezipf gives the probability mass function, pmoezipf gives the cumulative distribution function, qmoezipf gives the quantile function, and rmoezipf generates random values from a MOEZipf distribution.

References

Casellas, A. (2013) La distribució Zipf Estesa segons la transformació Marshall-Olkin. Universitat Politécnica de Catalunya.

Devroye L. (1986) Non-Uniform Random Variate Generation. Springer, New York, NY.

Duarte-López, A., Prat-Pérez, A., & Pérez-Casany, M. (2015). Using the Marshall-Olkin Extended Zipf Distribution in Graph Generation. European Conference on Parallel Processing, pp. 493-502, Springer International Publishing.

Pérez-Casany, M. and Casellas, A. (2013) Marshall-Olkin Extended Zipf Distribution. arXiv preprint arXiv:1304.4540.

Young, D. S. (2010). Tolerance: an R package for estimating tolerance intervals. Journal of Statistical Software, 36(5), 1-39.

Examples

dmoezipf(1:10, 2.5, 1.3)
pmoezipf(1:10, 2.5, 1.3)
qmoezipf(0.56, 2.5, 1.3)
rmoezipf(10, 2.5, 1.3)

zipfextR documentation built on July 8, 2020, 6:23 p.m.