Skip to contents

Density, distribution function, quantile function and random generation for a Bimodal Weibull distribution with parameters shape and scale.

Usage

dbwd(x, alpha, beta = 1, sigma, log = FALSE)

pbwd(q, alpha, beta = 1, sigma, lower.tail = TRUE, log.p = FALSE)

qbwd(p, alpha, beta = 1, sigma, lower.tail = TRUE)

rbwd(n, alpha, beta = 1, sigma)

Arguments

x, q

vector of quantiles.

alpha

a shape parameter.

beta

a scale parameter.

sigma

a control parameter that controls the uni- or bimodality of the distribution.

log, log.p

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

lower.tail

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

p

vector of probabilities.

n

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

Value

dbwd gives the density, pbwd gives the distribution function, qbwd gives the quantile function and rbwd generates random deviates.

Details

A Bimodal Weibull distribution with shape parameter \(\alpha\), scale parameter \(\beta\),and the control parameter \(\sigma\) that determines the uni- or bimodality of the distribution, has density $$f\left( x\right) =\frac{\alpha }{\beta Z_{\theta }} \left[ 1+\left( 1-\sigma~x\right) ^{2}\right] \left( \frac{x}{\beta } \right) ^{\alpha -1}\exp \left( -\left( \frac{x}{\beta }\right) ^{\alpha } \right),$$ where $$Z_{\theta }=2+\sigma ^{2}\beta ^{2}\Gamma \left( 1+\left( 2/\alpha \right)\right) -2\sigma \beta \Gamma \left( 1+\left( 1/\alpha \right) \right) $$ and $$x\geq 0,~\alpha ,\beta >0~ and ~\sigma \in\mathbb{R}.$$

References

Vila, R. ve Niyazi Çankaya, M., 2022, A bimodal Weibull distribution: properties and inference, Journal of Applied Statistics, 49 (12), 3044-3062.

Examples

library(new.dist)
dbwd(1,alpha=2,beta=3,sigma=4)
#> [1] 0.01594262
pbwd(1,alpha=2,beta=3,sigma=4)
#> [1] 0.003859685
qbwd(.7,alpha=2,beta=3,sigma=4)
#> [1] 4.759942
rbwd(10,alpha=2,beta=3,sigma=4)
#>  [1] 5.549071 4.330503 4.975350 5.421088 4.053200 3.965505 5.833866 3.590104
#>  [9] 3.555097 2.385201