bernexp: The Bernoulli-Exponential distribution

In qmap: Statistical Transformations for Post-Processing Climate Model Output

Description

Density, distribution function, quantile function and random generation for the Bernoulli-Exponential distribution with parameters prob, and rate.

Usage

dbernexp(x, prob, rate)
pbernexp(q, prob, rate)
qbernexp(p, prob, rate)
rbernexp(n, prob, rate)

Arguments

x,q	vector of quantiles.
p	vector of probabilities.
prob	probability of non-zero event.
n	number of random samples.
rate	rate parameter of the Exponential distribution.

Details

Mixture of the Bernoulli and the Exponential distribution. The mixture is analogue to the one described for the berngamma distribution.

Value

dbernexp gives the density (pdf), pbernexp gives the distribution function (cdf), qbernexp gives the quantile function (inverse cdf), and rbernexp generates random numbers.

Author(s)

Lukas Gudmundsson

See Also

Exponential, berngamma

Examples

data(obsprecip)

(ts <- startbernexp(obsprecip[,1]))
hist(obsprecip[,1],freq=FALSE) 
lines(seq(0,max(obsprecip[,1])),
      dbernexp(seq(0,max(obsprecip[,1])),
               prob=ts$prob,
               rate=ts$rate),
      col="red")

pp <- seq(0.01,0.99,by=0.01)
qq <-quantile(obsprecip[,1],probs=pp) 

plot(qq,pp)
lines(qbernexp(pp,
               prob=ts$prob,
               rate=ts$rate),
      pp,col="red")

plot(qq,pp)
lines(qq,
      pbernexp(qq,
               prob=ts$prob,
               rate=ts$rate),
      col="red")

hist(rbernexp(1000,prob=ts$prob,
              rate=ts$rate),freq=FALSE)

qmap documentation built on May 1, 2019, 7:31 p.m.