# NOT RUN {
n = 500
eps = generative_process(n,"AR1",c(0.7))
X = as.matrix(generative_model(n,"mod2"))
Y = 3 + 2*X[,2] + eps
reg = slm(Y ~ X, method_cov_st = "fitAR", model_selec = -1)
vcov(reg)
# }
Run the code above in your browser using DataLab