Полная парная сводная таблица для GLM в R

Question

Я пытаюсь запустить логистическую модель, используя 4 различных переменных предиктора: Глубина (непрерывная), DPM (2 уровня), Близость (3 уровня) и Площадь (2 уровня). Я знаю, что многие взаимодействия усложняют ситуацию, но мне было интересно, возможно ли создать сводную таблицу со всеми парными взаимодействиями (и глубиной). Большое спасибо за любую помощь. Код является:

mod <- glm(cbind(Towards, Total-Towards) ~ DPM*Proximity*Depth*Area,data=Data, family=binomial('logit'))

Я очень ценю любую помощь.

Некоторые данные, чтобы исправить идеи:

Data <- structure(list(Towards = c(4, 7, 9, 0, 15, 10, 11, 23, 1, 4), Total = c(6, 14,
    10, 7, 15, 12, 20, 41, 5, 8), Depth = c(-0.3, -0.25, -0.21, -0.17, -0.05, 0,
    0, 0.25, 0.5, 0.56), DPM = structure(c(2L, 1L, 2L, 2L, 1L, 2L, 1L, 1L, 1L,
    2L), .Label = c("Post", "Pre" ), class = "factor"), Proximity =
    structure(c(2L, 3L, 1L, 1L, 1L, 3L, 2L, 2L, 3L, 2L), .Label = c("East",
    "Far", "Near" ), class = "factor"), Area = structure(c(3L, 1L, 1L, 3L, 1L,
    1L, 3L, 3L, 2L, 3L), .Label = c("Control", "Control‌​", "DPM"), class =
    "factor")), .Names = c("Towards", "Total", "Depth", "DPM", "Proximity",
    "Area"), row.names = c(NA, -10L), class = "data.frame")