临床试验数据的SAS模块化分析课件

临床试验资料的SAS模块化分析我所理解的临床试验数据编程原则尽可能利用SAS固有的过程，尽量少编写SAS没有的算法；充分利用DATA步来实现数据的输出格式控制，以减少填表过程中出现的错误；编写适用于报告中的各种表格的宏程序库。2尽可能利用SAS固有的过程要用SAS，就必须采用SAS的固有输出结果，你所要做的仅仅是把你要的结果放在该放的地方。这样才能保证SAS固有的算法上的可靠性，才能保证统计分析结果的可靠性、可再现性和试验间的可比性。夏结来3技巧：频数表的输出结果控制直接用Freq命令的缺点：用Freq过程直接输出难以控制输出格式多余的”-””+”等等存在着某些类别没有频数的情形，导致Freq输出信息不完整。5datatest;inputfab@;cards;20111512022022;procfreqdata=test;title"带有频数为0的四格表";tablesa*b/nocolnopercent;weightf;run;

带有频数为0的四格表FREQ过程a*b表ab频数|行百分比|1|2|合计--------+--------+--------+A|20|15|35|57.14|42.86|--------+--------+--------+合计2015356解决方法：生成一个包含了全部可能取值的数据表；再将Freq的输出数据表和该数据表合并；利用Do循环分别从合并后的数据表中取出每一行的频数及频率，最后组合在一起。7ODSOUTPUTCrossTabFreqs=temp2;proc

freqdata=test;tablesa*b/NOCOLNOROWNOPERCENTNOCUM;weightf;run;ODSOUTPUTCLOSE;RUN;PROC

SORTDATA=temp2;BYab;RUN;PROC

SORTDATA=APPEND;BYab;RUN;DATAtemp3;MERGEtemp2(keep=abfrequency)APPEND;BYab;IFFREQUENCY=.THENFREQUENCY=0;ifa~=.andb~=.;RUN;带有频数为0的四格表abFrequency112012152102209%MACRODIST();DATAALLROW;SET_NULL_;RUN;%DOI=1%TO2;DATAROW&I;SET_NULL_;RUN;DATAROW&I;SETTEMP3;IFA=&I;KEEPFREQUENCY;RUN;PROCTRANSPOSEDATA=ROW&IOUT=ROW&I;RUN;DATAROW&I;SETROW&I;DROP_NAME__LABEL_;A=&I;OUTPUT;RUN;DATAALLROW;SETALLROWROW&I;RUN;%END;DATAALLROW;SETALLROW;COL3=COL1+COL2;COL4=PUT(100*COL1/COL3,8.2);COL5=PUT(100*COL2/COL3,8.2);IFCOL4=""THENCOL4="-";IFCOL5=""THENCOL5="-";LABELCOL1="阳性"COL2="阴性"COL3="合计"COL4="阳性率"A="分组";RUN;%MEND;带有频数为0的四格表

分组阳性阴性合计阳性率120153557.142000.10编写适用于报告中的各种表格的宏程序库从临床试验规范化的角度看，程序员所编写的程序应当由专人复核，这就要求程序员所撰写的程序必须尽可能清晰易懂；最好是统计分析报告中的每一个输出表格都对应一个宏。在宏程序中增加选项的设置，强化某一个宏程序的功能，减少宏的个数。这样程序审核者只要检查宏调用的参数是否正确便可。11%LABAUTO(TABNUM,TABTITLE,DATANAME,V1,V2,PJ1,PJ2,VLABEL,HOP,L,U,COND);用于输出试验室指标疗前、疗后正异常改变情况。%labauto(6.2.1,试验室指标改变情况,SA,HB0,HB1,YHB0,YHB1,血红蛋白,3,110

110

115,160

160

160,NONE);表6.2.1:试验室指标改变情况(SA)指标分组正/正正/异异/正异/异合计血红蛋白A6514171血红蛋白B6008270红细胞A6911071红细胞B6603170白细胞A6118171白细胞B5929070血小板A6361070血小板B6532070ALTA6144271ALTB6405170ASTA6532171ASTB6601168BUNA6731071BUNB6522170CrA378131371CrB4078146913%FFTAB(TABNUM,TABTITLE,DATANAME,VARX,VLABEL,CODE,POS,COND,TYPE);输出四格表及检验结果%FFTAB(4.2,疗前厌氧菌检查的比较,FAS,XJ,细菌名,12,2,,1);%FFTAB(4.2,疗前厌氧菌检查的比较,FAS,XJ,细菌名,12,2,,2);表4.2:疗前厌氧菌检查的比较(FAS)阴性/阳性/ID正常异常合计发生率发生率比较细菌名43297240.280.3858细菌名48237132.39表4.2:疗前厌氧菌检查的比较(FAS)ID阴性/正常阳性/异常合计发生率比较细菌名43(59.72%)29(40.28%)720.3858细菌名48(67.61%)23(32.39%)7114%QUANTEST(TABNUM,TABTITLE,DATANAME,V,VLABEL,K,METHOD,COND);用于定量资料的组间比较(包括t检验和秩和检验)%OUTDEF();%QUANTEST(5.5.1,疗后总分比较,FAS,ZF0ZF1ZF2,疗前总分用药后3天结束后1－2天,3,T,);%PRINTOUT;表5.5.1:疗后总分比较IDgroupNMean±SDM±Q最小值最大值组间比较疗前总分A729.90±2.9510.00±4.004.0017.00t=0.51，P=0.6103疗前总分B719.66±2.6810.00±4.003.0014.00用药后3天A723.53±2.363.00±3.000.0011.00t=-0.62，P=0.5346用药后3天B713.76±2.114.00±3.000.009.00结束后1－2天A721.01±1.970.00±1.000.008.00t=1.12，P=0.2635结束后1－2天B710.72±1.030.00±1.000.004.0015％SYMDIS(TABNUM,TABTITLE,DATANAME,V1,ITEM,NEG,TIMEV,COND);用于症状消失率的比较，指定时间变量后，也可进行Log-Rank检验。%SYMDIS(5.3.1,第三天症状、体征消失率比较,FAS,B1,发热,0,,A1~=0);%SYMDIS(5.8,痊愈天数比较,FAS,LX3,临床疗效,1,DAY,);表5.3.1:第三天症状、体征消失率比较(FAS)ID疗前阳性消失消失率消失率比较发热382257.890.6237.322165.63表5.8:痊愈天数比较(FAS)

LogRankID疗前阳性消失消失率消失率比较Med±QuarX2P临床疗效726083.331.00005.00±0.002.480.1155.705882.865.00±0.0017%QUAN_BET(TABNUM,TABTITLE,DATANAME,V,VLABEL,K,CI,CT,SEQ,COV,FACT,COND);用于定量资料组内及组间比较，可以控制协变量和混杂因素，可以进行中心间一致性检验，估计可信区间等。%OUTDEF();%QUAN_BET(5.6.2,总分下降情况(PP),PP,ZF0ZF1ZF2,疗前用药后3天结束后1－2天,3,1,1,1,ZF0,SEX,);%PRINTOUT;表5.6.2:总分下降情况(PP)ID分组N组内差值组内比较疗前-用药后3天A726.38±2.78(5.72，7.03)t=19.49，P<0.0001疗前-用药后3天B705.94±2.53(5.34，6.55)t=19.65，P<0.0001疗前-结束后1－2天A728.89±2.99(8.19，9.59)t=25.26，P<0.0001疗前-结束后1－2天B709.03±2.62(8.40，9.65)t=28.82，P<0.0001组间差值组间比较疗前-用药后3天A0.32±0.29(-0.26，0.91)F=1.21，P=0.2723疗前-用药后3天B中心间一致性检验F=1.23，P=0.2965疗前-结束后1－2天A-0.28±0.24(-0.76，0.19)F=1.37，P=0.2431疗前-结束后1－2天B中心间一致性检验F=0.86，P=0.427518%OC_BET(TABNUM,TABTITLE,DATANAME,V,VLABEL,K,CODE,COND);用于等级资料的列表及组间比较%OUTDEF();%OC_BET(5.2.1,停药后1－2天症状、体征评分,FAS,C1C2C3C4C5C6,发热牙痛叩痛开口度肿胀溢脓,6,0123,);%PRINTOUT();表5.2.1:停药后1－2天症状、体征评分指标分组0123合计MSD组间比较发热A601200