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.