面板的有序分析

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

求教各位，如何用stata做面板的有序概率模型？本人急用，所以请各位帮忙

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：请各位帮忙 Stata tata 如何用 模型 如何

stata13才有
xtologit . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Random-effects ordered logistic models 263
xtoprobit . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Random-effects ordered probit models 275

蓝色 发表于 2013-7-3 06:59
stata13才有
xtologit . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Random ...

太感谢了，谢谢谢谢，帮了我大忙了

看看

或者   regoprob2  命令

Estimating ordered categorical variables using panel data: a generalized ordered probit model with an autofit procedure
http://mpra.ub.uni-muenchen.de/23203/

小撂加油 发表于 2013-7-3 07:10
太感谢了，谢谢谢谢，帮了我大忙了

还有个问题想请教下，我用的是Stata11,现在是不是无法升级到13呀？

蓝色 发表于 2013-7-3 07:34
或者   regoprob2  命令

Estimating ordered categorical variables using panel data: a generalized o ...

本人为初学者，不是很懂，还想请教一下：我用的是Stata11,现在能升级到13吗？还有regoprob2是stata中的命令吗？非常感谢您的帮助

stata13需要购买
那个命令是可以下载的，是别人自己写的

蓝色 发表于 2013-7-3 09:50
stata13需要购买
那个命令是可以下载的，是别人自己写的

谢谢指教
我运行了regoprob2，但为什么总出现no observations？是数据输入有问题吗？

估计是数据的问题
clear
set obs 1000
forv i=0/4{
                    gen x`i'=uniform()
}

gen     y=1   if x0<=0.3
replace y=2   if x0>0.3 & x0<=0.6
replace y=3   if x0>0.6

egen id=seq(),b(200)
bysort id: gen t=_n


regoprob2 y x1-x4,i(id)


根据产生的模拟数据，是没有问题的


. clear

. set obs 1000
obs was 0, now 1000

. forv i=0/4{
  2.                     gen x`i'=uniform()
  3. }

.
. gen     y=1   if x0<=0.3
(694 missing values generated)

. replace y=2   if x0>0.3 & x0<=0.6
(321 real changes made)

. replace y=3   if x0>0.6
(373 real changes made)

.
. egen id=seq(),b(200)

. bysort id: gen t=_n

.
.
. regoprob2 y x1-x4,i(id)

Fitting constant-only model:

Fitting full model:

Iteration 0:   log likelihood = -1158.4897  (not concave)
Iteration 1:   log likelihood = -1125.6207  
Iteration 2:   log likelihood = -1099.2043  
Iteration 3:   log likelihood = -1097.3984  (not concave)
Iteration 4:   log likelihood = -1097.3983  
Iteration 5:   log likelihood = -1097.3965  
Iteration 6:   log likelihood = -1097.3965  

Random Effects Generalized Ordered Probit         Number of obs   =       1000
                                                  LR chi2(8)      =     -11.34
Log likelihood = -1097.3965                       Prob > chi2     =     1.0000

------------------------------------------------------------------------------
           y |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
mleq1        |
          x1 |  -.0702588   .1442088    -0.49   0.626    -.3529029    .2123854
          x2 |  -.0362581   .1449297    -0.25   0.802    -.3203152    .2477989
          x3 |   .0689639   .1482779     0.47   0.642    -.2216554    .3595831
          x4 |   .1475481   .1406295     1.05   0.294    -.1280807    .4231768
       _cons |   .8875506   266.6399     0.00   0.997     -521.717    523.4921
-------------+----------------------------------------------------------------
mleq2        |
          x1 |  -.1283095   .1392873    -0.92   0.357    -.4013076    .1446887
          x2 |   .0093806   .1420836     0.07   0.947    -.2690982    .2878593
          x3 |   .2402977   .1445976     1.66   0.097    -.0431084    .5237037
          x4 |  -.0676559   .1381038    -0.49   0.624    -.3383344    .2030226
       _cons |   .0828413   266.6399     0.00   1.000    -522.5217    522.6874
-------------+----------------------------------------------------------------
rho          |
       _cons |   .4871581   307.6014     0.00   0.999    -602.4005    603.3748
------------------------------------------------------------------------------

.
end of do-file