Resampling methods for randomly censored survival data
|
|
- Joleen Johnston
- 5 years ago
- Views:
Transcription
1 Resampling methods for randomly censored survival data Markus Pauly Lehrstuhl für Mathematische Statistik und Wahrscheinlichkeitstheorie Mathematisches Institut Heinrich-Heine-Universität Düsseldorf Wien, 5 June 2014
2 Thanks to J. Beyersmann (Ulm) M. Brendel (Boehringer Ingelheim, Biberach) D. Dobler (Düsseldorf) A. Janssen (Düsseldorf) C.-D. Mayer (BioSS, Aberdeen) S. di Termini (Freiburg)
3 Outline 1 Introduction and Motivation 2 Classical tests 3 Weighted projection-type permutation tests 4 Simulations and Data Analysis 5 More complex Multi-State Models
4 Introduction and Motivation Typcial study in medicine: patients Sample sizes are only sufficient if efficiently analyzed! Problems in practice Incompleteness of observations; e.g. right-censored data Observed are 2-types of data Complete observations, so called survival times Incomplete observations; last observed survival time or drop out of study for other reasons
5 Introduction and Motivation Censoring Source: Janssen
6 Introduction and Motivation Typcial study in medicine: patients Sample sizes are only sufficient if efficiently analyzed! Problems in practice Incompleteness of observations; e.g. right-censored data Observed are 2-types of data Complete observations, so called survival times Incomplete observations; last observed survival time or drop out of study for other reasons Example ( kidney data ): Infection times in two groups of dialysis patients
7 Introduction and Motivation Different catheterization procedures: percutaneous and surgical placements Source: Klein/Moeschberger
8 Introduction and Motivation Infection times of dialysis patients
9 Introduction and Motivation Basic Survival Analysis 1-4 T 0 survival time S(t) = P(T > t) survival function with ν-density f Λ(t) = t 0 f (t)/s(t )dν(t) cumulative hazard function λ = f /S hazard rate Interpretation (for smooth f ) P(T (t, t + ɛ] T t) ɛλ(t).
10 Introduction and Motivation Basic Survival Analysis 2-4 Random Censoring by C G; C and T independent We only observe X := min(t, C) and = 1{T C} Estimator for S in this case: Kaplan-Meier estimator Ŝ Estimator for Λ: Nelson-Aalen estimator ˆΛ
11 Introduction and Motivation Basic Survival Analysis 3-4 Easiest case: Let (T i ) i be i.i.d. discrete r.v.s; indep. from (C i ) i (also i.i.d.) Define (e.g. failure rates in actuarial sciences) 1 Theorem: Λ(t) = Λ T (t) = r(t) = P(T 1 = t T 1 t) = P(T 1 = t) P(T 1 t) 0 s t r(s) and S(t) = S T (t) = 0 s t (1 r(s)) Kaplan-Meier and Nelson-Aalen estimators by plug-in: ˆr(s) = number of failures at time s number under risk at s 1 r(t) = 0 for P(T 1 t) = 0
12 Introduction and Motivation Basic Survival Analysis 4-4 In counting process notation: Define counting processes N(t) = n 1{X i t, i = 1}, Y (t) = i=1 n 1{X i t} i=1 N(t) = N(t) N(t ) Nelson-Aalen and Kaplan-Meier are given by ˆΛ(t) = 0 s t N(s) Y (s) Uncensored case: 1 Ŝ = e.d.f. and Ŝ(t) = 0 s t ( 1 N(S) Y (s) Properties (under reg): Consistent and asymptotic Gaussian )
13 Introduction and Motivation Plots of Kaplan Meier estimators for the two sample kidney data solid line = percutaneous
14 Introduction and Motivation Here: Randomly right censored data Patients that survive without infection are censored Standard model: Cox proportional hazard model 2 2 λ ϑ (t) = e ϑ λ 0 (t)
15 Introduction and Motivation Questions: Procedures significant different? Is one survival time (or hazard rate) stochastic greater? Answer of classical log-rank test (for proportional hazards) p-value > 0.05 Reason: time depending hazard ratios! Log rank test has problems detecting them!
16 Classical tests Classical two sample problem Observe independent random variables X i := min(t i, C i ) and i = 1{T i C i } for 1 i n within two groups i.i.d. i.i.d. T i F 1, 1 i n 1, T n1 +i F 2, 1 i n 2, (continuous) i.i.d. i.i.d. C i = G 1, 1 i n 1, C n1 +i G 2, 1 i n 2. (continuous) Null hypothesis H 0 : {T 1,..., T n i.i.d.} = {F 1 = F 2 } = {Λ 1 = Λ 2 } against 2-sided or 1-sided alternatives H 0 1 : {Λ 1 Λ 2 }, H 1 1 : {λ 1 λ 2 }, H 2 1 : {Λ 1 Λ 2 }. Unknown nuisance param: G 1, G 2 and under the null L(T 1 )
17 Classical tests Class of weighted logrank test statistics w n (t) := w( F n (t )), F n Kaplan Meier estimator of the pooled sample w : [0, 1] R weight function Test statistics (ABGK 1993) are Tn(wn) σ(w n) n T n (w n ) = n 1 n 2 0 w n (s) Y 1(s)Y 2 (s) Y (s) σ 2 (w n ) adequate variance estimator 3 w = 1 for classical logrank statistic w(u) = 1 u for Prentice-Wilcoxon statistic Tn(wn) or σ(w 1{T n) n(w n ) > 0} with { d ˆΛ 1n (s) d ˆΛ } 2n (s) 3 Gill (1980): σ 2 (w n) = n w 2 n 1 n 2 0 n Y 1Y 2 d ˆΛ Y n
18 Classical tests Remarks Intuitively: w determines power behaviour Example: Prentice-Wilcoxon with decreasing w(u) = 1 u More weight on early times! Should have good power against early hazard differences! In comparison: Classical logrank: Equal weight on all time points! Now: ARE comparison via reparametrization!
19 Classical tests Local parametrization (scores) Classical approach: classical likelihood based approach parametric path ϑ P ϑ in a nonparametric model P of distributions score functions d dϑ log dp ϑ dµ = g(ϑ) score test
20 Classical tests Local parametrization (hazards) approach based on hazards: F = 1 S Λ cumulative hazard function Now: hazard rates λ(u) = f (u) S(u) rather than densities f (u) parameters: relative risk = hazard rate ratio represented by γ : [0, 1] R under ass Rem: Holds for (ϑ 1) Λ ϑ (t) := t 0 1 lim ϑ 0 ϑ (dλ ϑ 1) = γ F 0. dλ 0 (1 + ϑγ F)dΛ 0, Λ 0 baseline hazard.
21 Classical tests Examples of semiparametric models γ 1 (u) = c γ 2 (u) = cu u 1 proport. hazards γ 3 (u) = c(1 u) u 1 late hazards γ 4 (u) = cu(1 u) u 1 early hazards u 1 central hazards
22 Classical tests Consider model given by γ F 0 Recall w n (t) := w( F n (t )), w : [0, 1] R weight function Janssen (1991) and Neuhaus (2000): ARE of Tn(wn) σ(w n) 1{T n(w n ) > 0} for local alternatives 4 in direction of γ F 0 given by ARE = w, γ 2 µ w 2 µ γ 2 µ = cos 2 (β) w, γ µ = 1 0 wγdµ, β angle between w and γ n µ measure on R + depending on G i, F 0 and lim 1 n n Bad ARE for some directions γ! 4 L(T 1,..., T n) = i Pc ni
23 Classical tests On the ARE ARE (ψ): asymptotic relative efficiency of a test ψ ARE(ψ) N opt N(ψ) (in the limit) N(ψ) = No. of needed obs. for ψ to achieve a given power N opt = Minimum no. of needed obs. Pitman s interpretation: 100(1 ARE)% of observations are wasted by using ψ Use of the full data information only for ARE = 1!
24 Classical tests Example: ARE(ψ)= 1 2 twice as many obs. are necessary for ψ Clearly: Optimal procedures depend on the model which is in general unknown! Idea: Estimate the model! (locally)
25 Classical tests Closer look to weighted logrank statistics Consider emp. subspace V = {βw n : β R} and cone V + 1 = {βw n : β [0, )}, generated by w n. Lemma (Brendel et al., 2014) We have where T n (w n ) σ(w n ) = Π V1 (ˆγ n ) ˆµn ; T n(w n ) σ(w n ) 1{T n(w n ) > 0} = Π V + 1 (ˆγ n ) ˆµn, ˆµ n emp. estimator of µ, ˆγ n = emp. estimator of γ F 0 (locally a ) Π V (+) 1 = L 2 (ˆµ n )-projection into V (+) 1 a n all estimates depend on NA (and Y ); e.g. ˆγ n 1 n 2 (d Λ n 1 d Λ 2 )/d Λ n
26 Classical tests Interpretation ˆγ n estimates the γ F 0 -model (locally) Logrank stat measures distance of its projection into space/cone generated by w n Idea: Use larger spaces/cones to cover larger classes of alternatives!
27 Weighted projection-type permutation tests Projection approach Idea: Scientist chooses relevant weights 5 w in ( ) := w i ( F n ( )) 1 i r, e.g. to discover differences of the relative risk for proportional hazards (constant over time) early survival times central survival times late survival times. These generate larger linear space/cone: { r } { r } V := β i w in, β i R, and V + := β i w in, β i 0, i=1 i=1 5 Fn Kaplan Meier estimator of the pooled sample
28 Weighted projection-type permutation tests Projection statistics As in the one-dimensional case consider test statistics S n0 := Π V (ˆγ n ) 2ˆµ n and S n1 := Π V +(ˆγ n ) 2ˆµ n Theorem (Brendel et al., 2014) Under regularity conditions we have convergence in distribution under the null Π V (ˆγ n ) 2ˆµ n d χ 2 rank(σ) and Π V +(ˆγ n) 2ˆµ d n F Σ where F Σ is continuous and Σ = ( w i, w j µ ) i,j r.
29 Weighted projection-type permutation tests Projection-type tests Estimate Σ consistently by Σ Results in tests Theorem (Brendel et al., 2014) Both tests are φ 0n = 1{ Π V (ˆγ n ) 2ˆµ n χ 2 rank( Σ),1 α > 0}, φ 1n = 1{ Π V +(ˆγ n ) 2ˆµ n F 1 (1 α) > 0}, Σ asymptotical level α (even for G 1 G 2 ) and consistent for fixed alternatives if at least one of their associated weighted logrank tests is consistent. Posses broader power functions!
30 Weighted projection-type permutation tests Projection-type permutation tests Further modification: Use permutation version of the tests! Advantage: Finitely exact under {F 1 = F 2, G 1 = G 2 }. For c n0 = χ 2 and c rank( Σ),1 α n1 = F 1 (1 α): with Σ φ nk := 1{S nk c nk > cnk } + k n 1{S nk c nk = cnk } k = 0, 1, c nk = cond (1 α)-quantile of the perm dist of S nk c nk Role of S nk = S nk c nk ( Σ): Studentized-type statistics!
31 Weighted projection-type permutation tests Properties Theorem (Brendel et al., 2014) The permutation tests are asymptotically equivalent to their corresponding projection tests, i.e. lim E H n 0 ( φ 0,n φ 0,n ) = 0, lim E H n 0 ( φ 1,n φ 1,n ) = 0. Implies same power for contigouos alternatives! More math. details (as asymptotic admissability) in Brendel, Janssen, Mayer & Pauly (2014, SJS).
32 Simulations and Data Analysis Set-up 1-2 Group I: T 1,..., T n1 i.i.d. F 0 (x) = (1 exp( x))1 (0, ) (x) For group II: 3 different scenarios; given by directions γ i = w i w 1 (u) = 1, w 2 (u) = 1 2u, w 3 (u) = u(1 u), 0 u 1, corresponding to proportional, crossing and central hazards Group II: T n1 +1,..., T n i.i.d. with Λ ϑ,i (t) = t ϑ w i (F 0 (x))dx
33 Simulations and Data Analysis Set-up 2-2 u.i.v. Censoring: C i Simulated power of for Exp(0.2) in both groups. 2-sided weighted logrank test in T n ( w i ) (optimal for direction w i ) Projection test φ 0,n (covering all 3 directions.) n 1 = n 2 = 50, α = 5% and different values of ϑ Realizations of F ϑ,i by von Neuman s Acceptance-Rejection -procedure
34 Simulations and Data Analysis Results for proportional hazards
35 Simulations and Data Analysis Results for crossing hazards
36 Simulations and Data Analysis Resultats for central hazards
37 Simulations and Data Analysis Analysis of the data example
38 Simulations and Data Analysis Analysis of the data example Choose from class of weights w r,g (x) = x r (1 x) g, 0 x 1 Fleming/Harrington (1991)
39 Simulations and Data Analysis Results Kidney Data test with weight p-value (r,g) = (1,5) (r,g) = (2,4) (r,g) = (4,2) (r,g) = (5,1) (r,g) = (0,0) (logrank test) projection test φ 1,n 0.02
40 Simulations and Data Analysis One more example: Survival times of patients with tongue cancer Differences: Group 1: euploide cells Group 2: aneuploide cells (i.e. cells with abnormal number of chromosomes) Question: Can aneuploide cells be used as a prognostic parameter for survival time?
41 Simulations and Data Analysis Source: Klein/Moeschberger
42 Simulations and Data Analysis Results Tongue Data test with weight p-value (r,g) = (1,5) (r,g) = (2,4) (r,g) = (4,2) (r,g) = (5,1) (r,g) = (0,0) (logrank test) projection test φ n 0.09
43 More complex Multi-State Models Competing Risks Model (easiest case) Often: More than one event of interest! (X t ) t 0 càdlàg, X t : Ω {0, 1, 2}. 0 = initial state, P(X 0 = 0) = 1, 1 and 2 absorbing states (competing events). T = inf{t > 0 X t 0} (event time) X T {1, 2} (event) Regulated by cause-specific hazard intensities α j (t) P( T [t, t + ), X T = j T t) α j (t) = α 0j (t) = lim, j = 1, 2. 0
44 More complex Multi-State Models Cumulative Incidence Functions Aim: Statistical inference for CIFs F j (t) = P( T t, X T = j) = t 0 P( T > u )α j (u)du, j = 1, 2,
45 More complex Multi-State Models Study with 1 i n independent patients, (X (i) t ) t 0 independently right-censored and left-truncated 6 Counting processes n Y (t) = i=1 Y i(t) = No. under risk at t N j (t) = n i=1 N j,i(t) = Observed j-events in [0, t] Aalen-Johansen estimator for the CIFs: F j (t) = t 0 P( T > u )dn j (u), j = 1, 2. Y (u) Remark: are local L 2 -martingales! s M j,i (s) = N j,i (s) Y i (u)α j (u) du 0 6 can be relaxed as explained in Andersen et al. (1993, Chapter III); only multiplicative intensity model needed
46 More complex Multi-State Models Martingale Representation Let t < τ. Under sup u [0,t] with inf u y(u) > 0, it follows W n (t) = n 1/2 { F 1 (t) F 1 (t)} n { t = n 1/2 F 1 (t) i=1 t 0 Y (u)/n y(u) 0 0 p S 2 (u)dm 1,i (u) + Y (u) d(m 1,i + M 2,i )(u) Y (u) t 0 } + o P (1) for t < τ, where S 2 = 1 F 2, see Andersen et al. (1993). D Consequence: W n U on D[0, t] by Rebolledo F 1 (u)dm 2,i (u) Y (u)
47 More complex Multi-State Models Problems Covariance function ζ of the Gaussian process U unknown and lacks independent increments Solution: Apply Resampling version of W n 1. Possibility: Wild Bootstrap: Concrete: G j,i, 1 i n, 1 j 2, i.i.d. and data with (µ, σ 2 ) = (0, 1). Approx L(W n ) by cond dist of n { t Ŵ n (t) = n 1/2 F 1 (t) i=1 t 0 0 Ŝ 2 (u)g 1,i dn 1,i (u) Y (u) t + G 1,i dn 1,i (u) + G 2,i dn 2,i (u) Y (u) 0 F 1 (u)g 2,i dn 2,i (u) }. Y (u) Lin s resampling technique as special case for G j,i i.i.d. N(0, 1)
48 More complex Multi-State Models CCLT for Ŵn Theorem (Beyersmann et al., 2013) We have cond conv on D[0, t] given data Ŵ n D U in probability Results in various inference proc based on W n and crit. values from Ŵn as simultaneous CBs for Fj, see e.g. Beyersmann et al. (2013) or differences of CIFs from 2 ind groups tests for or = in the 2-sample case, see Bajrounaite and Klein (2007, 2008) and later... Q: Can we also apply other resampling techniques, e.g. the classical or even weighted bootstrap?
49 More complex Multi-State Models Answers A1: Application of Efron s or even the weighted bootstrap depends on the inference problem! Ex1: Testing for ordered CIFs in unpaired 2-sample problem works via studentization! (details in Dobler and Pauly, 2013) Ex2: Testing for equality does not work in general! Reason: Centering and too complicated limit distribution (e.g., weighted χ 2 ). A2: Yes, e.g. the Weird Bootstrap or the more general, data-dependent Wild Bootstrap (DDWB) Details:...
50 More complex Multi-State Models Rewrite n t Ŵ n (t) = n {G 1/2 1,i F 1 (t) = 2n i=1 t 0 2n i=1 0 G 1,i dn 1,i (u) Y (u) G 2n,i Z 2n,i (t). Ŝ 2 (u)dn 1,i (u) Y (u) + G 2,i t + G 2,i dn 2,i (u) Y (u) for i.i.d. (G 2n,i ) i with (µ, σ 2 ) = (0, 1). DDWB-weights D 2n,i conditionally independent given Z DDWB vers of AJE: Ŵ n D = 2n 2n i=1 D 2n,iZ 2n,i. Theorem (Dobler and Pauly) } 0 F 1 (u)dn 2,i (u) Y (u) If the weights satisfy a cond Lindeberg ass (and some regularity conditions), then we have cond conv on D[0, t] Ŵ D n D U in probability.
51 More complex Multi-State Models Examples Lin s Resampling technique and the (independent) Wild Bootstrap. The Weird Bootstrap of Andersen et al. (1993) corresponds to independent weights (D 2n,i ) i 2n = ((D i,j ) j=1,2 ) i n with 7 D i,j = ( B ( Y ( T i ), J( T i )/Y ( T i ) ) ) 1. Close in spirit to Wild Bootstrap with Poisson weights... 7 Ti = inf{s t : N i (s) > 0} t, where N i = N 1,i + N 2,i
52 More complex Multi-State Models Applications: 2-sample tests for CIFs 2 independent groups k = 1, 2, each with competing risks j = 1, 2. Null hypotheses 8 : H : {F (1) 1 st F (2) 1 } or H = : {F (1) 1 = F (2) 1 } Typical test statistic: Functional of n1 n 2 (1) (2) W n1 n 2 (t) = { F n 1 (t) F 1 (t)} Special case: T 1,n = I W n 1 n 2 (t)dt or T 2,n = I W 2 n 1 n 2 (t)dt 8 on interval I [0, τ)
53 More complex Multi-State Models Applications: 2-sample tests for ordered CIFs We have T 1,n = I W D n 1 n 2 (t)dt T 1 N(0, σζ 2), σ2 ζ unknown9 Solution: Direct resampling with DDWB or Resampling of a studentized version with T 1,stud := T 1,n /V 1,n T1 /σ ζ N(0, 1) For the general weighted bootstrap (and also for the DDWB) it can be shown that D T 1,stud := T 1,n / V 1,n D N(0, 1) in probability. Gives plenty of consistent resampling tests for H. Example: Bootstrap or permutation test 9 σζ 2 = ζ(s, t)dsdt I I
54 More complex Multi-State Models Applications: 2-sample tests for equality of CIFs We have T 2,n = Wn 2 D 1 n 2 (t)dt I j=1 λ j Z 2 j with Z j i.i.d. N(0, 1), λ j unknown Solution: Direct resampling with DDWB works due to Theorem 4. Additional possibilities/extensions: DDWB of standardized test statistic T2,n stan = (T 2,n Ê(T2,n))/ŜD(T 2,n) works as well! Other approximation techniques also! Gives plenty of consistent tests for H =.
55 More complex Multi-State Models Remarks and Outlook Done: Weighted Bootstrap and DDWB for the AJE in CR Generalizing the wild bootstrap and Lin s resampling technique In addition: Comparison of the different testing procedures (For H ) DDWB for more complex Multi-State models (theory )...
56 Literature Literature Andersen, P.K., Borgan, O., Gill, R.D. & Keiding, N. (1993). Statistical Models Based on Counting Processes. Springer, New York. Behnen, K. & Neuhaus, G. (1989). Rank Tests with Estimated Scores and Their Application. B.G. Teubner, Stuttgart. Beyersmann, J., Di Termini, S. & Pauly, M. (2013). Weak convergence of the wild bootstrap for the Aalen-Johansen estimator of the cumulative incidence function of a competing risk. Scand. J. Stat. 40, Brendel, M. & Janssen, A. & Mayer, C.-D. & Pauly, M. (2014). Weighted logrank permutation tests for randomly right censored life science data. Scand. J. Stat., Early View. Dobler, D. & Pauly, M. (2013). How to Bootstrap Aalen-Johansen Processes for Competing Risks? Handicaps, Solutions and Limitations. Submitted. Ehm, W., Mammen, E. & Mueller, D.W. (1995). Power robustification of approximately linear tests. J.A.S.A Fleming, T.R. & Harrington, D.P. (1991). Counting Processes and Survival Analysis. Wiley, New York.
57 Literature Literature Janssen, A. & Mayer, C.-D. (2001). Conditional studentized survival tests for randomly censored models. Scand.J.Stat. 28, Klein, J.P. & Moeschberger, M.L. (2003). Survival analysis: Techniques for censored and truncated data. Springer, New York. Neuhaus, G. (1993). Conditional rank tests for the two-sample problem under random censorship. Ann.Stat. 21, Neuhaus, G. (2000). A method of constructing rank tests in Survival Analysis. J.Stat.Plan.Inf. 91, Pauly, M. (2011). Weighted resampling of martingale difference arrays with applications, Electronic Journal of Statistics, 5, Strasser, H. and Weber, C. (1999). The asymptotic theory of permutation statistics. Math. Methods Stat. 8, No. 2,
58 Thank you for your attention!
STAT 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 informationHow to Bootstrap Aalen-Johansen Processes for Competing Risks? Handicaps, Solutions and Limitations.
How to Bootstrap Aalen-Johansen Processes for Competing Risks? Handicaps, Solutions and Limitations. Dennis Dobler and Markus Pauly November 13, 213 Heinrich-Heine University of Duesseldorf, Mathematical
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 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 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 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 informationA COMPARISON OF POISSON AND BINOMIAL EMPIRICAL LIKELIHOOD Mai Zhou and Hui Fang University of Kentucky
A COMPARISON OF POISSON AND BINOMIAL EMPIRICAL LIKELIHOOD Mai Zhou and Hui Fang University of Kentucky Empirical likelihood with right censored data were studied by Thomas and Grunkmier (1975), Li (1995),
More informationarxiv: v3 [math.st] 12 Oct 2015
Approximative Tests for the Equality of Two Cumulative Incidence Functions of a Competing Risk arxiv:142.229v3 [math.st] 12 Oct 215 Dennis Dobler and Markus Pauly September 28, 218 University of Ulm, Institute
More informationTESTS FOR LOCATION WITH K SAMPLES UNDER THE KOZIOL-GREEN MODEL OF RANDOM CENSORSHIP Key Words: Ke Wu Department of Mathematics University of Mississip
TESTS FOR LOCATION WITH K SAMPLES UNDER THE KOIOL-GREEN MODEL OF RANDOM CENSORSHIP Key Words: Ke Wu Department of Mathematics University of Mississippi University, MS38677 K-sample location test, Koziol-Green
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 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 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 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 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 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 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 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 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 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 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 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 informationUNIVERSITÄT POTSDAM Institut für Mathematik
UNIVERSITÄT POTSDAM Institut für Mathematik Testing the Acceleration Function in Life Time Models Hannelore Liero Matthias Liero Mathematische Statistik und Wahrscheinlichkeitstheorie Universität Potsdam
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 informationEfficient Semiparametric Estimators via Modified Profile Likelihood in Frailty & Accelerated-Failure Models
NIH Talk, September 03 Efficient Semiparametric Estimators via Modified Profile Likelihood in Frailty & Accelerated-Failure Models Eric Slud, Math Dept, Univ of Maryland Ongoing joint project with Ilia
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 informationAnalysis of competing risks data and simulation of data following predened subdistribution hazards
Analysis of competing risks data and simulation of data following predened subdistribution hazards Bernhard Haller Institut für Medizinische Statistik und Epidemiologie Technische Universität München 27.05.2013
More informationMulti-state models: prediction
Department of Medical Statistics and Bioinformatics Leiden University Medical Center Course on advanced survival analysis, Copenhagen Outline Prediction Theory Aalen-Johansen Computational aspects Applications
More informationChapter 7 Fall Chapter 7 Hypothesis testing Hypotheses of interest: (A) 1-sample
Bios 323: Applied Survival Analysis Qingxia (Cindy) Chen Chapter 7 Fall 2012 Chapter 7 Hypothesis testing Hypotheses of interest: (A) 1-sample H 0 : S(t) = S 0 (t), where S 0 ( ) is known survival function,
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 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 informationSTAT331. 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 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 informationarxiv: v1 [math.st] 15 Jul 2018
More powerful logrank permutation tests for two-sample survival data arxiv:1807.05504v1 [math.st] 15 Jul 2018 Marc Ditzhaus and Sarah Friedrich Address of the authors Ulm University Institute of Statistics
More informationUnderstanding product integration. A talk about teaching survival analysis.
Understanding product integration. A talk about teaching survival analysis. Jan Beyersmann, Arthur Allignol, Martin Schumacher. Freiburg, Germany DFG Research Unit FOR 534 jan@fdm.uni-freiburg.de It is
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 informationChapter 4 Fall Notations: t 1 < t 2 < < t D, D unique death times. d j = # deaths at t j = n. Y j = # at risk /alive at t j = n
Bios 323: Applied Survival Analysis Qingxia (Cindy) Chen Chapter 4 Fall 2012 4.2 Estimators of the survival and cumulative hazard functions for RC data Suppose X is a continuous random failure time with
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 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 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 informationGoodness-of-Fit Tests With Right-Censored Data by Edsel A. Pe~na Department of Statistics University of South Carolina Colloquium Talk August 31, 2 Research supported by an NIH Grant 1 1. Practical Problem
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 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 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 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 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 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 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 informationEfficiency Comparison Between Mean and Log-rank Tests for. Recurrent Event Time Data
Efficiency Comparison Between Mean and Log-rank Tests for Recurrent Event Time Data Wenbin Lu Department of Statistics, North Carolina State University, Raleigh, NC 27695 Email: lu@stat.ncsu.edu Summary.
More informationPhilosophy and Features of the mstate package
Introduction Mathematical theory Practice Discussion Philosophy and Features of the mstate package Liesbeth de Wreede, Hein Putter Department of Medical Statistics and Bioinformatics Leiden University
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 informationOn the Goodness-of-Fit Tests for Some Continuous Time Processes
On the Goodness-of-Fit Tests for Some Continuous Time Processes Sergueï Dachian and Yury A. Kutoyants Laboratoire de Mathématiques, Université Blaise Pascal Laboratoire de Statistique et Processus, Université
More informationOn the Breslow estimator
Lifetime Data Anal (27) 13:471 48 DOI 1.17/s1985-7-948-y On the Breslow estimator D. Y. Lin Received: 5 April 27 / Accepted: 16 July 27 / Published online: 2 September 27 Springer Science+Business Media,
More informationSurvival Times (in months) Survival Times (in months) Relative Frequency. Relative Frequency
Smooth Goodness-of-Fit Tests in Hazard-Based Models by Edsel A. Pe~na Department of Statistics University of South Carolina at Columbia E-Mail: pena@stat.sc.edu Univ. of Virginia Colloquium Talk Department
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 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 informationComparing Distribution Functions via Empirical Likelihood
Georgia State University ScholarWorks @ Georgia State University Mathematics and Statistics Faculty Publications Department of Mathematics and Statistics 25 Comparing Distribution Functions via Empirical
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 informationConsistency of bootstrap procedures for the nonparametric assessment of noninferiority with random censorship
Consistency of bootstrap procedures for the nonparametric assessment of noninferiority with random censorship Gudrun Freitag 1 and Axel Mun Institut für Mathematische Stochasti Georg-August-Universität
More informationSurvival Times (in months) Survival Times (in months) Relative Frequency. Relative Frequency
Smooth Goodness-of-Fit Tests in Hazard-Based Models by Edsel A. Pe~na Department of Statistics University of South Carolina at Columbia E-Mail: pena@stat.sc.edu Univ. of Georgia Colloquium Talk Department
More informationSurvival Regression Models
Survival Regression Models David M. Rocke May 18, 2017 David M. Rocke Survival Regression Models May 18, 2017 1 / 32 Background on the Proportional Hazards Model The exponential distribution has constant
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 informationThe Wild Bootstrap for Multivariate Nelson-Aalen Estimators
arxiv:1602.02071v2 [stat.me] 3 Feb 2017 The Wild Bootstrap for Multivariate Nelson-Aalen Estimators Tobias Bluhmki, Dennis Dobler, Jan Beyersmann, Markus Pauly Ulm University, Institute of Statistics,
More informationAnalysis of Time-to-Event Data: Chapter 6 - Regression diagnostics
Analysis of Time-to-Event Data: Chapter 6 - Regression diagnostics Steffen Unkel Department of Medical Statistics University Medical Center Göttingen, Germany Winter term 2018/19 1/25 Residuals for the
More informationEMPIRICAL ENVELOPE MLE AND LR TESTS. Mai Zhou University of Kentucky
EMPIRICAL ENVELOPE MLE AND LR TESTS Mai Zhou University of Kentucky Summary We study in this paper some nonparametric inference problems where the nonparametric maximum likelihood estimator (NPMLE) are
More informationMultistate Modeling and Applications
Multistate Modeling and Applications Yang Yang Department of Statistics University of Michigan, Ann Arbor IBM Research Graduate Student Workshop: Statistics for a Smarter Planet Yang Yang (UM, Ann Arbor)
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 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 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 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 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 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 informationModelling geoadditive survival data
Modelling geoadditive survival data Thomas Kneib & Ludwig Fahrmeir Department of Statistics, Ludwig-Maximilians-University Munich 1. Leukemia survival data 2. Structured hazard regression 3. Mixed model
More informationConcepts and Tests for Trend in Recurrent Event Processes
JIRSS (2013) Vol. 12, No. 1, pp 35-69 Concepts and Tests for Trend in Recurrent Event Processes R. J. Cook, J. F. Lawless Department of Statistics and Actuarial Science, University of Waterloo, Ontario,
More informationNonparametric rank based estimation of bivariate densities given censored data conditional on marginal probabilities
Hutson Journal of Statistical Distributions and Applications (26 3:9 DOI.86/s4488-6-47-y RESEARCH Open Access Nonparametric rank based estimation of bivariate densities given censored data conditional
More informationfor Time-to-event Data Mei-Ling Ting Lee University of Maryland, College Park
Threshold Regression for Time-to-event Data Mei-Ling Ting Lee University of Maryland, College Park MLTLEE@UMD.EDU Outline The proportional hazards (PH) model is widely used in analyzing time-to-event data.
More informationOn the generalized maximum likelihood estimator of survival function under Koziol Green model
On the generalized maximum likelihood estimator of survival function under Koziol Green model By: Haimeng Zhang, M. Bhaskara Rao, Rupa C. Mitra Zhang, H., Rao, M.B., and Mitra, R.C. (2006). On the generalized
More informationA CHANGE-POINT MODEL FOR SURVIVAL DATA WITH LONG-TERM SURVIVORS
Statistica Sinica 19 (2009), 377-390 A CHANGE-POINT MODEL FOR SURVIVAL DATA WITH LONG-TERM SURVIVORS Xiaobing Zhao 1, Xianyi Wu 2 and Xian Zhou 3 1 Jiang Nan University, 2 East China Normal University
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 informationPower and Sample Size Calculations with the Additive Hazards Model
Journal of Data Science 10(2012), 143-155 Power and Sample Size Calculations with the Additive Hazards Model Ling Chen, Chengjie Xiong, J. Philip Miller and Feng Gao Washington University School of Medicine
More informationPackage Rsurrogate. October 20, 2016
Type Package Package Rsurrogate October 20, 2016 Title Robust Estimation of the Proportion of Treatment Effect Explained by Surrogate Marker Information Version 2.0 Date 2016-10-19 Author Layla Parast
More informationLikelihood ratio confidence bands in nonparametric regression with censored data
Likelihood ratio confidence bands in nonparametric regression with censored data Gang Li University of California at Los Angeles Department of Biostatistics Ingrid Van Keilegom Eindhoven University of
More informationHarvard University. Harvard University Biostatistics Working Paper Series. A New Class of Rank Tests for Interval-censored Data
Harvard University Harvard University Biostatistics Working Paper Series Year 2008 Paper 93 A New Class of Rank Tests for Interval-censored Data Guadalupe Gomez Ramon Oller Pique Harvard School of Public
More informationFULL LIKELIHOOD INFERENCES IN THE COX MODEL: AN EMPIRICAL LIKELIHOOD APPROACH
FULL LIKELIHOOD INFERENCES IN THE COX MODEL: AN EMPIRICAL LIKELIHOOD APPROACH Jian-Jian Ren 1 and Mai Zhou 2 University of Central Florida and University of Kentucky Abstract: For the regression parameter
More informationSurvival Distributions, Hazard Functions, Cumulative Hazards
BIO 244: Unit 1 Survival Distributions, Hazard Functions, Cumulative Hazards 1.1 Definitions: The goals of this unit are to introduce notation, discuss ways of probabilistically describing the distribution
More informationIntroduction to Empirical Processes and Semiparametric Inference Lecture 01: Introduction and Overview
Introduction to Empirical Processes and Semiparametric Inference Lecture 01: Introduction and Overview Michael R. Kosorok, Ph.D. Professor and Chair of Biostatistics Professor of Statistics and Operations
More informationLecture 2: Martingale theory for univariate survival analysis
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
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 informationFrailty Models and Copulas: Similarities and Differences
Frailty Models and Copulas: Similarities and Differences KLARA GOETHALS, PAUL JANSSEN & LUC DUCHATEAU Department of Physiology and Biometrics, Ghent University, Belgium; Center for Statistics, Hasselt
More informationEmpirical Likelihood in Survival Analysis
Empirical Likelihood in Survival Analysis Gang Li 1, Runze Li 2, and Mai Zhou 3 1 Department of Biostatistics, University of California, Los Angeles, CA 90095 vli@ucla.edu 2 Department of Statistics, The
More informationPENALIZED LIKELIHOOD PARAMETER ESTIMATION FOR ADDITIVE HAZARD MODELS WITH INTERVAL CENSORED DATA
PENALIZED LIKELIHOOD PARAMETER ESTIMATION FOR ADDITIVE HAZARD MODELS WITH INTERVAL CENSORED DATA Kasun Rathnayake ; A/Prof Jun Ma Department of Statistics Faculty of Science and Engineering Macquarie University
More informationEffects of a Misattributed Cause of Death on Cancer Mortality
Effects of a Misattributed Cause of Death on Cancer Mortality by Jinkyung Ha A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy (Biostatistics) in
More informationPractice Exam 1. (A) (B) (C) (D) (E) You are given the following data on loss sizes:
Practice Exam 1 1. Losses for an insurance coverage have the following cumulative distribution function: F(0) = 0 F(1,000) = 0.2 F(5,000) = 0.4 F(10,000) = 0.9 F(100,000) = 1 with linear interpolation
More informationRobust estimates of state occupancy and transition probabilities for Non-Markov multi-state models
Robust estimates of state occupancy and transition probabilities for Non-Markov multi-state models 26 March 2014 Overview Continuously observed data Three-state illness-death General robust estimator Interval
More information8. Parametric models in survival analysis General accelerated failure time models for parametric regression
8. Parametric models in survival analysis 8.1. General accelerated failure time models for parametric regression The accelerated failure time model Let T be the time to event and x be a vector of covariates.
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 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 informationSurvival Analysis. Lu Tian and Richard Olshen Stanford University
1 Survival Analysis Lu Tian and Richard Olshen Stanford University 2 Survival Time/ Failure Time/Event Time We will introduce various statistical methods for analyzing survival outcomes What is the survival
More informationStat 710: Mathematical Statistics Lecture 31
Stat 710: Mathematical Statistics Lecture 31 Jun Shao Department of Statistics University of Wisconsin Madison, WI 53706, USA Jun Shao (UW-Madison) Stat 710, Lecture 31 April 13, 2009 1 / 13 Lecture 31:
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 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 information