现在的位置: 首页 > 综合 > 正文

常用的试验设计方法和SAS程序

2018年10月24日 ⁄ 综合 ⁄ 共 1950字 ⁄ 字号 评论关闭

摘自:

http://www.mysas.net/forum/viewtopic.php?f=4&t=1617

1.完全随机设计
例1:将10例实验对象随机地分配到A、B两组,先将实验对象按某项特征编号,然后从随机分配(数字)表中任意指定一行,舍掉10―19的数字,规定随机数字为偶数者分到A组,奇数者分到B组。

编 号:1 2 3 4 5 6 7 8 9 10
随机数:6 1 5 4 0 7 8 3 9 2
组 别:A B B A A B A B B A
结果将1、4、5、7、10号实验对象分到A组,
2、3、6、8、9号实验对象分到B组。

用SAS PLAN过程产生完全随机设计方案:
PROC PLAN SEED=385679;
FACTOR random=10;
OUTPUT OUT=ran_1;
DATA ran_2;
SET ran_1;
IF MOD(random,2)=0
THEN treat=“a”;
ELSE treat=“b”;
PROC PRINT;
RUN;

例2:将15例实验对象随机地分配到A、B、C三组,先将实验对象按某项特征编号,然后从随机分配(数字)表中任意指定一行,舍掉15―19的数字,随机数字为0―4分到A组,5―9分到B组,10―14分到C组。
编 号:1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
随机数:2 8 14 6 11 1 9 0 10 4 7 13 3 12 5
组 别:A B C B C A B A C A B C A C B
结果将1、6、8、10、13号实验对象分到A组,
2、4、7、11、15号实验对象分到B组,
3、5、9、12、14号实验对象分到C组。

用SAS PLAN过程产生完全随机设计方案:
PROC PLAN SEED=43697;
FACTOR random=15;
OUTPUT OUT=ran_3;
DATA ran_33;
SET ran_3;
IF random <6 THEN treat=“a”;
ELSE IF 5<random <11 THEN treat=“b”;
ELSE treat=“c”;
PROC PRINT;
RUN;

2.配对设计
将实验对象按一定的条件配成一对,如将年龄、性别、体重等条件一致的实验对象配成对子,再随机分配表每一对中的两个对象接受不同的处理,配对设计可增强处理间的均衡性,提高实验效率。
例:将10对实验对象随机地分配到两种处理。
编 号:1.1 2.1 3.1 4.1 5.1 6.1 7.1 8.1 9.1 10.1
1.2 2.2 3.2 4.2 5.2 6.2 7.2 8.2 9.2 10.2
随机数: 9 3 0 2 1 5 8 6 4 7
处 理:BA BA AB AB BA BA AB AB AB BA
规定随机数为偶数者按A、B分配,随机数为奇数者按B、A分配

用SAS PLAN过程产生配对设计方案:
PROC PLAN SEED=976521;
FACTOR random=10;
OUTPUT OUT=ran_11;
DATA ran_21;
SET ran_11;
IF MOD(random,2)=0
THEN treat=“ab”;
ELSE treat=“ba”;
PROC PRINT;
RUN;

3.配伍组设计(随机区组设计)
配伍组设计可看成是配对设计的扩大,将条件相同的实验对象配成一组,将每一配伍组内的受试对象随机地分配到不同的处理中去。
例:将20个受试对象分成五个配伍组,每一配伍组分别接受A、B、C、D四种处理。
配伍组: 1 2 3 4 5
编 号:1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
随机数:3 2 1 4 4 3 2 1 3 1 2 4 3 1 2 4 3 2 1 4
处理组:CBAD DCBA C A B D C A B D C B A D

用SAS PLAN过程产生配伍组设计方案:
PROC PLAN SEED=123456;
FACTOR block=5 ORDERED treat=4;
OUTPUT OUT=block;
DATA b_d;
SET block;
IF treat=1 THEN GROUP=“A”;
IF treat=2 THEN GROUP=“B”;
IF treat=3 THEN GROUP=“C”;
IF treat=4 THEN GROUP=“D”;
PROC PRINT;
RUN;

用SAS PLAN过程产生拉丁方(Latin Square)设计方案
例:一个4×4拉丁方设计的方案。
研究不同的蛇毒成分、不同的蛇毒剂量对不同的肿瘤的抑瘤效果。
A因素(瘤株): A B C D
B因素(成分): I II III IV
C因素(剂量): 0 1.5 2.0 2.5(mg/g)
PROC PLAN SEED=987654;
FACTOR a=4 ORDERED
b=4 ORDERED;
TRETMENTS c=4;
OUTPUT OUT=latin;
RUN;

用SAS FACTEX过程产生析因设计方案
例:一个2×2析因设计方案
A药:a1(用) a2(不用)
B药:b1(用) b2(不用)
共产生四种组合:
a1b1, a1b2, a2b1, a2b2

PROC FACTEX;
FACTORS A B/NLEV=2;
EXAMIN DESIGN;
RUN;
QUIT;

抱歉!评论已关闭.