1765 0

[统计软件] R语言竞争风险模型的生存分析:按照时依变量模型进行分析时的R语句是? [推广有奖]

  • 0关注
  • 0粉丝

小学生

0%

还不是VIP/贵宾

-

威望
0
论坛币
9 个
通用积分
1.0000
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
363 点
帖子
2
精华
0
在线时间
5 小时
注册时间
2012-11-20
最后登录
2021-7-1

楼主
扛着刀的诗人 发表于 2020-2-25 00:35:27 |AI写论文
5论坛币
采用r软件包timereg自带的骨髓移植数据(bmt)进行分析,在bmt数据中有408个样本,5个变量,其中,time表示生存时间,cause是数值编码变量(1:死亡,结局时间;0:截尾事件;2:竞争风险(复发或其他原因)),platelet是数值编码变量,其中1表示超过100*10^9per L,0表示低于上述值。tcell表示T-cell depleted BMT,也为数值编码变量(1:yes,0:no),age表示中心化年龄变量,所做的变换是(真实的年龄-35)/15.
注:竞争风险事件的生存分析模型中,截尾数据(研究中止)编码为0,结局时间编码1,其他原因引起的编码为2.
本例采用竞争风险模型进行分析。

install.packages(pkgs="timereg")

install.packages(pkgs="cmprsk")

library(timereg)

library(cmprsk)

data(bmt)

bmt

library(timereg)

library(cmprsk)

先按照时依变量模型进行分析。输入下面的语句,但怎么也运算不出来,各位大神能帮忙看一下吗,应该怎样编辑,谢谢!:

outf <- comp.risk(Surv(time,cause == 0)~platelet+age+tcell,data=bmt,bmt$cause,cause=1,n.sim=5000,cens.code=0,model="prop",cens.model="cox")

提示参数没有用(cens.code=0


关键词:生存分析 风险模型 R语句 R语言 Packages

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
jg-xs1
拉您进交流群
GMT+8, 2025-12-30 12:12