powered by
This function simulates a dynamic mixture. Currently only implemented for the lognormal - generalized Pareto case, with Cauchy or exponential weight.
rDynMix(nreps, x, weight)
ysim (nreps x 1) vector: nreps random numbers from the lognormal-GPD dynamic mixture.
integer: number of observations sampled from the mixture.
numerical vector: if weight = 'cau', values of μc, \(\tau\), \(\mu\), \(\sigma\), \(\xi\), \(\beta\); if weight = 'exp', values of \(\lambda\), \(\mu\), \(\sigma\), \(\xi\), \(\beta\).
'cau' or 'exp': name of weight distribution.
This function simulates a dynamic lognormal-GPD mixture using the algorithm of Frigessi et al. (2002, p. 221).
fri02FitDynMix
ysim <- rDynMix(100,c(1,2,0,0.5,0.25,3),'cau')
Run the code above in your browser using DataLab