通过使用R,人们如何能够开发出指数评分来预测患者的总体生存(OS)?
我有4个候选预测变量的候选列表,这些候选预测变量显示与OS相关。它们来自Cox多元回归(使用coxph()
运行)。预测因子是蛋白质水平,因此它们都是连续变量。
数据表如下所示(此处仅显示n = 10):
days Status Prot1 Prot13 Prot7 Prot21
Subj_1 115.69 0 2.284498 6.319168 6.070115 8.457412
Subj_2 72.30 1 2.473034 6.066573 6.140178 8.225987
Subj_3 1.08 1 2.662481 6.212845 6.971018 8.128949
Subj_4 69.63 1 2.761391 5.902610 6.433883 7.876319
Subj_5 78.41 1 3.038122 6.355257 6.852981 7.500973
Subj_6 42.90 1 2.058549 6.020681 7.231307 8.164025
Subj_7 31.00 1 2.305096 5.415107 8.126941 8.566320
Subj_8 51.12 1 2.931978 5.574601 7.503275 7.529957
Subj_9 11.01 1 2.218814 6.270222 6.710297 8.193895
Subj_10 27.68 1 2.821947 6.132379 6.911071 8.428218
问题是:我该如何创建一个公式,将这些患者分为两类:一组在一年内的估计生存率<60%,另一组将包括在同一时间段内的估计生存率> 60%?
R中是否有任何function()
可以解决这个问题?
非常感谢。
答案 0 :(得分:0)
我认为您应该在此处发布此问题
https://stats.stackexchange.com
因为这是统计问题。无论如何,您可以尝试从二项式回归开始,但是您可以尝试许多其他模型。你有几门科目?