各位好。“高级医学统计学_万崇华2014”一书中,在关于平衡不完全区组设计的SAS程序里,有一小段宏命令,我仿照着输入,却无法得到结果。在此请教大家,麻烦帮我看一看,不胜感激!
%let block = 5;
%let treat = 5;
%let subject = 3;
%macro aa( treat);;
data a; %do i = 1 %to & treat;
treatment = compress('疗法 '||& i);;
output; %end;
run;
%mend aa;
%aa( &treat);;
proc optex data = a seed = 10000 coding = orth;
class treatment; model treatment; blocks structure = ( &block)&subject;
output out = b;
run;
data c;
set b; number = _n_;
run;
ods html style analysis;
proc tabulate format = 2.0;
class block treatment; var number;
table block = ", number = '受试对象编号 '*treatment = "* sum = " /box = [label = '区组'];
run;
ods html close;


雷达卡






京公网安备 11010802022788号







