Объяснение значения p в функции ols_step_forward_p пакета olsrr

Question

Рассмотрим следующий пример :-( взят из, например, документации R)

model <- lm(y ~ ., data = surgical)
ols_step_forward_p(model)

В документации сказано, что для выбора конечной модели будет использован метод «Вперед пошагово». Если я не ошибаюсь, разве этот прямой выбор не основан на наблюдении значений p следующих тестовых статистик c: -

где мы рассматриваем модель с K всего предикторов. RSS обозначает остаточную сумму квадратов, а нижние индексы p и K обозначают RSS модели с предикторами p и K. соответственно.

Если я не дам команду penter в ols_step_forward_p, функция по умолчанию p -значение выбрано?

Answer 1

ols_step_forward_p имеет документацию здесь (также через ?olsrr::ols_step_forward_p, если у вас установлен пакет)

Вторая запись («Метод S3 по умолчанию») говорит

ols_step_forward_p (модель, penter = 0.3 , прогресс = FALSE, подробности = FALSE, ...)

(выделение добавлено).

Это означает, что значение по умолчанию penter равно 0,3.

Чтобы увидеть более подробно, что происходит, вы можете посмотреть код для olsrr::ols_step_forward_p.default. Вот соответствующий экспорт кода: *

for (i in seq_len(len_p)) {
            predictors <- c(preds, all_pred[i])
            m <- lm(paste(response, "~", paste(predictors, collapse = " + ")), 
                l)
            m_sum <- Anova(m)
            pvals[i] <- m_sum$`Pr(>F)`[ppos]
        }

Итак, вы можете видеть, что значение p основано на F-тесте из car::Anova. Подробнее см. ?car::Anova.lm ...