R/se.Cenarma.R

Defines functions se.Cenarma

Documented in se.Cenarma

#se.Cenarma.R
se.Cenarma <- function(obj, nBoot=250){
  phi <- numeric(nBoot)
  p <- obj$p
  for (iBoot in 1:nBoot){
    y <- boot.Cenarma(obj)$y
    outi <- try(cenarma(y, p=p), silent=TRUE)
    if (identical(class(outi), "try-error")) phi[iBoot] <- NA
    else phi[iBoot] <- coef(outi$outarima)["ar1"]
  }
  sd(phi, na.rm=TRUE)
}

Try the cents package in your browser

Any scripts or data that you put into this service are public.

cents documentation built on May 29, 2017, 11:04 a.m.