关于SAS中做蒙特卡洛模拟方面的问题

看到一篇文献用SAS编程做蒙特卡洛模拟进行率的多重比较方法研究，在下比较感兴趣，现发帖求SAS程序，麻烦高人予以解答，在此谢过。麻烦高人帮忙做SNK法，Bootstrap和Permutation抽样方法三种方法，同时附上所有过程和结果，鄙人不清楚之处，请多多给予解答。

方法

运用Monte Carlo方法，以二项分布的随机数函数进行重复抽样，重复次数50000次，用SAS 9.13软件编写程序，模拟计算各种方法的总I型错误率(family—wise error rate，FWER)、检验效能、判对率。

1．多重比较方法

样本率多重比较的方法

SNK法、(4)基于重复抽样的方法：Bootstrap和Permutation抽样方法。文中模拟Bootstrap抽样中采用了Freeman—Tukey Double Arcsine Test算法，简记为Boot；Permutation中运用了Fisher's exact test算法，简记为Perm。

2．评价指标

评价指标主要包括FWER、检验效能、判对率以及稳健性 。

3．模拟数据

在多重比较中，总组数为k，多重比较的次数为k(k一1)／2次，模拟数据主要考虑以下三种情形：

(1)多个总体率全部相等，考察各方法控制FWER的能力。能否将FWER的值控制在检验水准内，这是评价多重比较方法的一项重要指标。

(2)多个总体率全部不等，考察各方法的检验效能：计算其完全效能(all—pairs power)和比例效能(per—pair power)的大小。完全效能指在多重比较中发现了所有差异的能力，比例效能指差异被发现的平均概率。

(3)部分总体率相等，综合考察各方法的FWER、检验效能以及判对率的高低。判对率指对多重比较任务的推断完全正确的概率。

4．参数设置

检验水准d：0．05；率差δ=0.09；组数k取3～6、9、12；总体率取值0．09的整数倍：0．09、0.18…0.9；样本量n为5、10、20、40、80、160；两组间样本量比例为1：1、1：2、1：3、1：4。

各位大侠帮帮忙，感激不尽！！

自己顶一下，各位大侠帮帮忙，感激不尽！！

你可以在金融计量学基于SAS的金融实证研究里找到答案！

文献在哪儿啊？