par_logn_mix_known {LNPar}R Documentation

Estimate the parameters of a lognormal-Pareto density, assuming a known threshold

Description

This function estimates the parameters of a Pareto and a lognormal density, assuming a known threshold.

Usage

par_logn_mix_known(y, prior1, th, alpha, mu, sigma)

Arguments

y

non-negative numerical vector: random sample from the mixture.

prior1

scalar (0<prior1<1): starting value of the prior probability.

th

positive scalar: threshold.

alpha

non-negative scalar: starting value of the Pareto shape parameter.

mu

scalar: starting value of the lognormal parameter mu.

sigma

positive scalar: starting value of the lognormal parameter sigma.

Value

A list with the following elements:

xmin: estimated threshold.

prior: estimated mixing weight.

post: matrix of posterior probabilities.

alpha: estimated Pareto shape parameter.

mu: estimated expectation of the lognormal distribution on the lognormal scale.

sigma: estimated standard deviation of the lognormal distribution on the lognormal scale.

loglik: maximized log-likelihood.

nit: number of iterations.

Examples

mixFit <- par_logn_mix_known(TN2016, .5, 4700, 3, 7, 1.2)

[Package LNPar version 0.1.0 Index]