我想对我的数据进行一次定期生存分析,这是关于一家公司遭受网络破坏的信息。数据包括IT预算,一般财务变量,事件时间,以10k报告形式披露的风险因素。我有兴趣调查风险因素披露和IT预算如何影响被破坏的生存时间。所以这是我的代码:
在这里,start
是t,stop
是t + 1,status
是事件,gvkey
是索引,CTL_sic
是行业代码
survival1 <- coxph(
Surv(start, stop, status) ~ numwords + PE_it_budget + CTL_maketcap + CTL_roe + CTL_liq + CTL_mb + cluster(gvkey) + strata(CTL_sic),
method = "breslow",
data = combinedata)
我正确设置吗?
另外,如何在模型中指定numwords
和IT预算之间的交互项?
数据结构大致如此。