SAS随机分组方法及实现.docxVIP

SAS随机分组方法及实现

编号：研究对象（动物按体重大小、患者按预计样本量编号）从1到N编号；

获取随机数字：从随机数字表中任意一个数开始，沿同一方向顺序每个研究对象对应取一个随机数字；

求余数：随机数除以组数求余数。若整除，则取组数作为余数；

分组：按余数数值分组；

调整：假如某组待调整，该组共有n例。从中抽取1例，就取下一个随机数，随机数除以n。

???????除以n的余数(若整除则余数为n)作为在该组中所抽研究对象的序号，调整到其他组。

?

例1-1：两组

对心脑病区观察20例(编号1～20)心血管病患者分为2组，一组以灯盏花注射液为治疗组，另一组给予瓜蒌薤白汤。

从随机数字表任一行开始(以第11行第1个数(57)计)，按序查找，凡小于或等于20的数标记，查够10个数；

将与这10个数对应编号患者列为一组，余下患者为另一组。

57

35

27

33

72

24

53

63

94

09

.

41

10

.

76

47

91

44

04

.

95

49

66

39

60

04

.

59

81

48

50

86

54

48

22

06

.

34

72

52

82

21

15

.

65

20

.

33

29

94

71

11

.

15

.

91

29

12

.

03

.

61

96

48

95

03

.

07

第一组：9，10，4，6，15，20，11，12，3，7；第二组：1，2，5，8，13，14，16，17，18，19。

?

例1-2：多组(≥3组)

将15名血栓性血瘀证患者分为3组。

病人编号

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

随机数字

28

26

08

73

37

32

04

05

69

30

16

09

05

88

69

余数

1

2

2

1

1

2

1

2

3

3

1

3

2

1

3

分组

甲

乙

乙

甲

甲

乙

甲

乙

丙

丙

甲

丙

乙

甲

丙

调整

?

?

?

?

?

?

丙

?

?

?

?

?

?

?

?

第一次分组后，甲组6例，乙组5例，丙组4例。由于各组例数不等，须将甲组调整1例到丙组。

因此，继续查随机数字表，下一个随机数字为58。由于58/6=9……4，甲组中第4个研究对象调整到丙组。

?

5、SAS实现

对20例病人随机分成两个等比例组，使每组为10例。

方法一：PROC?PLAN?SEED=n过程。

PROCPLANSEED=210000;

FACTORSn=20;

OUTPUTOUT=patient;

RUN;

DATAresult;

SETpatient;

number=_n_;

IFn=10THENgroup=A;

ELSEgroup=B;

RUN;

PROCPRINTdata=resultNOOBS;

VARnumbergroup;

RUN;

?

方法二：UNIFORM(n)函数。在完全随机化时，UNIFORM函数法结果的平衡性较差。

DATApatient;

DOnumber=1to20;

r=UNIFORM(210000);

OUTPUT;

END;

RUN;

PROCRANKdata=patientOUT=rank;

RANKSr_rank;

VARr;

RUN;

DATAresult;

SETrank;

IFr_rank=10THENgroup=A;

ELSEgroup=B;

RUN;

PROCPRINTdata=resultNOOBS;

VARnumbergroup;

RUN;

?

二、区组随机化，又称均衡随机化、限制性随机化

1、定义：将随机加以约束，使各处理组的分配更加平衡，满足研究要求。在一个区间内包含一个预定的处理分组数目和比例。

??????区组：由若干特征相似的试验对象组成。如同一窝的动物、批号相同的试剂、体重相近的受试者等。

??????区组的长度：区组中对象的数目。

2、优点：区组随机化分组，避免简单随机化分组可能产生的不平衡现象，不仅提高统计学效率，而且保证分配率不存在时间趋势。

3、适用条件：区组的长度不宜太小，太小则形成不随机。一般区组的长度至少要求为组数的2倍以上。

??????????????????区组的长度也不宜太大，太大易使分段内不均衡，

?????????????????如果只有两个组别(试验组和对照组)，区组的长度一般可取?4~8，如果有4个组别则区组的长度至少为8。

??????????