Description Usage Arguments Details Value Author(s) References See Also Examples
This function implements an algortihm to estimate a linear AR(1) model without intercept and explosion by simplicial depth or one of the proposed simplified notions in Kustosz, Mueeller and Wendler (2016).
1 |
y |
An observed series form an linear AR(1) process without intercept. |
maxit |
A value for the maximal number of iterations of the optimisation algorithm. |
candy |
A switch deciding, if just edges of simplices defined by the residuals shall be evaluated |
acc |
A value in (0,1) definig the size of the search regions in each iteration step. |
plots |
A switch (TRUE/FALSE) which enables plots of the iterations steps. |
eps |
An epsilon defining the range around candidate parameters to inflate the candidate set. See |
unique |
A switch (TRUE/FALSE) defining, if a unique maximum shall be the result if multiple points with maximal depth exist. In this case the median point is selected. |
notion |
Here the function which is used for depth calculation is defined. The following notions are allowed: "dS1", "dS2", "dS3", "dS". |
The theoretical details can be found in Kustosz, Mueller and Wendler (2016) and Kustosz and Mueller (2014). The computational details are in Kustosz (2016).
The function returns a list with
estimate |
Parameter maximising the selected depth notion. |
value |
Depth at the resulting maximum. |
numit |
Number of iterations. |
Kustosz, Christoph
Kustosz, C. (2016). Depth based estimators and tests for
autoregressive processes with application. Ph. D. thesis. TU Dortmund.
Kustosz C., Mueller Ch. H. and Wendler M. (2016). Simplified Simplicial Depth for Regression and
Autoregressive Growth Processes. Journal of Statistical Planning and Inference. In press.
Kustosz C. and Mueller Ch. H. (2014). Analysis of crack growth with robust distribution-
free estimators and tests for nonstationary autoregressive processes. Statistical
Papers 55, 125-140.
dS1_lin2
,dS2_lin2
,dS3_lin2
,dS_lin1
,Tri_Mid
,lin1_theta
,Ele_Norm
1 2 | y <- RandomARMod_lin2(250, 0, 1.003, 15, 1)
est_lin1(y, eps = 1e-9)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.