我想进行治疗分析，当我运行下面的代码时：

报此错误：

rep(1, n) 中的错误：无效的“次”参数

我正在尝试使用rms包绘制受限三次样条模型。但是我找不到任何方法来调整我的 cox 比例风险模型，我只能得到未调整的拟合。这是我的代码：

通过这个编码，我得到了未调整的样条模型。我想知道如何添加混杂变量来调整模型。我试过了：

但这给了我每个变量的样条模型，任何人都知道如何调整我的模型？

我试图找出，我是否有足够的能力使用 R 进行 Cox 比例风险模型。我的预测指标是智力，我的结果是癌症。我的样本量是 22,000，我预测 25 年内有 600 例癌症病例。此外，我想知道，我是否有足够的能力来测试性别相互作用，如果我的样本中有一半是女性，而 600 例癌症病例中有 200 例是。我预测风险比至少为 1.3。我在 R 中找到的函数都需要二进制预测器，但智能是度量标准，我找不到任何函数。

我正在使用 Survival 包进行 cox 回归（coxph），但我无法弄清楚我做错了什么。

我的代码是

但我不断收到错误消息： eval 中的错误（predvars，data，env）：'character' 类型的无效 'envir' 参数

我已经检查了所有变量的类。age_ent 和 outage 是数字，a_menopause 和糖尿病是整数。

非常感谢您对此的任何帮助！

谢谢！

我试图通过在使用变量交互时绘制 cox 回归来绘制调整后的生存曲线。

阅读 survfit.coxph 页面https://stat.ethz.ch/R-manual/R-devel/library/survival/html/survfit.coxph.html我看到参数“newdata”

我想在我的 cox 输出中绘制作为交互的线条。即，如果我的 cox 输出看起来像：

我想在与Sex.

这让我想到了这个newdata参数。

与调用 newdata 相比，不包括 newdata 和仅使用协变量的平均值有什么区别。在这一点上，我什至不知道如何正确构建 newdata。

如果有人可以给我任何关于如何newdata根据我的 cox 模型构建的指示，以及与仅使用平均值相比有什么意义。在绘制基于 cox 数据的新生存图时，我应该期望原始生存曲线中的线条数量相同。

我是R中的新手，我想进行cox回归分析，如图所示，如何设置参考组并将HR更改为1，其他是参考组的具体值？

非常感谢，圣诞快乐！

我想为以下 Cox 比例风险模型绘制 ROC 曲线。

Cox 比例风险模型 - 使用时变协变量

我尝试了survivalROC 包。

有什么建议吗？你知道其他的包吗

coxph在我的数据上运行时出现错误 。你们中有人见过吗？

我的数据：

我有兴趣从 clogit 模型中获取设定位置的拟合值。这包括总体水平响应及其周围的置信区间。例如，我的数据大致如下所示：

在 Clogit 模型中，活动在层内没有变化，因此没有活动主效应。

我想要做的是建立一个新的数据框架，我最终可以在 Open 的指定位置绘制边缘拟合值，类似于传统 glm 模型中的新数据设计：例如，

但是，当我尝试在 clogit 上运行预测函数时，出现以下错误：

Surv (rep(1, 5000L), Used) 中的错误：找不到对象“Used”

我意识到这是因为 r 中的 clogit 正在通过 Cox.ph 运行，因此，预测功能试图预测同一层内的成对受试者之间的相对风险（在这种情况下 = 已使用）。

然而，我的问题是是否有办法解决这个问题。这很容易在 Stata 中完成（使用 Margins 命令），并在 Excel 中手动完成，但是我想在 R 中自动化，因为其他所有东西都是在那里编程的。我也在 R 中手动构建了这个（下面的示例代码），但是我一直在我的真实数据中发现似乎不正确的 CI，因此我希望尽可能依赖预测函数。我的手动预测代码是：

理想情况下，我的最终产品看起来像这样，但它是 predict 函数的产品......

拟合值的示例输出。

我有这个生存数据，描述了十年时间研究中三种服务（沙龙、餐厅和快递）的死亡率。

数据包含三个变量：服务类型（1=saloon、2=restaurant 和 3=express）、年数（1 到 11 的整数，其中 11 表示大于 10 年）和审查员。

我有两个问题：

1）我已经拟合了 Cox 比例风险模型，但是有什么方法可以检查比例风险假设。也就是说，我们假设每个人的风险比和基线风险与时间无关。

2）如何在 R 中拟合时间相关的 Cox 模型？

这是我的代码：