我正在使用标准的lm模型,并希望将一个或多个变量的系数设置为特定的整数。
例如,我希望我的天气和价格变量系数分别为647和15。
我正在使用带有标准公式的lm函数。
到目前为止,我发现的最接近的是glm中的偏移函数,或者systemfit中的restrict.rhs。
我还研究了用它们的系数集减去这些变量的总贡献,但这不是很可扩展。
我知道设置系数的所有问题,但是只想从探索的角度来看。
答案 0 :(得分:12)
您可以在公式中使用偏移项,并在其中包含所需的系数和变量:
df<-data.frame(aa=1:6,bb=2:7,cc=c(4,2,7,5,8,3))
lm(cc ~ aa + offset(647*bb), data = df)
所以这会在cc
加上固定字词aa
后退bb * 647
。对于多个给定系数,请添加适当的其他offset()
项。