.
多元logistic回归
1. 下面是子宫内膜癌的病例对照研究数据，暴露因素是雌激素。
分 组 病例组 对照组
使用过雌激素 55（a） 19（c）
未使用过雌激素 128（b） 164（d）
问题：使用过雌激素是否是子宫内膜癌的危险因素？危险强度为多少？
2. 为了探讨糖尿病与血压、血脂等因素的关系，研究者对56例糖尿病病人和65例对照者进行病例?对照研究，收集了性别、年龄、学历、体重指数、家族史、吸烟、血压、总胆固醇、甘油三脂、高密度脂蛋白、低密度脂蛋白11个因素的资料，各因素的观察结果见下表。 问题：糖尿病的相关因素有哪些？如何解释相关因素的作用大小？如何评价模型优劣？
因 素 变量名 赋 值 性别 X1 男=1，女=2 年龄 X2
小学以下=1，小学=2，初中=3，
学历 X3
高中=4，大专及以上=5
体重指数 X4 <24=1，24?<26=2， 26?=3 家族史 X5 无=1，有=2 吸烟 X6 不吸=1 吸=2 血压 X7 正常=1, 高=2 总胆固醇 X8 甘油三脂 X9 高密度脂蛋白 X10 低密度脂蛋白 X11 糖尿病 Y 对照=0，病例=1
体重总胆 甘油 高密度 低密度
编号 性别 年龄 学历 家族史 吸烟 血压 糖尿病
指数 固醇 三脂 脂蛋白 脂蛋白
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 .
1 1 2 1 1 1 1 1 2 2 1 1 2 2 2
60 48 63 68 45 45 59 68 63 58 44 46 62 65 58
2 3 2 3 2 3 2 3 2 2 2 3 1 1 2
2 2 1 2 1 3 1 3 2 2 2 1 2 2 1
1 1 1 2 2 2 1 1 1 1 2 1 1 1 1
1 1 1 1 1 1 1 1 1 1 1 2 1 1 1
1 1 2 1 1 1 1 1 1 1 2 1 2 1 2
4.30 4.60 4.60 4.15 3.42 4.16 4.32 3.80 3.87 5.42 4.35 3.42 3.18 3.30 4.41
1.50 1.32 1.15 1.43 1.22 0.96 1.02 1.42 1.55 0.87 1.01 1.26 1.38 0.85 1.05
1.24 1.15 1.15 1.07 0.63 0.98 1.05 2.86 2.44 4.46 5.13 1.40 1.67 1.92 2.97
2.30 2.30 2.30 3.21 2.30 2.65 3.49 0.85 0.81 3.14 2.20 0.28 0.48 0.69 1.79
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
16 1 17 1 18 2 19 2 20 1 21 2 22 2 23 2 24 1 25 1 26 1 27 2 28 1 29 2 30 2 31 1 32 1 33 1 34 2 35 2 36 1 37 1 38 1 39 2 40 2 41 1 42 1 43 2 44 2 45 2 46 2 47 2 48 1 49 2 50 1 51 1 52 2 53 1 54 1 55 2 56 2 57 2 58 2 59 1 60 1 61 2
.
22 5 69 1 55 2 66 1 47 2 47 3 64 2 63 1 64 1 70 2 48 1 64 2 47 2 54 4 57 1 45 2 62 2 52 1 50 2 50 2 58 3 54 3 60 1 60 2 63 3 65 2 37 3 52 2 46 1 63 1 63 2 60 3 58 2 65 3 43 2 58 1 61 4 46 4 58 2 54 4 45 1 38 2 22 5 44 1 49 3 49
2
1 2 1 1 1 1 2 1 1 2 2 2 2 2 2 1 1 1 2 1 2 1 1 1 2 2 2 1 2 1 2 1 2 1 2 2 2 1 3 2 1 1 2 1 3 2 1 1 1 1 2 2 2 1 2 1 3 1 2 2 2 2 2 2 2 1 1 1 2 2 2 2 2 1 2 1 2 2 2 1 2 2 2 1 3 2 3 1 3 1 3
1
.
2 2 5.10 2 1 5.09 1 2 3.98 1 1 4.27 1 1 4.19 1 1 4.19 1 1 4.47 1 1 4.23 2 1 4.16 2 1 4.77 1 1 3.78 1 1 5.20 1 2 4.25 1 2 3.56 1 1 4.80 1 1 4.65 1 2 4.23 1 2 3.70 1 1 4.11 2 1 3.85 1 1 3.01 2 1 3.02 1 1 4.16 1 2 5.12 1 1 4.30 1 1 3.67 1 1 5.46 1 2 4.35 1 1 3.60 1 2 3.16 1 2 3.82 1 1 4.63 1 2 4.05 1 2 3.80 1 1 4.98 2 1 5.04 1 1 3.53 1 2 3.50 1 1 4.83 1 1 2.88 1 1 3.31 1 2 4.13 2 1 3.66 2 2 4.15 1 2 4.21 1
2
3.97
1.63 3.24 1.50 3.29 1.41 0.66 1.35 2.30 1.24 2.49 1.04 2.43 1.60 3.01 1.50 4.04 1.44 1.63 1.08 1.79 1.00 1.19 1.56 1.21 1.06 3.34 1.35 1.24 1.26 3.19 0.92 2.99 1.15 3.28 0.85 2.55 1.23 2.54 1.35 2.34 1.08 1.86 1.23 1.55 1.02 2.46 1.31 3.50 1.40 2.80 0.72 2.29 1.45 3.56 1.13 3.86 0.87 2.27 0.68 2.30 1.36 1.98 1.28 2.46 0.85 2.60 1.42 1.85 1.48 5.07 1.13 3.77 0.74 1.66 1.53 1.94 0.83 2.82 1.04 1.58 1.69 1.63 1.38 3.92 0.73 0.71 1.31 1.90 1.59 2.49 1.34
4.10
0.96 0 0.75 0 4.25 0 1.05 0 1.65 0 1.02 0 0.95 0 1.14 0 1.95 0 2.69 0 2.63 0 1.28 0 0.94 0 0.88 0 0.93 0 0.79 0 1.07 0 2.97 0 1.67 0 1.09 0 1.14 0 0.40 0 0.80 0 1.11 0 0.97 0 0.66 0 0.94 0 1.11 0 1.00 0 1.23 0 1.18 0 0.33 0 0.77 0 0.79 0 1.03 0 0.77 0 0.55 0 0.44 0 0.96 0 0.55 0 0.72 0 1.26 0 0.53 0 0.73 0 1.02 0 1.10
0
62 2 63 1 64 1 65 1 66 2 67 2 68 1 69 1 70 2 71 2 72 1 73 2 74 1 75 1 76 2 77 1 78 2 79 1 80 1 81 2 82 2 83 2 84 2 85 2 86 2 87 2 88 1 89 1 90 2 91 1 92 2 93 2 94 2 95 2 96 2 97 2 98 1 99 1 100 2 101 2 102 2 103 1 104 1 105 2 106 2 107 1
.
41 3 44 3 48 1 48 2 51 2 60 1 67 1 48 1 58 2 41 1 52 1 54 1 60 2 67 2 43 1 61 1 52 1 60 1 50 1 45 1 45 3 60 1 34 1 54 1 58 2 61 1 53 1 46 1 72 1 58 2 50 1 60 2 55 1 34 1 64 2 62 1 65 3 43 2 52 1 61 1 53 1 64 1 38 1 60 1 44 1 62
1
2 1 3 2 1 1 2 1 3 1 3 1 2 1 3 1 2 2 2 2 3 1 2 1 2 2 3 1 3 2 2 1 2 2 2 2 2 1 2 1 3 1 3 2 3 1 2 2 3 1 2 2 3 2 2 2 1 1 2 1 3 2 2 2 3 2 3 1 2 2 2 1 3 2 1 2 3 1 2 2 3 2 3 2 2 2 3 1 2 2 2
1
.
1 1 3.13 1 1 4.46 2 1 4.06 1 1 4.33 1 1 2.75 1 2 4.50 2 1 2.70 2 1 3.54 2 1 3.09 2 1 5.12 2 1 4.81 1 1 4.93 2 2 4.38 1 2 3.61 2 2 4.23 1 2 4.97 1 2 4.42 2 2 4.80 2 2 4.44 2 2 4.83 2 1 5.68 1 1 5.85 2 2 4.69 2 1 4.30 2 2 4.67 2 2 4.91 2 2 3.44 1 2 3.95 2 2 3.86 2 2 3.68 2 2 5.17 1 2 4.07 2 2 3.38 2 2 4.75 2 2 3.89 2 1 5.73 2 1 3.68 2 2 3.57 2 2 5.50 2 1 3.82 2 2 3.08 1 2 3.87 1 1 9.47 1 1 3.06 2 2 4.17 2
1
6.79
0.75 1.62 1.33 2.86 1.54 2.35 1.60 4.81 1.10 0.84 1.22 0.74 1.00 0.96 0.63 1.21 0.82 2.05 0.71 1.38 1.66 0.75 1.43 0.61 1.56 1.93 1.14 1.17 1.34 1.26 1.32 1.05 1.31 1.24 1.65 0.85 1.49 1.20 1.66 1.04 1.43 1.44 2.57 1.39 2.04 1.12 3.10 0.75 1.76 2.29 2.56 2.37 0.74 2.77 1.65 1.18 1.32 1.19 0.97 1.31 4.31 2.54 2.45 2.42 1.15 2.68 2.45 2.51 1.29 1.93 1.72 4.41 0.97 1.97 1.17 2.82 1.79 3.95 6.47 2.67 1.39 2.28 1.19 1.72 2.38 2.26 1.03 2.07 0.94 2.21 2.65
3.82
0.57 0 0.27 0 0.80 0 1.17 0 1.59 1 3.27 1 0.80 1 1.90 1 2.50 1 2.61 1 3.29 1 3.12 1 3.71 1 2.60 1 3.06 1 0.45 1 2.30 1 2.92 1 2.91 1 0.96 1 1.44 1 3.03 1 2.86 1 2.05 1 0.84 1 0.86 1 0.78 1 2.04 1 2.30 1 1.97 1 1.63 1 1.19 1 0.49 1 0.32 1 0.61 1 0.66 1 0.61 1 1.08 1 1.40 1 1.11 1 1.12 1 0.35 1 0.86 1 0.67 1 0.83 1 0.68
1
.
108 109 110 111 112 113 114 115 116 117 118 119 120 121
2 2 1 2 1 1 1 1 1 1 1 1 1 1
60 55 60 65 63 67 42 60 57 69 62 55 67 65
1 4 2 1 1 1 1 1 1 3 1 4 5 1
2 2 2 2 2 1 2 3 2 2 2 2 2 2
2 2 2 1 2 2 2 1 2 1 2 2 2 2
2 2 2 1 1 1 2 2 2 1 2 2 2 2
2 2 1 1 1 2 1 1 2 1 1 2 1 1
3.03 4.51 4.98 5.01 4.28 4.02 5.12 3.58 5.48 4.81 5.72 4.71 5.41 4.10
0.89 1.24 1.97 2.59 1.69 0.93 1.58 1.51 1.77 1.73 1.24 1.20 1.30 1.10
2.14 2.98 3.31 3.28 3.01 2.76 4.05 2.12 3.76 2.95 3.50 2.92 2.99 1.72
0.90 0.87 0.99 1.02 0.60 0.74 1.36 1.43 1.02 0.86 0.96 0.88 1.08 0.73
1 1 1 1 1 1 1 1 1 1 1 1 1 1
3. 研究人员对使用雌激素与子宫内膜癌发病间的关系进行了1 ：1配对的病例-对照研究。病例与对照按年龄相近、婚姻状况相同、生活的社区相同进行了配对。收集了年龄、雌激素药使用、胆囊病史、高血压和非雌激素药的使用的数据。变量定义及具体数据如下。对使用雌激素与子宫内膜癌发病间的关系作1ogistic回归分析。 1: 2: 3: 4: 5: 6 7: match case age est gall hyper nonest match case age est gall hyper nonest 1 1 2 2 3 3 4 4 5 5 .
Match Case(即y) age est gall hyper nonest
配比组（匹配条件：年龄相近、婚姻状况相同、同一生活社区） y=1：病例， y=0：对照 _____岁
est=0：未使用过雌激素， est=1：使用过雌激素 gall=0：无胆囊病史， gall=1：有胆囊病史 hyper=0: 不是高血压， hyper=1: 是高血压
nonest=0：未使用过非雌激素， nonest=1：使用过非雌激素
1 0 1 0 1 0 1 0 1 0 62 62 79 79 71 71 72 71 83 82 1 0 1 0 1 1 1 1 1 0 1 1 0 0 1 0 0 0 0 0 0 0 1 0 0 1 0 1 1 0 1 0 1 1 1 1 0 1 1 0 31 31 32 32 33 33 34 34 35 35 1 0 1 0 1 0 1 0 1 0 65 65 80 79 68 68 64 64 61 61 1 0 1 0 1 0 1 1 0 0 1 0 0 0 1 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 1 0 1 0 1 1 1 1 1 1 .
6 6 7 7 8 8 9 9 10 10 11 11 12 12 13 13 14 14 15 15 16 16 17 17 18 18 19 19 20 20 21 21 22 22 23 23 24 24 25 25 26 26 27 .
1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 0 72 72 70 70 65 67 74 75 67 66 70 70 64 63 66 66 62 63 77 77 83 83 66 67 81 79 74 75 67 66 67 67 73 72 76 76 67 67 71 70 74 75 69 1 1 0 0 1 0 1 0 1 1 1 0 1 0 1 1 1 0 1 1 0 0 1 0 1 0 1 0 1 1 1 0 1 1 1 1 1 1 1 1 1 0 1 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 1 1 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 1 1 0 1 0 0 0 1 1 0 1 1 0 0 0 1 0 0 0 0 0 1 0 0 0 0 1 1 0 1 1 0 0 0 0 0 1 0 1 1 0 1 1 0 0 0 0 1 0 1 1 1 1 1 1 0 1 0 1 1 1 0 1 1 1 1 0 0 1 1 1 0 1 1 1 0 1 0 1 1 1 1 1 1 1 1 1 1 0 1 1 1 1 36 36 37 37 38 38 39 39 40 40 41 41 42 42 43 43 44 44 45 45 46 46 47 47 48 48 49 49 50 50 51 51 52 52 53 54 54 55 55 56 56 57 57 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 0 1 0 1 0 1 0 1 0 63 63 64 65 72 72 68 69 57 57 74 74 74 74 67 68 62 62 69 69 73 72 77 76 71 71 64 64 64 65 79 82 63 64 72 79 78 82 81 80 81 73 74 1 1 1 0 0 0 1 1 1 1 1 1 0 0 1 1 1 1 1 0 1 1 1 1 1 1 1 1 0 1 0 1 1 0 1 1 1 1 0 1 0 1 1 0 0 0 0 1 0 1 1 0 0 0 0 1 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 1 1 1 0 0 0 0 0 0 1 0 0 0 1 0 0 0 1 0 0 0 0 0 1 1 0 1 0 1 0 0 1 1 1 0 0 1 0 0 0 0 0 0 1 1 1 0 0 1 1 0 1 1 1 0 1 0 1 0 0 1 1 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 1 1 0 1 1 1 1 1 1 1 0 1 1 1 1
多元logistic回归实习题
