巧用SAS软件制作临床试验随机盲表
巧用S AS 软件制作临床试验随机盲表
长沙市中医院(410002) 邓巧君 复旦大学药学院(200032) 魏志军
摘要 临床试验通常需要制作各种随机盲表,如何准确、科学地制成符合临床试验需要的随机盲表是临床试验设计阶段急需解决的问题。本文拟通过一个临床试验的具体例子详细介绍利用S AS 软件的P LAN 过程及其分栏显示功能如何达到这一目的。
主题词 临床试验 软件设计
S AS 软件是全球久负盛名的最权威的数据分析软件之一,具有完备的数据访问、数据管理、数据分析和数据呈现功能。目前已广泛应用于生物制药行业的生物统计分析及临床试验设计中[1]。现笔者对临床试验设计中如何应用
S AS 软件情况介绍如下。
随机原则是临床试验设计的三原则之一(另两个原则
为对照、重复),因而在具体的临床试验设计中,我们会经常遇到随机盲表的制作问题。一个准确、科学的随机盲表是保证临床试验成功的关键之一,而传统上通过查找随机数字表来进行随机编码,显然已不能满足临床试验,尤其是大型临床试验的要求。笔者将通过一个多中心的随机、双盲、对照的新药Ⅱ期临床等效试验的随机盲表制作例子,具体介绍如何利用S AS 软件的PROCP LAN 过程及其分栏显示功能快速、准确地制成临床试验所需的随机盲表。根据临床试验有关要求,新药的Ⅱ期临床等效试验通常安排在多个中心进行(本例假定3个),样本含量要求为至少60对即
120例(假定为1个适应症),处理组分别为治疗组与对照
组。在进行此类临床试验设计时,为了让不同组别的受试者分配均衡及减少季节等时间因素的可能影响,通常进行区组的随机化设计(又称分段随机),即先将受试者分成若干例数相等的区组,然后再在每区组内按简单随机化原则进行分配[2]。通常将受试者分成4人一区组(也可为2,6,8…但需为处理组的倍数),每组随机分配两人入A 组、两人入B 组,共有AABB 、ABAB 、BBAA 、BAAB 、ABBA 、BABA 等6种组合形式。利用S A 1S 软件的PROC P LAN 过程即可轻松实现这样的临床试验设计。程序如下:/3-2
-Pr ogra m f or
Random ized Comp lete B l ock Design -3
/op ti on nodate nonu m 2
ber p s =40ls =100;title’xxx 临床试验随机盲表’;title2’======;%let seedno =20050101;/3-利用宏设定随机数
的初值(种子数)-3/3-p r oc P LAN 过程-3/p r oc p lan
seed =&seedno;/3调用宏“seedno ”-3/3-本试验有三个
因子:中心、区组及组别,其中中心、区组水平数依顺序产生,组别则随机安排-3/fact ors center =3ordered bl ock =10
ordered gr oup =4random /nop rint;/3-将结果输出至临时
文件“OUT ”中-3/out put out =out gr oup cvals =(’A ’’A ’’B ’’B ’);/3-变量GROUP 表示处理组别-3/Run 。以上程序运行后将依据3个中心,每中心10个区组,每区组
4个病例,产生120个随机数,其中1～40号分配给中心1、
41～80号分配给中心2、81～120号分配给中心3。由于加
了NOPR I N T 及OUTP UT 选项,结果将不在OUTP UT 窗口中显示,而是直接输出至临时数据文件OUT 中。但此时生成的结果不能以我们实际工作中需要的样式显示。为了让生成的结果更符合实际使用的需要,我们需要结合运用S AS 程序的P UT 语句及其分栏显示功能。具体程序如下[3];data
gr oup;set out;code +1;/3-变量CODE 表示随机编码-
3/data .null;file p rint n =pagesize;/3-指定结果显示在
OUTP UT 窗口-3/do r ow =1t o 12;/3-设定随机表的行
数-3/do col =1,8,15,22,29,36,43,50,57,64;/3-设定随机表的栏数及显示位置-3/set gr oup;put #r ow @col
code 5.+1gr oup ＄4.;/3-设定显示内容及格式-3/end;end;put;put @4“种子数=&seedno ”;/3-显示随机
数的种子数-3/put @4“-”;put page;run .
程序运行结果如下表显示。用户可根据实际使用习惯,对上述程序稍作修改即可生成其他式样的随机盲表,同样其他临床试验的随机盲表亦可照此方法生成(以上程序在
S AS V8.2以上通过)。
XXX 临床试验随机盲表
1A 2B 3A 4B 5B 6A 7B 8A 9A 10A 11B 12B 13A 14B 15A 16B 17B 18A 19B 20A 21A 22A 23B 24B 25B 26A 27B 28A 29A 30A 31B 32B 33B 34B 35A 36A 37A 38B 39A 40B 41A 42B 43B 44A 45A 46A 47B 48A 49A 50A 51B 52B 53A 54B 55A 56B 57A 58B 59A 60B 61A 62A 63B 64B 65B 66A 67A 68B 69B 70A 71A 72B 73B 74B 75A 76A 77B 78A 79A 80B 81A 82B 83A 84B 85A 86A 87B 88B 89A 90A 91B 92B 93A 94B 95A 96B 97B 98A 99B 100A 101B 102A 103A 104B 105B 106A 107A 108B 109A 110B 111A 112B 113A 114B 115B 116A 117B 118A 119B 120A
注:种子数=20050101
参考文献
[1] 高惠旋.S AS 系统S AS/ST AT 软件使用手册[M ].北京:中国统
计出版社,2001.
[2] 王瑞莲.新药临床研究实用手册[M ].北京:化学工业出版社,
2003.
[3] htt p://support .sas .com /rnd /sa mp les/index .ht m l .
(收稿日期:2005-03-23)
?
621?第21卷第3期2005年5月湖南中医杂志HUNAN JOURNAL OF T RAD I TI O NAL CH I N ESE M E D I C I N E Vol 121 No 13
May 2005