Генерация переменных взаимодействия в кадрах данных R

Существует ли путь - кроме для цикла - для генерации новых переменных в кадре данных R, который будет всеми возможными взаимодействиями с 2 путями между существующими? т.е. предполагая кадр данных с тремя числовыми переменными V1, V2, V3, я хотел бы генерировать следующие новые переменные:

Inter.V1V2 (= V1 * V2) 
Inter.V1V3 (= V1 * V3)
Inter.V2V3 (= V2 * V3)

Использование в качестве примера для цикла:

x <- read.table(textConnection('
   V1 V2 V3 V4
1  9   25   18
2  5   20   10
3  4   30   12
4  4   34   16'
), header=TRUE)

dim.init <- dim(x)[2]
for (i in 1: (dim.init - 1) ) {
        for (j in (i + 1) : (dim.init) ) {
                x[dim(x)[2] + 1]    <- x[i] * x[j]
                names(x)[dim(x)[2]] <- paste("Inter.V",i,"V",j,sep="")

        }
}