使用模型矩阵而不是公式调用lm

Question

我想在R中使用lm来拟合线性模型以获得系数估计和p值+用于总模型拟合的p值（类似ANOVA），因此基本上是来自summary.lm的输出。

问题是我想使用自己的模型矩阵，而不是在调用lm时使用公式指定它。

我已经尝试使用底层的lm.fit函数来实现这一点，但后来我lm.fit了summary.lm函数的便利性，我不想再自己重新计算所有的测试统计数据。

有没有办法'欺骗' lm并给它一个模型矩阵而不是公式？

谢谢！

Answer 1

以下是使用内置BOD数据框的示例：

# inputs
demand <- BOD$demand
mm <- model.matrix(~ Time, BOD) # model matrix

summary(lm(demand ~ mm + 0))

要么

summary(lm(demand ~. + 0, as.data.frame(mm)))