Jak wykonać wszystkie interakcje przed użyciem glmnet

Question

Mam macierz x z 8 kolumnami. Chcę uruchomić glmnet, aby przeprowadzić regresję lasso. Wiem, że muszę zadzwonić:

glmnet(x, y, family = "binomial", ...). 

Jednak, jak mogę dostać x rozważyć wszystkie jeden sposób interakcji, jak również? Czy muszę ręcznie przerobić ramkę danych: jeśli tak, czy istnieje łatwiejszy sposób? Przypuszczam, że miałem nadzieję, że coś zrobię, używając formuły R.

Answer 1

Tak, istnieje wygodny sposób na to. Ważne są dwa kroki.

library(glmnet) 
# Sample data 
data <- data.frame(matrix(rnorm(9 * 10), ncol = 9)) 
names(data) <- c(paste0("x", 1:8), "y") 
# First step: using .*. for all interactions 
f <- as.formula(y ~ .*.) 
y <- data$y 
# Second step: using model.matrix to take advantage of f 
x <- model.matrix(f, data)[, -1] 
glmnet(x, y) 

Answer 2

f <- as.formula(~ .^2) powinny również działać w tym Główne wpływy i interakcje wszystkich parach