我正在使用R中的mboost
包进行仿真。我想使用下面的代码来拟合逐分量线性最小二乘提升算法。但是,我想使用校正后的AIC确定最佳的提升迭代次数。
我是否必须先使用默认的迭代次数mstop = 100
(默认值/用户定义的)运行增强算法,还是可以直接在代码中应用停止准则(即,使代码减少到boost_fit_aic
中的单线)?我还需要修改predict
方法吗?
boost_fit_initial <- glmboost(y ~ ., data = data.frame(y,X), control =
boost_control(mstop = 500))
m_stop <- mstop(aic <- AIC(boost_fit_initial))
boost_fit_aic <- glmboost(y ~ ., data = data.frame(y,X), control =
boost_control(mstop = m_stop))
boost_pred <- predict(boost_fit_aic, data.frame(X_test))