| simACD-class | R Documentation |
A specification class to create an object of a simulated piecewise constant conditional duration model of order (1,1).
x_t / \psi_t = \varepsilon_t \; \sim \mathcal{G}(\theta_2)
\psi_t = \omega(t) + \sum_{j=1}^p \alpha_{j}(t)x_{t-j} + \sum_{k=1}^q \beta_{k}(t)\psi_{t-k}.
where \psi_{t} = \mathcal{E} [x_t | x_t,\ldots,x_1| \theta_1] is the conditional mean duration of the t-th event with parameter vector \theta_1 and \mathcal{G}(.)
is a general distribution over (0,+\infty) with mean equal to 1 and parameter vector
\theta_2. In this work we assume that \varepsilon_t \; \sim \exp(1).
Returns an object of simACD class.
xThe durational time series.
psiThe psi time series.
NSample sze of the time series.
cp.locThe vector with the location of the changepoints. Takes values from 0 to 1 or NULL. Default is NULL.
lambda_0The vector of the parameters lambda_0 in the ACD series as in the above formula.
alphaThe vector of the parameters alpha in the ACD series as in the above formula.
betaThe vector of the parameters beta in the ACD series as in the above formula.
BurnInThe size of the burn-in sample. Note that this only applies at the first simulated segment. Default is 500.
Korkas, K.K., 2022. Ensemble binary segmentation for irregularly spaced data with change-points. Journal of the Korean Statistical Society, 51(1), pp.65-86.
pw.acd.obj <- new("simACD")
pw.acd.obj@cp.loc <- c(0.25,0.75)
pw.acd.obj@lambda_0 <- c(1,2,1)
pw.acd.obj@alpha <- rep(0.2,3)
pw.acd.obj@beta <- rep(0.7,3)
pw.acd.obj@N <- 3000
pw.acd.obj <- pc_acdsim(pw.acd.obj)
ts.plot(pw.acd.obj@x)
ts.plot(pw.acd.obj@psi)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.