0

我目前正在使用两组虚拟变量进行线性回归 -

model1 <- lm(y ~ factor1 + factor2, data=mydata)

其中 factor1 包含 30 个分类变量,factor2 包含 5 个分类变量。

我有两个约束 - 总和(每个因子 1 的系数 * 每个因子 1 的某个常数项)= 0 和总和(每个因子 2 的系数 * 每个因子 2 的某个常数项)= 0。

我意识到该解决方案在某种程度上与向 lm 函数添加“对比度”有关,但我不知道如何正确编码。任何人都可以请你启发我。或者请随时留下您对解决此类问题的替代方法的建议

非常感谢。

4

1 回答 1

0

你可以使用 glmc 包吗?但我不知道如何设置约束。

于 2012-05-06T17:27:20.887 回答