dlindley_mix: Lindley Mixture Density

In matheuscastro43/finiteMix: Finite Mixture Models

Lindley Mixture Density

Description

Gives the density value of a univariate Lindley Mixture.

Usage

dlindley_mix(x, pi, beta, log = FALSE)

Arguments

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

Details

Each component of the mixture has density

f_j(x) = (x + 1)/(β_j * (β_j + 1)) * exp(-x/β_j)

where β_j is the parameter of the jth component. j = 1, 2, ..., g.

Author(s)

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

Examples

f <- function(x){dlindley_mix(x, pi = c(0.1, 0.9), beta = c(5, 2))}
## Plot using graphics package
require(graphics)
curve(f, 0, 15)

require(ggplot2)
## Plot using ggplot2 package
ggplot(data.frame(x=c(0, 15)), 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.