我进行了 Cox 回归分析，包括四个变量（性别、年龄和两个二元解释变量），它们都与结果有显着关联。我使用了 R 中“survival”包中的 coxph 函数：

我想制作一个预测得分表，其中包括四个变量和 4 个年龄分层组（30、40、50、60 岁）。此表中的所有危害必须除以一个预定义的危害，以获得每个年龄组的 HR。

如何计算 R 中每个特定年龄组的 95% CI 的 HR？

我想在表中估算一些缺失的数据，并在估算的表上运行 Cox 模型。

我可以在我的数据上运行插补，并在插补数据上运行 cox 模型，但我不明白如何查看数据集中的 cox 输出，其中一些值是插补的（即我特别需要我输出中的风险比和 P 值）。

命令是：

然后我确保我的相关变量是因素（例如，队列可以是 0 或 1，以确保这些被视为不同的类别）。

然后我重新调整因素以使 Cox 模型在以后更容易解释：

然后我对数据进行了估算：polyreg 用于具有两个以上级别的分类数据，logreg 用于具有 3 个级别的因子。

然后，我运行 Cox 模型以在估算数据集上运行：

输出是 5 个 cox 模型分析。我无法将信息集中在一起。当我输入“pool(cox_with_imp)”时，它会给我一些统计数据。但我想要一个包含 HR 和 P 值的“汇总”表。

有谁知道我输入的命令，将 5 个估算的 Cox 模型汇集到一个具有 HR 和 P 值的共识 Cox 模型中。

谢谢。

在我的生存分析任务中，我使用 cox 比例模型来计算数据集不同组中的一致性指数 (c-index) 值。我想知道如何将 p 值添加到我的 c-index 图中以比较不同的组看起来像这个图？

这是我的代码：

提前致谢，

在 R 中运行 Cox PH 回归后，我需要在数据框中添加预测风险比列。数据框是面板数据，其中 numgvkey 如果公司标识符和年龄是时间标识符。您可以从此链接下载一小部分日期： https ://drive.google.com/file/d/0B8usDJAPeV85VFRWd01pb0h1MDA/view?usp=sharing

我做了以下事情：

如何在我的 'sme' 数据框中添加第 5 列“危险比”（hr）？另外，有什么方法可以预测 EVENT2 的概率而不是“hr”？

我正在尝试使用 cox 模型来预测一段时间后的失败概率（称为停止）3。

但是，预测函数的输出都不在 0-1 范围内。是否有任何函数或如何使用 lp 预测和基线风险函数来计算概率？

每次我运行 model<– coxph(enter,exit, event)~X1+X2+ GDP时，会话突然崩溃而没有错误消息。只有当我添加 gap 变量时它才会崩溃。

我想知道这是否与GDP变量有关？我尝试了不同的来源，它的会话再次崩溃。似乎具有更高值的变量（例如“223945604003”或“39455684e-5”）是问题所在？

我在 Mac Yosemite 版本 10.10.3 上使用最新版本的 R-Studio MAX OS 10.6（64 位）

感谢您的各种评论。

coxph对于具有基于样条项的时间相关系数的模型，我想将估计的风险比绘制为时间的函数。我使用函数创建了时间相关系数tt，类似于直接来自的示例?coxph：

调用survfit(cox)会导致survfit无法理解带有tt术语的模型的错误（如 Terry Therneau 在 2011 年所述）。

您可以使用 提取线性预测器cox$linear.predictors，但我需要以某种方式提取年龄，并且不那么琐碎地提取每个时间。因为tt在事件时间上拆分数据集，所以我不能只将输入数据框的列与coxph输出匹配。此外，我真的很想绘制估计函数本身，而不仅仅是对观察到的数据点的预测。

这里有一个涉及样条的相关问题，但不涉及tt.

编辑 (7/7)

我仍然坚持这一点。我一直在深入研究这个对象：

所以，我有结，但我仍然不确定如何将coxph系数与结结合起来以实际绘制函数。任何线索都非常感谢。

我正在使用 cox 模型进行预测。

在 predict.coxph 文档中，它说预期类型是“给定协变量和后续时间的预期事件数”。这里的跟进时间是什么意思？我怎样才能改变它？

1) 在以下示例中，如何将 y 轴更改为“适合”的“优势比”、“死亡率”和“死亡率”？

2) 在以下示例中，如何将“fit2”的 y 轴更改为“风险比”？

我期待着您的帮助！

更新 1 在 BondedDust 的回答之后，我写了以下内容：

由于总的死亡概率是 44.5%，所以预测概率的计算和结果图对我作为非统计学家来说似乎是正确的，不是吗？

我使用 MICE 估算数据并使用生存运行 Cox 模型。

输出是：

我知道 SE 与 est 相关，所以我知道要获得 CI 来进行估计：

（这就是 lo 95 和 hi 95 列的含义）。我可以做哪些具体的计算来获得 HR 的 CI（即 exp(est)）？

我可以这样做：

(exp(est) -1.96 X se, exp(est) + 1.96 X se)为 HR 获取 CI 而不是 est？

说白了，如果我想用这个表来发布结果（HR=X, 95% CI=X, P Val=X），我要怎么从上表中去呢？

谢谢