glmnet을 사용하기 전에 모든 상호 작용을하는 방법

Question

나는 x- 행렬이 8 열입니다. 올가미 회귀를 수행하려면 glmnet을 실행하고 싶습니다. 내가 호출 할 필요가 알고

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

그러나, 어떻게 나뿐만 아니라 모두 하나의 방법으로 상호 작용을 고려 x을받을 수 있나요? 수동으로 데이터 프레임을 다시 만들어야합니까? 그렇다면 더 쉬운 방법이 있습니까? 나는 R 공식을 사용하여 무언가를하기를 희망한다고 가정합니다.

Answer 1

예, 편리한 방법이 있습니다. 두 단계가 중요합니다.

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)는 주 효과 모든 페어의 상호 작용을 포함하여 작동한다