# NOT RUN {
Xcat <- c(rep(0, 5), rep(1, 5))
th <- rvmc(10, 0, 4) + Xcat
# Compare a model that includes group differences with a model that does not.
IC_compare.circGLM(circGLM(th = th), circGLM(th = th, X = Xcat))
# }
Run the code above in your browser using DataLab