如何解决如何计算R中具有时间相互作用的危险比?
我有一个生存分析模型,该模型具有时间相关的协变量。我可以使用以下代码在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 举报,一经查实,本站将立刻删除。