Включение перекрестной проверки в пошаговую регрессию в R

У меня есть набор данных из 162 наблюдений со 151 различной переменной, и я хотел бы выполнить для него пошаговую регрессию, но также выполнить для него 10-кратную перекрестную проверку. Я использовал пакет DAAG раньше, чтобы выполнить 10-кратную перекрестную проверку с множественной линейной регрессией, и смог использовать одну из его формул: -

CVlm(df = data, seed=1500, m = 10, form.lm = formula(RT..seconds.~.,), printit=TRUE)

Мне было интересно, поддерживает ли пакет то же самое, но с пошаговой регрессией? Я просмотрел его pdf, но ничего не нашел.

Я знаю, что могу выполнить пошаговую регрессию, используя пакет MASS, выполнив

step <- stepAIC(fit, direction="both")

но не понимаю, как я могу включить в уравнение 10-кратное CV.

Спасибо


person user2062207    schedule 05.12.2013    source источник
comment
Взгляните на пакет caret.   -  person tcash21    schedule 06.12.2013


Ответы (1)


Пакет SuperLearner поддерживает v-кратную перекрестную проверку пошаговой регрессии.

person Michael    schedule 05.12.2013
comment
Привет, спасибо за это. Я только что взглянул на пакет SuperLearner, но мало что могу сказать о том, как его построить для пошагового выполнения. Я тоже немного запутался с PDF. Я понимаю, что вам нужно сделать: - 'CV.SuperLearner(Y, X, V = 20, family = gaussian(), SL.library, method = method.NNLS, id = NULL, verbose = FALSE, control = list( saveFitLibrary = FALSE), cvControl = list(), obsWeights = NULL, saveAll = TRUE, parallel = seq)', но где и как вы указываете, что хотите использовать здесь пошаговую регрессию? - person user2062207; 06.12.2013
comment
Кроме того, как мне указать значения Y и X, как и во всех других алгоритмах регрессии, я бы просто сделал RT..seconds.~. и указал, какой у меня набор данных. - person user2062207; 06.12.2013
comment
Вам нужно определить SL.library. Если вы используете listWrappers(), вы увидите все доступные алгоритмы, включая SL.step. Так что, если вы просто хотите, чтобы SL.library <- c("SL.step"). Исходный код находится здесь. Вы также можете найти код из этого вопроса. - person Michael; 06.12.2013