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

实习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=WEIGHT0WEIGHT1，表示治疗后体重减轻量。(2) S0_1=S0S1，表示治疗后食欲减退情况。(3) EFFECT，表示治疗是否有效如果RANK=1,那么EFFECT=1，表示治疗无效；如果RANK=2或RANK=3,那么EFFECT=2,表示治疗有效。建立sas数据集：data a;set trial;w0_1=weight0-weight1;s0_1=s0-s1;if rank=1 then effect=1;if rank=2 or rank=3 then effect=2; run;2进行如下统计分析：（1）性别和年龄的两组比较。（2）两组治疗前后体重的变化情况和比较。（3）3个医院间试验组体重减轻情况的比较。（4）两组治疗前后食欲情况的变化和比较。（5）两组临床综合疗效的情况和比较。（6）两组有效率的比较。表1 两组一般情况的比较 试验组对照组统计量（卡方值）P值性别男38430.75160.3860女8475年龄（岁）37.02511.62536.8910.590.09 0.9254注：性别中填入例数（构成比），年龄中填入均数标准差性别构成比的比较为计数资料的2X2卡方检验proc freq data=a; table sex*group/chisq nopercent norow expected;run;年龄的均数比较用计量资料的成组设计的ttestproc ttest data=a;class group;var age;run;统计结论：性别：P0.05,按照a=0.05的水准，不拒绝H0，差异无统计学意义，可认为试验组和对照组男女构成比相同。年龄：P0.05,按照a=0.05的水准，不拒绝H0，差异无统计学意义，可认为试验组和对照组年龄均数相同。表2 两组治疗前后体重的变化（均数标准差）治疗前治疗后治疗前后差值试验组 82.8213.3378.6613.184.162.16对照组84.5312.5383.4412.571.091.56Means过程用于描述 均数 标准差 标准误 单样本t检验 或者配对t检验proc means data=a n mean std stderr t prt;class group;var weight0 weight1 w0_1;run;表3 两组治疗前后体重的比较治疗前比较（成组）两组分别治疗前后比较（配对）治疗前后差值比较（成组）统计量tP值统计量他P值统计量他P值试验组-1.02 0.309521.28.000112.68（方差不齐）.0001对照组7.550.05,按照a=0.05的水准，不拒绝H0，差异无统计学意义，可认为试验组和对照组治疗前体重均数相同。2. P0.05,按照a=0.05的水准，拒绝H0，差异有统计学意义，可认为试验组治疗前后体重不相同，对照组治疗前后体重不相同。3.P0.052中山414.062.720.180.8369中山和长海0.053长海394.331.84长海和瑞金0.05计量资料的多组设计用单因素方差分析glm程序data b;set a;if group=1 then output;run;proc sort data=b;by hospital;run;proc univariate data=b normal;var w0_1;by hospital;run;proc glm data=b;class hospital;model w0_1=hospital;means hospital/hovtest;means hospital/snk bon dunnett(1);3种方法means hospital/snk ALPHA=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等级资料 率或构成比的差异检验用四格表或者R x C列表的卡方检验 用freq 程序proc freq data=a; tables group*(s0 s1)/nopercent nocol; run; 表6 两组治疗前后食欲情况的比较治疗前比较成组两组分别治疗前后比较配对治疗前后差值比较成组统计量ZP值统计量SP值统计量ZP值试验组-0.1760.861983 .0001 -4.1346.0001对照组927.5 0.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等级资料 率或构成比的差异检验用四格表或者R x C列表的卡方检验 用freq 程序proc freq data=a; tables rank*group/nopercent norow; run; 等级资料两独立样本wilcoxon秩和检验proc npar1way data=a wilcoxon; class group; var rank; run;统计结论：P0.05,按照a=0.05的水准，拒绝H0，差异有统计学意义，可认为试验组和对照组临床综合疗效有差异。表8 两组有效率的比较总例数有效例数有效率（%）统计量P值试验组12211594.26