看到一篇文献用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。 各位大侠帮帮忙,感激不尽!! |