Simulate a synthetic Markov chain from a fitted
simmcpot(object, plot = TRUE, ...)
An object of class
Other optional arguments to be passed to the
The simulated Markov chain is computed as follows:
Simulate a Markov chain
prob with uniform margins on
(0,1) and with the fixed extreme value dependence given by
prob such as prob <= 1
- pat, set mc = NA (where
pat is given by
prob such as prob >= 1
- pat, set prob2 = (prob
- 1 + pat) / pat. Thus,
prob2 are uniformly distributed on
mc = qgpd(prob2, thresh,
scale, shape), where
thresh, scale, shape are given by the
object$threshold, object$param["scale"] and
A Markov chain which has the same features as the fitted object. If
plot = TRUE, the Markov chain is plotted.
data(ardieres) flows <- ardieres[,"obs"] Mclog <- fitmcgpd(flows, 5) par(mfrow = c(1,2)) idx <- which(flows <= 5) flows[idx] <- NA plot(flows, main = "Ardieres Data") flowsSynth <- simmcpot(Mclog, main = "Simulated Data")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.