使用R做plm回归后的结果有错

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

panel <- summary(plm(y~x1+x2, data, model="within"))
错误于model.frame.default(terms(formula, lhs = lhs, rhs = rhs, data = data),  :
变数的长度不一样 (發現於 'gdp')
上面红色部分是运行结果出现的错误,
我检查了我的数据,这是一个panel data包含了29个国家,24年的数据,数据并没有unbalanced的部份,
我该如何解决这个问题？非常感谢。

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：plm Panel pane within terms

如果方便，可以把数据上传
这样才能帮忙调试

数据如附件

其中
> sy<-data[ ,3]#ssy
> y<-log(1+sy/100)
然后变数x1即为gdp
谢谢

99258001hong 发表于 2011-11-4 10:50
数据如附件

其中

1. > library(plm)
   
2. > aa=read.table("data.txt",header=TRUE)
   
3. >
   
4. > attach(aa)
   
5. The following object(s) are masked from 'aa (position 3)':
   
6. 
   
7.     i, ssy, x1, x2, Year
   
8. > panel <- plm(ssy~x1+x2, data=aa,index=c("i","Year"))
   
9. > panel
   
10. 
    
11. Model Formula: ssy ~ x1 + x2
    
12. <environment: 03153804>
    
13. 
    
14. Coefficients:
    
15.        x1        x2
    
16. 38.844187  0.020689

复制代码

1. > summary(panel)
   
2. Oneway (individual) effect Within Model
   
3. 
   
4. Call:
   
5. plm(formula = ssy ~ x1 + x2, data = aa, index = c("i", "Year"))
   
6. 
   
7. Balanced Panel: n=29, T=24, N=696
   
8. 
   
9. Residuals :
   
10.    Min. 1st Qu.  Median 3rd Qu.    Max.
    
11. -18.700  -2.340   0.725   3.110  12.500
    
12. 
    
13. Coefficients :
    
14.     Estimate Std. Error t-value Pr(>|t|)   
    
15. x1 38.844187   3.001903 12.9399   <2e-16 ***
    
16. x2  0.020689   0.044522  0.4647   0.6423   
    
17. ---
    
18. Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
    
19. 
    
20. Total Sum of Squares:    20792
    
21. Residual Sum of Squares: 16607
    
22. R-Squared      :  0.2013
    
23.       Adj. R-Squared :  0.19234
    
24. F-statistic: 83.8032 on 2 and 665 DF, p-value: < 2.22e-16

复制代码

谢谢!: D

ywh19860616 发表于 2011-11-4 11:08

请问如果我希望这个within model的参数估计中包含截距项,那麽我该如何修改指令呢?谢谢。

99258001hong 发表于 2011-11-4 23:12
请问如果我希望这个within model的参数估计中包含截距项,那麽我该如何修改指令呢?谢谢。

加上一个变量，取值为
1
1
1
。。。
1
类型的

ywh19860616 发表于 2011-11-5 13:48
加上一个变量，取值为
1
1

我在资料中加入变数x3取值为
1
1
1
．．．
1
再跑一次plm则显示
serie x3  is constant and has been removed
估计结果同未包含截距项
它会自动删除constant
这样的话我该如何处理
谢谢

99258001hong 发表于 2011-11-5 17:15
我在资料中加入变数x3取值为
1
1

根据固定效应的估计方法，截距项会被drop
若要得到截距项，可以参考
https://bbs.pinggu.org/forum.php? ... ;page=1#pid10041384