r - MuMIn::dredge()。仅排除包含主效应的模型，仅保留具有交互作用的模型

Question

是否可以为子集化（排除）所有仅包括例如 A 的主要影响的潜在模型做出一般逻辑陈述？

y ~ A + B + C + A:C + A:B

对于包括 A 的模型，只包括那些 A 是交互的一部分的模型（因为关系 y~A 本身没有意义）。

Answer 1

要么 要么!A || {A:B} || {A:C}但dc(A:B, A) && dc(A:C, A)没有用于指定“任何包含 A 的交互”的符号。

更新：在 MuMIn >= 1.42.3（目前在 R-Forge 上）中，您可以使用 dot.伪函数来指定包含特定一阶项的任何或所有交互。例如 .(A, 2:100)，如果模型包含任何 2 到 100 阶交互项，则包括模型A。它仍然是实验性功能，因此请检查您是否得到了预期的结果。

Answer 2

我得出的一个不优雅的解决方案（在注意到 dredge() 如何与 poly() 项一起工作；包括或排除整个 nx 幂矩阵，而不是各个列）是创建一个名为“int”的小函数来创建一个类似的矩阵两个预测变量及其相互作用。根据 poly() 的约定，我将产生的元素称为“a”、“b”和“a:b”，poly() 将其两个元素称为“1”和“2”。

int = function(a,b){
      matOut = matrix(c(a,b,a*b),nrow=length(a))
      dimnames(matOut) =list(NULL,c("a","b","a:b"))
      return(matOut)}

我们以mtcars数据集为例，用“carb”代替A，用“gear”代替B，用“am”代替C来预测“mpg”。从标准方法开始，我们可以试试这个。

require(MuMIn)
data(mtcars)

fm = lm(mpg ~ carb*gear + carb*am, mtcars,na.action=na.fail)
length(dredge(fm))

这会产生 11 个可能的模型，其中包含三个线性项和两个交互项的所有组合。

现在，如果我们使用 int() 而不是 *，我们会看到我们只得到四个模型；每个交互、它们的组合和空模型。排除包含没有交互作用的线性项的所有模型。

fm = lm(mpg ~ int(carb,gear) + int(carb,am), mtcars,na.action=na.fail)
length(dredge(fm))