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