ExtDist (version 0.7-2)

Uniform: The Uniform Distribution.

Description

Density, distribution, quantile, random number generation and parameter estimation functions for the uniform distribution on the interval \([a,b]\). Parameter estimation can be based on an unweighted i.i.d. sample only and can be performed analytically or numerically.

Usage

dUniform(x, a = 0, b = 1, params = list(a, b), ...)
pUniform(q, a = 0, b = 1, params = list(a, b), ...)
qUniform(p, a = 0, b = 1, params = list(a, b), ...)
rUniform(n, a = 0, b = 1, params = list(a, b), ...)
eUniform(X, w, method = c("analytic.MLE", "moments", "numerical.MLE"), ...)
lUniform(X, w, a = 0, b = 1, params = list(a, b), logL = TRUE, ...)

Value

dUniform gives the density, pUniform the distribution function, qUniform the quantile function, rUniform generates random deviates, and eUniform estimates the parameters. lUniform provides the log-likelihood function.

Arguments

x, q

A vector of quantiles.

a, b

Boundary parameters.

params

A list that includes all named parameters.

...

Additional parameters.

p

A vector of probabilities.

n

Number of observations.

X

Sample observations.

w

An optional vector of sample weights.

method

Parameter estimation method.

logL

logical;if TRUE, lUniform gives the log-likelihood, otherwise the likelihood is given.

Author

Haizhen Wu and A. Jonathan R. Godfrey.
Updates and bugfixes by Sarah Pirikahu.

Details

If a or b are not specified they assume the default values of 0 and 1, respectively.

The dUniform(), pUniform(), qUniform(),and rUniform() functions serve as wrappers of the standard dunif, punif, qunif, and runif functions in the stats package. They allow for the parameters to be declared not only as individual numerical values, but also as a list so parameter estimation can be carried out.

The uniform distribution has probability density function $$p_x(x) = 1/(b-a)$$ for \(a \le x \le b\). The analytic maximum likelihood parameter estimates are as given by Engineering Statistics Handbook. The method of moments parameter estimation option is also avaliable and the estimates are as given by Forbes et.al (2011), p.179.

The log-likelihood function for the uniform distribution is given by $$l(a,b|x) = -n log(b-a)$$

References

Johnson, N. L., Kotz, S. and Balakrishnan, N. (1995) Continuous Univariate Distributions, volume 2, chapter 26, Wiley, New York.

Engineering Statistics Handbook

Forbes, C. Evans, M. Hastings, N. & Peacock, B. (2011) Statistical Distributions, 4th Ed, chapter 40, Wiley, New Jersey.

See Also

ExtDist for other standard distributions.

Examples

# Parameter estimation for a distribution with known shape parameters
X <- rUniform(n=500, a=0, b=1)
est.par <- eUniform(X, method="analytic.MLE"); est.par
plot(est.par)

# Histogram and fitted density
den.x <- seq(min(X),max(X),length=100)
den.y <- dUniform(den.x,a=est.par$a,b=est.par$b)
hist(X, breaks=10, probability=TRUE, ylim = c(0,1.2*max(den.y)))
lines(den.x, den.y, col="blue")  # Original data
lines(density(X), lty=2)         # Fitted curve

# Extracting boundary parameters
est.par[attributes(est.par)$par.type=="boundary"]

# log-likelihood
lUniform(X,param = est.par)

# Example of parameter estimation for a distribution with
# unknown parameters currently been sought after.