Lecture 2: Martingale theory for univariate survival analysis
|
|
- Joella Bradley
- 5 years ago
- Views:
Transcription
1 Lecture 2: Martingale theory for univariate survival analysis In this lecture T is assumed to be a continuous failure time. A core question in this lecture is how to develop asymptotic properties when studying statistical methods for univariate survival data? Empirical process approach - a general tool for asymptotic theory Martingale theory - enjoy some advantages in variance simplifications for right-censored data; widely used!
2 2.1 Notation. f(t): density function of T F (t) = P (T t) = t f(u)du: cumulative distribution function S(t) = 1 F (t) = f(u)du: survival function of T t S(t) = exp{ Λ(t)}, Λ(t): cumulative hazard function Λ(t) = logs(t) λ(t) = Λ (t): hazard function C: censoring time
3 more notation: i = 1, 2,..., n: index for subjects X i = min(t i, C i ): observed survival time (possibly censored) i = I(T i < C i ): censoring indicator X (1), X (2),..., X (k) : ordered uncensored times R(t) = {(X j, j ) : X j t, j = 1, 2,..., n}: risk set at t Y i (t) = I(X i t): at-risk indicator Y (t) = Y i(t): total no. of subjects at risk at t Z i or Z i (t): (time-varying) covariates N i (t) = I(X i t, i = 1): count subject i s failure event prior to time t N(t) = N i(t): total number of observed events prior to time t {(X i, i, Z i ) : i = 1, 2,..., n}: collected data
4 2.2 Martingale theory: Initial setting Probability space (Ω, F, P ) Ω: sample space F : a class of subsets of Ω; the class is a σ-algebra P : a probability measure on (Ω, F ) Conditional expectation E[Y X] 1. E[Y X] is σ(x)-measurable 2. E[Y X] = E[Y σ(x)] - average prediction of Y given all the information on X 3. E[Y ] = E[E[Y X]]
5 Stochastic process {X(t); t τ} on (Ω, F, P ) 1. a collection of random variables indexed by t τ 2. for a fixed sample ω Ω, X(t; ω) is a sample path as a function of t; for convenience, we usually write X(t) instead of X(t; ω)
6 - Define F t as the σ-field generated by {X(u) : u t}: information about X( ) from to t - {W (t); t } is adapted to F t if W (t) is F t -measureable - Define F as the collection of F t, t τ. It is called the filtration of the underlying probability space A stochastic process M(t) is a martingale with respect to the stochastic process X(t), t τ, if - M(t) is adapted to filtration F : for each t, the random variable M(t) is a F t -measurable function - E[M(t)] < - for s, E[M(t + s) F t ] = M(t) (fair game)
7 Under the same setting, for s - M(t) is submartingale if E[M(t + s) F t ] > M(t) (winning) - M(t) is supermartingale if E[M(t + s) F t ] < M(t) (losing) If M(t) is a martingale, then - E[dM(t)] = - if E[M()] =, then E[M(s)] =, s Predictable process W (t) - W (t) is F t -measurable - E[W (t) F t ] = W (t) - left-continuous process adapted to F t is predictable - example: at-risk process E[Y (t) F t ] = Y (t)
8 2.3 Doob-Meyer Decomposition Theorem Doob-Meyer Decomposition Theorem If X(t) is an adapted, right-continuous, non-negative submartingale, then there exists a unique left-continuous, increasing, predictable process A(t) such that A() =, E[A(t)] < and Q(t) = X(t) A(t) is a martingale. A(t) is called a compensator. da(t) = E[dX(t) Ft ]
9 Example. - X(t) = N i (t) is a sub-martingale. - Y i (t) = I(X i t) is left-continuous with left-hand limit and is a predictable process. - Define a predictable process A(t) = t Y i(u)λ(u)du. - Then, da(t) = Y i (t)λ(t)dt = E[dN i (t) F t ] - Q(t) = M i (t) = N i (t) t Y i(u)λ(u)du is a martingale. - We can show by a more direct proof that M i (t) is a martingale in next few slides.
10 Consider a Hilbert space H which is a complete metric space with respect to the distance function induced by the inner product < X, Y >= E[XY ]. The norm of the inner product for X H based on the inner product <, > is defined as X = < X, X > = (E[X 2 ]) 1/2. martingale: M(t) = N(t) t Y (u)λ(u)du variance: < M, M > (t) = var[m(t)] = E[M(t) 2 ] Define: d c < M, M > (t) = E[d(M(t)) 2 F t ]
11 2.4 variance and covariance of martingales Calculation of E[dN i (t) F t ], E[dM i (t) F t ] dn i (t) = or 1 E[dN i (t) F t ] = P r(dn i (t) = 1 F t ) - given F t, Y i (t) is known (predictable) - if Y i (t) = (i.e., subject i has failed or been censored before t), then P r(dn i (t) = 1 F t ) = - if Y i (t) = 1 (i.e., subject i at risk at t), then P r(dn i (t) = 1 F t ) = λ(t)dt under independent censoring. Exercise. Prove this result. E[dN i (t) F t ] = Y i (t)λ(t)dt E[{dN i (t) Y i (t)λ(t)} F t ] = E[dM i (t) F t ] =
12 M i (t) = N i (t) t Y i(u)λ(u)du is a martingale. proof. - M i (t) is adapted to F - E[M i (t)] = < - Note that for s, [{ t+s E[M i (t + s) F t ] = E dn i (u) t t+s = = = t+s t } ] Y i (u)λ(u)du F t + M i (t) E [dn i (u) Y i (u)λ(u)du F t ] + M i (t) t t+s [ ] E[E dn i (u) Y i (u)λ(u)du F u F t ] + M i (t) t t+s E[ F t ] + M i (t) = M i (t) t M(t) = N(t) t Y (u)λ(u)du is a martingale.
13 From now on, use martingales: M i (t) = N i (t) t Y i(u)λ(u)du M(t) = N(t) t Y (u)λ(u)du Note that - E[M i (t)] =, t τ; E[dM i (t) F t ] =, < t τ - E[M(t)] =, t τ; E[dM(t) F t ] =, < t τ
14 Calculation of variance and covariance: d c < M i, M j > (t) d c < M i, M j > (t) = E[d(M i (t)m j (t)) F t ] d[m i (t)m j (t)] = M i (t)m j (t) M i (t )M j (t ) Thus, = [M i (t ) + dm i (t)][m j (t ) + dm j (t)] M i (t )M j (t ) = M j (t )dm i (t) + M i (t )dm j (t) + dm i (t)dm j (t) E[d(M i (t)m j (t)) F t ] = E[dM i (t) dm j (t) F t ] = cov[dm i (t), dm j (t) F t ] Note: M i (t) = M j (t) d c < M i, M i > (t) = E[d(M i (t)) 2 F t ] = var[dm i (t) F t ]
15 d c < M i, M i > (t) = Y i (t)λ(t)dt + o p (dt) - dm i (t) = dn i (t) Y i (t)λ(t)dt - Given F t, Y i (t)λ(t)dt is a constant term - var[dm i (t) F t ] = var[dn i (t) F t ] = E[(dN i (t)) 2 F t ] E[dN i (t) F t ] 2 = E[dN i (t) F t ] E[dN i (t) F t ] 2 = Y i (t)λ(t)dt Y i (t)[λ(t)dt] 2 = Y i (t)λ(t)dt + o p (dt) Thus, d c < M i, M i > (t) = var[dm i (t) F t ] = Y i (t)λ(t)dt + o p (dt), this is a predictable variation. (Remark: o p (dt)/dt as dt o p (dt) )
16 if j k, d c < M j, M k > (t) = o p (dt) M j (t) = N j (t) t Y j(u)λ(u)du M k (t) = N k (t) t Y k(u)λ(u)du d c < M j, M k > (t) = cov[dm j (t), dm k (t) F t ] = E[(dN j (t) Y j (t)λ(t)dt)(dn k (t) Y k (t)λ(t)dt) F t ] = E[dN j (t)dn k (t) F t ] Y j (t)λ(t)dt E[dN k (t) F t ] E[dN j (t) F t ] Y k (t)λ(t)dt + Y j (t)y k (t)(λ(t)dt) 2 = E[dN j (t)dn k (t) F t ] Y j (t)y k (t)(λ(t)dt) 2 = E[dN j (t)dn k (t) F t ] + o p(dt) Note that E[dN j (t)dn k (t) F t ] = as long as N j (t) and N k (t) do not have jump at the same t (with positive probability). Thus, for continuous failure time models, E[dN j (t)dn k (t) F t ] =.
17 2.5 Martingale Central Limit Theorem (MCLT) Property. Suppose M(t) is a martingale and H(t) is predictable. Then L(t) = t H(u)dM(u) is martingale subject to F t. Recall ˆΛ(t) = t dn(u) Y (u), U n (t) = n[ˆλ(t) Λ(t)] = t 1 n n [dn i (u) Y i (u)λ(u)du] Y (u) n t n t n = Y (u) dm i(u) = Y (u) dm(u) for t < τ = sup t {t : P r(x t) > }.
18 What is the limiting distribution of U n (t)? An application of MCLT Assumptions: 1. convergent variance: < U n, U n > (t) v(t) 2. U n (t) is smooth (condition skipped) Results: 1. U n (t) converges weakly to U(t) 2. U(t) is a Brownian motion E[U(t)] = var[u(t)] = lim n < U n, U n > (t) independent increment: U(s) and U(t) U(s) are independent for s t
19 Use the results of the variance and covariance calculations, - if u < v, E[ n Y (u) dm i(u) n Y (v) dm j(v) F v ] = n Y (u) dm i(u) E[ n Y (v) dm j(v) F v ] = - E[ n Y (u) dm i(u) n Y (u) dm i(u) F u ] = n Y (u) var[dm 2 i (u) F u ] = n Y (u) Y 2 i (u)λ(u)du [ ] t n Also, < U n, U n > (t) = E Y (u) 2 Y i(u)λ(u)du [ t ] n = E Y (u) 2 Y (u)λ(u)du [ t ] λ(u)du t λ(u)du = E Y (u)/n E[Y 1 (u)] = v(t) where v(t) = t λ(u)du E[Y du = t 1(u)] λ(u)du S c(u)s(u)
20 Thus, n[ˆλ(t) Λ(t)] converges weakly to Brownian motion U(t), where E[U(t)] = and var[u(t)] = t λ(u)du S. c(u)s(u) Note that S(t) = e Λ(t), and the Kaplan-Meier estimator satisfies Ŝ(t) = e t dn(u) Y (u) + o p (n 1/2 ) (reference: Breslow and Crowley, 1974). By the functional Delta method, n[ Ŝ(t) S(t)] = n[e ˆΛ(t) e Λ(t) ] + o p (1) d S(t) U(t) (weak convergence)
21 Results: 1. n[ŝ(t) S(t)] converges weakly to V (t) = S(t) U(t) 2. V (t) is a mean zero Gaussian process E[V (t)] = var[v (t)] = [S(t)] 2 t λ(u)du S c(u)s(u)
22 Remark: 1. Estimation of var[v (t)]: df t U (v) t var[v (t)] = [S(t)] 2 S X (v) S c (v)s(v) = df [S(t)]2 U (v) S 2 (v) X Thus, by plugging in the KM and empirical distribution estimators, var[v (t)] can be estimated by t d var[v (t)] = [Ŝ(t)]2 ˆF U (v) Ŝ 2 (v), X which is approximately the Greenwood s formula: ṽar(v (t)) = [Ŝ(t)]2 t dn(u) t n 1 Y (u) (Y (u) dn(u)) = [Ŝ(t)]2 d ˆF U (u) Ŝ (u)ŝx (u+ X )
23 2. Both survival and censoring distributions plays roles in var[v (t)]. Note: The censoring distribution does not play a role in the likelihood function, but does play a role implicitly in estimation inference.
24 2.6 Proportional hazards model Now consider the proportional hazards model, λ(t Z) = λ (t)exp(βz), where Z is a 1-dimensional covariate. For simplicity of notation, we consider time-independent covariates only. The results can also be extended to time-dependent-covariates cases. L = (i) p(z (i) H(x (i) ), x (i)) {a residual likelihood}, where H(t ) = data history prior to t. Partial likelihood L p = (i) p(z (i) H(x (i) ), x (i)) = n ( e βz i j=1 yj(ti)eβz j Partial score equation: S n (β) = {Z i Z(u; β)}dn i (u) where ZiYi(u) exp(βzi) Z(u; β) = Yi(u) exp(βzi) (expected value of Z i in the risk set at u). The partial likelihood estimator ˆβ is derived by solving S n (β) =. ) δi
25 A representation of S n (β) S n(t; β) = = + = + = n t {Z i Z(u; β)}dn i (u) n t {Z i Z(u; β)}[dn i (u) Y i (u) exp(βz i )λ (u)du] n t {Z i Z(u; β)}y i (u) exp(βz i )λ (u)du n t {Z i Z(u; β)}dm i (u; β) { t n } n Z i Y i (u) exp(βz i ) Y i (u) exp(βz i ) Z(u; β) λ (u)du n t {Z i Z(u; β)}dm i (u; β) ( )
26 Let U n (t) = n 1/2 S n (t; β). Then under proper regularity conditions < U n, U n > (t) = 1 n n t t {Z i Z(u)} 2 Y i (u) exp(βz i )λ (u)du E[{Z µ Z (u)} 2 Y i (u) exp(βz i )]λ (u)du Properties of S n (t; β) 1. E[S n (t; β)] = 2. using representation (*) to derive independent increment property: for s t; cov[s n (s; β){s n (t; β) S n (s; β)}] = 3. var[n 1/2 S n (t; β)] =< U n, U n > (t)
27 The partial likelihood estimator ˆβ n is consistent (proof skipped) Asymptotic normality partial score equation S n (t; β) = t {Z i Z(u; β)}dn i (u) Taylor expansion (mean value expansion): for β n lying between ˆβ n and β, Thus, = S n (t; ˆβ n ) = S n (t; β) + S n(β) β=β β n ( ˆβ n β) { n 1/2 ( ˆβ n β) = n 1 S } 1 n(β) β=β β n n 1/2 S n (β) n 1/2 S n (β) d N(, σ 2 )
28 slope: S n (β) β = since β n β, n j=1 {Z j Z(u; β)}y j (u) exp(βz j ) j=1 Y dn i (u) j(u) exp(βz j ) n 1 S n(β) β=β β n E[{Z j Z(u; β)}y j (u) exp(βz j )] dp (X j u, j = 1) = σ 2 E[Y j (u) exp(βz j )] thus, by likelihood theory, n( ˆβ n β) d N(, σ 2 )
29 therefore, standardized version [ n ˆβ n β ] 1/2 D N (, 1) j=1 {Zj Z(u)} 2 exp( ˆβ nz j)y j(u) j=1 exp( ˆβ dn i (u) nz j)y j(u) estimation of σ 2 : n 1 n i j=1 {Z j Z(X i ; ˆβ n )}Y j (X i ) exp( ˆβ n Z j ) j=1 Y j(x i ) exp( ˆβ n Z j )
30 Denote S Z (u; β) = 95% confidence interval for β is ˆβ n ± 1.96 j=1 {Zj Z(u)} 2 exp(βz j)y j(u) j=1 exp(βzj)yj(u) [ n S Z (u; ˆβ n )dn i (u) ] 1/2
31 Hypothesis testing on H : β = β : rejects H at 5% type-i error Wald s test: ˆβ n β [ n S Z (u; ˆβ n )dn i (u) The partial likelihood score test: ] 1/2 n 1 2 S n (β ) var[s n (β )] i.e., {Z i Z(u; β )}dn i (u) [ S Z (u; β )dn i (u) ] 1/2 1.96
32 what if β = and Z = /1 in the partial likelihood score test: {Z i Z(u; β )}dn i (u) [ S Z (u; β )dn i (u) ] 1/2 where S Z (u; β ) = j=1 {Z j Z(u)} 2 exp(β Z)Y j (u) j=1 exp(β Z)Y j (u) = Z(u){1 Z(u)} and Z(u) is the fraction of treated subjects in the risk set at u. Note: The partial likelihood score test is the Log-rank test with a more general variance estimator. This variance estimator estimator coincides with that of the Log-rank test when failure time T is continuous.
33 Proportional hazards model for multivariate covariates Z R p λ(t Z) = λ (t)exp(β T Z) β = (β 1,..., β p ) T is p-vector partial score function where Z p 1 i = S n (β) = Z i1 Z i2... Z ip n p 1 {Z p 1 i Z p 1 (u, β)}dn i (u), Zp 1 (u, β) = Z 1 (u, β) Z 2 (u, β)... Z p (u, β) p 1
34 n 1/2 S n (β) D N (, Σ(β)), where Σ(β) is the limit of n 1 = n 1 n j=1 {Z j Z} p 1 {Z j Z} T 1 pexp(βz j )Y i (u) j=1 exp(βz dn i (u) j)y i (u) where v = 1, v 1 = v, v 2 = vv T. partial derivative Sn(β n ) β Σ(β) p p n 1/2 ( ˆβ n β) D N (, Σ 1 (β)) j=1 {Z 2 j Z} p pexp(βz j )Y i (u) j=1 exp(βz dn i (u) j)y i (u)
35 Estimation of baseline hazard function (Breslow estimator) ˆΛ (t; ˆβ n ) = n t dn i (u) t Y i(u)exp( ˆβ n T Z i ) = dn i(u) Y i(u)exp( ˆβ n T Z i ) consistent asymptotic normality: n 1/2 [ˆΛ (t; ˆβ n ) Λ (t)] equals n 1/2 [ˆΛ (t; β) Λ (t)] + n 1/2 [ˆΛ (t; ˆβ n ) ˆΛ (t; β)] = Term I + Term II
36 Term I: n 1/2 n t dm i (t; β) exp(βt Z i )Y i (u), which is a sum of martingale residuals with predictable variation: n n t Y i (u)exp(β T Z i )λ (u)du t [ exp(βt Z i )Y i (u)] λ (u)du 2 E[exp(β T Z)Y (u)] = σ2 1(t)
37 Term II: mean value expansion t = n 1/2 ( ˆβ n β) t n 1/2 [ˆΛ (t; ˆβ n ) ˆΛ (t; β)] exp(β T n Z i )Z i Y i (u)dn(u) [ exp(β T n Z i )Y i (u)] 2 exp(β T n Z i )Z i Y i (u)dn(u) [ µ(t) exp(β T n Z i )Y i (u)] 2
38 n 1/2 ( ˆβ n β) = n 1/2 {Z i Z}dM i (u; β) S 2 Z (u; β)dn(u) D N (, σ 2 ) covariance between Term I and Term II n t dm i (t; β) exp(βt Z i )Y i (u), which is a predictable covariation n t n {Z i Z(u)}dM i (u; β), {Z i Z(u)}exp(β T Z i )Y i (u)λ (u)du exp(βt Z i )Y i (u) =
39 Term I and II are uncorrelated Thus, n 1/2 [ˆΛ (t; ˆβ n ) Λ (t)] N (, σ 2 1(t) + µ(t)2 σ 2 )
40 Appendix Consistency of K-M estimator We consider the general case that the failure time T could be partly continuous and partly discrete. Recall S(t) = P {T t}. Define S(t) = 1 S(t) = P {T < t} and Define the subsurvival functions S (t) = S Y (t) = P {Y t} = [S(t)][S C (t)] Then, SU (t) = P {Y t, δ = 1} = SC(t) = P {Y t, δ = } = t S (t) = S U (t) + S C(t). t [S C (u)]d S(u) [S(u)]d S C (u)
41 We will show that S(t) can be expressed as a function of SU (t) and S C (t). (i) Suppose SU (t) is continuous at u < t. t ds U (u) S U (u) + S C (u) = = t t = log S(t). S C (u)df (u) S(u)S C (u) ds(u) S(u) = log[s(u)] t
42 Therefore, [ t dsu S(t) = exp (u) ] SU (u) + S C (u). (ii) Suppose SU has a jump at t, but S C is continuous at t. log S U (t+) + S C (t+) S U (t ) + S C (t ) = log [S(t+)][S C(t+)] [S(t )][S C (t )], = log [S(t+)] [S(t )] (The second equality follows from the fact that S c is continuous at t so S C (t+) = S C (t ).)
43 Therefore, S(t+) S(t) = S(t+) { S(t ) = exp log [ S U (t+) + SC (t+) ]} SU (t ) + S C (t ). (The first equality is due to the left continuity of survival function.) If the underlying distributions S and S C have no common jumps, then from (i) and (ii) t S(t) = exp{c ds U (u) S U (u) + S C (u) + d u<t log [ S U (u+) + SC (u+) ] SU (u ) + S C (u ) }, where c denotes integration over the continuity intervals of Su and d denotes summation over the discrete jumps of Su. The above expression is called Peterson s representation, showing that S(t) can be represented as a function of SU, S C, and t, that is, S(t) = ψ(s U, S C; t).
44 Peterson s representation gives us a proof that the KM-estimator Ŝ(t) is consistent. The proof proceeds as follows. Define the empirical sub-survival functions Ŝ U (t) = 1 n Ŝ C(t) = 1 n n I(y i > t, δ i = 1), n I(y i > t, δ i = ). It can be seen that the PL estimator is Ŝ(t) = ψ(ŝ U, Ŝ C; t), provided any ties between uncensored and censored observations are interpreted as uncensored observations preceding censored.
45 Notice that since Ŝ U is discrete, ψ(ŝ U, Ŝ C ; t) involves only the summation over the discrete jumps of Ŝ U. By the Glivenko-Cantelli theorem, Ŝ U (t) a.s. S U (t), Ŝ C(t) a.s. S C(t), uniformly in t. (The notation a.s. denotes converges almost surely to. ) Also, ψ is a continuous function of SU, S C in the sup norm. That is, if SU S U = sup t SU (t) S U (t) and S C S C, then ψ(su, S C ; t) ψ(s U, S c ; t). Therefore, Ŝ(t) = ψ(ŝ U, Ŝ C; t) a.s. ψ(s U, S C; t) = S(t).
46 REFERENCE Peterson, JASA (1977).
47 Functional Delta method Suppose that n 1/2 (ˆΛ( ) Λ( )) d U( ). If φ( ) is compact differentiable, then that is, n 1/2 [φ(ˆλ( )) φ(λ( ))]) d dφ(λ( ); U( )) n 1/2 [φ(λ n ( )) φ(λ( ))]) d dφ(λ( ); n 1/2 (ˆΛ(t) Λ(t))) Note: dφ(f ; G) means differentiability at F in the direction of G
STAT331. Combining Martingales, Stochastic Integrals, and Applications to Logrank Test & Cox s Model
STAT331 Combining Martingales, Stochastic Integrals, and Applications to Logrank Test & Cox s Model Because of Theorem 2.5.1 in Fleming and Harrington, see Unit 11: For counting process martingales with
More informationSurvival Analysis: Counting Process and Martingale. Lu Tian and Richard Olshen Stanford University
Survival Analysis: Counting Process and Martingale Lu Tian and Richard Olshen Stanford University 1 Lebesgue-Stieltjes Integrals G( ) is a right-continuous step function having jumps at x 1, x 2,.. b f(x)dg(x)
More informationSTAT331 Lebesgue-Stieltjes Integrals, Martingales, Counting Processes
STAT331 Lebesgue-Stieltjes Integrals, Martingales, Counting Processes This section introduces Lebesgue-Stieltjes integrals, and defines two important stochastic processes: a martingale process and a counting
More informationLinear rank statistics
Linear rank statistics Comparison of two groups. Consider the failure time T ij of j-th subject in the i-th group for i = 1 or ; the first group is often called control, and the second treatment. Let n
More informationSTAT Sample Problem: General Asymptotic Results
STAT331 1-Sample Problem: General Asymptotic Results In this unit we will consider the 1-sample problem and prove the consistency and asymptotic normality of the Nelson-Aalen estimator of the cumulative
More informationExercises. (a) Prove that m(t) =
Exercises 1. Lack of memory. Verify that the exponential distribution has the lack of memory property, that is, if T is exponentially distributed with parameter λ > then so is T t given that T > t for
More informationSTAT 331. Martingale Central Limit Theorem and Related Results
STAT 331 Martingale Central Limit Theorem and Related Results In this unit we discuss a version of the martingale central limit theorem, which states that under certain conditions, a sum of orthogonal
More informationUNIVERSITY OF CALIFORNIA, SAN DIEGO
UNIVERSITY OF CALIFORNIA, SAN DIEGO Estimation of the primary hazard ratio in the presence of a secondary covariate with non-proportional hazards An undergraduate honors thesis submitted to the Department
More informationDAGStat Event History Analysis.
DAGStat 2016 Event History Analysis Robin.Henderson@ncl.ac.uk 1 / 75 Schedule 9.00 Introduction 10.30 Break 11.00 Regression Models, Frailty and Multivariate Survival 12.30 Lunch 13.30 Time-Variation and
More informationStatistics 262: Intermediate Biostatistics Non-parametric Survival Analysis
Statistics 262: Intermediate Biostatistics Non-parametric Survival Analysis Jonathan Taylor & Kristin Cobb Statistics 262: Intermediate Biostatistics p.1/?? Overview of today s class Kaplan-Meier Curve
More informationOther Survival Models. (1) Non-PH models. We briefly discussed the non-proportional hazards (non-ph) model
Other Survival Models (1) Non-PH models We briefly discussed the non-proportional hazards (non-ph) model λ(t Z) = λ 0 (t) exp{β(t) Z}, where β(t) can be estimated by: piecewise constants (recall how);
More informationUniversity of California, Berkeley
University of California, Berkeley U.C. Berkeley Division of Biostatistics Working Paper Series Year 24 Paper 153 A Note on Empirical Likelihood Inference of Residual Life Regression Ying Qing Chen Yichuan
More informationPhD course in Advanced survival analysis. One-sample tests. Properties. Idea: (ABGK, sect. V.1.1) Counting process N(t)
PhD course in Advanced survival analysis. (ABGK, sect. V.1.1) One-sample tests. Counting process N(t) Non-parametric hypothesis tests. Parametric models. Intensity process λ(t) = α(t)y (t) satisfying Aalen
More informationSTAT331. Cox s Proportional Hazards Model
STAT331 Cox s Proportional Hazards Model In this unit we introduce Cox s proportional hazards (Cox s PH) model, give a heuristic development of the partial likelihood function, and discuss adaptations
More informationAsymptotic Distributions for the Nelson-Aalen and Kaplan-Meier estimators and for test statistics.
Asymptotic Distributions for the Nelson-Aalen and Kaplan-Meier estimators and for test statistics. Dragi Anevski Mathematical Sciences und University November 25, 21 1 Asymptotic distributions for statistical
More informationLecture 5 Models and methods for recurrent event data
Lecture 5 Models and methods for recurrent event data Recurrent and multiple events are commonly encountered in longitudinal studies. In this chapter we consider ordered recurrent and multiple events.
More informationPart III. Hypothesis Testing. III.1. Log-rank Test for Right-censored Failure Time Data
1 Part III. Hypothesis Testing III.1. Log-rank Test for Right-censored Failure Time Data Consider a survival study consisting of n independent subjects from p different populations with survival functions
More informationSurvival Analysis Math 434 Fall 2011
Survival Analysis Math 434 Fall 2011 Part IV: Chap. 8,9.2,9.3,11: Semiparametric Proportional Hazards Regression Jimin Ding Math Dept. www.math.wustl.edu/ jmding/math434/fall09/index.html Basic Model Setup
More informationLecture 22 Survival Analysis: An Introduction
University of Illinois Department of Economics Spring 2017 Econ 574 Roger Koenker Lecture 22 Survival Analysis: An Introduction There is considerable interest among economists in models of durations, which
More informationEfficiency of Profile/Partial Likelihood in the Cox Model
Efficiency of Profile/Partial Likelihood in the Cox Model Yuichi Hirose School of Mathematics, Statistics and Operations Research, Victoria University of Wellington, New Zealand Summary. This paper shows
More informationLecture 6 PREDICTING SURVIVAL UNDER THE PH MODEL
Lecture 6 PREDICTING SURVIVAL UNDER THE PH MODEL The Cox PH model: λ(t Z) = λ 0 (t) exp(β Z). How do we estimate the survival probability, S z (t) = S(t Z) = P (T > t Z), for an individual with covariates
More information9 Estimating the Underlying Survival Distribution for a
9 Estimating the Underlying Survival Distribution for a Proportional Hazards Model So far the focus has been on the regression parameters in the proportional hazards model. These parameters describe the
More information1. Stochastic Processes and filtrations
1. Stochastic Processes and 1. Stoch. pr., A stochastic process (X t ) t T is a collection of random variables on (Ω, F) with values in a measurable space (S, S), i.e., for all t, In our case X t : Ω S
More informationGoodness-of-fit test for the Cox Proportional Hazard Model
Goodness-of-fit test for the Cox Proportional Hazard Model Rui Cui rcui@eco.uc3m.es Department of Economics, UC3M Abstract In this paper, we develop new goodness-of-fit tests for the Cox proportional hazard
More information1 Glivenko-Cantelli type theorems
STA79 Lecture Spring Semester Glivenko-Cantelli type theorems Given i.i.d. observations X,..., X n with unknown distribution function F (t, consider the empirical (sample CDF ˆF n (t = I [Xi t]. n Then
More informationIntroduction to Empirical Processes and Semiparametric Inference Lecture 25: Semiparametric Models
Introduction to Empirical Processes and Semiparametric Inference Lecture 25: Semiparametric Models Michael R. Kosorok, Ph.D. Professor and Chair of Biostatistics Professor of Statistics and Operations
More informationFollow this and additional works at: Part of the Applied Mathematics Commons
Louisiana State University LSU Digital Commons LSU Master's Theses Graduate School 2004 Cox Regression Model Lindsay Sarah Smith Louisiana State University and Agricultural and Mechanical College, lsmit51@lsu.edu
More informationCox s proportional hazards model and Cox s partial likelihood
Cox s proportional hazards model and Cox s partial likelihood Rasmus Waagepetersen October 12, 2018 1 / 27 Non-parametric vs. parametric Suppose we want to estimate unknown function, e.g. survival function.
More informationMAS3301 / MAS8311 Biostatistics Part II: Survival
MAS3301 / MAS8311 Biostatistics Part II: Survival M. Farrow School of Mathematics and Statistics Newcastle University Semester 2, 2009-10 1 13 The Cox proportional hazards model 13.1 Introduction In the
More informationTMA 4275 Lifetime Analysis June 2004 Solution
TMA 4275 Lifetime Analysis June 2004 Solution Problem 1 a) Observation of the outcome is censored, if the time of the outcome is not known exactly and only the last time when it was observed being intact,
More informationHypothesis Testing Based on the Maximum of Two Statistics from Weighted and Unweighted Estimating Equations
Hypothesis Testing Based on the Maximum of Two Statistics from Weighted and Unweighted Estimating Equations Takeshi Emura and Hisayuki Tsukuma Abstract For testing the regression parameter in multivariate
More informationLecture 3. Truncation, length-bias and prevalence sampling
Lecture 3. Truncation, length-bias and prevalence sampling 3.1 Prevalent sampling Statistical techniques for truncated data have been integrated into survival analysis in last two decades. Truncation in
More informationTentamentsskrivning: Stochastic processes 1
Tentamentsskrivning: Stochastic processes 1 Tentamentsskrivning i MSF2/MVE33, 7.5 hp. Tid: fredagen den 31 maj 213 kl 8.3-12.3 Examinator och jour: Serik Sagitov, tel. 772-5351, mob. 736 97 613, rum H326
More informationProduct-limit estimators of the survival function with left or right censored data
Product-limit estimators of the survival function with left or right censored data 1 CREST-ENSAI Campus de Ker-Lann Rue Blaise Pascal - BP 37203 35172 Bruz cedex, France (e-mail: patilea@ensai.fr) 2 Institut
More informationBrownian motion. Samy Tindel. Purdue University. Probability Theory 2 - MA 539
Brownian motion Samy Tindel Purdue University Probability Theory 2 - MA 539 Mostly taken from Brownian Motion and Stochastic Calculus by I. Karatzas and S. Shreve Samy T. Brownian motion Probability Theory
More informationFULL LIKELIHOOD INFERENCES IN THE COX MODEL
October 20, 2007 FULL LIKELIHOOD INFERENCES IN THE COX MODEL BY JIAN-JIAN REN 1 AND MAI ZHOU 2 University of Central Florida and University of Kentucky Abstract We use the empirical likelihood approach
More informationA Comparison of Different Approaches to Nonparametric Inference for Subdistributions
A Comparison of Different Approaches to Nonparametric Inference for Subdistributions Johannes Mertsching Johannes.Mertsching@gmail.com Master Thesis Supervision: dr. Ronald B. Geskus prof. Chris A.J. Klaassen
More informationEstimation and Inference of Quantile Regression. for Survival Data under Biased Sampling
Estimation and Inference of Quantile Regression for Survival Data under Biased Sampling Supplementary Materials: Proofs of the Main Results S1 Verification of the weight function v i (t) for the lengthbiased
More informationSurvival Analysis for Case-Cohort Studies
Survival Analysis for ase-ohort Studies Petr Klášterecký Dept. of Probability and Mathematical Statistics, Faculty of Mathematics and Physics, harles University, Prague, zech Republic e-mail: petr.klasterecky@matfyz.cz
More informationMartingale Theory for Finance
Martingale Theory for Finance Tusheng Zhang October 27, 2015 1 Introduction 2 Probability spaces and σ-fields 3 Integration with respect to a probability measure. 4 Conditional expectation. 5 Martingales.
More informationSurvival Analysis: Weeks 2-3. Lu Tian and Richard Olshen Stanford University
Survival Analysis: Weeks 2-3 Lu Tian and Richard Olshen Stanford University 2 Kaplan-Meier(KM) Estimator Nonparametric estimation of the survival function S(t) = pr(t > t) The nonparametric estimation
More informationYou know I m not goin diss you on the internet Cause my mama taught me better than that I m a survivor (What?) I m not goin give up (What?
You know I m not goin diss you on the internet Cause my mama taught me better than that I m a survivor (What?) I m not goin give up (What?) I m not goin stop (What?) I m goin work harder (What?) Sir David
More informationSurvival Analysis. Stat 526. April 13, 2018
Survival Analysis Stat 526 April 13, 2018 1 Functions of Survival Time Let T be the survival time for a subject Then P [T < 0] = 0 and T is a continuous random variable The Survival function is defined
More informationMATH 6605: SUMMARY LECTURE NOTES
MATH 6605: SUMMARY LECTURE NOTES These notes summarize the lectures on weak convergence of stochastic processes. If you see any typos, please let me know. 1. Construction of Stochastic rocesses A stochastic
More information1 Introduction. 2 Residuals in PH model
Supplementary Material for Diagnostic Plotting Methods for Proportional Hazards Models With Time-dependent Covariates or Time-varying Regression Coefficients BY QIQING YU, JUNYI DONG Department of Mathematical
More informationAttributable Risk Function in the Proportional Hazards Model
UW Biostatistics Working Paper Series 5-31-2005 Attributable Risk Function in the Proportional Hazards Model Ying Qing Chen Fred Hutchinson Cancer Research Center, yqchen@u.washington.edu Chengcheng Hu
More informationEMPIRICAL LIKELIHOOD AND DIFFERENTIABLE FUNCTIONALS
University of Kentucky UKnowledge Theses and Dissertations--Statistics Statistics 2016 EMPIRICAL LIKELIHOOD AND DIFFERENTIABLE FUNCTIONALS Zhiyuan Shen University of Kentucky, alanshenpku10@gmail.com Digital
More informationCox s proportional hazards/regression model - model assessment
Cox s proportional hazards/regression model - model assessment Rasmus Waagepetersen September 27, 2017 Topics: Plots based on estimated cumulative hazards Cox-Snell residuals: overall check of fit Martingale
More informationSTAT 331. Accelerated Failure Time Models. Previously, we have focused on multiplicative intensity models, where
STAT 331 Accelerated Failure Time Models Previously, we have focused on multiplicative intensity models, where h t z) = h 0 t) g z). These can also be expressed as H t z) = H 0 t) g z) or S t z) = e Ht
More informationlog T = β T Z + ɛ Zi Z(u; β) } dn i (ue βzi ) = 0,
Accelerated failure time model: log T = β T Z + ɛ β estimation: solve where S n ( β) = n i=1 { Zi Z(u; β) } dn i (ue βzi ) = 0, Z(u; β) = j Z j Y j (ue βz j) j Y j (ue βz j) How do we show the asymptotics
More information(A n + B n + 1) A n + B n
344 Problem Hints and Solutions Solution for Problem 2.10. To calculate E(M n+1 F n ), first note that M n+1 is equal to (A n +1)/(A n +B n +1) with probability M n = A n /(A n +B n ) and M n+1 equals
More informationStatistical Methods for Alzheimer s Disease Studies
Statistical Methods for Alzheimer s Disease Studies Rebecca A. Betensky, Ph.D. Department of Biostatistics, Harvard T.H. Chan School of Public Health July 19, 2016 1/37 OUTLINE 1 Statistical collaborations
More informationDefinitions and examples Simple estimation and testing Regression models Goodness of fit for the Cox model. Recap of Part 1. Per Kragh Andersen
Recap of Part 1 Per Kragh Andersen Section of Biostatistics, University of Copenhagen DSBS Course Survival Analysis in Clinical Trials January 2018 1 / 65 Overview Definitions and examples Simple estimation
More information11 Survival Analysis and Empirical Likelihood
11 Survival Analysis and Empirical Likelihood The first paper of empirical likelihood is actually about confidence intervals with the Kaplan-Meier estimator (Thomas and Grunkmeier 1979), i.e. deals with
More informationStatistical Inference and Methods
Department of Mathematics Imperial College London d.stephens@imperial.ac.uk http://stats.ma.ic.ac.uk/ das01/ 31st January 2006 Part VI Session 6: Filtering and Time to Event Data Session 6: Filtering and
More informationFiltrations, Markov Processes and Martingales. Lectures on Lévy Processes and Stochastic Calculus, Braunschweig, Lecture 3: The Lévy-Itô Decomposition
Filtrations, Markov Processes and Martingales Lectures on Lévy Processes and Stochastic Calculus, Braunschweig, Lecture 3: The Lévy-Itô Decomposition David pplebaum Probability and Statistics Department,
More informationSurvival Analysis. STAT 526 Professor Olga Vitek
Survival Analysis STAT 526 Professor Olga Vitek May 4, 2011 9 Survival Data and Survival Functions Statistical analysis of time-to-event data Lifetime of machines and/or parts (called failure time analysis
More information6. Brownian Motion. Q(A) = P [ ω : x(, ω) A )
6. Brownian Motion. stochastic process can be thought of in one of many equivalent ways. We can begin with an underlying probability space (Ω, Σ, P) and a real valued stochastic process can be defined
More informationAsymptotic Properties of Kaplan-Meier Estimator. for Censored Dependent Data. Zongwu Cai. Department of Mathematics
To appear in Statist. Probab. Letters, 997 Asymptotic Properties of Kaplan-Meier Estimator for Censored Dependent Data by Zongwu Cai Department of Mathematics Southwest Missouri State University Springeld,
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More information9 Brownian Motion: Construction
9 Brownian Motion: Construction 9.1 Definition and Heuristics The central limit theorem states that the standard Gaussian distribution arises as the weak limit of the rescaled partial sums S n / p n of
More informationEmpirical Processes & Survival Analysis. The Functional Delta Method
STAT/BMI 741 University of Wisconsin-Madison Empirical Processes & Survival Analysis Lecture 3 The Functional Delta Method Lu Mao lmao@biostat.wisc.edu 3-1 Objectives By the end of this lecture, you will
More informationSTAT 200C: High-dimensional Statistics
STAT 200C: High-dimensional Statistics Arash A. Amini May 30, 2018 1 / 59 Classical case: n d. Asymptotic assumption: d is fixed and n. Basic tools: LLN and CLT. High-dimensional setting: n d, e.g. n/d
More informationUniversal examples. Chapter The Bernoulli process
Chapter 1 Universal examples 1.1 The Bernoulli process First description: Bernoulli random variables Y i for i = 1, 2, 3,... independent with P [Y i = 1] = p and P [Y i = ] = 1 p. Second description: Binomial
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 15. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationThe instantaneous and forward default intensity of structural models
The instantaneous and forward default intensity of structural models Cho-Jieh Chen Department of Mathematical and Statistical Sciences University of Alberta Edmonton, AB E-mail: cjchen@math.ualberta.ca
More informationIn contrast, parametric techniques (fitting exponential or Weibull, for example) are more focussed, can handle general covariates, but require
Chapter 5 modelling Semi parametric We have considered parametric and nonparametric techniques for comparing survival distributions between different treatment groups. Nonparametric techniques, such as
More informationMultivariate Survival Data With Censoring.
1 Multivariate Survival Data With Censoring. Shulamith Gross and Catherine Huber-Carol Baruch College of the City University of New York, Dept of Statistics and CIS, Box 11-220, 1 Baruch way, 10010 NY.
More informationECE534, Spring 2018: Solutions for Problem Set #4 Due Friday April 6, 2018
ECE534, Spring 2018: s for Problem Set #4 Due Friday April 6, 2018 1. MMSE Estimation, Data Processing and Innovations The random variables X, Y, Z on a common probability space (Ω, F, P ) are said to
More informationSolution for Problem 7.1. We argue by contradiction. If the limit were not infinite, then since τ M (ω) is nondecreasing we would have
362 Problem Hints and Solutions sup g n (ω, t) g(ω, t) sup g(ω, s) g(ω, t) µ n (ω). t T s,t: s t 1/n By the uniform continuity of t g(ω, t) on [, T], one has for each ω that µ n (ω) as n. Two applications
More informationContinuous case Discrete case General case. Hazard functions. Patrick Breheny. August 27. Patrick Breheny Survival Data Analysis (BIOS 7210) 1/21
Hazard functions Patrick Breheny August 27 Patrick Breheny Survival Data Analysis (BIOS 7210) 1/21 Introduction Continuous case Let T be a nonnegative random variable representing the time to an event
More informationPOWER AND SAMPLE SIZE DETERMINATIONS IN DYNAMIC RISK PREDICTION. by Zhaowen Sun M.S., University of Pittsburgh, 2012
POWER AND SAMPLE SIZE DETERMINATIONS IN DYNAMIC RISK PREDICTION by Zhaowen Sun M.S., University of Pittsburgh, 2012 B.S.N., Wuhan University, China, 2010 Submitted to the Graduate Faculty of the Graduate
More informationSurvival Analysis I (CHL5209H)
Survival Analysis Dalla Lana School of Public Health University of Toronto olli.saarela@utoronto.ca January 7, 2015 31-1 Literature Clayton D & Hills M (1993): Statistical Models in Epidemiology. Not really
More informationST745: Survival Analysis: Nonparametric methods
ST745: Survival Analysis: Nonparametric methods Eric B. Laber Department of Statistics, North Carolina State University February 5, 2015 The KM estimator is used ubiquitously in medical studies to estimate
More informationApplications of Ito s Formula
CHAPTER 4 Applications of Ito s Formula In this chapter, we discuss several basic theorems in stochastic analysis. Their proofs are good examples of applications of Itô s formula. 1. Lévy s martingale
More informationCIMPA SCHOOL, 2007 Jump Processes and Applications to Finance Monique Jeanblanc
CIMPA SCHOOL, 27 Jump Processes and Applications to Finance Monique Jeanblanc 1 Jump Processes I. Poisson Processes II. Lévy Processes III. Jump-Diffusion Processes IV. Point Processes 2 I. Poisson Processes
More informationExercise Exercise Homework #6 Solutions Thursday 6 April 2006
Unless otherwise stated, for the remainder of the solutions, define F m = σy 0,..., Y m We will show EY m = EY 0 using induction. m = 0 is obviously true. For base case m = : EY = EEY Y 0 = EY 0. Now assume
More information9. Estimating Survival Distribution for a PH Model
9. Estimating Survival Distribution for a PH Model Objective: Another Goal of the COX model Estimating the survival distribution for individuals with a certain combination of covariates. PH model assumption:
More informationChapter 2 Inference on Mean Residual Life-Overview
Chapter 2 Inference on Mean Residual Life-Overview Statistical inference based on the remaining lifetimes would be intuitively more appealing than the popular hazard function defined as the risk of immediate
More informationStatistical Analysis of Competing Risks With Missing Causes of Failure
Proceedings 59th ISI World Statistics Congress, 25-3 August 213, Hong Kong (Session STS9) p.1223 Statistical Analysis of Competing Risks With Missing Causes of Failure Isha Dewan 1,3 and Uttara V. Naik-Nimbalkar
More informationA Very Brief Summary of Statistical Inference, and Examples
A Very Brief Summary of Statistical Inference, and Examples Trinity Term 2009 Prof. Gesine Reinert Our standard situation is that we have data x = x 1, x 2,..., x n, which we view as realisations of random
More informationTests of independence for censored bivariate failure time data
Tests of independence for censored bivariate failure time data Abstract Bivariate failure time data is widely used in survival analysis, for example, in twins study. This article presents a class of χ
More informationPoint Process Control
Point Process Control The following note is based on Chapters I, II and VII in Brémaud s book Point Processes and Queues (1981). 1 Basic Definitions Consider some probability space (Ω, F, P). A real-valued
More informationTyler Hofmeister. University of Calgary Mathematical and Computational Finance Laboratory
JUMP PROCESSES GENERALIZING STOCHASTIC INTEGRALS WITH JUMPS Tyler Hofmeister University of Calgary Mathematical and Computational Finance Laboratory Overview 1. General Method 2. Poisson Processes 3. Diffusion
More informationEstimation of the Bivariate and Marginal Distributions with Censored Data
Estimation of the Bivariate and Marginal Distributions with Censored Data Michael Akritas and Ingrid Van Keilegom Penn State University and Eindhoven University of Technology May 22, 2 Abstract Two new
More informationSTAT 6350 Analysis of Lifetime Data. Failure-time Regression Analysis
STAT 6350 Analysis of Lifetime Data Failure-time Regression Analysis Explanatory Variables for Failure Times Usually explanatory variables explain/predict why some units fail quickly and some units survive
More informationQuick Review on Linear Multiple Regression
Quick Review on Linear Multiple Regression Mei-Yuan Chen Department of Finance National Chung Hsing University March 6, 2007 Introduction for Conditional Mean Modeling Suppose random variables Y, X 1,
More informationn E(X t T n = lim X s Tn = X s
Stochastic Calculus Example sheet - Lent 15 Michael Tehranchi Problem 1. Let X be a local martingale. Prove that X is a uniformly integrable martingale if and only X is of class D. Solution 1. If If direction:
More informationResiduals and model diagnostics
Residuals and model diagnostics Patrick Breheny November 10 Patrick Breheny Survival Data Analysis (BIOS 7210) 1/42 Introduction Residuals Many assumptions go into regression models, and the Cox proportional
More informationLikelihood Construction, Inference for Parametric Survival Distributions
Week 1 Likelihood Construction, Inference for Parametric Survival Distributions In this section we obtain the likelihood function for noninformatively rightcensored survival data and indicate how to make
More informationLecture 7. Proportional Hazards Model - Handling Ties and Survival Estimation Statistics Survival Analysis. Presented February 4, 2016
Proportional Hazards Model - Handling Ties and Survival Estimation Statistics 255 - Survival Analysis Presented February 4, 2016 likelihood - Discrete Dan Gillen Department of Statistics University of
More informationREGRESSION ANALYSIS FOR TIME-TO-EVENT DATA THE PROPORTIONAL HAZARDS (COX) MODEL ST520
REGRESSION ANALYSIS FOR TIME-TO-EVENT DATA THE PROPORTIONAL HAZARDS (COX) MODEL ST520 Department of Statistics North Carolina State University Presented by: Butch Tsiatis, Department of Statistics, NCSU
More informationEconomics 508 Lecture 22 Duration Models
University of Illinois Fall 2012 Department of Economics Roger Koenker Economics 508 Lecture 22 Duration Models There is considerable interest, especially among labor-economists in models of duration.
More informationBRANCHING PROCESSES AND THEIR APPLICATIONS: Lecture 15: Crump-Mode-Jagers processes and queueing systems with processor sharing
BRANCHING PROCESSES AND THEIR APPLICATIONS: Lecture 5: Crump-Mode-Jagers processes and queueing systems with processor sharing June 7, 5 Crump-Mode-Jagers process counted by random characteristics We give
More information4. Comparison of Two (K) Samples
4. Comparison of Two (K) Samples K=2 Problem: compare the survival distributions between two groups. E: comparing treatments on patients with a particular disease. Z: Treatment indicator, i.e. Z = 1 for
More informationModule 9: Stationary Processes
Module 9: Stationary Processes Lecture 1 Stationary Processes 1 Introduction A stationary process is a stochastic process whose joint probability distribution does not change when shifted in time or space.
More informationA Concise Course on Stochastic Partial Differential Equations
A Concise Course on Stochastic Partial Differential Equations Michael Röckner Reference: C. Prevot, M. Röckner: Springer LN in Math. 1905, Berlin (2007) And see the references therein for the original
More informationProblem Sheet 1. You may assume that both F and F are σ-fields. (a) Show that F F is not a σ-field. (b) Let X : Ω R be defined by 1 if n = 1
Problem Sheet 1 1. Let Ω = {1, 2, 3}. Let F = {, {1}, {2, 3}, {1, 2, 3}}, F = {, {2}, {1, 3}, {1, 2, 3}}. You may assume that both F and F are σ-fields. (a) Show that F F is not a σ-field. (b) Let X :
More informationAnalysis of Time-to-Event Data: Chapter 2 - Nonparametric estimation of functions of survival time
Analysis of Time-to-Event Data: Chapter 2 - Nonparametric estimation of functions of survival time Steffen Unkel Department of Medical Statistics University Medical Center Göttingen, Germany Winter term
More informationTESTINGGOODNESSOFFITINTHECOX AALEN MODEL
ROBUST 24 c JČMF 24 TESTINGGOODNESSOFFITINTHECOX AALEN MODEL David Kraus Keywords: Counting process, Cox Aalen model, goodness-of-fit, martingale, residual, survival analysis. Abstract: The Cox Aalen regression
More information