Come rendere tutte le interazioni prima di utilizzare glmnet

Question

Ho una matrice x di 8 colonne. Voglio eseguire glmnet per eseguire una regressione lazo. So che ho bisogno di chiamare:

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

Tuttavia, come faccio ad ottenere x di prendere in considerazione tutte le interazioni in un modo così? Devo rifare manualmente il frame dei dati: se sì, c'è un modo più semplice? Suppongo che speravo di fare qualcosa usando una formula R.

Answer 1

Sì, c'è un modo conveniente per quello. Due passaggi in esso sono importanti.

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) dovrebbe funzionare anche per l'inclusione di effetti principali e tutte le interazioni a coppie