按孙振球《医学统计学》197页表11-29模式设计(程序已经安排好将观测随机引入分组)
title 'wxq1638提供赞助!!!';
%macro wxq(a,n,b);/*完全随机3*2裂区设计分组*/
data a;do no=1 to &a*&b;/**/x=substr(normal(0),8,1);output;end;proc sort;by x;run;
data b;set a(keep=x);call symput('seed',x);run;%let wxq168=%sysevalf(&seed);
proc plan seed=&wxq168; factors a=&a ordered n=&n b=&b; output out=wxq88;run;
data wxq1688;set a;do i=1 to &n;keep no;output;end;
data wxq1638; merge wxq88 wxq1688;proc print;run;quit; %mend wxq(a,n,b);%wxq(3,2,2);
结果
wxq1638提供赞助!!! 2009年01月01日 星期四 下午08时28分33秒 1
The PLAN Procedure
Factor Select Levels Order
a 3 3 Ordered
n 2 2 Random
b 2 2 Random
a n -b-
1 2 2 1
1 1 2
2 2 2 1
1 1 2
3 2 2 1
1 1 2
wxq1638提供赞助!!! 2009年01月01日 星期四 下午08时28分33秒 2
Obs a n b no
1 1 2 2 1
2 1 2 1 1
3 1 1 1 5
4 1 1 2 5
5 2 2 2 2
6 2 2 1 2
7 2 1 1 4
8 2 1 2 4
9 3 2 2 3
10 3 2 1 3
11 3 1 1 6
12 3 1 2 6