rgamma_mix: Random Generator of a Gamma Mixture

View source: R/rgamma_mix.R

Random Gamma MixtureR Documentation

Random Generator of a Gamma Mixture

Description

Generates random samples of a univariate Gamma Mixture.

Usage

 rgamma_mix(n, pi, alpha, beta, plot.it = TRUE, empirical = FALSE, 
            col.pop = "red3", col.empirical = "navy", ...) 

Arguments

n

number of observations.

pi

vector of mixture weights.

alpha

shape parameters.

beta

scale parameters.

plot.it

logical; if TRUE (default), plot the histogram with estimated distribution curve.

empirical

logical; if TRUE, add the empirical curve ("Kernel Density Estimation") in the plot.

col.pop

a colour to be used in the curve of population density.

col.empirical

a colour to be used in the curve of empirical density.

...

further arguments and graphical parameters passed to hist.

Author(s)

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

Examples

## Generate a sample.
data = rgamma_mix(n = 1000, pi = c(0.6, 0.4), alpha = c(2, 4), beta = c(1, 2))

## Sample vector.
data$sample
## The classification of each observation.
data$classification
## The histogram of the sample with population density curve.
data$plot

## Not plotting the graphic.
rgamma_mix(n = 1000, pi = c(0.6, 0.4), alpha = c(2, 4), beta = c(1, 2), plot.it = FALSE)

## Adding the empirical curve to the graphic.
rgamma_mix(n = 1000, pi = c(0.6, 0.4), alpha = c(2, 4), beta = c(1, 2), plot.it = TRUE,
          empirical = TRUE)

## Changing the color of the curves.
rgamma_mix(n = 1000, pi = c(0.6, 0.4), alpha = c(2, 4), beta = c(1, 2), plot.it = TRUE,
          empirical = TRUE, col.pop = "blue", col.empirical = "green")

## Using "...".
rgamma_mix(n = 1000, pi = c(0.6, 0.4), alpha = c(2, 4), beta = c(1, 2), breaks = 300)

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