http://www.pinggu.org/bbs/thread-547419-1-1.html

可能先要调入stata中然后画图。

e,搞不定，好像不知道怎么将模拟过程得到的结果赋值给某一特定的矩阵变量，比如对某一结果模拟了1000次得到一个结果，再I++了，这样又得到一个结果，这个怎么赋值给一特定矩阵呢？
现在唯一的解决办法就是对着屏幕复制黏贴，不过数据量大的话要搞死人啊，比如1000次。。。靠，一屏幕才多少数据啊。。。

可否把你的模拟过程贴出来？

哦,前几天出差,我的代码是下面的:
mata
     A=uniform(1000,1):>=0.95
     B=uniform(1000,1):>=0.9
     C=uniform(1000,1):>=0.85
for (q=1;q<=200;q++){
n=uniform(1000,q)  
k=J(1000,1,0)
for (j=1;j<=rows(n);j++) {
     for (i=1;i<=cols(n);i++) {
     if (n[j,i]>=0&n[j,i]<0.5){
      k[j,1]=k[j,1]+A[j,1]
      }
     else if (n[j,i]>=0.5&n[j,i]<0.85){
        k[j,1]=k[j,1]+B[j,1]
       }
    else if (n[j,i]>=0.85&n[j,i]<1){
       k[j,1]=k[j,1]+C[j,1]
       }
      }
   }
mean(k/cols(n))
}
end

各位同学请勿乱用......谢谢!

那么，你想画出什么样的图呢？