lagged <- function(x, k) {
L <- length(x+k)
z <- matrix(NA, nrow = L + k - 1, ncol = k)
for (i in 1:k) {
z[(1:L) + (i - 1), i] <- x
}
return(z)
}
x <- rep(1, 30)
x <- lagged(x, 5)
xx <- rowSums(x)
wf <- function(x, r) {
w <- exp(x * r)
w/sum(w)
}
z <- t(t(x) * wf(1:5, -0.1))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.