SAS随机分组方法和实现.docVIP

专业资料 word完美格式 随机分组方法包括: 简单随机化(simple randomization) 区组随机化(block randomization) 分段(或分层)随机化(stratified randomization) 分层区组随机化(stratified block randomization) 动态随机化(dynamic randomization) ? 一、简单随机化，又称完全随机化 1、定义：在事先或者实施过程中不作任何限制和干预或调整，对研究对象直接进行随机分组。 ???????????? 通常，通过掷硬币、随机数字表、计算机产生随机数来进行随机化。 2、适用条件：在研究例数较少、总体中个体差异较小时，采用此法。 3、缺点：在研究对象例数较少时，由于随机误差难以保证组间病例数的均衡，各组例数可能会出现不平衡现象。 4、解决办法：随机数表法、随机数余数分组法。 随机数余数分组法的具体操作： 编号：研究对象（动物按体重大小、患者按预计样本量编号）从 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 过程。 PROC PLAN SEED=210000; FACTORS n=20; OUTPUT OUT=patient; RUN; DATA result; SET patient; number=_n_; IF n=10 THEN group=A; ELSE group=B; RUN; PROC PRINT data=result NOOBS; VAR number group; RUN; ? 方法二：UNIFORM(n)函数。在完全随机化时，UNIFORM函数法结果的平衡性较差。 DATA patient; DO number=1 to 20; r=UNIFORM(210000); OUTPUT; END; RUN; PROC RANK data=patient OUT=rank; RANKS r_rank; VAR r; RUN; DATA result; SET rank; IF r_rank=10 THEN group=A; ELSE group=B; RUN; PRO