source("example_IPprobit.R")
library(mvtnorm, lib.loc="/scratch/ufhpc/viknesh/Rlibs")
library(MASS)
library(varSelectIP, lib.loc="/scratch/ufhpc/viknesh/Rlibs")
varSelectIP(a=0.2, keep=5, covariates.retain=mydata2[,4], model.type="probit", 
            q=5, covariates.test=mydata2[,5:10], response=mydata2[,1], 
            nsim=32, out.fname="probit.models_02.csv")
varSelectIP(a=0.2, keep=5, covariates.retain=mydata2[,4], model.type="reg", 
            q=5, covariates.test=mydata2[,5:10], response=mydata2[,2], nsim=32,
            out.fname="linReg.models_02.csv")