dgsm: Computing probability density function of the gamma shape...

In ForestFit: Statistical Modelling for Plant Size Distributions

Computing probability density function of the gamma shape mixture model

Description

Computes probability density function (pdf) of the gamma shape mixture (GSM) model. The general form for the pdf of the GSM model is given by

f(x,{Θ}) = ∑_{j=1}^{K}ω_j \frac{β^j}{Γ(j)} x^{j-1} \exp\bigl( -β x\bigr),

where Θ=(ω_1,…,ω_K, β)^T is the parameter vector and known constant K is the number of components. The vector of mixing parameters is given by ω=(ω_1,…,ω_K)^T where ω_js sum to one, i.e., ∑_{j=1}^{K}ω_j=1. Here β is the rate parameter that is equal for all components.

Usage

dgsm(data, omega, beta, log = FALSE)

Arguments

data	Vector of observations.
omega	Vector of the mixing parameters.
beta	The rate parameter.
log	If TRUE, then log(pdf) is returned.

Value

A vector of the same length as data, giving the pdf of the GSM model.

Author(s)

Mahdi Teimouri

References

S. Venturini, F. Dominici, and G. Parmigiani, 2008. Gamma shape mixtures for heavy-tailed distributions, The Annals of Applied Statistics, 2(2), 756–776.

Examples

data<-seq(0,20,0.1)
omega<-c(0.05, 0.1, 0.15, 0.2, 0.25, 0.25)
beta<-2
dgsm(data, omega, beta)

ForestFit documentation built on March 7, 2023, 8:27 p.m.