dexp_mix: Exponential Mixture Density

In matheuscastro43/finiteMix: Finite Mixture Models

Exponential Mixture Density

Description

Gives the density value of a univariate Exponential Mixture.

Usage

dexp_mix(x, pi, rate, log = FALSE)

Arguments

x	value or vector of values, to apply in the mixture probability density function.
pi	vector of mixture weights.
rate	vector of rate parameters.
log	logical; if TRUE, probabilities p are given as log(p).

Details

Each component of the mixture has density

f_j(x) = λ_j exp(-λ_j * x)

where λ_j is the rate paramter of the jth component. j = 1, 2, ..., g.

Author(s)

CASTRO, M. O.; MONTALVO, G. S. A.

Examples

f <- function(x){dexp_mix(x, pi = c(0.1, 0.9), rate = c(1, 1/30))}
## Plot using graphics package
require(graphics)
curve(f, 0, 40)

require(ggplot2)
## Plot using ggplot2 package
ggplot(data.frame(x=c(0, 40)), aes(x)) +
  stat_function(fun=f, size = 1L, colour = "#0c4c8a") +
  theme_minimal()

## Integrate
integrate(f, 0, "Inf")$value

matheuscastro43/finiteMix documentation built on March 30, 2022, 12:49 p.m.