SAS中文论坛

标题: 常用的试验设计方法与SAS程序 [打印本页]

作者: shiyiming    时间: 2005-1-30 16:00
标题: 常用的试验设计方法与SAS程序
[b:3b5ed]常用的试验设计方法与SAS程序[/b:3b5ed]

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过程产生完全随机设计方案:

[quote:3b5ed]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;[/quote:3b5ed]
   

例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过程产生完全随机设计方案:

[quote:3b5ed]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;[/quote:3b5ed]

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过程产生配对设计方案:

[quote:3b5ed]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;[/quote:3b5ed]

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过程产生配伍组设计方案:

[quote:3b5ed]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;[/quote:3b5ed]

用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)

[quote:3b5ed]PROC PLAN SEED=987654;
   FACTOR  a=4 ORDERED   
                     b=4  ORDERED;
   TRETMENTS  c=4;
OUTPUT  OUT=latin;
RUN;[/quote:3b5ed]

用SAS  FACTEX过程产生析因设计方案

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

[quote:3b5ed]PROC  FACTEX;
  FACTORS   A   B/NLEV=2;
  EXAMIN   DESIGN;
RUN;
QUIT;[/quote:3b5ed]




欢迎光临 SAS中文论坛 (http://mysas.net/forum/) Powered by Discuz! X3.2