/* 研究对象随机分组 */ proc plan seed =091015; /* proc plan必有,seed参数可选,默认以及其时间为种子 */ factors rand =384 random; /* factors为关键字不可省略,rand为变量名可以任取。384表示抽样例数*/ output out =outA; /* 将随机结果输出值数据集outa */ data outB; set outA; /*链接数据集outA outB 。此时数据B中包含rand变量 */ /* 使用if语句分组 */ if rand <=192 then group =‘A‘; else group =‘B‘; /* 192以下分入A组,其余B组 */ num =_n_; /* 产生新变量num,值为顺序号 */ proc print; run;
/* 随机区组设计样本含量估计 */ proc power; onewayanova /* 表示为单因素方差分析 */ test =overall /* 三组间比较而非两两比较 */ groupmeans =61.4 |72.2 |68.7 /*一下可理解为赋值语句,等号不能省略 切记*/ stddev =11.6 power =0.8 ntotal =.; run;
原文:http://blog.csdn.net/lelieven/article/details/20032557