Nothing
library(midasr)
theta_h0 <- function(p, dk) {
i <- (1:dk-1)/100
pol <- p[3]*i + p[4]*i^2
(p[1] + p[2]*i)*exp(pol)
}
##Generate coefficients
theta0 <- theta_h0(c(-0.1,10,-10,-10),4*12)
xx <- ts(arima.sim(model = list(ar = 0.6), 3000*12), frequency = 12)
aa <- lapply(c(50,100,200,500,1000,1500,2000), function(n) {
y <- midas_auto_sim(n, 0.5, xx, theta0, n_start = 200)
x <- window(xx,start=start(y))
midas_r(y~mls(y,1,1)+fmls(x,4*12-1,12,theta_h0),start=list(x=c(-0.1,10,-10,-10)))
})
sapply(aa,function(x)c(nrow(x$model),coef(x)))
bb <- lapply(c(50,100,200,500,1000,1500,2000), function(n) {
y <- midas_auto_sim(n, c(0.5, 0.1), xx, theta0, n_start = 200)
x <- window(xx, start=start(y))
midas_r(y~mls(y,1:2,1)+fmls(x,4*12-1,12,theta_h0),start=list(x=c(-0.1,10,-10,-10)))
})
sapply(bb,function(x)c(nrow(x$model),coef(x)))
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.