logistic 临床预测模型分析

R语言临床logistics模型预测

数据:
CHD Age FCHD
1 0 23 No
2 0 24 No
3 0 25 No
4 0 26 No
5 0 26 No
6 0 28 No
7 0 28 Yes
8 0 29 No
9 0 30 No
10 0 30 No
11 0 30 Yes
12 0 30 No
13 0 30 No
14 0 32 No
15 0 32 No
16 0 33 No
17 0 33 Yes
18 0 34 No
19 0 34 No
20 0 34 No
21 0 34 No
22 0 35 No
23 0 35 Yes
24 0 36 No
25 0 36 No
26 0 37 No
27 0 37 Yes
28 0 38 No
29 0 38 Yes
30 0 39 No
31 0 40 No
32 0 41 No
33 0 41 Yes
34 0 42 No
35 0 42 Yes
36 0 42 No
37 0 43 No
38 0 43 No
39 0 44 No
40 0 44 Yes
41 0 45 No
42 0 46 No
43 0 47 No
44 0 47 Yes
45 0 48 No
46 0 49 Yes
47 0 49 No
48 0 50 No
49 0 51 No
50 0 52 Yes
51 0 55 Yes
52 0 57 No
53 0 57 Yes
54 0 58 No
55 0 60 Yes
56 0 64 No
57 1 25 Yes
58 1 30 Yes
59 1 34 No
60 1 36 No
61 1 37 Yes
62 1 39 No
63 1 40 No
64 1 42 Yes
65 1 43 No
66 1 44 No
67 1 44 Yes
68 1 45 No
69 1 46 No
70 1 47 Yes
71 1 48 Yes
72 1 48 No
73 1 49 No
74 1 50 Yes
75 1 52 Yes
76 1 53 No
77 1 53 Yes
78 1 54 No
79 1 55 Yes
80 1 55 No
81 1 56 Yes
82 1 56 No
83 1 56 Yes
84 1 57 No
85 1 57 No
86 1 57 Yes
87 1 57 Yes
88 1 58 No
89 1 58 Yes
90 1 59 Yes
91 1 59 No
92 1 60 No
93 1 61 Yes
94 1 62 Yes
95 1 62 No
96 1 63 Yes
97 1 64 No
98 1 65 No
99 1 69 Yes

问题1)估计如下 CHD 对 FCHD 和 Age 的逻辑回归模型(只需估计参数 β_0,〖 β〗_1 和 β_2):

log⁡〖(p/(1-p))〗=β_0+β_1 FCHDYes+β_2 Age ,
其中 p 为CHD=1的概率,FCHDYes为FCHD=“Yes”的示性(哑)变量。

问题2) :利用上面的模型,估计某个无家族冠心病史、年龄35岁的个体得冠心病的风险。
我不会构建Logistics模型,同时不会对构建的模型进行预测。
我的代码:
C$CHD<- as.factor(C$CHD)
class(C$CHD)
C$FCHD<- as.factor(C$FCHD)
class(C$FCHD)
CHD1<-glm(CHD~., data = C,family = binomial())
summary(CHD1)



指导如何构建模型病指导预测。

谢谢