arxiv: v3 [stat.me] 6 Feb 2018

Transcription

1 TRANSFORMING CUMULATIVE HAZARD ESTIMATES PÅL C RYALEN, MATS J STENSRUD, AND KJETIL RØYSLAND arxiv: v3 [same] 6 Feb 218 Deparme of Biosaisics, Uiversiy of Oslo, Domus Medica Gausad, Sogsvasveie 9, 372 Oslo, Norway Absrac Time o eve oucomes are ofe evaluaed o he hazard scale, bu ierpreig hazards may be difficul Recely, here has bee cocer i he causal iferece lieraure ha hazards acually have a buil i selecio-effec ha preves simple causal ierpreaios This is eve a problem i radomized corolled rials, where hazard raios have become a sadard measure of reame effecs Modelig o he hazard scale is everheless coveie, eg o adjus for covariaes Usig hazards for iermediae calculaios may herefore be desirable Here, we provide a geeric mehod for rasformig hazard esimaes cosisely o oher scales a which hese buil i selecio effecs are avoided The mehod is based o differeial equaios, ad geeralize a well kow relaio bewee he Nelso-Aale ad Kapla-Meier esimaors Usig he marigale ceral limi heorem we also fid ha covariaces ca be esimaed cosisely for a large class of esimaors, hus allowig for rapid calculaios of cofidece iervals Hece, give cumulaive hazard esimaes based o eg Aale s addiive hazard model, we ca obai may oher parameers wihou much more effor We prese several examples ad associaed esimaors Coverage ad covergece speed are explored usig simulaios, suggesig ha reliable esimaes ca be obaied i real-life scearios 1 Oulie Applied researchers are ofe faced wih ime o eve oucomes Causal iferece for such oucomes is challegig due o he dyamic aspec of ime, wih possibly ime-varyig effecs ad biases I paricular, he proporioal hazard model is he sadard approach for aalyzig survival daa i medicie, bu assumig proporioal hazards cao be jusified i may real-life scearios Furhermore, sice esimaig hazards ivolves codiioig o rece survival, he presece of a seemigly ioce uobserved heerogeeiy meas ha we Dae: February 7, 218 1

2 2 TRANSFORMING CUMULATIVE HAZARD ESTIMATES will codiio o a so-called collider, ad will herefore acivae ocausal pahways from he exposure o he risk of a eve a shor erm [1] These commoly used effec esimaes ca herefore ofe o be ierpreed causally as shor-erm risks [1, 2, 3, 4, 5, 6] However, o evaluae ime o eve oucomes, we do o eed o assess hazards per se We may raher be ieresed i effec measures o eg he survival scale Such measures could be easier o ierpre, ad allow for causal evaluaios [6, 7] Sill, modelig he hazard scale may be useful as a iermediae sep, eg o adjus for covariaes Addiive hazard models seem o become more popular, a leas i he causal iferece lieraure [8, 9] I coras o he proporioal hazards model, he addiive models are collapsible ad hey easily allow for ime-depede effecs The addiive models are also less proe o selecio biases whe umeasured facors follow liear srucural models [1,11] O he oher had, addiive models have bee criicized because he effec esimaes may be harder o udersad: Ofe effec esimaes are direcly ploed as cumulaive hazard curves, which seems o be usaisfacory o applied researchers [12] Ideed, he cumulaive hazard curves may eiher have a immediae causal ierpreaio Wih his i mid, we aim o rasform cumulaive hazard esimaes o oher scales by exploiig srucure imposed by ordiary differeial equaios We will esimae parameers ha solve such equaios, which are ypically drive by cumulaive hazards The esimaors we cosider are soluios of aurally associaed sochasic differeial equaios ha are sraigh forward o solve umerically o a compuer We are mosly cocered wih such equaios ha are drive by cumulaive hazards, bu our mehod is o a all limied o his seig Our approach eables us o apply a simple varia of he dela-mehod o prove ha applyig our rasformaios o cosise cumulaive hazard esimaors also yield cosise esimaors Furhermore, we are able o provide esimaes of he asympoic variace I fac, his is much simpler ad allows us o go furher ha he fucioal delamehod suggesed by Richard D Gill e al [13, 14], sice he laer ivolves opologies ha someimes eed o be carefully desiged for each example, which may be hard o hadle i pracice; see eg [15] 2 Parameers The parameers we will cosider are fucios o a fiie ierval [, T ] ha are soluios o a ordiary differeial equaio sysem o he form 1 X = X + F X s da s Here A is a coiuous k-dimesioal fucio of bouded variaio, ad F : R M,k R is sufficiely smooh ad saisfies a liear

3 TRANSFORMING CUMULATIVE HAZARD ESTIMATES 3 growh boud There are several examples of parameers o his form i he survival aalysis lieraure, bu 1 has o received much aeio I he followig we show examples of parameers ha ca be wrie i his way 21 Survival curves Le T be a survival ime ad le S := P T, ie S is he survival fucio If A deoes he correspodig cumulaive hazard, he we have he relaio S = S S s da s 22 Relaive survival Whe assessig effecs of exposures, i is iriguig o cosider he probabiliy of survival direcly, raher ha esimaig eg hazard raios The relaive survival is rouiely used i cliical medicie [16] I paricular, five year relaive survival raes are sadard measures of he survival i cacer paies vs he geeral populaio, ad such umbers are coveioally repored from cacer regisries Suppose ha we are give wo groups 1 ad, ad we wa o ideify he relaive survival RS := S 1 /S of he posiive group, relaively o he egaive group a ime If we le A 1 ad A deoe he correspodig cumulaive hazards, he a simple calculaio of derivaives shows ha S1 S = RS S1 S RS + S1 s Ss RS s RS s d A 1 s A s 23 Mea ad resriced mea survival Isead of sudyig he relaive survival as a fucio of ime, we may be ieresed i he mea survival i differe populaios Ideed, he mea survival is simply he area uder he survival curve I pracice, however, he ail of he esimaed survival fucio will srogly ifluece he iegral of he esimaed survival curve, ad he mea survival esimaes may be ureliable As a remedy, i is possible o sudy he resriced mea survival E [ T IT τ ], ha is he area uder he survival curve up o some resriced ime τ, R τ I may be ierpreed as he life expecacy bewee = ad τ Advocaes claim ha he resriced mea survival o a prespecified τ should be repored i cliical rials [17, 18, 19], i paricular whe he proporioal hazards assumpios is ivalid Ideed, R τ is readily foud by solvig he sysem Rτ = + S τ 1 τ Ss d S s where A is he cumulaive hazard for deah s, As

4 4 TRANSFORMING CUMULATIVE HAZARD ESTIMATES 24 Life expecacy differece ad life expecacy raio The life expecacy differece LED ad life expecacy raio LER ca be useful o compare resriced mea survivals i differe groups [2] The LED bewee wo groups 1 ad 2 is he differece bewee heir resriced mea survival imes I solves he sysem LED12 S 1 = 1 + S 2 1 S1 s Ss 2 Ss 1 d s A 1 Ss 2 s A 2 s Here, S 1, S 2 ad A 1, A 2 are he survival, ad cumulaive hazards i group 1 ad 2 respecively The life expecacy raio is he raio of he resriced mea survival imes I solves he sysem LER 12 1 S 1 S 2 R 1 = R 2 where R 2 s cao be oo small S 1 s R2 s S2 s R1 s R 2 s 2 S 1 s S 2 s S 1 s S 2 s d s A 1 s, A 2 s 25 Cumulaive icidece ad compeig risks I pracice we are ofe ieresed i he ime o a paricular eve, bu our sudy may be disruped by oher, compeig eves I medicie, we may eg be ieresed i he ime o he ose of a disease D, bu subjecs may die before hey develop D If here are compeig risks, i is icorrec o assume a oe-o-oe relaio bewee cause-specific hazards ad he cumulaive icidece [21]; if we eg rea deah as cesorig, we obai esimaes from a hypoheical populaio i which subjecs cao die wihou D Aleraively, we may iclude he hazards of he compeig eves io our model Cosider a siuaio wih k compeig causes of deah, wih cause-specific cumulaive hazards A 1,, A k Followig Aderse eal [22], we see ha 2 S C 1 C k = 1 S s S s S s + S s d S s describes all he cause-specific cumulaive icideces C 1,, C k, subjec o he remaiig compeig risks A 1 s A k s

5 TRANSFORMING CUMULATIVE HAZARD ESTIMATES 5 26 Mea frequecy fucio for recurre eves Cosider idividuals ha may experiece a umber of recurre eves, eg hospial visis H, while beig a risk of experiecig some ermiaig eve, eg deah D We le N H be he couig process ha cous he umber of recurre eves up o ime, ad A H, A D be he cumulaive hazards for he recurre ad ermiaig eve, respecively The mea frequecy fucio K = E[N D H ] he solves he sysem 3 K = + S 1 Ss A H d s S s A D s A H ad A D ca be esimaed usig Nelso-Aale esimaors, ad we ca esimae K by pluggig hem io 3 Asympoic resuls for his esimaor have previously bee described [23] The cosisecy resul will be a cosequece from a more geeral resuls i our paper 27 Cumulaive sesiiviy ad specificiy afer screeig The effeciveess of disease screeig is frequely sudied i medicie Ofe he screeig is performed a =, bu he disease may advace or develop a imes > Hece, subjecs may be followed over ime o assess if eves occur, eg he ose of disease I such scearios, coveioal mehods o assess sesiiviy ad specificiy are iadequae, because hey do o accou for differeces i follow-up imes We may raher express sesiiviy ad specificiy as fucios of ime, as suggesed i [24] Le A 1 be he hazard of he disease eve for he group wih posiive screeig resuls, le A be he hazard of he egaive oes, le β deoe he prevalece of he disease ad le Z deoe he screeig oucome We le U deoe he cumulaive posiive predicive value P N > Z = 1, le V deoe he cumulaive egaive predicive value P N = Z =, le W deoe he cumulaive sesiiviy P Z = 1 N > ad le X deoe he cumulaive specificiy P Z = N = Applyig Bayes rule gives he equaio 4 U V W X = U V W X + 1 U s V s Ws 2 1 V s 1 U s γus 2 γxs 2 1 U s V s W s 2 V s U s γus 2 γx2 s 1 U s V s A d 1 s A s, where we have iroduced he odds γ = cao be oo small β 1 β, ad where V s ad U s

6 6 TRANSFORMING CUMULATIVE HAZARD ESTIMATES 3 Plugi esimaors Our suggesed esimaors will build o a cumulaive hazard esimaor ÂÂÂ for A ha is give by couig process iegrals, ie such ha 5 ÂÂÂ = G s dn s, where G = g 1 g l is a predicable k l dimesioal marixvalued process, ad N is a adaped l-dimesioal couig process I should be oed ha he resuls here is applicable for more geeral siuaios ha 5 Suppose, however, ha his esimaor is cosise i he sese ha 6 lim P ÂÂÂ s A s ɛ =, sup s T for every ɛ >, or such ha 7 W := ÂÂÂ A coverges i law o a idepede icremes mea zero Gaussia local marigale W relaive o he Skorohod-meric The laer is ypically a cosequece of he marigale ceral limi heorem, as i [14, Theorem II51] or [25, Theorem VIII322] Two sadard examples i eve hisory aalysis where his holds are he Nelso-Aale cumulaive hazard esimaor [14, Theorem IV11 ad Theorem IV12] ad Aale s addiive hazard regressio [14, Theorem VII41] I addiio o cosisecy, here is aoher propery of ÂÂÂ ha will be crucial i his seig This is predicably uiformly ighess, abbreviaed P-UT The exac defiiio of P-UT is give i [25, VI6a] However, as we will o eed he full geeraliy, we ca use he followig Lemma o deermie if processes are P-UT i mos siuaios we ecouer Lemma 1 Le {Z } be a sequece of semi-marigales o [, T ], ad le {ρ } be predicable processes such ha every M lim sup J := Z ρ s ds defies a square iegrable local marigale, ad suppose ha i P ρ s 1 J = ad ii 8 lim sup J The {Z } is P-UT sup s T P Tr M T J =

7 TRANSFORMING CUMULATIVE HAZARD ESTIMATES 7 We are ow ready o defie he esimaors for parameers ha are defied hrough 1 Defiiio 1 Le ˆXˆXˆX be a esimaor for X The plugi-esimaor ˆXˆXˆX for X ha saisfies 1 is he soluio of he followig sochasic differeial equaio: 9 ˆXˆXˆX = ˆXˆXˆX + F ˆXˆXˆX s dâââ s Plugi esimaors are relaively easy o impleme o a compuer due o heir recursive form If τ 1, τ 2, deoe he jump imes of ÂÂÂ, he we have ha ha 1 ˆXˆXˆX = ˆXˆXˆX 1 + F ˆXˆXˆX 1 ÂÂÂ, wheever < +1 These esimaors are also cosise i may siuaios This is formulaed as a heorem, ad is proved i he Appedix Theorem 1 Suppose ha i X saisfies he ordiary differeial equaio 1, where F is locally Lipschiz-coiuous, ad saisfies he liear growh codiio o a domai ha coais {X [, T ]}, ii {ÂÂÂ } ad { ˆXˆXˆX } are cosise, ie lim P sup s T ÂÂÂ s A s ɛ = ad lim P ˆXˆXˆX X ɛ =, for every ɛ >, iii ÂÂÂ is P-UT see Proposiio 1 for addiive hazards, The ˆXˆXˆX 11 lim P sup s T s X s ɛ = for every ɛ >, ie { ˆXˆXˆX } defies a cosise esimaor of X Remark 1 We ca also hadle he siuaio where A = Cosider a oe-dimesioal deermiisic couig process N wih jumpig imes τ1 < τ2 < such ha 12 lim sup k Le h := ad oe ha, due o 12, 13 Â := τ k τk 1 = {, = sup{ɛ N ɛ = N ad ɛ < }, >, h s dn s = 1 = max{ }

8 8 TRANSFORMING CUMULATIVE HAZARD ESTIMATES defies a process ha saisfies 6, ie lim sup T Â =, ad Theorem 1 holds If ÂÂÂ A defies a local marigale ha saisfies he marigale ceral limi heorem, he he sequece { ˆXˆXˆX X } coverges o a soluio of a give SDE ha is drive by a Gaussia process This is he coe of he followig heorem A proof ca be foud i he appedix Theorem 2 Suppose ha i F := F 1,, F k have bouded ad coiuous firs ad secod derivaives o a domai ha coais {X [, T ]}, ii Z := ˆXˆXˆX X coverges i law o a mea zero Gaussia radom vecor Z as, iii W := ÂÂÂ A coverges i law relaive o he Skorohodmeric owards a mea zero Gaussia marigale wih idepede icremes ad is P-UT see Proposiio 1 for addiive hazards, iv {[W ]} is P-UT see Proposiio 1 for addiive hazards, Le Z deoe he uique soluio o he sochasic differeial equaio 14 Z = Z + k j=1 F j X s Z s da j s + Now, Z is a mea zero Gaussia process, ad Z := ˆXˆXˆX X F X s dw s coverges i law relaive o he Skorohod-meric o Z as Moreover, le V deoe he covariace for Z, suppose ha {ˆVˆVˆV } is a cosise esimaor for V, ad le ˆVˆVˆV deoe he soluio of he sochasic differeial equaio: 15 ˆVˆVˆV =ˆVˆVˆV + + k j=1 ˆVˆVˆV s F j ˆXˆXˆX s + F j ˆXˆXˆX s ˆVˆVˆV s dâj s F ˆXˆXˆX s d[âââ ] s F ˆXˆXˆX s Now, ˆVˆVˆV defies a cosise esimaor of V, ie 16 lim P ˆVˆVˆV s V s ɛ = sup s T

9 TRANSFORMING CUMULATIVE HAZARD ESTIMATES 9 for every ɛ >, ad V solves he followig ordiary differeial equaio: k V =V + V s F j X s + F j X s V s da j s j= F X s d[w ] s F X s We ca compue covariaces ad herefore also poi-wise cofidece iervals for our esimaes based o he ODEs from Secio 1, by pluggig i ˆXˆXˆX for X ad ÂÂÂ for A Noe ha, aalogously o 9, i is sraigh forward o solve he sochasic differeial equaio 15, sice i is drive by jump-processes If τ 1, τ 2, deoe he jump imes of ÂÂÂ, he we have ha 18 ˆVˆVˆV =ˆVˆVˆV 1 + k j=1 ˆVˆVˆV 1 F j ˆXˆXˆX 1 + F j ˆXˆXˆX 1 ˆVˆVˆV τ Âj k 1 + F ˆXˆXˆX 1 ÂÂÂ ÂÂÂ F ˆXˆXˆX 1 wheever < +1 The properies iii of Theorem 1 ad iii ad iv of Theorem 2 are ideed saisfied for Aales addiive regressio uder relaively weak assumpios Proposiio 1 Assume ha we have iid idividuals ad ha he iesiy of each N i, a ime, is o he form λ i = Y i α U, i where α is bouded ad coiuous Wrie U for he marix whose i h row is equal o U i, ad Y for he diagoal marix wih i h diagoal eleme is equal o Y i Suppose ha 1 E[sup T U i 3 3] < for each i, 2 U lim if P Y U 1 sup Tr J = J T The 19 A := U s Y s U s 1 U s Y s dn s saisfies iii of Theorem 1 ad iii ad iv of Theorem 2 4 Plugi esimaor examples I he followig we wrie ou he plugi esimaors for some of he examples we showed i Secio 2, usig he oaio iroduced i Defiiio 1 The mos ivolved variace expressios are moved o

10 1 TRANSFORMING CUMULATIVE HAZARD ESTIMATES he Appedix We will wrie := 1, ad use e i,j o deoe marices wih eries equal o zero i all posiios apar from he ih row ad jh colum where he value is 1, i order o keep he oaio simple Simulaios of seleced parameers are displayed i Figures 1 ad 2 41 Survival curves The survival plugi esimaor is Ŝ = Ŝ 1 Ŝ 1 Â, which is ohig bu he Kapla-Meier esimaor expressed as a differece equaio The variace plugi esimaor is ˆV = ˆV τk 1 2 ˆV τk 1  + Ŝ 1  2 42 Relaive survival The relaive survival differece equaio is Ŝ 1 Ŝτ 1 Ŝ k 1 = Ŝτ Ŝ1 1 Â1 k 1 + Ŝ 1 τk RS RS  τk 1 RSτk 1 RS τk 1 The plugi variace solves ˆVˆVˆV = ˆVˆVˆV τk 1 ˆVˆVˆV e τk 1 1,1 + e 3,3 + e 1,1 + e 3,3 ˆVˆVˆV τk 1 Â1 + ˆVˆVˆV e τk 1 3,3 e 2,2 + e 3,3 e 2,2 ˆVˆVˆV τk 1  Ŝ1 1 Â1 + Ŝ 1 τk RS   1  Ŝ1 1 RS τk 1 Ŝ τ RSτk 1 k 1 τk 1 RSτk 1 43 Mea ad resriced mea survival The resriced mea survival plugi esimaor solves ˆR ˆRτk 1 Ŝτk 1 τk = + Ŝ Ŝ τk 1 Ŝ 1  τk The variace expressio reads ˆVˆVˆV = ˆVˆVˆV τk 1 + ˆVˆVˆV e τk 1 2,1 + e 1,2 ˆVˆVˆV τk 1 + ˆVˆVˆV e τk 1 2,2 + e 2,2 ˆVˆVˆV τk 1  Ŝτk 1 τk + Ŝ 1  τk  τk Ŝ τk 1 Ŝ 1

11 TRANSFORMING CUMULATIVE HAZARD ESTIMATES Life expecacy differece ad life expecacy raio The LED plugi esimaor is LED 12 LED 12 τ k 1 Ŝ 1 τ Ŝ 1 = Ŝτ 1 k 1 Ŝ2 1 k 1 + Ŝ1 1  1, Ŝ 2 Ŝ2 1  2 Ŝ 2 1 wih variace equaio ˆVˆVˆV = ˆVˆVˆV τk 1 + ˆVˆVˆV e τk 1 2,1 e 3,1 + e 1,2 e 1,3 ˆVˆVˆV τk 1 ˆVˆVˆV e τk 1 2,2 + e 2,2 ˆVˆVˆV τk 1 Â1 ˆVˆVˆV e τk 1 3,3 + e 3,3 ˆVˆVˆV τk 1 Â1 + F τk 1  1  1  2  2 F τk 1 Here F is he LED iegrad fucio, ie Ŝτ 1 k 1 Ŝ2 1 F τk 1 = Ŝτ 1 k 1 The LER plugi esimaor reads LER 12 LER 12 Ŝ 1 1 Ŝ 1 1 Ŝ 2 = Ŝ 2 τ ˆR 1 k 1 + ˆR 1 ˆR 2 1 ˆR τ 2 k 1 Ŝ 2 1 Ŝ 1 1 ˆR2 τk 1 Ŝ2 1 ˆR1 τk 1 ˆR Ŝ1 1 Ŝ2 1 Ŝτ 1 k 1 Ŝτ 2 k 1 The variace expressio is ivolved, ad ca be foud i he Appedix  1  2 45 Cumulaive icidece ad compeig risks The plugi esimaor for he cumulaive icidece example is Ŝ Ĉ 1 Ĉ m = Ŝ τk 1 Ĉτ 1 k 1 Ĉ m 1 + Ŝ 1 Ŝ 1 Ŝ 1 Ŝ τk 1 Ŝτk 1  1  m,

12 12 TRANSFORMING CUMULATIVE HAZARD ESTIMATES The plugi variace equaio is ˆVˆVˆV = ˆVˆVˆV τk 1 + ˆVˆVˆV e τk 1 1,2 e 1,1 + e 2,1 e 1,1 ˆVˆVˆV τk 1 Â1 + ˆVˆVˆV e τk 1 1,3 e 1,1 + e 3,1 e 1,1 ˆVˆVˆV τk 1 Â2 + ˆVˆVˆV e τk 1 1,m+1 e 1,1 + e m+1,1 e 1,1 ˆVˆVˆV τk 1 + F τk 1  1  m  1  m F τk 1 Âm 46 Mea frequecy fucio for recurre eves We ge a plugi esimaor equaio ha reads ÂH ˆK ˆKτk 1 Ŝτk 1 = + τk Ŝ Ŝ τk 1 Ŝ 1  D The variace is ˆVˆVˆV = ˆVˆVˆV τk 1 + ˆVˆVˆV e τk 1 2,1 + e 1,2 ˆVˆVˆV τk 1 ˆVˆVˆV e τk 1 2,2 + e 2,2 ˆVˆVˆV τk 1 ÂH Ŝτk 1 + τk Ŝ 1  D  H  D Ŝ τk 1 Ŝ 1 47 Cumulaive sesiiviy ad specificiy afer screeig We have 1 Û Û Û 1 τk 1 ˆV Ŵ = ˆV τk 1 ˆV τk 1 Ŵ τk 1 + Ŵτ 2 1 ˆV τk 1 1 Û 1 k 1 Ŵ 2 τk 1 ˆVτk 1 Û τk 1 γû τ 2 γû 2 k 1 1 ˆX ˆX τk 1 γ ˆX Û 1 ˆV τk 1 γ ˆX Û 1 ˆV τk 1 The variace expressio ca be foud i he Appedix 5 Performace 51 Covergece Our covergece assessme ivolves wo seps I sep oe we fid close approximaios o he parameers, usig he oaio X I sep wo we calculae plugi esimaors for a rage of sample sizes We he check how fas he plugi esimaors coverge o X as he sample size icreases For performig sep oe we have developed code ha geeraes survival daa for give ipu hazards Usig he hazards we ca calculae Â1 τk Â,

13 TRANSFORMING CUMULATIVE HAZARD ESTIMATES 13 X o desired precisio For isace, kowig he hazard α we ca obai he survival e α sds by umerical iegraio For performig sep wo we simulae daa, eablig us o esimae cumulaive hazards, ad herefore parameers specified by 9 This sep is performed for a rage of sample sizes For each sample size we simulae k populaios, such ha we obai plugi esimaes ˆX,j, for j = 1,, k We check covergece usig he followig L 2 crierio: 2 L = 1 k k T j=1 X s,j ˆX s 2 ds For assessig covergece of he plugi variace esimaors we firs obai a variace esimae from a large boosrap sample As i sep wo above, we simulae daa ad obai k plugi variace esimaes for a rage of sample sizes I his way he same crierio 2 ca be used Covergece plos are show i Figure 3 Figure 3 suggess ha he LED esimaor performs worse ha he oher esimaors This happes because he variace depeds o he ime scale, sice we are approximaig Lebesgue iegrals The LER plugi variace also depeds o he ime scale, bu here he ime coribuio is smaller 52 Coverage To esimae coverage we firs obai close approximaios o he parameers, X, as described i sep oe i Secio 51 Nex, we simulaed daa o calculae plugi esimaes, which were used o obai cofidece iervals This sep is repeaed uil we have a large collecio of cofidece iervals A ime, a Beroulli rial decides wheher a cofidece ierval covers X We esimae he expeced coverage a ime by his Beroulli probabiliy, ie by he average umber of cofidece iervals ha cover X We calculae coverage for wo scearios; cosa hazards, ad liear crossig hazards over a ime period [, T ] We seleced liear hazards ha were large iiially bu liearly decreasig, or small iiially bu liearly icreasig, such ha hey crossed each oher a he halfway poi T /2 Coverage is show i Figure 4, suggesig ha he coverage drops below he cofidece level i some scearios This behavior is due o he plugi variace esimaors The Survival plugi variace, for isace, will drop drasically if here are eves whe few people are a risk, causig he Survival cofidece ierval o be arrow hereafer Icideally, he Greewood esimaor would give similar performace, sice i is small wheever he Kapla Meier curve is small Overall he plugi esimaors behave saisfacory, o oly whe he hazards are cosa bu also whe he hazards are liearly crossig

14 14 TRANSFORMING CUMULATIVE HAZARD ESTIMATES 6 Sofware R sofware for calculaig plugi esimaors, illusraed o simulaed daa examples, ca be foud o hps://gihubcom/palryale/rasformhazards For quesios or commes regardig he shared code, coac pcryale@medisiuioo 7 Fudig The auhors were all suppored by he research gra NFR239956/F2 - Aalyzig cliical healh regisries: Improved sofware ad mahemaics of ideifiabiliy 71 Proof of Lemma 1 Appedix: Proofs Proof We sar by provig ha {M } is P-UT Suppose ha H is predicable, marix-valued ad such ha TrH sh s bouded by 1 Noe ha for every opioal soppig ime τ: E [ =E [ Tr τ ] [ τ H s dm = E Tr τ s 2 2 H s dm s τ H s dm ] s ] [ τ ] H s d M s H s E TrH sh s dtr M s E [ τ ] [ ] dtr M s = E Tr M τ The Leglar-iequaliy [25, I-331] herefore implies ha {M } is P-UT For he firs erm i we have ha Var ρ s ds = ρ s 1 ds is uiformly igh for every, which by [25, VI 612] implies ha {ρ } is P-UT Fially, he claim follows sice {Z } is he sum of wo sequeces ha are P-UT 72 Proof of Theorem 1 Proof We are give a series of SDEs ha coverges o a SDE acually jus a ODE, ad wa o assure ha his meas ha heir soluios also coverge o he soluio of he limi equaio This is a claim abou sabiliy of sochasic differeial equaios We wa o prove our claim wih a applicaio of [25, Theorem IX69] Traslaed io our siuaio, i says ha if ˆX,ÂÂÂ coverges i law o X, A ad he series {ÂÂÂ } is predicably uiformly igh P-UT, he X ˆ, ÂÂÂ, ˆXˆXˆX also coverges i law o X, A, X Moreover, [25, VI333] esures covergece of ˆXˆXˆX,ÂÂÂ The adverised resul follows from he coiuiy of he projecio mappig

15 TRANSFORMING CUMULATIVE HAZARD ESTIMATES Proof of Theorem 2 Proof We will firs prove ha {ẐẐẐ } coverges i law o Z by aoher applicaio of he sabiliy resul for SDEs foud i [25, Theorem IX69] The we will prove ha he covariace of Z is give by he ODE 17 Fially, we will apply he SDE-sabiliy resul oce more o see ha he soluios of he SDE 15 coverge i law o he covariace of Z To prove ha {Z } coverges i law, we sar wih he followig approximaio of he derivaive of F j i he direcio z: 21 DF j x, z := F j x + 1/2 z F j x, le DF x, z := DF 1 x, z,, DF l x, z ad oe ha we have lim DF j x, z = F j xz Now, Z =Z + =Z + =Z + = ˆXˆXˆX X k j=1 k j=1 F j ˆXˆXˆX s F j X s DF j X s, Z s dâj s + DF X s, Z s dâââ s + dâj s + F X s dw s F X s dw s F X s dw s Now [25, Theorem VI622] implies ha {Z + F X s dw s } coverges is law o Z + F X s dw s, ad sice A is coiuous, [25, Corollary VI333] implies ha {Z + F X s dw s,âââ} coverges i law o Z + F X s dw s, A Sice F has bouded ad coiuous firs as secod derivaives, a calculus argume shows ha {DF } have uiformly bouded Lipschiz cosas by he Lipcschiz cosa of F This eables us o apply he sabiliy resul for SDEs foud i [25, Theorem IX69] agai I says ha sice {ÂÂÂ } is uiformly predicably igh, we also have ha {Z + F X s dw s,âââ, Z } coverges i law o Z + F X s dw s, A, Z To see ha he covariace of Z is give by he ODE 17, we oe ha [W ] is deermiisic sice W is a coiuous Gaussia marigale

16 16 TRANSFORMING CUMULATIVE HAZARD ESTIMATES wih idepede icremes, ad ha E [ [ ] [ ] ] Z Z =E Z Z + E Z s dz s [ ] + E dz s Z s + E [ ] [Z, Z] =E [ [ ] k ] Z Z + E Z s da j sz s F j X s j=1 [ k + E j=1 F j X s Z s da j sz s [ ] + E F X s d[w ] s F X s =E [ ] k Z Z k j=1 j=1 E [ Z s Z s ] Fj X s da j s F j X s E [ ] Z s Z s da j F X s d[w ] s F X s We ow show ha he soluios of he SDE 15 coverge i law o he covariace of Z We firs oe ha [25, Theorem VI626] implies ha {[W ]} coverges i law o [W ] sice {W } is predicably uiformly igh Moreover, sice A is coiuous, [25, Corollary VI 333] implies ha {ÂÂÂ, [W ]} coverges i law o A, [W ] If his sequece was predicably uiformly igh, he we could apply [25, Theorem IX69] oe more ime o esure ha he soluios of he SDEs 15 coverge i law o V Fially, oe ha {ÂÂÂ, [W ]} is P-UT sice boh {ÂÂÂ } ad {[W ]} are P-UT 74 Proof of Proposiio 1 Proof I is well kow ha W := A A coverges weakly o a Gaussia process wih idepede icremes, see [14] To see ha {W } is P-UT, le 22 Γ ad oe ha Tr W T = T := U Tr Γ 1 Y U U λ ], Y U Γ 1 d

17 TRANSFORMING CUMULATIVE HAZARD ESTIMATES 17 From he submuliplicaiviy of he race orm we have Tr Γ 1 Tr Γ 1 U λ Y U Tr O he oher had, we have ha 23 Tr U λ Y U U λ Y Γ 1 U = 1 k,i Tr Γ 1 U k,i 2 Y i λ i, ad E P [sup T U k,i 2 Y i λ i ] <, so Markovs iequaliy implies ha 24 lim J if U P Tr λ Y U J = Fially, [14, Proposiio II53] implies ha 8 of Lemma 1 is saisfied Noe ha {A } is P-UT sice A = 1 W + A To see ha {[W ]} is P-UT, oe ha we have already see ha he compesaor of [W ], W, saisfies 8 of Lemma 1 Moreover, oe ha Tr [W ] W T = Tr W T, so 8 of Lemma 1 is saisfied for {[W ] W } as well, which herefore implies ha {[W ]} is P-UT Appedix: Oher plugi variace expressios Variace for LER Here, he plugi variace reads ˆVˆVˆV = ˆVˆVˆV τk 1 + ˆVˆVˆV τk 1 G 1 + G τk 1 ˆVˆVˆV τk 1 ˆVˆVˆV e τk 1 2,2 + e 2,2 ˆVˆVˆV τk 1 Â1 ˆVˆVˆV e τk 1 3,3 + e 3,3 ˆVˆVˆV τk 1 Â2 + F τk 1 Â 1 Â 1 Â 2 Â 2 F τk 1,

18 18 TRANSFORMING CUMULATIVE HAZARD ESTIMATES where G ad F is give by 1 ˆR τ 1 k 1 ˆR τ 2 ˆR Ŝ2 τk 1 k 1 τ 2 2 k 1 G τk 1 = 1 1 Ŝτ 1 ˆR2 k 1 τk 1 Ŝ2 τ ˆR1 k 1 τk 1 ˆR τ 2 2 k 1 F τk 1 = Ŝ1 1 Ŝ2 1 Ŝτ 1 k 1 Ŝτ 2 k 1 2Ŝ2 τk 1 ˆR1 τk 1 Ŝ1 τk 1 ˆR2 τk 1 ˆR τ 2 2 ˆR 2 k 1 τ 3 k 1 Variace for cumulaive sesiiviy ad specificiy I his case he plugi variace is Here we have defied G 1 1 = G 2 1 = F τk 1 = ˆVˆVˆV = ˆVˆVˆV τk 1 + ˆVˆVˆV τk 1 G1 1 + G 1 1 ˆVˆVˆV Â1 τk 1 + ˆVˆVˆV τk 1 G2 1 + G 2 1 ˆVˆVˆV τk 1  τk + F τk 1 Â1 τk   1  F τk 1 1 Ŵ ˆV τk 1 Û 1 2 γû 3 1 γ ˆX 2 1 ˆV τk 1 Ŵ 2 1 Û 1 1 γû 2 1 2Ŵ 1 1 Û 1 1 ˆV τk 1 γû 2 1 γ ˆX 2 τ 1 Û 1 k 1 ˆV τ 2 k 1 1 Ŵ 2 1 ˆVτk 1 γû 2 1 γ ˆX 2 1 ˆV τk 1 Ŵ 1 γû 1 2Ŵ 1 ˆV τk 1 γû 1 γ ˆX 2 τ 1 Û 1 k 1 2γ ˆX τk 1 1 Û 1 ˆV τ 2 ˆV τk 1 k 1 1 Û 1 ˆV τk 1 Ŵ ˆV τk 1 1 Û 1 γû τ 2 k 1 γ ˆX τ 2 1 Û 1 k 1 ˆV τk 1 Ŵ 2 1 ˆVτk 1 Û τk 1 γû τ 2 k 1 γ ˆX 2 τ 1 Û 1 k 1 ˆV τk 1 Refereces [1] O O Aale, R J Cook, ad K Røyslad, Does cox aalysis of a radomized survival sudy yield a causal reame effec?, Lifeime daa aalysis, vol 21, o 4, pp , 215, 2γ ˆX τk 1 1 Û 1 ˆV τk 1,

19 TRANSFORMING CUMULATIVE HAZARD ESTIMATES fietimes fietimes RelSurv S_KM = 2 fietimes fietimes fietimes W fietimes RelSurv S_KM = 1 W Cumulaive sesiiviy Relaive Survival Survival = 5 19 Figure 1 Plugi esimaes for seleced parameers wih populaio sizes 5, 1, ad 2, alog wih 95 % poiwise cofidece iervals Exac parameers is show i doed lies [2] J Robis ad S Greelad, The probabiliy of causaio uder a sochasic model for idividual risk, Biomerics, pp , 1989 [3] S Greelad, Absece of cofoudig does o correspod o collapsibiliy of he rae raio or rae differece, Epidemiology, pp , 1996 [4] M J Sesrud, M Valberg, K Røyslad, ad O O Aale, Explorig selecio bias by causal fraily models: The magiude maers, Epidemiology, vol 28, o 3, pp , 217 [5] O Aale, O Borga, ad H Gjessig, Survival ad eve hisory aalysis: a process poi of view Spriger Sciece & Busiess Media, 28 [6] M A Hera, The hazards of hazard raios, Epidemiology Cambridge, Mass, vol 21, o 1, p 13, 21 [7] M A Hera ad J M Robis, Causal iferece CRC Boca Rao, FL:, 217 [8] T Mariusse, S Vaseelad, M Gerser, ad J v B Hjelmborg, Esimaio of direc effecs for survival daa by usig he aale addiive hazards model, Joural of he Royal Saisical Sociey: Series B Saisical Mehodology, vol 73, o 5, pp , 211 [9] T Lage ad J V Hase, Direc ad idirec effecs i a survival coex, Epidemiology, vol 22, o 4, pp , 211 [1] S Srohmaier, K Røyslad, R Hoff, Ø Borga, T R Pederse, ad O O Aale, Dyamic pah aalysis a useful ool o ivesigae mediaio processes i cliical survival rials, Saisics i medicie, vol 34, o 29, pp , 215 [11] T Mariusse ad S Vaseelad, O collapsibiliy ad cofoudig bias i cox ad aale regressio models, Lifeime daa aalysis, vol 19, o 3, pp , 213

20 TRANSFORMING CUMULATIVE HAZARD ESTIMATES RMS_diff = 2 fietimes RMS_raio RMS_raio fietimes fietimes C_cumic LED LER fietimes Cumulaive icidece = 1 RMS_diff = 5 fietimes C_cumic 2 fietimes Figure 2 Plugi esimaes for seleced parameers wih populaio sizes 5, 1, ad 2, alog wih 95 % poiwise cofidece iervals Exac parameers is show i doed lies LED ad LER are abbreviaios for life expecacy differece ad life expecacy raio, respecively [12] M J Bradbur, T G Clark, S Love, ad D Alma, Survival aalysis par ii: mulivariae daa aalysis a iroducio o coceps ad mehods, Briish joural of cacer, vol 89, o 3, p 431, 23 [13] R D Gill, J A Weller, ad J Præsgaard, No-ad semi-parameric maximum likelihood esimaors ad he vo mises mehod par 1[wih discussio ad reply], Scadiavia Joural of Saisics, pp , 1989 [14] P K Aderse, Ø Borga, R D Gill, ad N Keidig, Saisical models based o couig processes Spriger Series i Saisics, New York: Spriger-Verlag, 1993 [15] A W Va der Vaar, Asympoic saisics, vol 3 Cambridge uiversiy press, 1998 [16] M P Perme, J Sare, ad J Ese ve, O esimaio i relaive survival, Biomerics, vol 68, o 1, pp , 212 [17] P Royso ad M K Parmar, The use of resriced mea survival ime o esimae he reame effec i radomized cliical rials whe he proporioal hazards assumpio is i doub, Saisics i medicie, vol 3, o 19, pp , 211 [18] P Royso ad M K Parmar, Resriced mea survival ime: a aleraive o he hazard raio for he desig ad aalysis of radomized rials wih a ime-o-eve oucome, BMC medical research mehodology, vol 13, o 1, p 152, 213 [19] L Triquar, J Jaco, S C Coer, ad R Porcher, Compariso of reame effecs measured by he hazard raio ad by he raio of resriced mea

21 TRANSFORMING CUMULATIVE HAZARD ESTIMATES 21 Parameer Variace 6 LED LER Survival CumIc L L log log Figure 3 Covergece of seleced plugi esimaors; parameer o he lef ad variaces o he righ We used dashed lies covergece order 1 i he lef pael ad 23 i he righ pael LED ad LER are abbreviaios for life expecacy differece ad life expecacy raio, respecively survival imes i ocology radomized corolled rials, Joural of Cliical Ocology, vol 34, o 15, pp , 216 [2] H-M Dehbi, P Royso, ad A Hackshaw, Life expecacy differece ad life expecacy raio: wo measures of reame effecs i radomised rials wih o-proporioal hazards, BMJ, vol 357, p j225, 217 [21] P K Aderse, R B Geskus, T de Wie, ad H Puer, Compeig risks i epidemiology: possibiliies ad pifalls, Ieraioal joural of epidemiology, vol 41, o 3, pp , 212 [22] P K Aderse, S Z Abildsrom, ad S Roshøj, Compeig risks as a mulisae model, Saisical mehods i medical research, vol 11, o 2, pp , 22 [23] D Ghosh ad D Li, Noparameric aalysis of recurre eves ad deah, Biomerics, 2 [24] M Nygård, K Røyslad, S Campbell, ad J Diller, Comparaive effeciveess sudy o huma papillomavirus deecio mehods used i he cervical cacer screeig programme, BMJ ope, vol 4, o 1, p e346, 214 [25] J Jacod ad A N Shiryaev, Limi heorems for sochasic processes, vol 288 of Grudlehre der Mahemaische Wisseschafe [Fudameal Priciples of Mahemaical Scieces] Berli: Spriger-Verlag, secod ed, 23

22 22 TRANSFORMING CUMULATIVE HAZARD ESTIMATES Survival Relaive survival Cumulaive sesiiviy = 5 = 1 = LED LER Cumulaive icidece Figure 4 Esimaed mea coverage for seleced parameers simulaed wih cosaupper pael ad liearly crossiglower pael hazards The doed lie idicaes he cofidece level LED ad LER are abbreviaios for life expecacy differece ad life expecacy raio, respecively