# NOT RUN {
myRegressionData <- rmultreg(25, k=5, p=.15, sdnoise = .25)
pairs(myRegressionData$data)
out <- ices(y ~ ., data = myRegressionData$data) # fit model to simulated data
summary(out) # estimates and standard errors for all coefficients
myRegressionData$coefficients # compare with true coefficients
# }
Run the code above in your browser using DataLab