请教，哪位知道如何在模拟程序中输出多次模拟的相关统计值的均值及显著性值呢？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

capture prog drop mono
prog mono
  drop _all
  set obs 200
  gen y=uniform()
  gen for=1 if y>=0&y<=0.3333334
  gen dom=1 if y>0.33333334&y<=0.6666667
  gen joi=1 if y>0.6666667&y<=1
  replace for=0 if for==.
  replace dom=0 if dom==.
  replace joi=0 if joi==.
  gen overall=5.25+invnormal(uniform()) if for==1
  replace overall=5.34+invnormal(uniform()) if dom==1
  replace overall=5.24+invnormal(uniform()) if joi==1
  gen business=5.25+invnormal(uniform()) if for==1
  replace business=5.34+invnormal(uniform()) if dom==1
  replace business=5.24+invnormal(uniform()) if joi==1
  gen political=5.27+invnormal(uniform()) if for==1
  replace political=5.4+invnormal(uniform()) if dom==1
  replace political=5.17+invnormal(uniform()) if joi==1
  manova overall business political=for dom
  end
simulate _b,reps(1000):mono
sum
以这个为例，如何输出Manova的模拟相关统计指标呢？？

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：模拟程序 Uniform capture drop Prog 程序 如何 统计

发现程序模拟还是xxxxx说明模拟出问题了，不过还是不知道如何将结果结合起来？

simulate 命令有个 -saving()- 选项，可以把模拟的中间结果保存下来，你可以看看。

得到t值以后，根据相应的函数转换即可得到p值。

谢谢连老师啊。我试试。

咋还是调不出来呢？？郁闷咯。

capture prog drop mono
prog mono
clear
set obs 200
g y=uniform()
g for=(y>=0&y<=1/3)
g dom=(y>1/3&y<=2/3)
g joi=(y>2/3&y<=1)
g overall=5.25+invnormal(uniform()) if for==1
replace overall=5.34+invnormal(uniform()) if dom==1
replace overall=5.24+invnormal(uniform()) if joi==1
g business=5.25+invnormal(uniform()) if for==1
replace business=5.34+invnormal(uniform()) if dom==1
replace business=5.24+invnormal(uniform()) if joi==1
g political=5.27+invnormal(uniform()) if for==1
replace political=5.4+invnormal(uniform()) if dom==1
replace political=5.17+invnormal(uniform()) if joi==1
n manova overall business political=for dom
end

forv i=1(1)1000{
mono
}

这样是否可以达到目的？

谢谢啊。不过你的结果是独立的。我需要把所有的统计指标统一放在一起最后获得平均值，这样如：每次都有一个F值，这样算1000次，最后得到这1000次的均值，还可以得到显著性程度。有点类似meta-analysis的意思。有点麻烦。