基于SAS统计软件的多变量随机化分组及均衡性检验的方法

Ｊｏｕｒｎａ１　ｏｆ　Ｍａｔｈｅｍａｔｉｃａ１　Ｍｅｄｉｅｉｎｅ　Ｖｏ１．２６　Ｎｏ．４　●　２０１３　文章编号：１００４—４３３７（２０１３）０４—０３９８—０２　中图分类号：Ｒ３１１　文献标识码：Ａ　・方法评介・　基于ＳＡＳ统计软件的多变量随机化分组　及均衡性检验的方法　杨　勇　（山东中医药大学济南２５０３５５）　ｄｏｉ：１０．３９６９／ｊ．ｉｓｓｎ．１００４－４３３７．２０１３．０４．００７　在中医症候模型的实验研究中，一般采用较长时间的复　小鼠用雌二醇腹腔注射１５ｄ复制肾阳虚模型，再随机分组，然　合因素造模，动物体重等生理状态会出现较大的变化和差异。　后考察受试药对其治疗作用。课题设计检测指标包括模型动　因此，在进行药物干预之前，需要对模型动物再进行随机化分　物的体温、活动度、精神状态、免疫器官脏器系数的变化，这就　组。这就要求不能只按体重的波动进行单因素随机化样本分　需要重点考察症候模型动物的体重、体温等因素分组时是否　配，需要根据实验目的及症候的特点，选择几个重要的影响因　均衡。在考察变量的均衡性时，因这些变量都属于计量数据　素作为分组依据，比如对实验动物的体重、体温、血压、血糖等　类型，因此可直接考虑采用方差分析，ＬＳＤ检验方法进行均数　多变量指标进行重新完全随机化分组并进行均衡性检验，这　间多重比较ｌ＿２］。　才能为随后进行的研究提供平行、随机及可对照的实验设计。　本实验设计为：需要把２Ｏ只已复制成肾阳虚的小鼠再随　ＳＡＳ统计软件可较方便地实现多变量随机分组及均衡性检　机分配到模型组、阳性药对照组、中药治疗Ｉ组、中药治疗ＩＩ　验，完美地完成实验设计要求［１］。　组，动物的编号（ｂｈ）、体重（ｔｚ）、体温（ｔｗ）情况见表１。　１　多变量的实验分组设计　举例：为研究某中药复方治疗肾阳虚的效果，需先将雄性　表１　２０只肾阳虚小鼠的体重、体温数据　２　ＳＡＳ源程序　ｐｒｏｃ　ｐｌａｎ　ｓｅｅｄ＝０７５２；／＊随机种子＊／　ｆａｃｔｏｒｓ　ｕｎｉｔ＝２０；　ｄａｔａ　ａ；　ｏｕｔｐｕｔ　ｄａｔａ＝Ｍ１　ｏｕｔ＝ｂ；　ｄｏ　ｕｎｉｔ＝１　ｔｏ　２Ｏ；／＊待分组的样本数＊／　ｄａｔａ　ｃ；　ｏｕｔｐｕｔ；　ｓｅｔ　ｂ；　ｅｎｄ；　ｉｆ（ｕｎｉｔ＜＝５）ｔｈｅｎ　ｇｒｏｕｐ＝１：　ｒｕｎ；　ｅｌｓｅ　ｉｆ（ｕｎｉｔ＜一１０）ｔｈｅｎ　ｇｒｏｕｐ＝２　ｄａｔａ　Ｍ１；　ｅｌｓｅ　ｉｆ（ｕｎｉｔ＜＝１５）ｔｈｅｎ　ｇｒｏｕｐ＝３：　ｓｅｔ　ａ；　ｅｌｓｅ　ｇｒｏｕｐ＝４；／＊总共组别数＊／　ｉｎｆｉｌｅ　Ｃ：ｔｅｓｔ．ｔｘｆ；／＊导人体重、体温的纯文本数据＊／　ｏｕｔｐｕｔ；　ｉｎｐｕｔ　ｂｈ　ｔｚ　ｔｗ；　ｒｕｎ；　ｐｒｏｃ　ｓｏｒｔ；ｂｙ　ｇｒｏｕｐ；　收稿日期：２０１３—０２—０８　・　３９８　・　数理医药学杂志　ｐｒｏｃ　ｐｒｉｎｔ；　ｖａｒ　ｇｒｏｕｐ　ｂｈ　ｕｎｉｔｔｚｔｗ；／＊ｂｈ编号ｕｎｉｔ随机分配号ｔｚ体重ｔｗ体温　｜　ｐｒｏｃ　ａｎｏｖａ；　ｃｌａｓｓ　ｇｒｏｕｐ；　ｍｏｄｅｌ　ｔｚ＝ｇｒｏｕｐ；／＊方差分析，计算体重变量＊／　ｍｅａｎｓ　ｇｒｏｕｐ／ｌｓｄ；／＊采用ＬＳＤ检验，进行均数间多重比较＊／　ｐｒｏｃ　ａｎｏｖａ；　ｃｌａｓｓ　ｇｒｏｕｐ；　ｍｏｄｅｌ　ｔｗ＝ｇｒｏｕｐ；／＊方差分析，计算体温变量＊／　ｍｅａｎｓ　ｇｒｏｕｐ／ｌｓｄ；／＊采用ＬＳＤ检验，进行均数间多重比较＊／　ｑｕｉｔ；　３结果　表２随机分组结果　以上结果显示将２Ｏ个数据进行随机化分为４组。　Ｄｅｐｅｎｄｅｎｔ　Ｖａｒｉａｂｌｅ：ｔｚ　以ｔｚ（体重）为变量，结果输出包括模型（Ｍｏｄｅ１）的自由度（ＤＦ）、　平方和（Ｓｕｍ　ｏｆ　Ｓｑｕａｒｅｓ）、均方（Ｍｅａｎ　ｑｕａｒｅ）、Ｆ值和Ｐ值（Ｐｒ＞Ｆ）。　２０１３年第２６卷第４期　Ｍｅａｎｓ　ｗｉｔｈ　ｔｈｅ　ｓａｍｅ　ｌｅｔｔｅｒ　ａｒｅ　ｎｏｔ　ｓｉｇｎｉｆｉｃａｎｔｌｙ　ｄｉｆｆｅｒｅｎｔ．　ｔ分组的英文字母是用来说明各组之间差别有无统计意义，如两　组或多组字母相同，则差别无意义。本例四个组字母都是Ａ，说明组　间两两比较差别都无统计意义。这里均数已按大小重新排列。　图１　随机分组后体重ｔｚ的均衡性检验　Ｄｅｐｅｎｄｅｎｔ　Ｖａｒｉａｂｌｅ：ｔｗ　以ｔｗ（体温）为变量，结果输出包括模型（Ｍｏｄｅ１）的自由度（ＤＦ）、　平方和（Ｓｕｍ　ｏｆ　Ｓｑｕａｒｅｓ）、均方（Ｍｅａｎ　ｑｕａｒｅ）、Ｆ值和Ｐ值（Ｐｒ＞Ｆ）。　Ｍｅａｎｓ　ｗｉｔｈ　ｔｈｅ　ｓａｍｅ　ｌｅｔｔｅｒ　ａｒｅ　ｎｏｔ　ｓｉｇｎｉｆｉｃａｎｔｌｙ　ｄｉｆｆｅｒｅｎｔ．　本例四个组字母都是Ａ，说明组间两两比较差别也都无统计意义。　图２　随机分组后体温ｔｗ的均衡性检验　４结论　以上程序可方便地对多变量数据进行随机化分组，并对　每个变量进行均衡性检验，ｍｅａｎｓ　ｇｒｏｕｐ／ｌｓｄ语句要求对变量　ｇｒｏｕｐ的各水平用ｌｓｄ法作均数间的多重比较，还可根据要求　选用ｓｎｋ、ｄｕｎｃａｎ、ｓｃｈｅｆｆｅ等方法＿３　］。一般需要进行多次不同　的随机种子计算，选择Ｐ值较大，均衡性更好的数据作为随机　分组的依据。　参考文献　１周仁郁．ＳＡＳ的分层设计及均衡性检验．数理医药学杂志，２００７，　２：２２６￣２２７．　２丁元林，孔丹莉．多个样本及其两两比较的秩和检验ＳＡＳ程序．　中国卫生统计，２００２，５：５８￣５９．　３金丕焕，主编．医用ＳＡＳ统计分析．上海：上海医科大学出版社，　２０００，３３～４２．　４薛富波，等编著．ＳＡＳ　８．２统计应用教程．北京：兵器工业出版社，　２００４．】２３～１４６．　・　３９９　・