# NOT RUN {
myRegressionData <- rmultreg(100, k=20, p=.1, sdnoise = 1)
pairs(myRegressionData$data)
out <- ices(y ~ ., data = myRegressionData$data) # fit model to simulated data
confint(out) # calculate 95% confidence intervals for all coefficients
myRegressionData$coefficients # compare with true coefficients
# }
Run the code above in your browser using DataLab