1  | 
x | 
|
y | 
|
regfun | 
|
varfun | 
|
... | 
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21  | ##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.
## The function is currently defined as
function (x, y, regfun = tsreg, varfun = pbvar, ...) 
{
    xy = elimna(cbind(x, y))
    x = as.matrix(x)
    px = ncol(x)
    px1 = px + 1
    n = nrow(xy)
    vals = NA
    for (i in 1:n) {
        est = regfun(xy[-i, 1:px], xy[-i, px1])$coef
        vals[i] = xy[i, px1] - (est[1] + sum(est[2:px1] * xy[i, 
            1:px]))
    }
    pe = varfun(vals)
    pe
  }
 | 
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.