Я работаю со следующим набором данных и разделяю его на набор обучающих данных и набор тестовых данных
library(mvdalav)
data(Penta)
set.seed(321)
> idx <- sample(1:30, size = 20, replace = F)
> ptrain <- Penta[idx,]
> ptest <- Penta[-idx,]
Я использую модель множественной регрессии и хочу выбрать лучшую модель по функции шага и R2. Итак, я выполнил шаг:
>selectAIC = step(lm(log.RAI~.,data = ptrain[,-1]),direction="backward")
> selectAIC$call
lm(formula = log.RAI ~ S1 + P1 + S2 + L2 + S3 + L3 + L4 + P4,
data = ptrain[, -1])
Теперь я хочу сделать то же самое, но используя R ^ 2 модели в качестве критерия выбора.
Как я могу это сделать?
С уважением