rDynMix: Simulating a dynamic lognormal-Pareto mixture
Description
This function fits a dynamic mixture by Approximate Maximum Likelihood and by standard maximum likelihood.
Currently only implemented for the lognormal - generalized Pareto case,
with Cauchy or exponential weight.
Usage
rDynMix(nreps, x, weight)
Value
ysim (nreps x 1) vector: nreps random numbers from the lognormal-GPD dynamic mixture.
Arguments
nreps
integer: number of observations sampled from the mixture.
x
numerical vector: if weight = 'cau', values of \(mu_c\), \(\tau\), \(\mu\),
\(\sigma\), \(\xi\), \(\beta\); if weight = 'exp', values of \(\lambda\), \(\mu\), \(\sigma\), \(\xi\), \(\beta\).
weight
'cau' or 'exp': name of weight distribution.
Details
This function simulates a dynamic lognormal-GPD mixture using
the algorithm of Frigessi et al. (2002, p. 221).