2.风险函数回归(多因素分析)：风险函数[hazardfunction，h(t)]即t时刻存活的观察个体在t时死亡瞬时的概率强度，公式为
CS公式为
其中CS表示t时刻未发生任何事件的观察个体发生第k类事件的瞬时概率强度。
SD公式为
其中SD表示t时刻未发生第k类事件的观察个体发生第k类事件的瞬时概率强度。
上述两式中λk(t)称为粗危险率，存在竞争事件时粗危险率≠净危险率，故采用经典Cox分析HR值及其95%CI存在偏倚。
本研究以非参数CIF及右删失半参数法CS、SD为基础，并提供了SAS9.4分析程序，临床生存数据复杂多样，总结右删失数据竞争风险模型具有良好的临床指导意义，罕见的左删失、区间删失数据、参数法、非等比例风险数据还有待进一步研究。
附录：SAS9.4竞争风险模型程序
/*cov1-covn为协变量名*/
odsgraphicson；odshtmldpi＝300；
%cif(data＝cohort，out＝cif1_data，time＝time，status＝y，event＝1，censored＝0)；/*不分组A总的CIF*/
%cif(data＝cohort，out＝cif1_data，time＝time，status＝y，event＝1，censored＝0，group＝treat)；/*手术组treat组间gray’s比较A的CIF*/
%cif(data＝cohort，out＝cif2_data，time＝time，status＝y，event＝2，censored＝0，group＝treat)；/*手术组treat组间gray’s比较B的CIF*/
/*A校正CIF，SD模型*/
procphregdata＝cohortplots(overlay＝stratum)＝cif；
modeltime*y(0)＝cov1~covn/eventcode＝1rl；run；
/*A的CS模型，y(0，2)表示删失事件，zph为Schoenfeld残差图、assess为鞅式残差图，均可验证PH比例风险假设*/
procphregdata＝cohortzph；
modeltime*y(0，2)＝cov1~covn/ties＝efronrl；assessph/resampleseed＝1234；run；
procphregdata＝cohortplots(overlay＝stratum)＝cif；/*B校正CIF，SD模型*/
modeltime*y(0)＝cov1~covn/eventcode＝2rl；run；
procphregdata＝cohortzph；/*B的CS模型*/
modeltime*y(0，1)＝cov1~covn/ties＝efronrl；assessph/resampleseed＝1234；run；