问题描述
我有一个生存分析模型,该模型具有时间相关的协变量。我可以使用以下代码在R中对此建模:
coxph(Surv(FUTIME,status) ~ exposure + tt(exposure),tt = function(x,t,...) x * t,data = dataset)
但是,如果我发现有相互作用,我希望能够在不同时间计算出HR和95%CI。
为进行比较,这里是SAS:
proc phreg data=dataset;
model FUTIME*status(0)=exposure time1/risklimits ties=efron;
time1=FUTIME*exposure;
contrast 'expo 20' exposure 1 time1 20/estimate=exp;
run;
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)