上海交通大学医学统计学上机实习14答案

-1-实习14（上机）在一个减肥药治疗单纯性肥胖的多中心临床试验中，入组病例随机分为两组共240例合格病例；试验组122例给以减肥药治疗，对照组118例给以安慰剂。源数据在SASEXER文件夹的Excel数据文件trial.xls中，各指标及其意义如下：HOSPITAL医院1=瑞金，2=中山，3=长海GROUP组别1=试验组，2=对照组SEX性别1=男性，2=女性AGE年龄（岁）WEIGHT0治疗前体重（kg）WEIGHT1治疗后体重（kg）S0治疗前食欲1=差，2=一般，3=强烈，4=很强烈S1治疗后食欲1=差，2=一般，3=强烈，4=很强烈RANK临床综合疗效1=无效，2=有效，3=显效1．读入数据建立SAS数据集，同时产生3个新变量(1)W0_1=WEIGHT0－WEIGHT1，表示治疗后体重减轻量。(2)S0_1=S0－S1，表示治疗后食欲减退情况。(3)EFFECT，表示治疗是否有效如果RANK=1,那么EFFECT=1，表示治疗无效；如果RANK=2或RANK=3,那么EFFECT=2,表示治疗有效。建立sas数据集：dataa;settrial;w0_1=weight0-weight1;s0_1=s0-s1;ifrank=1theneffect=1;ifrank=2orrank=3theneffect=2;run;2．进行如下统计分析：（1）性别和年龄的两组比较。（2）两组治疗前后体重的变化情况和比较。（3）3个医院间试验组体重减轻情况的比较。（4）两组治疗前后食欲情况的变化和比较。（5）两组临床综合疗效的情况和比较。（6）两组有效率的比较。-2-表1两组一般情况的比较试验组对照组统计量（卡方值）P值性别男38430.75160.3860女8475年龄（岁）37.025±11.62536.89±10.590.090.9254注：性别中填入例数（构成比），年龄中填入均数±标准差性别构成比的比较为计数资料的2X2卡方检验procfreqdata=a;tablesex*group/chisqnopercentnorowexpected;run;年龄的均数比较用计量资料的成组设计的ttestprocttestdata=a;classgroup;varage;run;统计结论：性别：P0.05,按照a=0.05的水准，不拒绝H0，差异无统计学意义，可认为试验组和对照组男女构成比相同。年龄：P0.05,按照a=0.05的水准，不拒绝H0，差异无统计学意义，可认为试验组和对照组年龄均数相同。表2两组治疗前后体重的变化（均数±标准差）治疗前治疗后治疗前后差值试验组82.82±13.3378.66±13.184.16±2.16对照组84.53±12.5383.44±12.571.09±1.56Means过程用于描述均数标准差标准误单样本t检验或者配对t检验procmeansdata=anmeanstdstderrtprt;classgroup;varweight0weight1w0_1;run;-3-表3两组治疗前后体重的比较治疗前比较（成组）两组分别治疗前后比较（配对）治疗前后差值比较（成组）统计量tP值统计量他P值统计量他P值试验组-1.020.309521.28.000112.68（方差不齐）.0001对照组7.55.0001计量资料的t检验公式成组（ttest）不同于配对(means)★注意成组设计的方差齐性检验procttestdata=a;classgroup;varweight0;run;procmeansdata=ameanstdstderrtprt;classgroup;varw0_1;run;procttestdata=a;classgroup;varw0_1;run;统计结论：1.P0.05,按照a=0.05的水准，不拒绝H0，差异无统计学意义，可认为试验组和对照组治疗前体重均数相同。2.P0.05,按照a=0.05的水准，拒绝H0，差异有统计学意义，可认为试验组治疗前后体重不相同，对照组治疗前后体重不相同。3.P0.05,按照a=0.05的水准，拒绝H0，差异有统计学意义，可认为试验组和对照组治疗前后体重差值不相同。表4三个医院间试验组体重减轻情况的比较治疗后体重的减轻三个医院间比较医院间两两比较医院例数均数标准差统计量FP值比较的医院P值1瑞金424.111.84瑞金和中山0.052中山414.062.720.180.8369中山和长海0.053长海394.331.84长海和瑞金0.05-4-计量资料的多组设计用单因素方差分析glm程序datab;seta;ifgroup=1thenoutput;run;procsortdata=b;byhospital;run;procunivariatedata=bnormal;varw0_1;byhospital;run;procglmdata=b;classhospital;modelw0_1=hospital;meanshospital/hovtest;meanshospital/snkbondunnett('1');3种方法meanshospital/snkALPHA=0.05;run;统计结论：1.三个医院之间的比较P0.05,按照a=0.05的水准，不拒绝H0，差异无统计学意义，尚不能认为三个医院间试验组体重减轻情况有差异。2.由第一步已知三个医院之间治疗前后体重减轻情况相同，即不需进一步进行三组之间的两两差异比较?。表5两组治疗前后食欲情况的变化食欲情况治疗前s0治疗后s1例数构成比（%）例数构成比（%）1试验组1差10.823024.592一般3528.698267.213强烈6351.64108.204很强烈2318.8500.002对照组差21.6943.39一般3428.817866.10强烈6050.853025.42很强烈2218.6465.08等级资料率或构成比的差异检验用四格表或者RxC列表的卡方检验用freq程序procfreqdata=a;tablesgroup*(s0s1)/nopercentnocol;run;-5-表6两组治疗前后食欲情况的比较治疗前比较成组两组分别治疗前后比较配对治疗前后差值比较成组统计量ZP值统计量SP值统计量ZP值试验组-0.1760.861983.0001-4.1346.0001对照组927.5.0001等级资料两独立样本wilcoxon秩和检验配对资料的符号秩和检验procnpar1waydata=awilcoxon;（成组秩和）classgroup;vars0;run;procunivariatedata=anormal;（配对秩和）classgroup;vars0_1;run;procnpar1waydata=awilcoxon;（成组秩和）classgroup;vars0_1;run;统计结论：1.P0.05,按照a=0.05的水准，不拒绝H0，差异无统计学意义，可认为试验组和对照组治疗前食欲相同。2.两组均P0.05,按照a=0.05的水准，拒绝H0，差异有统计学意义，可认为试验组和对照组分别治疗前后食欲不相同。3.P0.05,按照a=0.05的水准，拒绝H0，差异有统计学意义，可认为试验组和对照组治疗前后食欲差值不相同。-1-表7两组临床综合疗效的比较试验组对照组两组比较综合疗效例数构成比（%）例数构成比（%）统计量P值无效75.747361.86-10.8435.0001有效4032.794134.75显效7561.4843.39等级资料率或构成比的差异检验用四格表或者RxC列表的卡方检验用freq程序procfreqdata=a;tablesrank*group/nopercentnorow;run;等级资料两独立样本wilcoxon秩和检验procnpar1waydata=awilcoxon;classgroup;varrank;run;统计结论：P0.05,按照a=0.05的水准，拒绝H0，差异有统计学意义，可认为试验组和对照组临床综合疗效有差异。表8两组有效率的比较总例数有效例数有效率（%）统计量P值试验组12211594.2685.0320.0001对照组1184538.14两样本率的比较用四格表卡方检验procfreqdata=a;tablesgroup*effect/chisqnopctnocol;run;统计结论：P0.05,按照a=0.05的水准，拒绝H0，差异有统计学意义，可认为试验组和对照组治疗有效率不同。