第四章 定量资料的统计描述
图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