Prova della relazione tra tasso di pericolosità, densità di probabilità, funzione di sopravvivenza


11

Sto leggendo un po 'sulle analisi di sopravvivenza e la maggior parte dei libri di testo lo afferma

h(t)=limΔt0P(t<Tt+Δt|Tt)Δt=f(t)1F(t)(1)

dove è la percentuale di rischio,h(t)

f(t)=limΔt0P(t<Tt+Δt)Δt(2) la funzione di densità,

F(t)=Pr(T<t)(3) e

S(t)=Pr(T>t)=1F(t)(4)

Inoltre lo affermano

S(t)=e0th(s)ds(5)

La maggior parte dei libri di testo (almeno quelli che ho) non fornisce prove per (1) o (5). Penso di essere riuscito a superare (1) come segue

limΔt0P(Tt|t<Tt+Δt)P(t<Tt+Δt)h(t)=limΔt0P(t<Tt+Δt|Tt)Δt= limΔt0P(Tt|t<Tt+Δt)P(t<Tt+Δt)P(Tt)Δt che a causa di (2) e (4) diventa limΔt0P(Tt|t<Tt+Δt)f(t)S(t)Δt ma P(Tt|t<Tt+Δt)=1 quindi h(t)=f(t)1F(t)

Come si dimostra (5)?


5
Hai notato che è la derivata di ? - log S ( t )h(t)logS(t)
Stéphane Laurent,

Sì, non capisco neanche questo ...
nostock

Nella tua prova di (1), dovresti prima sostenere che la seconda probabilità nel numeratore è 1, quindi applicare (2) e (4).
Ocram,

Perché l'ordine è importante?
nostock,

1
Se mantieni l'ordine, dovresti sostenere che il limite come (piuttosto che la proba stessa) è uguale a . Comunque, questo è un dettaglio ...1Δt01
Ocram,

Risposte:


15

La derivata di è Pertanto, come menzionato da @ StéphaneLaurent, abbiamo dove l'ultima uguaglianza segue da (1).d S ( t )S-dlog(S(t))

dS(t)dt=d(1F(t))dt=dF(t)dt=f(t)
dlog(S(t))dt=dS(t)dtS(t)=f(t)S(t)=h(t)

Prendendo l'integrale su entrambi i lati della relazione precedente, otteniamo modo che S ( t ) = exp { - t 0 h ( s )

log(S(t))=0th(s)ds
S(t)=exp{0th(s)ds}

Questa è la tua equazione (5). La parte integrale dell'esponenziale è il rischio integrato, chiamato anche rischio cumulativo [in modo che ].S ( t ) = exp (H(t)S(t)=exp(H(t))


Potresti essere un po 'più esplicito su
dlog(S(t))dt=dS(t)dtS(t)
nostock,

1
Questa è la regola della catena. Abbiamo modo chedlog(x)dx=1x
dlog(f(x))dx=df(x)dxx
ocram

La x nella parte destra dell'ultima equazione dovrebbe essere f (x) ?, cioè Per differenziare y = log S (t). Sia u = S (t) quindi . Inoltre, abbiamo e quindi . Secondo la regola della catena, quindi
dudt=dS(t)/dt=S(t)
y=logS(t)=log(u)
dydu=1u=1S(t)
dydt=dydududt=1S(t)S(t)=S(t)S(t)
user1420372

@ user1420372: Sì, hai ragione. Avrebbe dovuto essere f (x).
Ocram,

3

h(t)=f(t)S(t) 
=f(t)1F(t)
=f(t)10tf(s)ds

Integra entrambi i lati: Differenzia entrambi i lati:

0th(s)ds=0tf(s)10tf(s)dsds
=ln[10tf(s)ds]0t+c
10tf(s)ds=exp[0th(s)ds]
f(t)=h(t)exp[0th(s)ds]
f(t)=h(t)exp[0th(s)ds]

Poiché

h(t)=f(t)S(t)

S(t)=f(t)h(t)

Sostituisci con , Pertanto, f(t)h(t)exp[0th(s)ds]

S(t)=h(t)exp[0th(s)ds]h(t)
S(t)=exp[0th(s)ds]

3

Dimostriamo la seguente equazione: proof:

S(t)=exp{0th(u)du}

Dimostriamo innanzitutto proof:

f(t)=dS(t)dt

f(t)=dF(t)dt=dP(T<t)dt=d(1S(t))dt=dS(t)dt 
E sappiamo che Sostituisci in otteniamo quindi continua la nostra prova principale. Integrando entrambi i lati dell'equazione precedente, abbiamo Quindi otteniamo il risultato
h(t)=f(t)S(t)
f(t)h(t)
h(t)=dS(t)dtS(t)
S ( t ) = exp { - t 0 h ( u ) d u }
0th(u)du=0tdS(t)dtS(t)dt=0tS(t)1dS(t)=[logS(t)logS(0)]=logS(t)
S(t)=exp{0th(u)du} 
Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.