stata 用极大似然估计求pin，报错说variable __00000D not found，求问程序错在哪里

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

1. use BS.dta,clear
   
2. cap program drop poisson
   
3. program define poisson
   
4.   version 1.0
   
5.   args lnf Alpha Epsilon Delta Mu
   
6.   tempvar lnf P M N R Q S T G H I
   
7.   quietly{
   
8.   gen double `P'=exp(-`Epsilon')
   
9.   gen double `M'=`Epsilon'^$ML_y1
   
10.   gen double `N'=`Epsilon'^$ML_y2
    
11.   gen double `S'=exp(lnfactorial($ML_y1))
    
12.   gen double `T'=exp(lnfactorial($ML_y2))
    
13.   gen double `G'=exp(-(`Epsilon'+`Mu'))
    
14.   gen double `H'=(`Mu'+`Epsilon')^$ML_y1
    
15.   gen double `I'=(`Mu'+`Epsilon')^$ML_y2
    
16.   replace `lnf'=ln((1-`Alpha')*`P'*`M'*`P'*`N'/(`S'*`T')+`Alpha'*`Delta'*`P'*`M'`G'*`I'/(`S'*`T')+`Alpha'*(1-`Delta')*`G'*`H'*`P'*`I'/(`S'*`T'))
    
17.   }
    
18.   end
    
19.   ml model lf poisson (Alpha:B S=)(Epsilon:B S=)(Delta:B S=)(Mu:B S=)
    
20.   ml maximize

复制代码

  

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Variable 极大似然估计 found Stata tata version double

数据呢？

夏目贵志 发表于 2017-4-18 23:08
数据呢？

我知道我的程序里少打了一个乘号，改了之后可以运行，但是找不到初始值。ml search也用了，没有用。然后，给四个估计值写循环，来赋值初始值也没有用。我想知道是不是ml model lf poisson (Alpha:B S=)(Epsilon:B S=)(Delta:B S=)(Mu:B S=)这句话写错了？？？

B        S
1127        950
1668        1235
1217        1322
1161        1059
1166        1100
1409        1563
1932        1253
1291        1311
1147        1154
870        1083
1093        1135
1392        1025
1412        1479
1306        1251
1246        1507
1207        1304
1610        1436
1277        1232
1344        1335
1011        1025
1490        1400
1210        1090
1079        988
884        1217
1007        1071
1068        981
1296        977
938        1110
1329        1465
935        923
1102        1117
918        1233
908        1090
1096        962
1107        1235
1159        1287
1056        1056
970        1117
1104        1005
1239        1343
1221        1042
985        1185
1071        1182
1162        1208
1075        1140
897        1116
861        1028
937        938
1023        1244
1067        1017
1129        972
957        887
1065        816
1042        1157
1004        1407
1312        1119
1273        1376
1160        1170
1459        1387
1244        1231
1562        1408

求问有人知道吗？？

你的参数估计PIN做出来了吗？我现在也想用stata，就是不知道该怎么做

我刚才看了一下，在你的tempvar中好像多了一个Inf

hittilehua 发表于 2017-7-17 17:16
我刚才看了一下，在你的tempvar中好像多了一个Inf

改了也不行，你有没有做出来？？？

lilyyung 发表于 2017-4-18 16:53

哪位大神可以帮忙用var分析下数据吖?写论文，但是不会用

hittilehua 发表于 2017-7-17 17:16
我刚才看了一下，在你的tempvar中好像多了一个Inf

大神，可以帮忙用var分析下数据吗?毕业论文，不会用

懒吖头 发表于 2017-9-10 21:57
大神，可以帮忙用var分析下数据吗?毕业论文，不会用

用var分析数据时什么意思？什么问题？