dpareto_em: EM algorithm function for discrete Pareto distribution

In camponsah/BivMixDist: Fit bivariate continous and discrete distributions as well as provide EM algorithm for discrete Pareto

Description

This function computes the parameter estimates of discrete Pareto distribution using the EM algorithm.

Usage

dpareto_em(N, delta = 1, p = NULL, maxiter = 1000, tol = 1e-08, verb = FALSE)

Arguments

N	vector of random sample from discrete Pareto distribution
delta	shape parameter which must be numeric greater than or equal to 0
p	numeric parameter between 0 and 1
maxiter	maximum number of iterations
tol	tolerance value
verb	If TRUE, estimates are printed during each iteration.

Details

Takes initial guess for the parameters in discrete Pareto and the algorithm will estimate the MLE.

Value

list containg parameter estimates, Deviance and data frame of iteration

References

Amponsah, C. K., Kozubowski, T. J. and Panorska (2019). A computational approach to estimation of discrete Pareto parameters. Inprint.

Examples

N<-rdpareto(500, delta=0.2,p=0.6)
fit<-dpareto_em(N,maxiter = 1000)
fit$par

camponsah/BivMixDist documentation built on Nov. 15, 2021, 3:11 a.m.