- DATA f1;
- INPUT grp value @@;
- CARDS;
- 1 0.74 1 1.05 1 1.20 1 1.83 1 1.42 1 2.97 1 1.54
- 2 0.54 2 1.44 2 0.65 2 0.84 2 0.86 2 0.91 2 1.16
- ;
- RUN;
- proc multtest data=f1 nsample=1000 outperm=a1 permutation seed=112233 pvals;
- class grp;
- test mean(value);
- contrast 'a vs b' 0 1;
- run;
这个是关于permutation的一段程序,理论上重复1000次的P=0.0340.但是我用这段程序做出来的结果是不对的~