foldnormal: Folded Normal Distribution Family Function

Description

Fits a (generalized) folded (univariate) normal distribution.

Usage

foldnormal(lmean = "identitylink", lsd = "loglink", imean = NULL, isd = NULL,
           a1 = 1, a2 = 1, nsimEIM = 500, imethod = 1, zero = NULL)

Arguments

lmean, lsd

Link functions for the mean and standard deviation parameters of the usual univariate normal distribution. They are $μ$ and $σ$ respectively. See Links for more choices.

imean, isd

Optional initial values for $μ$ and $σ$ . A NULL means a value is computed internally. See CommonVGAMffArguments.

a1, a2

Positive weights, called $a_{1}$ and $a_{2}$ below. Each must be of length 1.

nsimEIM, imethod, zero

See CommonVGAMffArguments.

Value

An object of class "vglmff" (see vglmff-class). The object is used by modelling functions such as vglm and vgam.

Warning

Under- or over-flow may occur if the data is ill-conditioned. It is recommended that several different initial values be used to help avoid local solutions.

Details

If a random variable has an ordinary univariate normal distribution then the absolute value of that random variable has an ordinary folded normal distribution. That is, the sign has not been recorded; only the magnitude has been measured.

More generally, suppose $X$ is normal with mean mean and standard deviation sd. Let $Y = max (a_{1} X, - a_{2} X)$ where $a_{1}$ and $a_{2}$ are positive weights. This means that $Y = a_{1} X$ for $X > 0$ , and $Y = a_{2} X$ for $X < 0$ . Then $Y$ is said to have a generalized folded normal distribution. The ordinary folded normal distribution corresponds to the special case $a_{1} = a_{2} = 1$ .

The probability density function of the ordinary folded normal distribution can be written dnorm(y, mean, sd) + dnorm(y, -mean, sd) for $y \geq 0$ . By default, mean and log(sd) are the linear/additive predictors. Having mean=0 and sd=1 results in the half-normal distribution. The mean of an ordinary folded normal distribution is $E (Y) = σ \sqrt{2 / π} \exp (- μ^{2} / (2 σ^{2})) + μ [1 - 2 Φ (- μ / σ)]$ and these are returned as the fitted values. Here, $Φ ()$ is the cumulative distribution function of a standard normal (pnorm).

References

Lin, P. C. (2005). Application of the generalized folded-normal distribution to the process capability measures. International Journal of Advanced Manufacturing Technology, 26, 825--830.

Johnson, N. L. (1962). The folded normal distribution: accuracy of estimation by maximum likelihood. Technometrics, 4, 249--256.

Examples

Run this code

# NOT RUN {
 m <-  2; SD <- exp(1)
fdata <- data.frame(y = rfoldnorm(n <- 1000, m = m, sd = SD))
hist(with(fdata, y), prob = TRUE, main = paste("foldnormal(m = ", m,
     ", sd = ", round(SD, 2), ")"))
fit <- vglm(y ~ 1, foldnormal, data = fdata, trace = TRUE)
coef(fit, matrix = TRUE)
(Cfit <- Coef(fit))
# Add the fit to the histogram:
mygrid <- with(fdata, seq(min(y), max(y), len = 200))
lines(mygrid, dfoldnorm(mygrid, Cfit[1], Cfit[2]), col = "orange")
# }

Run the code above in your browser using DataLab

Get 50% off unlimited learning