# Generate interaction.ind
interaction.ind <- t(combn(4,2))
# Generate data
set.seed(0)
X <- matrix(rnorm(20), ncol = 4)
y <- X[, 2] + rnorm(5)
# Extract X1 and X1X2 from X1, ..., X4
Extract(X, varind = c(1,5), interaction.ind)
# Extract X5 from X1, ..., X4
Extract(X, varind = 5, interaction.ind)
# Extract using duplicated values
try(Extract(X, varind = c(1,1), interaction.ind)) # this will not run
Run the code above in your browser using DataLab