我目前正在使用两组虚拟变量进行线性回归 -
model1 <- lm(y ~ factor1 + factor2, data=mydata)
其中 factor1 包含 30 个分类变量,factor2 包含 5 个分类变量。
我有两个约束 - 总和(每个因子 1 的系数 * 每个因子 1 的某个常数项)= 0 和总和(每个因子 2 的系数 * 每个因子 2 的某个常数项)= 0。
我意识到该解决方案在某种程度上与向 lm 函数添加“对比度”有关,但我不知道如何正确编码。任何人都可以请你启发我。或者请随时留下您对解决此类问题的替代方法的建议
非常感谢。