## tests
# load data
data(glmsimdata)
dat <- data.frame(y=glmsimdata$y, glmsimdata$x) ####x is a covariate matrix.
models <- list(y~x1, y~x1+x2, y~x1+x2+x3)
output<-ELCICglm(models, dat, poisson())
output$model.selection
output$glm.output
Run the code above in your browser using DataLab