R与STATA中的COX比例风险模型
我正在尝试使用以下数据在R中复制来自STATA的Cox比例风险模型估计http://iojournal.org/wp-content/uploads/2015/05/FortnaReplicationData.dta Stata中的命令如下： stset enddate2009, id(VPFid) fail(warends) origin(time startdate) stcox
..
survival-analysis相关内容
如何绘制survreg生成的生存曲线(R的包生存)?
我正在尝试将 Weibull 模型拟合并绘制到生存数据中.数据只有一个协变量，即队列，运行时间为 2006 年到 2010 年.那么，对于绘制 2010 年队列生存曲线的以下两行代码有什么想法吗? 图书馆(生存)s
..
将函数与 R 中的向量参数集成
我对之前的帖子有类似的挑战:如何将向量传递给整合功能 我有一个函数，我想对曲线下的面积进行积分. 一、[生存]功能: surv score 来自风险计算器，它调整生存估计.患者有不同的分数，例如: score 计算所有 7 位患者的 surv 很容易，如果我们有一个特定的时间点 x: surv(5, score) # 生存到第 5 年[1] 0.7161497 0.691
..
使用 R 中“rpart"包中的生存树来预测新的观察结果
我正在尝试使用 R 中的“rpart"包来构建生存树，我希望使用这棵树然后对其他观察结果进行预测. 我知道有很多关于 rpart 和预测的 SO 问题；但是，我还没有找到任何可以解决(我认为)特定于将 rpart 与“Surv"对象一起使用的问题. 我的特殊问题涉及解释“预测"函数的结果.一个例子很有帮助: 库(rpart)图书馆(OIsurv)# 制作数据:set.seed(4)d
..
R - cox 风险模型不包括因子水平
我正在将 cox 模型拟合到一些结构如下的数据中: str(test)'data.frame':147 个观察.共 8 个变量:$ 年龄:int 71 69 90 78 61 74 78 78 81 45 ...$ Gender : Factor w/2 levels "F","M": 2 1 2 1 2 1 2 1 2 1 ...$ RACE :因子 w/5 个级别 "","BLACK","H
..
如何预测R中的生存概率?
我将名为 veteran 的数据存储在 R 中.我创建了一个生存模型，现在希望预测生存概率.例如，具有 80 karno 值、10diagtime、age 65 和 prior=10 和 trt = 2 寿命超过 100 天? 在这种情况下，设计矩阵是 x = (1,0,1,0,80,10,65,10,2) 这是我的代码: 图书馆(生存)附上(老手)weibull
..
在 `ggsurv` 图中(或在 `plot` 中)区分不同类型的每条线
我正在使用 Rstudio.我正在使用 GGally 包中的 ggsurv 函数为我的数据绘制 Kaplan-Meier 曲线(用于生存分析)，来自教程 此处.我使用它而不是 plot 因为 ggsurv 自己处理图例. 如链接所示，多条曲线按颜色区分.我想根据线型进行区分.本教程似乎没有任何选项.以下是我的命令: surv1 lty.est=3(或 2)为所有行提供相同的虚线.我希望
..
游泳者生存图
有没有一种简单的方法可以在 R 中生成游泳者图?与 KM 曲线中的数据相同，但每个个体的生存都用一条线表示.示例: 我搜索了 stackoverflow、R-help 邮件列表，并咨询了 Google 博士，但没有得到明显的答案，尽管我的搜索技术可能不是最理想的.谢谢！ **** 添加 ****为不恰当地提问而道歉 - 这是我第一次！玩玩，我已经能够做到以下几点: OS Death
..
根据Survreg预测绘制生存曲线
我在生存分析方面还比较陌生，并且已经使用了一些标准的电信用户流失数据示例，下面的示例称为“电信": 电信公司--read.csv(文本=“状态，帐户长度，区域代码，国际计划，Day_Mins，Day_Calls，Day_Charge，Eve_Mins，Eve_Calls，Eve_Charge，Night_Mins，Night_Calls，Night_Charge，Intl_Mins，Intl_M
..
将字符串数据包装在ggplot2的surv_miner的风险表中
对于在RStudio中用surv_miner构造的Kaplan Meier曲线，我具有以下风险表: Kaplan Meier曲线 我想在风险表中将文本字符串换行，以使数字位于顶部，而百分比(用括号括起来)位于数字下方的另一行，以便它们都适合表.我尝试使用stringr包，但我不太了解如何将其合并到主题中.任何帮助，将不胜感激.我只发布了我的代码，因为添加实际的生存数据有些棘手.但是，如
..
在调整后的生存曲线ggadjustedcurves(survminer，ggplot2)中自定义线型
我正在尝试绘制一条调整后的生存曲线，但是努力地按组更改线型.我可以使用典型的ggplot2语言自定义绘图的其他方面，但是我遇到了线型变化的问题. 示例: 图书馆(生存)图书馆(survminer)fit2
..
在R中找到函数的最大值
我具有以下功能. 并调整 xlim 以找到最大值 df = curve(surviveFunction，从= 0，到= 45，xlim = c(30,40))plot(df，type ='l')
..
如何从R中的Coxph模型摘要中提取公式?
我能够按照本教程中的概述复制coxph模型和提取组件: http://www.sthda.com/english/wiki/cox-proportional-hazards-model .但是，我正在努力从模型中提取公式对象. library("survival")库("survminer")数据(“肺")头(肺)协变量
..
R:绘制个人预测
我正在使用R编程语言.我正在尝试按照本教程进行操作: https://rdrr.io/cran/randomForestSRC/man/plot.competing.risk.rfsrc.html 本教程显示了如何使用“生存随机森林".algorithm-用于分析生存数据的算法.在该示例中，“叶面"为白色.使用数据集，使用生存随机森林算法来分析观察到的经历“状态1"的即时危险.对比“状态2"
..
删除缺失的数据值
我删除了原始帖子，以便能够发布更大版本的数据集.实际上总共有418行. 以下是我正在进行的生存分析的数据.第一列是ID号，其他列标记为V2-V20.有许多丢失的数据用“."表示. 我使用 coxph()函数获得以下信息: #将生存时间保存在向量"time"中.时间= surv.df [，"V2"]#将关于检查状态的信息保存在向量“状态"中.在考克斯的回归中；死亡代码应始终设置为1.状
..
拼写数据管理：过去24个月在给定状态下花费的月份
我正在使用具有以下形式的咒语数据集： 清除所有 输入persid start end t_start t_end spell_type年spell_number事件 1 8 9 44 45 1 1999 1 0 1 12 12 60 60 1 2000 1 0 1 1 1 1 61 61 1 2001 1 0 1 7 11 67 71 1 2001 2 0 1 1 4
..
R NeuroNet:“不一致的参数".
啊！尝试使用神经网络compute时，始终出现以下错误: > net.compute
..
WinBUGS Weibull网络元分析
我目前正在进行多项临床试验中生存数据的荟萃分析. 为此，我使用相同的方法从发布的分析中获取代码.但是，当使用已发布的分析中的数据运行此代码时，我无法复制其结果.实际上，结果无法收敛到任何合理的估计. 代码本身(不包括数据)应正确无误，因为它直接来自作者.我认为问题必须与初始值或 采样运行方式的参数，但在多次播放之后 初始值，老化时间，变薄等...我没有得到有意义的结果. 任何人都
..
如何在R中模拟左截断的Weibull故障时间数据
我想模拟来自Weibull分布的左截断的故障时间数据. 我的目标是通过拟合Weibull回归模型来模拟数据并检索系数(我用于模拟的x1，x2，x3，x4和x5).在这里，xt=runif(N, 30, 80)表示研究的开始，Tm
..
使用curve()绘制Survreg的生存率和危险函数
我有以下survreg模型: Call: survreg(formula = Surv(time = (ev.time), event = ev) ~ age, data = my.data, dist = "weib") Value Std. Error z p (Intercept) 4.0961 0.5566 7.36 1
..