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的模拟相关统计指标呢？？