# \donttest{
set.seed(1)
n = 1000
d = 2
X = cbind(rnorm(n), rnorm(n))
y = rbinom(n,1,.5)
fit = maximum.context(y, X)
coef(fit, c(0,0,1,0))
## context in the order: y_{t-1} = 0, y_{t-2} = 0, y_{t-3} = 1, y_{t-4} = 0
# }
Run the code above in your browser using DataLab