# \donttest{
set.seed(1)
n = 1000
X = cbind(rnorm(n))
y = rbinom(n,1,.5)
fit = maximum.context(y, X)
## using the original data
predict(fit)
## using new data
predict(fit, new.y = c(0,0,1,0,0), new.X = c(2.3, 1.1, -.2, -3,1))
# }
Run the code above in your browser using DataLab