Nothing
ShamanStine1T.PI <-
function(x,p,h,nboot,prob)
{
set.seed(12345)
n <- nrow(x)
BC <- Shaman.StineT(x,p,h)
BCB <- Shaman.StineBT(x,p,h)
bb <- BCB$coef
eb <- sqrt( (n-p) / ( (n-p)-length(bb)))*BCB$resid
ef <- sqrt( (n-p) / ( (n-p)-length(bb)))*BC$resid
fore <- matrix(NA,nrow=nboot,ncol=h)
for(i in 1:nboot)
{
index <- as.integer(runif(n-p, min=1, max=nrow(eb)))
es <- eb[index,1]
xs <- ysbT(x, bb, es)
bs <- Shaman.StineT(xs,p,h)$coef
fore[i,] <- ART.ForeB(xs,bs,h,ef,length(bs)-1)
}
Interval <- matrix(NA,nrow=h,ncol=length(prob),dimnames=list(1:h,prob))
for( i in 1:h)
Interval[i,] <- quantile(fore[,i],probs=prob)
return(list(PI=Interval,forecast=BC$forecast))
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.