(精选)上海交通大学医学统计学上机实习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