BMT.Phi {BMT}

The BMT-Phi Distribution.

Description

Density, distribution function, quantile function, random number generation for the BMT-Phi distribution with mean equal to mean and standard deviation equal to sd.

Usage

dBMT.Phi(x, mean = 0, sd = 1, log = FALSE)

pBMT.Phi(q, mean = 0, sd = 1, lower.tail = TRUE, log.p = FALSE)

qBMT.Phi(p, mean = 0, sd = 1, lower.tail = TRUE, log.p = FALSE)

rBMT.Phi(n, mean = 0, sd = 1)

Arguments

x, q

vector of quantiles.

mean

vector of means.

sd

vector of standard deviations.

log, log.p

logical; if TRUE, probabilities p are given as log(p).

lower.tail

logical; if TRUE (default), probabilities are P[X \le x], otherwise, P[X > x].

p

vector of probabilities.

n

number of observations. If length(n) > 1, the lenght is taken to be the number required

Details

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

The BMT-Phi distribution is the BMT distribution with \kappa_l = \kappa_r = 0.58029164978583758. The BMT-Phi cumulative distribution function (cdf) is the closest BMT cdf to the normal cdf with the same mean and standard deviation.

Value

dBMT.Phi gives the density, pBMT.Phi the distribution function, qBMT.Phi the quantile function, and rBMT.Phi generates random deviates.

The length of the result is determined by n for rBMT.Phi, and is the maximum of the lengths of the numerical arguments for the other functions.

The numerical arguments other than n are recycled to the length of the result. Only the first elements of the logical arguments are used.

sd <= 0 is an error and returns NaN.

Author(s)

Camilo Jose Torres-Jimenez [aut,cre] [e-mail]

References

Torres-Jimenez, C. J. (2018), The BMT Item Response Theory model: A new skewed distribution family with bounded domain and an IRT model based on it, PhD thesis, Doctorado en ciencias - Estadistica, Universidad Nacional de Colombia, Sede Bogota.

See Also

Distributions for other standard distributions. pBMT for the BMT distribution and pBMT.Psi for the BMT-Psi distribution.

Examples

layout(matrix(1:4,2,2))

curve(pnorm(x), -4, 4, col = "red", lty = 2, ylab = "cdf")
curve(pBMT.Phi(x), add = TRUE, col = "blue", lty = 3)
legend("topleft", legend = c("norm(0,1)","BMT-Phi(0,1)"), 
       bty = "n", col = c("red","blue"), lty = 2:3)

curve(pnorm(x)-pBMT.Phi(x), -4, 4)

curve(qnorm(x), col = "red", lty = 2, xlab = "p", ylab = "qf")
curve(qBMT.Phi(x), add = TRUE, col = "blue", lty = 3)

hist(rBMT.Phi(10000), freq = FALSE, breaks = seq(-4,4,0.25), border = "blue")
curve(dnorm(x), add = TRUE, col = "red", lty = 2)
curve(dBMT.Phi(x), add = TRUE, col = "blue", lty = 3)