求救，Stata 编程

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

nb1
. clear
. set obs 50000
. gen x1= runiform()
. gen x2= runiform()
. gen xb = .5 + 1.25*x1 -- 1.5*x2
. gen exb=exp(xb)
. gen idelta = 2*exb
. gen xg=rgamma(idelta, 1/idelta)
. gen xbg = exb*xg
. gen nb1y = rpoisson(xbg)
NB2
clear
set obs 50000
set seed 4321
gen x1 = rnormal()
gen x2 = rnormal()
gen xb = 2 + 0.75*x1 - 1.25*x2
gen a = .5
gen ia = 1/a
gen exb = exp(xb)
gen xg = rgamma(ia, a)
gen xbg = exb * xg
gen nby = rpoisson(xbg)
glm nby x1 x2, nolog fam(nb ml)
我用NB1模拟产生数据，如何再用NB2去拟合NB1产生的数据啊？上面程序是用NB1和NB2分别产生数据，然后拟合，我如何用一个模型去拟合另一个模型？
还有谁会NB-P模型模拟数据的编程？求帮忙

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Stata tata poisson Uniform normal

把变量名弄一致了就行了

davalic 发表于 2014-2-17 11:09
把变量名弄一致了就行了

还请你帮帮我，说详细点，谢谢啦。