(精选)上海交通大学医学统计学上机实习14答案

实习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数据集：
data a;
set trial;
w0_1=weight0-weight1;
s0_1=s0-s1;
if rank=1then effect=1;
if rank=2 or rank=3then effect=2;
run;
2．进行如下统计分析：
（1）性别和年龄的两组比较。

（2）两组治疗前后体重的变化情况和比较。

（3）3个医院间试验组体重减轻情况的比较。

（4）两组治疗前后食欲情况的变化和比较。

（5）两组临床综合疗效的情况和比较。

（6）两组有效率的比较。

表1 两组一般情况的比较
试验组对照组统计量
（卡方
值）
P值
性别男38 43
0.7516 0.3860 女84 75
年龄（岁）37.025±11.625 36.89±10.59 0.09 0.9254 注：性别中填入例数（构成比），年龄中填入均数±标准差
性别构成比的比较为计数资料的2X2卡方检验
proc freq data=a;
table sex*group/chisq nopercent norow expected;
run;
年龄的均数比较用计量资料的成组设计的ttest
proc ttest data=a;
class group;
var age;
run;
统计结论：
性别：P>0.05,按照a=0.05的水准，不拒绝H0，差异无统计学意义，可认为试验组和对照组男女构成比相同。

年龄：P>0.05,按照a=0.05的水准，不拒绝H0，差异无统计学意义，可认为试验组和对照组年龄均数相同。

表2 两组治疗前后体重的变化（均数±标准差）
治疗前治疗后治疗前后差值
试验组 82.82±13.3378.66±13.18 4.16±2.16对照组84.53±12.5383.44±12.57 1.09±1.56
M eans过程用于描述均数标准差标准误单样本t检验或者配对t检验
proc means data=a n mean std stderr t prt;
class group;
var weight0 weight1 w0_1;
run;
表3 两组治疗前后体重的比较
治疗前比较（成组）两组分别治疗前后比较
（配对）
治疗前后差值比较
（成组）
统计量
P值统计量他P值统计量他P值
（方差不齐）<.0001
对照组7.55<.0001
统计结论：
1.P>0.05,按照a=0.05的水准，不拒绝H0，差异无统计学意义，可认为试验组和对照组治疗前体重均数相同。

2.P<0.05,按照a=0.05的水准，拒绝H0，差异有统计学意义，可认为试验组治疗前后体重不相同，对照组治疗前后体重不相同。

3.P<0.05,按照a=0.05的水准，拒绝H0，差异有统计学意义，可认为试验组和对照组治疗前后体重差值不相同。

表4 三个医院间试验组体重减轻情况的比较
治疗后体重的减轻三个医院间比较医院间两两比较
医院例数均数标准差统计量F P值比较的
医院
P值
1瑞金42 4.11 1.84 瑞金和
中山
>0.05
2中山41 4.06 2.72
0.18 0.8369 中山和
长海
>0.05
3长海39 4.33 1.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. 三个医院之间的比较P>0.05,按照a=0.05的水准，不拒绝H0，差异无统计学意义，尚
不能认为三个医院间试验组体重减轻情况有差异。

2. 由第一步已知三个医院之间治疗前后体重减轻情况相同，即不需进一步进行三组之间的
两两差异比较?。

表5 两组治疗前后食欲情况的变化
食欲情况
治疗前s0
治疗后s1 例数
构成比（%）
例数 构成比（%）
1试验组
1差
1 0.8
2 30 24.59 2一般
35 28.69 82 67.21 3强烈 63 51.64 10 8.20 4很强烈 23 18.85 0 0.00 2对照组
差
2 1.69 4 3.39 一般 34 28.81 78 66.10 强烈 60 50.85 30 25.42 很强烈
22
18.64
6
5.08
表6 两组治疗前后食欲情况的比较
治疗前比较
成组两组分别治疗前后比较
配对
治疗前后差值比较
成组
统计量
Z
P值统计量S P值统计量Z P值
试验组
-0.1760.861983 <.0001
-4.1346
<.0001
对照组927.5 <.0001
统计结论：
1.P>0.05,按照a=0.05的水准，不拒绝H0，差异无统计学意义，可认为试验组和对照组治疗前食欲相同。

2.两组均P<0.05,按照a=0.05的水准，拒绝H0，差异有统计学意义，可认为试验组和对照组分别治疗前后食欲不相同。

3.P<0.05,按照a=0.05的水准，拒绝H0，差异有统计学意义，可认为试验组和对照组治疗前后食欲差值不相同。

表7 两组临床综合疗效的比较
试验组对照组两组比较
综合疗效例数构成比（%）例数构成比（%）统计量P值无效7 5.74 73 61.86
-10.843
5<.0001
有效40 32.79 41 34.75
显效75 61.48 4 3.39
统计结论：
P<0.05,按照a=0.05的水准，拒绝H0，差异有统计学意义，可认为试验组和对照组临床综合疗效有差异。

表8 两组有效率的比较
总例数有效例数有效率（%）统计量P值
试验组122 115 94.26
85.0320<.0001
对照组118 45 38.14
proc freq data=a;
tables group*effect/chisq nopct nocol;
run;
统计结论：
P<0.05,按照a=0.05的水准，拒绝H0，差异有统计学意义，可认为试验组和对照组治疗有效率不同。

（注：专业文档是经验性极强的领域，无法思考和涵盖全面，素材和资料部分来自网络，供参考。

可复制、编制，期待你的好评与关注）。