关于lm()拟合的一个问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

各位大神，我在用lm()拟合过程中，有一些系数我已经知道了不想拟合，只想拟合其他未知系数可以吗？
例如我想拟合y~x1+x2，他会拟合出3个参数：截距，x1的系数和x2的系数。但是我现在我想固定x1的系数只拟合x2的系数和截距可以吗？我把程序改成
y~1*x1+x2，这样可以吗？我试了试好像算出来结果不对，不知各位大神有何见解？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：何见解 程序

lm是不可以固定一个已经参数，然后进行其它约束估计的。

你所说的，已经知道某个变量的系数，这样的说法其实是有问题的。因为你所知道的值，也是估计值。唯一不同的是，你用的数据不同，得到的估计值也是不同的。

对于每一组样本，你得到的估计值都是固定的，这个是用样本去估计总体。而总体的值是很难知道的。这就是模型的模糊性所在。

jiangbeilu 发表于 2016-2-25 10:45
lm是不可以固定一个已经参数，然后进行其它约束估计的。

你所说的，已经知道某个变量的系数，这样的说法 ...

谢谢。我知道了。但是我在利用systemfit包做联立方程组求解过程中遇到了另一个问题，我的方程组如下
e1<-LV~LD+LH
e2<-LABOVE~LD+LH+LV
e3<-LROOT~LD+LH+LABOVE
其中外生变量是LD和LH，内生的是LV LAVOVE LROOT，在第二个方程中我需要LV的系数是1，第三个方程中LABOVE的系数也是1，我该如何设置呢？
注：这个方程组是通过幂函数取对数后转化成为线性的，转化过来以后第二个方程里的LV和第三个里的LABOVE的系数就是1，不需要估计。。该怎么做？万分感谢您的解答。。

zhao_ 发表于 2016-2-25 11:09
谢谢。我知道了。但是我在利用systemfit包做联立方程组求解过程中遇到了另一个问题，我的方程组如下
e1

结构方程，建议用stata做方便一些。R做结构方程我没有用过，不过可以用sem包。

jiangbeilu 发表于 2016-2-25 12:29
结构方程，建议用stata做方便一些。R做结构方程我没有用过，不过可以用sem包。

那我再试试你说的sem包

还有没有人啦[loveliness][loveliness]

把y2=y-a*x1的值算出来，lm(y2~x2)不就行了

czrdiao 发表于 2016-2-25 19:14
把y2=y-a*x1的值算出来，lm(y2~x2)不就行了

这样是可以拟合出系数，但是拟合之后得到的R平方和残差都不是y的了，而是y2的