新药临床试验数据中计数资料统计分析的SAS宏实现

新药临床试验数据

’#’

!

设计 统计 方法

!

中国临床药理学与治疗学

中国药理学会主办

!"#$%&’()*+，,--"&((.%’/(&

0%1234：5567.)8’&59:5;1

（’）：’(()<=>；&&’#’?’#)

新药临床试验数据中计数资料统计分析的!"!宏实现

李婵娟，薛富波，于莉莉，徐丽，王素珍，王陵，夏结来

第四军医大学卫生统计学教研室，西安A&((#’，陕西

摘要

基于-H-软件中的<+0I与!J++过程，设

计新药临床试验中计数资料统计分析的-H-宏，结合实例对其各项参数的含义及其可实现的功能做了详细说明。通过调用此-H-宏程序，能直接产生出各种类型计数资料的-H-统计分析报表，实现了分析内容的客户定制化以及结果输出的规范化，大大提高了新药临床试验中针对该类数据统计分析过程的工作效率和质量。关键词

临床试验；计数资料；-H-宏；

统计报表中图分类号：+.).文献标识码：H

文章编号：&((.%’/(&（’(()）(’%(’#’%(/

-H-（-7273G73524H924KG3G-KG7=1，统计分析系统）是当今国际上著名的数据分析软件系统，是目前国内外新药临床试验中数据管理和统计分析的首选。然而-H-软件输出结果项目繁多，内容复杂，不能满足新药临床试验对统计分析报表简洁直观的要求，而从JLBMLB中用人工的方法挑选所需要的结果再制表，工作十分烦琐而且容易产生错误。工作效率和工作质量的双重要求迫使越来越多的生物统计专业人员通过编制-H-宏来实现新药临床试验数据

处理工作的自动化、流程化和规范化［&?#］

。为适应

新药临床试验中对数据处理的特殊要求，我们应用-H-提供的宏功能，

编制了针对计数资料统计分析的宏程序，直接生成针对各种类型计数资料的统计分析报表，无需任何加工处理直接用于统计分析报告，大大提高了该类数据处理操作的工作效率。

’((/%&’%(/收稿’(()%(&%($修回国家自然科学基金资助项目（@#(’A&&/A）李婵娟，女，博士生，主要从事临床医学数据统计学处理方法研究。B=4：(’.%C$AA’&C(0%1234：5D29EF29438G;DF:5;1夏结来，通讯作者，男，博士，教授，博士生导师，主要从事新药临床试验与统计学方法研究。

B=4：(’.%C$AA$C)&万方数据0%1234：9;G28’&59:5;1

#

!"!宏设计方案

#$#

设计思想

在新药临床试验中，计数资料包括

双向无序、单向有序、双向有序属性相同和双向有序属性不同四种类型。针对不同的研究目的和不同的资料性质，统计分析应采用不同的方法。如果存在分层因素，那么计数资料的统计分析还要考虑对分层因素的校正。概括来讲，计数资料的统计分析方

法包括分层校正的双向有序表!NO卡方

［$］

、分层校正的单向有序表!NO卡方、分层校正的双向无序表!NO卡方、

双向有序表!NO卡方、单向有序表!NO卡方、双向无序表!NO卡方（四格表时的M=2PG;9卡方）、直接概率法、等级相关等C种方法。对计数资料的统计分析主要采用-H-软件的<+0I过程和!J++过程

［/?A］

，常规的方法是直接调用其过程，根

据不同的资料性质选择不同的参数，但这些过程十分复杂，参数很多，运行结果繁杂，给使用者造成了极大的不便。通过编制-H-宏，控制宏参数可以达到控制分析的数据集、分析的变量、采用的分析方法等目的，使得临床试验数据处理中相关的操作趋于灵活和高效，实现了分析内容的客户定制化以及结果输出的规范化。#$%

实现过程

针对-H-提供的<+0I和!J++

过程，将其中的关键选项或参数转化为-H-宏的输入参数，通过对宏参数的定制及控制来实现用户对<+0I和!J++过程调用，

采用不同的统计分析方法来实现不同的分析目的。利用-H-的JQ-功能和-H-宏工具提供的灵活的条件控制语句，提取统计分析报表中所需的<+0I和!J++过程的参数估计值和关键统计量以及其他指标值，并对其进行适当的裁剪和拼接，然后采用输出控制语句，输出用户所需的统计报表。以分析结果中的参数估计值和关键统计量为依据，自动生成并输出相应的文本，对数据处理结果的相应统计学结论进行概括的描述和说

新药临床试验数据

中国临床药理学与治疗学&??PS+T；（&）%%

&EE

；.+3812，8+/2）…

./0*193/3+.+:39；…

-.+:2)*+；

各参数的具体含义及其设置方式见表%。

明，输出结果一目了然。

!"#"宏参数设置

以)*+命名该,!,宏，宏的定义及其参数设置方

式如下：

（2/3/%，-./0*1)*+2/3/&，4/*，4/*).3，5*167，5).3，0，

表$

%&’宏参数设置参数

2/3/%2/3/&

意义

数据集%名称数据集&名称

设置方式说明

仅输出一个数据集的分析结果，反之并列输2/3/%和2/3/&相同时，

出两个数据集统计分析结果173=1:

?不打印行标目和行合计

不打印行合计%打印行标目，

（默认）@打印行标目和行合计

如性别

方式一：性别;9+A).3B

男’&D’女’；4/C6+9+A).3%D’

方式二：9+A).3B

性别’%D’男’&D’女’；4/C6+9+A).3’@@@D’173=1:

%不打印列百分比

（默认）?打印列百分比

例如4/C6+5).3%D’组间’’@@@D’分组’；!’&D’"’%&D’如中心0:分层变量格式0:).3B

中心’；4/C6+0:).3’@@@D’

173=1:%仅用于四格表F=9G%D$%;$&F=9G&D(%;(&!*%D$%’$&!*&D(%’(&

H73=1:&仅用于四格表当173=1:%指定时

（非劣效下限或等效下限）II@J-KL下限

H73=1:E仅用于四格表当173=1:&指定时

（等效上限）II@J-KL上限

.+3812

%分层校正的双向有序表KMN卡方&分层校正的单向有序表KMN卡方E分层校正的双向无序表KMN卡方O双向有序表KMN卡方J单向有序表KMN卡方P双向无序表KMN卡方（四格表时的$+/*91:卡方）Q直接概率R等级相关?不做检验

8+/2

?输出表头%正常输出

&输出表头和表尾’%输出表尾173=1:

对行变量的文字性描述，输出统计学结论。例如两组间有效率

4/*;<173=1:>行变量名及其控制选项

4/*).3;<173=1:>行变量描述或变量格式

5*167;<173=1:>5).3;<173=1:>

分组变量名及其控制选项分组变量描述或变量格式

!#’

$%(%

"$&(&

0;<173=1:%>

<173=1:&><173=1:E>

分层变量名及其控制选项

.+3812检验方法

8+/2;<173=1:>

统计表格控制及其输出结论控制选项

万方数据

新药临床试验数据

29M

（2）HGC8RH+C8EG(#’(%$+SG-#2;;?B-T；//

!

!"#

实例分析

实例#

一项!!期临床试验，对其一般情况部

；.#$,"，.&’)>，%8，?，/）

<&#-（=()(&(6，=()(""6，(8)，(8)&’)>，.#$,"，.&’)>，%8，A，；4/）

参数说明如下：

（B,++18(+:6C6=()(&(6表示该项临床试验的全分析集。D-)，B1D）

（E-#E#$)$F=()(""6表示该项临床试验的符合方案集。%$+D-)，EED）

6-7、%+(66、.’、(8)分别为待分析的变量名称，6-7&’)、%+(66&’)、:-68$、(8)&’)是其对应的格式。

.#$,"、%8分别为临床试验的分组变量名和中心变量名，.&’)是.#$,"对应的格式。

由于病情程度(8)是单向有序表资料，所以宏参数’-)G$=取A（单向有序表HIJ卡方），其他分析变量的宏参数’-)G$=取?（双向无序表HIJ卡方）。

宏参数G-(=分别从;到4/，表示拼接四个变量的分析结果输出一个完整的统计表。

提交运行上述程序，所得分析结果如下（表2、。9）

分指标进行分析的程序如下："#$%&$#’()；*(+,-.&’)分组’；*(+,-6-7&’)*(+,-%+(66&’)入院情况’；*(+,-(8)&’)*(+,-:-68$#,8；

<&#-（=()(&(6，=()(""6，6-7，6-7&’)>，.#$,"，.&’)>，%8，?，；;）

<&#-（=()(&(6，=()(""6，%+(66，%+(66&’)>，.#$,"，.&’)>，；%8，?，/）

青霉素类过敏@:-68$>，<&#-（=()(&(6，=()(""6，.’，

表$

一般情况项目性别

分类男!（<）女!（<）合计

入院情况

住院!（<）门诊!（<）

合计

青霉素类过敏

否!（<）是!（<）合计

病情程度

轻度!（<）中度!（<）（<）重度!

合计

B1D

1处理（A?>?L）?K

（M9>99）A2

/2;（92>A;）95

（?L>A;）K/

/2;（5K>99）//K（/>?L）2/2;（MA>;;）AM

（A2>A;）?9（2>A;）9/2;

/0’1处理’20’3处理’45550’男’20’女’45550’性别’；/0’

住院’20’门诊’/0’

45550’

轻度’20’中度’90’重度’/0’

是’;0’否’；/0’

病情程度’；45550’

EED

3处理（A2>A;）?9

（ML>A;）AL

/2;（9;>K9）9L

（?5>/L）K9

/2;（5K>99）//K（/>?L）2/2;（MA>K9）AA

（A/>?L）?2（2>A;）9/2;

1处理（A?>M/）??

（M9>A5）A/

//L（9/>?2）9L

（?K>9K）K;

//L（5K>25）//A（/>L/）2//L（M?>/A）AM

（A2>/M）?/（/>L/）2//L

3处理（A9>;M）?/

（M?>5?）AM

//A（25>AL）9M

（L;>M9）K/

//A（55>/9）//M（;>KL）///A（M?>5?）AM

（A/>9;）A5（/>LM）2//A

表!一般情况组间比较结果变量性别入院情况

检验方法HJ!FDNO1PQHJ!FDNO1PQHJ!FDNO1PQ单向有序HIJFHJ!

B1D

统计量;>M2;;>;LL;>;;;;>;/A

";>A/?K;>LK/M/>;;;;;>5;/;

统计量;>2?A;>//?;>92;;>;/M

""#

";>?;?A;>L99L;>AL/9;>5;L/

青霉素类过敏病情程度

以上程序运行结果均是由D1D系统自动生成万方数据

并保存在特定文件中，是新药临床试验统计报表的

新药临床试验数据

中国临床药理学与治疗学5""LF1W；（5）%%

5?E

；’=;>%"#$%#%，?，5<细菌培养转阴率）

参数说明：

符合:,-,*,;、:,-,&&;、2’(0&、)6分别是全分析集、方案集、分组和中心的变量名称。宏参数@1-A(:4?表示对两组进行分层校正的双向无序表B@C卡方检验。

宏参数A1,:45表示输出表头和表尾。率<’=;>%"#$%#%表明对资料进行等效性检验，比（试验组<对照组）等效区间为（"#$，。%#%）

<细菌培养转阴率表明对细菌培养转阴率直接输出统计分析结论。

提交运行上述程序，所得分析结果如下（表D、：E）

规范形式，无需进行任何手工编辑。!"#

实例#

一项多中心!!期临床试验，主要疗效

指标为细菌培养转阴率，方案要求进行等效性检验，，其分析的程序如两组率比的等效区间为（"#$，%#%）下：

&’()*(’+,-；

分组’%4’试验组’54’对.,/012*+-3$$$4’照组’；

中心’；.,/01)6*+-3$$$4’

细菌培养’%4’转阴’54’.,/0178&*3$$$4’未转阴’；

’06；

9*’1（:,-,*,;，:,-,&&;，78&，78&*#，2’(0&，2*+-#，)6<

表$

细菌培养转阴率细菌培养转阴未转阴

数量（9）!（9）!合计

FGH

试验组（JK#L%）$$

（%5#?$）%D

%%?

IIH

对照组（JJ#E"）%""（%%#E"）%?

%%?

试验组（JK#?$）$K

（%5#L%）%D

%%%

对照组（JK#$L）$E

（%5#"D）%?

%"J

表%细菌培养转阴率组间比较结果变量细菌培养O’1;/(PQR,S

检验方法（$E9B!）NMM

齐性检验

"#$

统计量（"#$"，"#$$%#"$）

%4"#$$JJ

%"#J"JJ

统计量（"#$"，"#$$%#%"）

%4"#$$L?

%%$

%"#$%KL

检验结果考虑了中心效应；等效性检验界值，方案中拟定的率比（试验组<对照组）等效区间为（"#$，N：MM4试验组率<对照组率；注：%#%）

主要统计学结论：组间细菌培养转阴率比较转阴率FGH（IIH）：试验组JK#L%9（JK#?$9）、对照组（JK#$L9），试验组<对照组率MM（$E9B!），JJ#E"9

（"#$"，FGH结果为"#$$%#"$），IIH结果为"#$$（"#$"，，两组间等效（表D、。%#%"）E）

以上程序运行结果及统计学结论均是由HGH系统自动生成并保存在特定文件中，无需进行任何手工编辑。

宏，很大程度上实现了新药临床试验数据处理工作的自动化、流程化和规范化，从而使得工作效率和工作质量得到双重的保证。

本文介绍的针对计数资料处理的宏*’1是我室针对新药临床试验数据管理与处理开发的一整套（打包形成7=,;+,)’(宏库）中的功能单元之HGH宏

一。宏*’1能实现计数资料的各种常用统计分析方法，适用于新药临床试验中一般情况、疗效评价及安全性评价中的各类计数资料的统计描述和分析，以及等效性检验、非劣效检验等。只需编制少数几行指定所需的相关参数，调用方式简单，可HGH程序，

以和本室编制的定量资料统计分析HGH宏+1,6交替拼接直接输出最终的统计分析报表，而不需要任何人为的中间加工处理，可以大大减少人为产生的错误，保证结果的准确性，并能减少分析过程占用时间，从而为整个临床试验过程节省宝贵的时间。

$讨论

在新药临床试验数据管理与分析中，HGH软件

功能虽然强大，但由于其统计输出结果过于繁杂，需要用人工的方法从TUVIUV中挑出所需要的结果，工作十分烦琐而且容易产生错误。因此如何兼顾工作效率和工作质量就成为数据管理和统计分析工作者必须考虑的关键问题之一。通过编制针对不同数

万方数据据处理功能单元（包括数据管理和统计分析）的HGH

新药临床试验数据

&(,

（&）0I;9%0J;9FI5AX5=EJ?I7A&’’,R7S；!!

*,’

史周华，汪涛，刘勤等"012统计分析方法在多中心试验［%］&3&表资料的应用"中国卫生统计，!...；!,：.,*+高慧璇，李贵斌，耿直，主编"#$#系统4567#$#软件使用手册［1］中国统计出版社，"北京：&’’!,+

［0F>6B<;C7D7A6;E9/"&#$#896:;:<:789="#$#>#?$?@67A’0G］"HEA:I05AEJ;95"!...［0F>0G］#$#896:;:<:789="#$#15=AEK7L7A79=7D7A6;E9/"&"HEA:I05AEJ;95"!...

参考文献

!&

童新元，张高魁，姚晨"定性指标的#$#统计分析报表

［%］"中国卫生统计，&’’(；&’：)’*!邹建东，熊宁宁，卜擎燕，等"四格表指标统计分析报表的#$#宏程序［%］"中国临床药理学与治疗学，&’’)；!’：(

()+*,’童新元，张高魁，姚晨"新药临床试验中协方差分析的［%］#$#统计报表"中国临床药理学与治疗学，&’’-；.：.)/

-)

!"!"#$%&’(’&)*$*+)*+%$,$-$,.)/)’(%’0-*1$*$+-%,+-+%$,*&+$,’(-/21&03

M80I59NO<59，P@QR<NSE，T@M;NJ;，P@M;，U$HB#<NVI79，U59WM;9W，P8$%;7NJ5;

!"#$%&’"(&)*+"$,&-.&$&/0&/10，2)3%&-4/,/&$%54"6/1$,7(/8"%0/&5，9/’$(:;<<=>，.-$$(?/，@-/($!4"56!75

$#$#X5=AELEA595JYV;9W=E<9:C5:5;6

=E<9:C5:5，=<6:EX;V7:I7=E9:79:ELC5:5595JY6;6，59C6:59C5AC;V7;:6E<:\<:59C;X\AE_7;:6595JYV;9W7LL;=;79=Y59C‘<5J;:Y;9=J;9;=5J:A;5JEL97^CA<W"

89:;<6="=J;9;=5J:A;5J；=E<9:C5:5；#$#X5=AE；6:5:;6:;=5JA7\EA:LEAX6

C76;W97CS567CE9RKQZ59C0[KK\AE=7C<A76"Q5=I\5A5X7:7A59CL<9=:;E9EL;:;67]\J5;97C;9C7:5;JC7\79CN;9WE96EX77]5X\J76"U;:I:I;6X5=AE，^7=59C;A7=:JY\AEC<=7:I7A7\EA:LEAX6EL#$#6:5:;6:;=5J595JY6;6EL

万方数据

新药临床试验数据

新药临床试验数据中计数资料统计分析的SAS宏实现

作者：作者单位：刊名：英文刊名：年，卷(期)：被引用次数：

李婵娟， 薛富波， 于莉莉， 徐丽， 王素珍， 王陵， 夏结来， LI Chan-juan， XUEFu-bo， YU Li-li， XU Li， WANG Su-zhen， Wang Ling， XIA Jie-lai第四军医大学卫生统计学教研室,西安,710032,陕西

中国临床药理学与治疗学

CHINESE JOURNAL OF CLINICAL PHARMACOLOGY AND THERAPEUTICS2006,11(2)1次

参考文献(7条)

1.童新元.张高魁.姚晨 定性指标的SAS统计分析报表[期刊论文]-中国卫生统计 2003(1)

2.邹建东.熊宁宁.卜擎燕 四格表指标统计分析报表的SAS宏程序[期刊论文]-中国临床药理学与治疗学 2005(10)3.童新元.张高魁.姚晨 新药临床试验中协方差分析的SAS统计报表[期刊论文]-中国临床药理学与治疗学 2004(0q)4.史周华.汪涛.刘勤 CMH统计分析方法在多中心试验2×2表资料的应用[期刊论文]-中国卫生统计 1999(2)5.高慧璇.李贵斌.耿直.主编 SAS系统Base SAS软件使用手册 20016.SAS Institute Inc SAS/STAT User's Guide Version 8.2 19997.SAS Institute Inc SAS Macro Reference Version 8.2 1999

本文读者也读过(10条)

1. 邹建东.熊宁宁.卜擎燕.蒋萌.刘芳.ZOU Jian-dong.XIONG Ning-ning.BO Qing-yan.JIANG Meng.LIU Fang 四格表指标统计分析报表的SAS宏程序[期刊论文]-中国临床药理学与治疗学2005,10(3)

2. 王杨.李卫.成小如.贾宣.胡泊.WANG Yang.LI Wei.CHENG Xiao-ru.JIA Xuan.HU Bo 临床试验统计分析报告通用的SAS生成方法[期刊论文]-中华疾病控制杂志2008,12(6)

3. 邹建东.熊宁宁.卜擎燕.蒋萌.刘芳 正态分布定量指标统计分析报表的SAS宏程序[期刊论文]-中国临床药理学与治疗学2004,9(7)

4. 童新元.张高魁.姚晨 新药临床试验中定量指标的SAS统计报表[期刊论文]-军医进修学院学报2005,26(1)5. 姚利香.王文辉.曲峰.刘玉娥.YAO Li-xiang.WANG Wen-hui.QU Feng.LIU Yu-e Proc tabulate过程在临床试验统计报告中的应用[期刊论文]-中国临床药理学与治疗学2005,10(12)

6. 于莉莉.薛富波.李婵娟.王素珍.徐丽.王陵.夏结来.YU Li-li.XUE Fu-bo.LI Chan-juan.WANG Su-zhen.XU Li.WANG Ling.XIA Jie-lai 新药临床试验中定量指标常用统计分析方法的SAS宏实现[期刊论文]-中国临床药理学与治疗学2006,11(1)

7. 童新元.张高魁.姚晨 定性指标的SAS统计分析报表--SAS软件在新药临床试验统计分析中的应用(一)[期刊论文]-中国卫生统计2003,20(1)

8. 薛富波.于莉莉.李婵娟.王素珍.徐丽.王陵.夏结来.XUE Fu-bo.YU Li-li.LI Chan-juan.WANG Su-zhen.XU Li.WANG Ling.XIA Jie-lai 新药临床试验中log-rank检验的SAS宏实现[期刊论文]-中国临床药理学与治疗学2006,11(2)

9. 田晓燕.金丕焕.张文彤 应用SAS软件的动态数据交换技术自动生成统计报告[期刊论文]-中华预防医学杂志2005,39(1)

10. 谢绍锋.陈炯华 精密度实验变异系数计算的SAS宏程序[期刊论文]-数理医药学杂志2008,21(1)

引证文献(1条)

1.许林勇.张伟.曾芳.付陈超.赵利.赵晓华.孙振球 动态血压监测评价降压药物疗效的谷峰比率的SAS统计自动报表[期刊论文]-中国临床药理学与治疗学 2010(7)

新药临床试验数据

本文链接：/Periodical_zglcylxyzlx200602023.aspx

本文来源：https://www.bwwdw.com/article/e0c1.html