医学统计学(李晓松主编 第2版 高等教育出版社)附录 思考与练习(3)

2019-03-29 14:25

第四章 定量资料的统计描述

图3-2 Data View窗口内录入数据

分析:

Transform→Random Number Generators … Active Generator Initialization

Set Starting Point Fixed Value

Value: 键入20071222 OK Transform→Compute… Target Valuable: 键入 random Numeric Expression: RV.Uniform(0,1) OK Transform→Rank Cases… Variables(s): random By: pair OK

设定随机种子为20071222

产生范围在0~1之间的伪随机数

对每个对子中的两个伪随机数进行排序

注:当样本量较大时,最好编程实现,以免除数据录入的不便。以下为SPSS的Syntax窗口中用编程来完成本题。

File→New→Syntax,打开Syntax窗口,键入如下程序:

input program.

numeric k b n number pair (F8.0)/random(F8.6). string treat(A1). compute k=2.

开始数据录入程序段

处理组为2

4-11

第四章 定量资料的统计描述

compute b=6. compute n=b*k. loop number=1 to n.

compute pair=rnd((number-1)/k+0.5). end case.

compute k=lag(k). end loop. end file.

对子数为6

结束数据录入程序段

设定随机种子为20071222 产生范围在0~1之间的伪随机数

对伪随机数排序

根据伪随机数的序号进行分组

end input program. set seed 20071222.

compute random=uniform(1). rank variables=random by pair. formats rrandom(F3.0). compute

treat=substr(\\list number pair treat.

结果中显示受试对象编号及被分到的处理组

在Syntax窗口中选择Run→All 提交运行。

注:程序中的b(即对子数)可根据实际情况赋予不同的值。随机区组设计受试对象的分配也可以用此程序来实现,仅需将k值设为处理组数,b值设为相应的区组数。

结果及解释 Data View窗口

4-12

第四章 定量资料的统计描述

图3-3 SPSS编程实现受试对象随机分配结果

Output窗口

List

number 1 2 3 4 5 6 7 8 9 10 11 12

pair 1 1 2 2 3 3 4 4 5 5 6 6

treat A B A B B A A B B A B A

编程实现随机分配的结果见data view窗口(图3-3)或output 窗口。可以看出12只小鼠被随机分配到A组和B组。在用此程序进行随机分组前,规定A组表示实验组,B组表示对照组。则1,3,6,7,10,12号小鼠被分到实验组,2,4,5,8,9,11号小鼠被分到对照组。

3. 解:

4-13

第四章 定量资料的统计描述

(1) 该研究最好采用2×2析因设计的实验方案,如表3-2所示。因为析因设计不仅可以检验两因素各水平之间的差异有无统计学意义,而且可以检验两因素间的交互作用。因此既可分析MEL和TXL对糖尿病模型大鼠血管紧张素的影响,又可分析两药合用时的交互效应。

表3-2 MEL和TXL对32只大鼠血管紧张素影响的析因设计

MEL (用)a1 (不用)a2

TXL

(用)b1 a1b1 a2b1

(不用)b2 a1b2 a2b2

(2) 先将32只大鼠按体重从小到大的顺序编号。再从随机数字表中任一行,如第2列最上端开始纵向连续取32个两位数字。对随机数字排序,事先规定,序号为1~8的豚鼠分为A组,序号为9~16的大鼠分为B组,序号为17~24的大鼠分为C组,序号为25~32的大鼠分为D组,分配结果见表3-3。

表3-3 32只大鼠随机分组的结果

编号

1

2

3 77 25 D 19 94 29 D

4 43 11 B 20 64 17 C

5 28 5 A 21 44 13 B

6 22 4 A 22 76 23 C

7 76 22 C 23 17 2 A

8 68 19 C 24 17 3 A

SPSS操作 数据录入:

打开SPSS Data Editor窗口,点击Variable View标签,定义要输入的变量number,见图3-4。再点击Data View标签,录入变量number的值,即从1~32。

9 39 9 B 25 76 24 C

10 71 20 C 26 29 6 A

11 35 7 A 27 80 26 D

12 50 14 B 28 40 10 B

13 96 30 D 29 56 16 B

14 15 16 93 87 56 28 27 15 D

D

B

随机数 17 36 序号 分组 编号

1 A

8 A

17 18 30 31 32 65 43 96 18 12 32 C

B

D

随机数 72 96 序号 分组

21 31 C

D

图3-4 Variable View窗口内定义要输入的变量number

4-14

第四章 定量资料的统计描述

分析:

Transform→Random Number Generators… Active Generator Initialization

Set Starting Point Fixed Value

Value: 键入20071222 OK Transform→Compute… Target Valuable: 键入 random Numeric Expression: RV.Uniform(0,1) OK Transform→Rank Cases… Variables(s): random OK Transform→Recode→Into Different Variables… Input VariableOutput Variable: Rrandom Output Variable Name: 键入treat 点击Change 点击Old And New Values… Old Value

Range: 键入1 through 8

Output variables are strings New Value

Value: 键入A OldNew: 点击Add Old Value

Range: 键入9 through 16 New Value

Value: 键入B OldNew: 点击Add Old Value

Range: 键入17 through 24 New Value

Value: 键入C

4-15

设定随机种子为20071212

产生范围在0~1之间的伪随机数

对伪随机数排序

原变量名为伪随机数的序号Rrandom

新变量名为treat

Rrandom值为1~8时,treat值为A

Rrandom值为9~16时,treat值为B

Rrandom值为17~24时,treat值为C


医学统计学(李晓松主编 第2版 高等教育出版社)附录 思考与练习(3).doc 将本文的Word文档下载到电脑 下载失败或者文档不完整,请联系客服人员解决!

下一篇:赤峰学院2011年迎新工作方案

相关阅读
本类排行
× 注册会员免费下载(下载后可以自由复制和排版)

马上注册会员

注:下载文档有可能“只有目录或者内容不全”等情况,请下载之前注意辨别,如果您已付费且无法下载或内容有问题,请联系我们协助你处理。
微信: QQ: