SAS中文论坛

标题: 跟crackman读SAS文献（2）--1《SAS程序在随机化实验设计中的应用》 [打印本页]

作者: shiyiming 时间: 2010-11-6 19:15
标题: 跟crackman读SAS文献（2）--1《SAS程序在随机化实验设计中的应用》
From crackman's blog on Baidu

《SAS程序在随机化实验设计中的应用》中药新药与临床药理2001年7月第12卷第4期
这篇文章主要讲述的是SAS PLAN过程在完全随机化设计、随机区组设计、系统(嵌套)设计、分割(裂区)设计、拉丁方设计的应用。
TITLE '完全随机设计'; 
/*建立一个数据集，产生两个变量，一个UNIT 一个GROUP*/ 
DATA CRD; 
DO UNIT=1 TO 12 ; 
IF(UNIT<=6) THEN GROUP=1; 
ELSE GROUP=2; 
OUTPUT; 
END; 
RUN; 
PROC PLAN SEED=27371; /*随机种子为27371*/ 
FACTORS UNIT=12; /*试验处理的因子，可以这样理解，就是你要随机分配的对象，例如，12只老鼠如何随机分配 
在两组里面，这里UNIT就是你要分配对象---老鼠的数量*/ 
OUTPUT DATA=CRD OUT=OUTDAT; /*这里如果去掉DATA=CRD和不除去DATA=CRD的差别，可以看看一下。DATA=CRD，是说明 
原始的UNIT对应GROUP，但是不符合随机分配的要求，所以要加入PLAN过程，对CRD数据集中的UNIT的分配进行随机化，将随机 
分配的结果输出保存在数据集OUTDAT中。等于说，你先用一个房间把兔子都放进去，那么UNIT有了，脑壳一亮靠近我的兔子1 
远离我的兔子是2，那么GROUP有了，所以这个房间相当于CRD，但是不是随机化，用PLAN对房间的兔子来个随机化吧*/ 
PROC PRINT; 
RUN; 
TITLE; 
以下程序也可以实现上述功能：
PROC PLAN SEED=27371; 
FACTORS UNIT=12; 
TREATMENTS GROUP=12 CYCLIC<1 1 1 1 1 1 2 2 2 2 2 2>; 
OUTPUT OUT=OUTDAT; 
PROC SORT; 
BY UNIT; 
PROC PRIN; 
RUN; 

关于PLAN的过程参数大家可以看看SAS MANUAL:
<a target="_blank" href="http://support.sas.com/documentation/cdl/en/statug/63347/HTML/default/viewer.htm#plan_toc.htm">support.sas.com/documentation/cdl/en/statug/63347/HTML/default/viewer.htm#plan_toc.htm</a> <a href="http://hi.baidu.com/crack%5Fman/blog/item/a0fdc039934c2dd79e3d62e8.html">阅读全文</a>

 类别：<a href="http://hi.baidu.com/crack%5Fman/blog/category/%B8%FAcrackman%B6%C1sas%CE%C4%CF%D7">跟crackman读sas文献</a> <a href="http://hi.baidu.com/crack%5Fman/blog/item/a0fdc039934c2dd79e3d62e8.html#comment">查看评论</a>

欢迎光临 SAS中文论坛 (https://mysas.net/forum/)