采用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）