问题标签 [survival-analysis]

下面是一个最小的例子：

在 for 循环中的两个图中，仅生成第二个（通常：仅生成最后一个）。问题似乎与 survplot 具体有关：如果我们在其中一个 survplot 之后添加一个普通图，例如：

或者：

所有地块都在您期望的地方生成。

有谁知道这里发生了什么？

提前致谢，文森特

我有以下 survreg 模型：

我想根据上述估计绘制危险函数和生存函数。
我不想使用predict()或pweibull()（如此处介绍的参数生存或此处的SO 问题。

我想使用该curve()功能。有什么想法可以做到这一点吗？似乎 survreg 的 Weibull 函数使用了与通常不同的比例和形状定义（并且不同于例如 rweibull）。

更新：我想我真正需要它来表达作为估计函数的危险/生存，Intercept而不使用任何现成的函数。age (+ other potential covariates)Scale*weilbull

我在分析我的数据集时考虑了正确的删失数据。我正在使用该survival软件包 - 给出了癌症治疗策略以及患者最后一次在我的客户诊所登记的时间。

标准包是否有建议的方法或操作来survival解释右删失数据？

我们的行是独特的个体患者...

以下是我们填写的栏目：

• 项目清单
• 我们的治疗类型（常数）
• 自最初诊断后的天数
• “censored”是这一天最后一次听到的患者人数。因此，我们现在不确定他们在停止就诊时是否还活着或死了。它们应该在未来的所有时间点从概率估计中删除。
• # 当天死亡的患者人数（来自原始诊断）

那么您是否建议对标准生存包进行操作？还是使用另一个包？我已经看到了survSNP，这可能会有所帮助。我想避免重新计算各个列/字段并创建 R 算法的新对象，因为这是一个非常大的数据集的一小部分。survPRESMOOTHsurvBIVAR

我认为这应该很容易，但我无法完全理解它。

我有以下代码：

但我希望将结果分解为各个年份。

以下是 SPSS 语法和解决方案的样子：

和 STATA 中的相同：

这是输出

我们实际上是在尝试在 R 中重现模型的结果，该模型已在 SAS 中编码。模型如下所示：ln(Duration)=X'B+S*e，其中 X 是 10 个自变量的矩阵，B 是系数向量，S 是尺度参数，e 是误差项。

我们使用的数据集在这里

在那里您也可以找到 SAS 代码。

第一次尝试如下所示：

但是这个模型中的系数是不正确的。在下图中，您可以看到左侧的 R 输出和右侧的正确 SAS 输出：

我们检测到平方项（Acq_Expense_SQ，Ret_Expense_SQ）存在问题，因为当我们排除这些项时，所有其他估计都更接近正确值。因此，我们尝试将平方项缩小 0.001 倍。

现在，系数更接近正确的系数，但我不知道为什么。这个问题有可能的解释吗？或者您是否发现我们的代码存在另一个问题？

我们正在尝试在 R 中重现加速故障时间 (aft) 模型的结果，该模型已在 SAS 中编码。

我们使用的数据集在这里

在那里您也可以找到 SAS 代码。

我们的问题是如何定义 err.mat 术语？err.mat 指定具有测量误差的变量。由于我们的数据集是正确删失的，我认为具有测量误差的变量可能是 Duration 和/或 Censor。但并不是那么简单，err.mat 必须是一个对称的平方数字矩阵。

可能是一个愚蠢的问题，但是如何在卡普兰迈耶图中更改符号又名标记颜色？我知道lines.survfit 可以选择更改实际标记符号，但显然不是符号的颜色（只有线条的颜色）

我尝试使用 line() 或 points() 手动添加标记，并使用 pch 和 col 更改这些标记的颜色和符号。但是标记不完整；并非所有被审查的事件都有使用该方法的标记。

有谁知道该怎么做？

如果我在 R 中有一个 survdiff 对象，是否有一种快速的方法来确定哪个组的生存率最高（如果有显着差异）？

基本上，我循环遍历数据帧的行并每次获取 survdiff p 值，但如果特定组显示出最佳生存率，我真的只对 p 值感兴趣。

我正在研究心脏移植数据，我必须在R. 由于它是为此目的而设计的，因此我选择使用该coxme软件包，但在选择模型时遇到了很多麻烦。我想使用转发过程，但它不适用于coxme对象。我在这里问了另一个问题，他的答案允许我访问step()'s 的代码来修改它，以便它适用于coxme模型。但是，这非常困难（至少对我而言），因为功能相当复杂（绝对不仅仅是修改exctractAIC()，还有许多其他功能需要适应）。

所以这是我的问题：

1) 有没有更简单的方法在模型上运行前锋coxme？

2）我应该使用另一个包吗？

我正在尝试打印包survfit创建的对象的受限平均存活率和标准错误。survival当我从我的完整数据中运行以下示例摘录时，我无法rmean在摘要中显示（在我的完整数据中的所有 1000 个 obs 上运行时也不会显示）。