# we adapted the function a little for the dodona exercise # for testing your code please use the predict function defined above # for submitting your code on dodona, use this function instead predict.regsubsets <- function(object, newdata, id, ...) { mat <- model.matrix(medv ~ ., data = newdata) coefi <- coef(object, id = id) xvars <- names(coefi) mat[, xvars] %*% coefi } k <- ... folds <- ... cv.errors <- ... for (j in 1:...) { best.fit <- ... for (i in 1:...) { pred <- ... cv.errors[...] <- ... } } mean.cv.errors <- ... reg.best <- ... reg.best.coef <- ...