### R code from vignette source 'IBGS.Rnw'

###################################################
### code chunk number 1: setup
###################################################
library(IBGS)
set.seed(1)
n <- 100
p <- 60
x <- matrix(rnorm(n * p), n, p)
y <- as.numeric(x[, 1:3] %*% c(3, -3, 3) + rnorm(n))


###################################################
### code chunk number 2: fit
###################################################
fit <- glmIBGS(y, x, criterion = "BIC")
fit


###################################################
### code chunk number 3: summary
###################################################
summary(fit)


###################################################
### code chunk number 4: coef
###################################################
head(coef(fit))


###################################################
### code chunk number 5: margprob
###################################################
plotMargProb(fit, n.vars = 15)


###################################################
### code chunk number 6: predict
###################################################
predict(fit, x[1:5, ])
head(fitted(fit))