cox-regression相关内容

从字符串创建公式调用

我使用最佳子集选择包来确定用于构建模型的最佳独立变量(我确实有这样做的特定原因,而不是直接使用最佳子集对象).我想以编程方式提取特征名称,并使用生成的字符串构建我的模型公式.结果将是这样的: x ..
发布时间:2021-05-09 19:13:42 其他开发

删除缺失的数据值

我删除了原始帖子,以便能够发布更大版本的数据集.实际上总共有418行. 以下是我正在进行的生存分析的数据.第一列是ID号,其他列标记为V2-V20.有许多丢失的数据用“."表示. 我使用 coxph()函数获得以下信息: #将生存时间保存在向量"time"中.时间= surv.df [,"V2"]#将关于检查状态的信息保存在向量“状态"中.在考克斯的回归中;死亡代码应始终设置为1.状 ..
发布时间:2021-04-28 20:44:53 其他开发

间隔检查数据:R的Cox比例风险和生存差异

关于如何使用 survival 软件包在R中进行生存分析的信息(互联网和教科书)很多。但是,当您左审查数据时,我没有找到有关如何执行此操作的任何信息。 问题背景: 我有一个带有公开生存数据的自建数据集。通常会给出事件的时间和最近一次随访的日期(正确的检查)。但是,有一项研究仅指出该事件发生在360天之前。所以我左删掉了此数据。 想要做的事情: 我想分析带有左截断,事件和右截 ..
发布时间:2020-06-08 18:55:13 其他开发

使用Sbrier(R)估计Cox生存模型的预测准确性

格拉夫(Graf)等人(1999)在一篇论文中提出了综合的Brier分数(IBS),作为衡量生存模型预测准确性的良好方法(请参见例如 Wiering等人的综述文件。,第23页)。 它在软件包 ipred 中作为函数 sbrier 。但是,虽然brier分数定义显然适用于Cox比例风险模型,但我无法让 sbrier 返回 coxph 模型。 这是问题所在。 图书馆(生存) 图书馆( ..
发布时间:2020-06-08 18:55:10 其他开发

使用与时间有关的系数和样条曲线绘制来自Coxph对象的估计HR

在 coxph 模型具有基于时间的系数(基于样条项)的模型的情况下,我想将估计的危险比作为时间的函数进行绘制。我使用函数 tt 创建了时间相关系数,类似于本示例直接来自?coxph : #使用当前年龄 拟合时间变换模型cox = coxph(Surv(time,status)〜ph.ecog + tt(age),data = lung, tt = function(x,t,...)pspl ..
发布时间:2020-06-08 18:55:03 其他开发

R coxph()警告:Loglik在变量前收敛

我在使用coxph()时遇到了一些麻烦。 我有两个分类变量:“性别”和“可能的原因”,我想用作预测变量。性别只是典型的男性/女性,但可能的原因有5种选择。 我不知道警告消息是什么问题。为什么置信区间是从0到Inf,而p值这么高? 以下是代码和输出: > my_coxph ..
发布时间:2020-06-08 18:54:55 其他开发

如何用样条曲线绘制Cox危害模型

我有以下模型: coxph(Surv(长度,死亡率== 1)〜pspline(预测变量)) 其中fulength是随访时间(包括死亡率),预测因子是死亡率的预测因子。 上面命令的输出是这样的: coef se(coef)se2 Chisq DF p pspline(预测变量),行0.174 0.0563 0.0562 9.52 1.00 0.002 ..
发布时间:2020-06-08 18:54:52 其他开发

Coxph预测与系数不匹配

下午好, 我可以发布可复制的代码,并且如果每个人都同意出了点问题,我当然会的,但是现在我想我的问题很简单,有人会为我指明正确的路径。 我正在处理以下数据集: created_as_free_user tc 1是36 0 2是36 0 3是0 1 4是28 0 5是9 0 6是0 1 7是13 0 8是 ..
发布时间:2020-06-08 18:54:47 其他开发

R中用于对过程表单数据进行计数的分层对数秩检验?

背景:在半年的随访时间(为期4年)中,患者可能会改用其他药物治疗。为了解决这个问题,我将生存数据转换为计数过程表格。我想比较药物组A,B和C的生存曲线。我正在使用扩展的Cox模型,但想对每个危险函数进行成对比较或进行分层对数秩检验。我认为 pairwise_survdiff 会由于我的数据形式而引发错误。 示例数据: x x 问题: 在 survdiff ..
发布时间:2020-06-08 18:54:45 其他开发

coxph()X矩阵被认为是奇异的;

我在使用coxph()时遇到了一些麻烦。我有两个类别变量:“ tecnologia”和“ pais”,我想评估“ pais”对“ tecnologia”的可能的交互作用。“ tecnologia”是一个具有2个级别的变量因子:gps和convencional。而“ pais”分为2个级别:PT和ES。我不知道为什么这个警告不断出现。 这是代码和输出: cox_AC ..
发布时间:2020-06-08 18:54:42 其他开发

在ggplot中绘制连续协变量的预测生存曲线

如何在cox比例风险模型中绘制连续协变量代表值的生存曲线?具体来说,我想在ggplot中使用“ survfit.cox”“ survfit”对象进行此操作。 这个问题似乎已经得到了解答,但是我用“ survfit”和“ newdata”(以及许多其他搜索)搜索了SO中的所有内容条款)。到目前为止,这是最接近回答我的问题的线索: Plot Kaplan -Meier进行Cox回归 与该帖 ..
发布时间:2020-06-07 18:42:42 其他开发