带约束的线性模型

Question

我对 R 很陌生，但遇到以下问题：

我有一个简单的 2 因子线性模型：

Rate~factor1 + factor2 //factor1 has 8 categorical values, factor2 has 6 categories;
model1 <- lm(Rate~factor1+factor2, data=myData)

并且想将factor1系数的约束SUM = 0，factor2也一样。

没有任何手册提供任何线索。

我找到了类似问题here 的链接，但它有所不同，我不知道如何修改它...

Answer 1

在 MASS (Modern Applied Statistics with S) 的第 6 章中有描述。使用lm 的contrasts arg（例如，查看?contr.sum 和?model.matrix.default）。