Faculty of Health Sciences. Cox regression. Torben Martinussen. Department of Biostatistics University of Copenhagen. 20. september 2012 Slide 1/51
|
|
- Thomasina Eaton
- 6 years ago
- Views:
Transcription
1 Faculty of Health Sciences Cox regression Torben Martinussen Department of Biostatistics University of Copenhagen 2. september 212 Slide 1/51
2 Survival analysis Standard setup for right-censored survival data. IID copies of (T, D) where T = T C D = I(T C) ubiquitous with T being the true survival time and C the (potential) censoring time and possibly covariates X i (t). Hazard-function: α(t) = lim h 1 h P(t T < t + h T t, F t ). Counting process: N i (t) = I(T i t, D i = 1) Martingale: M i (t) = N i (t) Λ i (t) where Λ i (t) = t Y i(s)α(s) ds (compensator), Y i (t) = I(t T i ) (at risk process). Intensity: Y (t)α(t). Slide 2/51
3 Modelling Survival Data To model and analyze survival data and deal with censorings in a convenient way one often specifies the model through the hazard rate: α(t) = risk of dying among people at risk at time t The hazard function is closely related to the survival probability in other words So t P(T > t) = S(t) = exp( α(s)ds) logs(t) = A(t) = t α(t) = d dt logs(t) α(s)ds Notation, may use both α(t) and λ(t) for the hazard function. Forgive me. But keep in mind that intensity function is Y (t)α(t). Slide 3/51
4 Parametric Models To give precise and clear answers. Kaplan-Meier curves and log-rank tests are not always sufficient. By making simplifying reasonable assumptions something more can be said through parametric models. The simplest possible survival model simply claims that the intensity is constant over time λ(t) = λ, Λ(t) = λt When the rate is constant the survival time is exponentially distributed. One may validate the model by considering the non-parametric estimator of the cumulative intensity (Nelson-Aalen) ˆΛ(t) Should be approximately linear if the model is correct Slide 4/51
5 Malignant melanoma In the period patients had their tumour removed and were followed until At the end of 1977: 57 died of mgl. mel. 14 died of non-related mgl. mel. 134 were still alive. Purpose: Study effect on survival of sex, age, thickness of tumour, ulceration,.. Slide 5/51
6 Malignant melanoma N time status sex age year thickness ulcer Slide 6/51
7 Parametric models - the Weibull model Let X, a (p-dimensional) covariate. The Weibull regression model has hazard function α θ (t) = λγ(λt) γ 1 exp(x T β), (1) where θ = (λ, γ, β); β denotes the regression parameters. The hazard (1) may be written as a so-called proportional hazards model λ (t) exp(x T β), Likelihood function based on T i = Ti C i D i = I(Ti C i ), X i (ignoring censoring part) is : L(θ) = i α θ (T i ) D i e T i αθ (t) dt Standard MLE applies, n(ˆθ θ) D N(, Σ), ˆΣ = I 1 (ˆθ) where I(θ) is minus the derivative of U(θ). How would you compute ˆθ if you were to implement it in R? Slide 7/51
8 Parametric models - the Weibull model Do exercise 3.7 (a) in MS, p. 76. Read in the melanoma-data, and fit the Weibull-model with only sex in the model, see p. 69 in MS Any suggestion on how to check the validity of the Weibull-model in this case? Slide 8/51
9 Cox-Regression Cox-regression is the regression technique for survival analysis, Cox (1972). Per 25. The two most cited statistical papers: (1) With 25,869 citations (currently cited 1,984 times per year): Kaplan, E. L. and Meier, P. (1958) Nonparametric estimation from incomplete observations, Journal of the American Statistical Association, 53, pp (2) With 18,193 citations (1,342 per year), Cox, D. R. (1972) Regression models and life tables, Journal of the Royal Statistical Society, Series B, 34, pp Regression techniques are very useful for dealing with many covariates Can be used to learn about treatment effect while correcting for other covariates. Cox model: λ i (t) = λ (t) exp(β 1 X i β px ip ), where λ (t) is baseline hazard for a subject with covariates. Note: λ (t) is not further specified! Slide 9/51
10 The Cox model The regression coefficients β 1,..., β p represent the effects of the covariates. β 1 is the effect of X i1 when we have corrected for the other covariates. β 1 may be interpreted in terms of the relative risk when the covariate X i1 is increased 1: λ (t) exp(β 1 (X i1 + 1) β px ip ) λ (t) exp(β 1 X i β px ip ) = exp (β 1 ) If β 1 > the risk of dying increases as X i1 increases, and if β 1 < the risk of dying decreases as X i1 increases. The quantity ˆβ 1 X i ˆβ px ip is called the prognostic index for the ith subject. Slide 1/51
11 Two-Sample Cox Model Consider a simple 2-sample situation where we wish to study effect of sex. Defining a covariate for the ith patient { 1 Male X i = Female The hazard can be written as λ i (t) = λ (t) exp(βx i ) giving the contributions λ (t) exp(β) when X i = 1 and λ (t) when X i =. Note again that λ (t) is unspecified meaning that which distribution is unspecified? In a moment we derive estimators for the Cox model, but let s already use them. Can be calculated in R (or SAS, Stata, SPSS). forskere-paaviser-klar-sammenhaeng-mellem-antistof-og-leddegigt/ BMJ-paper Slide 11/51
12 Example: Melanoma data Consider the Cox-model for the melanoma with the explanatory variables sex and log(thickness) (lthick): We get: λ i (t) = λ (t) exp (β 1 sex i + β 2 lt i ) > melanoma$lthick=log(melanoma$thick) > fit=coxph(surv(days,status==1)~ factor(sex)+lthick,data=melanoma) > fit Call: coxph(formula = Surv(days, status == 1) ~ factor(sex) + lthick, data = melanoma) coef exp(coef) se(coef) z p factor(sex) e-2 lthick e-7 Likelihood ratio test=33.5 on 2 df, p=5.45e-8 n= 25 What does the relative risks of 1.58 and 2.18 mean? Slide 12/51
13 PBC-data PBC data (primary biliary cirrhosis): 418 patients are followed until death or censoring. PBC is a fatal chronic liver disease. Important explanatory variables: Age: in years Albumin: serum albumin (mg/dl) Bilirubin: serum bilirunbin (mg/dl) Edema: no edema,.5 untreated or successfully treated 1 edema despite diuretic therapy Prothrombin time: standardised blood clotting time Slide 13/51
14 Fitting Cox s model in R > library(survival) > data(pbc) > with(pbc,cbind(time,status,age,edema,bili,protime,albumin)[1:5,]) time status age edema bili protime albumin [1,] [2,] [3,] [4,] [5,] > > fit.pbc<-coxph(surv(time, status==2) ~ age+factor(edema)+log(bili)+log(protime)+ log(albumin),data=pbc) > fit.pbc Call: coxph(formula = Surv(time, status == 2) ~ age + factor(edema) + log(bili) + log(protime) + log(albumin), data = pbc) coef exp(coef) se(coef) z p age e-7 factor(edema) e-1 factor(edema) e-4 log(bili) e+ log(protime) e-3 log(albumin) e-4 Likelihood ratio test=232 on 6 df, p= n=416 (2 observations deleted due to missingness) Slide 14/51
15 Likelihood construction Likelihood function based on T i = Ti C i D i = I(Ti C i ), X i is: i α(t i ) D i e T i α(t) dt = i { α(t) N i (t) } e τ t Y i (t)α(t) dt τ is end of observation period Y i (t) = I(t T i ) is the at risk indicator. Can also write it as { da(t) N i (t) } e τ i t Y i (t)da(t) where A(t) = t α(s) ds. Cox-model: A(t) = A (t)e X T β with A (t) = t α (s) ds. Slide 15/51
16 Likelihood construction Let S (t, β) = n i=1 Y i(t) exp(xi T β). Based on martingale equation can you then suggest an estimator of A (t) assuming β is known? Slide 16/51
17 Likelihood construction The "estimator" Â (t, β) = is called the Breslow-estimator. t 1 S (t, β) dn (s) Now plug this estimator into the likelihood function, and arrive at a function only depending on β! This function, L(β) is called Cox s partial likelihood function. Slide 17/51
18 Cox s proportional hazards model The Cox model is λ i (t) = Y i (t)λ (t) exp(x T i β), (2) where X = (X 1,..., X p) is a p-dimensional locally bounded predictable covariate and Y i (t) is an at risk indicator. The regression parameter β is estimated as the maximizer to Cox s partial likelihood function L(β) = ( ) exp (X T Ni (t) i β), (3) S (t, β) t where S (t, β) = i n i=1 Λ (t) is estimated by the Breslow estimator ˆΛ (t) = ˆΛ (t, ˆβ) = Y i (t) exp(xi T β). t 1 S (t, ˆβ) dn (t). (4) Slide 18/51
19 Cox s proportional hazards model We have presented theory with time-const. X s, but can easily handle time-var. covariates. Now let s look into the asymptotical properties of Cox s estimators. The first partial derivative of S (t, β) with respect to β is denoted : S 1 (t, β) = n i=1 Y i (t) exp(xi T β)x i Show that ˆβ is found as the solution to the score equation U( ˆβ) =, where U(β) = n i=1 τ (X i (t) E(t, β))dn i (t) E(t, β) = S 1(t, β) S (t, β). (5) Show also that U(β ) can be written as a martingale evaluated at τ. Here, β denotes the true β. Slide 19/51
20 Asymptotic properties Theorem Under the standard conditions, then, as n, n 1/2 U(β ) D N(, Σ), n 1/2 ( ˆβ β ) D N(, Σ 1 ), and Σ is estimated consistently by n 1 I( ˆβ), with I minus the derivative of U. Under the standard conditions, then, as n, n 1/2 (ˆΛ (t, ˆβ) Λ (t)) D U(t) where U(t) is a zero-mean Gaussian process with covariance function Φ(t). Not a martingale though, why? Get back to how Φ(t) can be estimated consistently Slide 2/51
21 Asymptotic properties Key to the proof is that the score evaluated in the true point β is a martingale (evaluated at τ): U(β ) = n τ i=1 (X i E(t, β ))dm i (t) (6) The predictable variation process of n 1/2 U(β ) is n 1/2 U(β ) = n 1 n where i=1 τ (X i (t) E(t, β )) 2 Y i (t) exp(xi T (t)β )dλ (t) τ = n 1 V (t, β )S (t, β )dλ (t) Σ. P V (t, β) = S 2(t, β) S (t, β) E(t, β) 2. (7) Lindeberg condition of martingale CLT is also fulfilled Slide 21/51
22 Asymptotic properties It follows that U(β ) converges in distribution to a normal variate with zero-mean and variance Σ. A Taylor series expansion of the score gives n 1/2 ( ˆβ β ) = (n 1 I(β )) 1 n 1/2 U(β ), where β is on the line segment between β and ˆβ. Consistency of ˆβ and the results above give that n 1/2 ( ˆβ β ) converges to the postulated normal distribution. Slide 22/51
23 IID decomposition The score process evaluated at β where n 1/2 U(β, t) = n 1/2 ɛ 1i (t) = t n ɛ 1i (t) + o p(1) i=1 ( X i s ) 1(β, t) dm i (s) s (β, t) with s j (β, t) the limit in prob. of n 1 S j (β, t). That is a sum of zero-mean iid terms! The variance of n 1/2 ( ˆβ β ) may be estimated consistently by { n } Σ = ni 1 ( ˆβ, τ) ˆɛ 2 1i (τ) I 1 ( ˆβ, τ). where ˆɛ 1i (τ) is given by ɛ 1i (τ) replacing unknowns with their empirical counterparts. i=1 Slide 23/51
24 IID decomposition More important! Try the same for the Breslow-estimator. n 1/2 (ˆΛ (t, ˆβ) Λ (t)) n 1/2 (ˆΛ (t, β ) Λ (t)) + D β ˆΛ (t, ˆβ)n 1/2 ( ˆβ β ) D β ˆΛ (t, ˆβ) conv. in probability so last term is essentially a sum of zero-mean iid s. But also n 1/2 (ˆΛ (t, β ) Λ (t)) = since n 1 S (s, β) P s (s, β). t n 1/2 t 1 S (s, β) dm (s) 1 s (s, β) dm i(s), Slide 24/51
25 IID decomposition It thus follows that n 1/2 (ˆΛ (t, ˆβ) Λ (t)) is asymptotically equivalent to where n 1/2 n ɛ 2i (t), (8) i=1 ɛ 2i (t) =ɛ 3i (t) + H T (β, t)i(β, τ) 1 ɛ 1i (τ), ɛ 3i (t) = t S 1 (β, s)dm i (s), and where H(β, t) is the derivative of t S 1 (β, s)dn.(s) wrt β. The variance of n 1/2 (ˆΛ (t, ˆβ) Λ (t)) thus may be estimated by n 1 n i=1 ɛ 2 2i (t) Slide 25/51
26 Resampling Now, n 1/2 (ˆΛ (t, ˆβ) Λ (t)) is asymptotically equivalent to n 1/2 where G 1,..., G n are iid N(, 1). n ˆɛ 2i (t)g i (9) This is very useful for constructing confidence bands for the survival predictions and to make tests about the baseline and survival function predictions. Taylor expansion: i=1 n 1/2 (S (ˆΛ, ˆβ, { t) S (Λ, β, t)) = S (Λ, β, t) n 1/2 exp(x T β)(ˆλ (t) Λ (t)) + Λ (t) exp(x T β)x T ( ˆβ } β). (1) Slide 26/51
27 Resampling Just saw that : n 1/2 ( ˆβ β, ˆΛ Λ ) was asymptotically equivalent to the processes 1 = n 1/2 2 (t) = n 1/2 n ˆɛ 1i G i, i=1 n ˆɛ 2i (t)g i, where G 1,..., G n are independent standard normals. It follows that n 1/2 (Ŝ S ) has the same asymptotic distribution as i=1 where n S (t) = S (t)n 1/2 ˆɛ 4i (t)g i, (11) i=1 ɛ 4i (t) = exp(z T β)ɛ 2i (t) + A(t) exp(x T β)x T ɛ 1i. Slide 27/51
28 Cox-Regression: checking assumptions Wish to study a treatment effect while correcting for other variables. May use Cox-model: λ i (t) = λ (t) exp(β 1 X i β px ip ) λ (t) is the baseline hazard for a subject with covariates. The Cox models ability to deal with many covariates comes from the regression structure. Some assumptions have been made The effects of covariates are additive and linear on the log risk scale. If covariates interact with each other the regression model should include interaction terms. The relative risk between the hazard rate for two subjects is constant over time c(β 1,..., β p) = λ i (t) λ i (t) We will try to check all these assumptions with the latter being the key assumption. Slide 28/51
29 Interaction between categorical and cont. variables Melanoma data: Consider the Cox-model for melanoma data with explanatory variables sex and log(thickness) (lt): λ i (t) = λ (t) exp (β 1 sex i + β 2 lt i ) How do we check for interaction in this situation? Group the continuous variable (log(thickness)) and proceed as before. Allow for two different β 2 s, one for each value of sex. Fit this latter model in R, and interpret output. Slide 29/51
30 Interaction between categorical and cont. variables > fit=coxph(surv(days,status==1)~ factor(sex)+factor(sex)*lthick,data=melanoma) > fit Call: coxph(formula = Surv(days, status == 1) ~ factor(sex) + factor(sex) * lthick, data = melanoma) coef exp(coef) se(coef) z p factor(sex) e-1 lthick e-5 factor(sex)1:lthick e-1 Cox-model: λ (t) exp (β 1 sex + β 2 lt + β 3 sexlt) with { β1 + (β β 1 sex + β 2 lt + β 3 sexlt = 2 + β 3 ) lt, sex = 1 β 2 lt, sex = The baseline λ (t) is the hazard for? The effect of lt is estimated to? for male and to? for female Is there an interaction between sex and lt? Slide 3/51
31 Checking proportionality Constant relative risk between the hazard rates for two subjects c(β 1,..., β p) = λ i (t) λ i (t) In the 2-sample case the assumption is equivalent to proportional intensities for the two groups, i.e., exp(β 1 )λ 1 (t) = λ 2 (t). Implies that the cumulative intensities are proportional Λ 2 (t) = t λ 2 (s)ds = exp(β 1 )Λ 1 (t) A plot of ˆΛ 2 (t) versus ˆΛ 1 (t) should give a straight line through (,) with slope exp(β 1 ). Similarly log(λ 2 (t)) log(λ 1 (t)) = β 1 so plotting log(ˆλ k (t)), k = 1, 2, versus t should give parallel curves. Slide 31/51
32 The Stratified Cox model The stratified Cox model contains different baselines for different strata : λ ik (t) = λ k (t) exp(β 1 X i β px ip ) k = 1,..., K λ k (t) is the baseline hazard for a subject in strata k. The regression coefficients β 1,..., β p represent the effects of the covariates as in the simple Cox regression model. Melanoma-data: we might stratify according to ulceration. The baselines give the mortality of the two defined by ulceration (yes/no) when all other covariates are. Slide 32/51
33 Log-cumulative hazard plot > fit=coxph(surv(days,status==1)~ factor(sex)+lthick+strata(ulc),data=melanoma) > fit.detail=coxph.detail(fit) > logcum1=log(cumsum(fit.detail$hazard[c(17:57)])) > logcum=log(cumsum(fit.detail$hazard[1:16])) > plot(c(fit.detail$time[17:57],37),c(logcum1,logcum1[41]),type= s,xlab="time (in days)",ylab="logcums") > lines(c(fit.detail$time[1:16],37),c(logcum,logcum[16]),type= s,lty=2) Logcums Time (in days) Parallel? Slide 33/51
34 Checking Proportionality: tests Graphical test may reveal serious violations, but they may be difficult to use in general. If effect is expected to wear off with time or to increase with time, one may try to describe this effect and then fit a Cox model with a time-varying explanatory variable f (t) exp(β + θ f (t)) = λ 1(t) λ 2 (t) For a given choice of f (t) we can then test if θ is, i.e, if the time-varying effect can be left out of the model. A common choice of f (t) is ln(t). More specialized methods and programs are available. Slide 34/51
35 Checking Proportionality: tests > fit=coxph(surv(days/365,status==1)~ factor(sex)+ulc+lthick,data=melanoma) > time.test=cox.zph(fit,transform="log") > time.test rho chisq p factor(sex) ulc lthick GLOBAL NA Based on this, the Cox model is rejected. Slide 35/51
36 Checking Proportionality: tests Another option is to have a time-changing effect of ulc. Imagine for instance that effect of ulc is most important in an initial phase. We try the first 4 years, approx 14 days. We hence replace β ulc ulc by β ulc (t) ulc with Note that with β ulc (t) = β ulc,1 + β ulc,2 I(t > 14) β ulc (t) ulc = (β ulc,1 + β ulc,2 I(t > 14)) ulc being a time-varying covariate. = β ulc,1 ulc + β ulc,2 X(t) X(t) = I(t > 14)) ulc So this is just a Cox-model again, but now with a time-varying covariate. Slide 36/51
37 Time-changing effect of covariates Such a model can easily be fitted in R using the survsplit-function. > melanoma1=survsplit(melanoma,cut=c(14),end="days",start="start", event="status") > melanoma1$ulcnew=melanoma1$ulc*as.numeric(melanoma1$days>14) > > fit1=coxph(surv(start,days,status==1)~ factor(sex)+lthick+ actor(ulc)+ factor(ulcnew), data=melanoma1) > summary(fit1) Call: coxph(formula = Surv(start, days, status == 1) ~ factor(sex) + lthick + factor(ulc) + factor(ulcnew), data = melanoma1) n= 367 coef exp(coef) se(coef) z Pr(> z ) factor(sex) lthick ** factor(ulc) *** factor(ulcnew) * --- Relative risk Ulc/no Ulc in the first 4 years, and after the first 4 years? Conclusion? Slide 37/51
38 Cox s proportional hazards model The procedures described above are the traditional goodness-of-fit tools. Make tests against specific deviations: Replace X 1 with (X 1, X 1 (log (t))), say (β 1 β 1 + β p+1 log (t)). Test the null β p+1 =. Graphical method Time-changing effect. These methods are quite useful but also have some limitations: Graphical method: Not parallel. What is acceptable? What if a given covariate is continuous? Test: Ad hoc method. Which transformation to use? Time-changing effect. Also ad-hoc in that the time where effects changes are rarely known in practice. All methods: They assume that model is ok for all the other covariates. Slide 38/51
39 Cumulative martingale residuals Alternative: Cumulative martingale residuals, Lin, Wein and Ying (1993). The martingales under the Cox regression model can be written as t M i (t) = N i (t) Y i (s) exp(x T β)dλ (s) ˆM i (t) = N i (t) = N i (t) t t Y i (s) exp(x T ˆβ)d ˆΛ (s) Y i (s) exp(x T i (s) ˆβ) 1 S (s, ˆβ) dn (s). One idea is now to look at different groupings of the these residuals and see if they behave as they should under the model. Slide 39/51
40 Cumulative martingale residuals The score function, evaluated in the estimate ˆβ, and seen as a function of time, can for example be written as U( ˆβ, t) = n i=1 t (X i S 1(s, ˆβ) S (s, ˆβ) )dn i(s) = n i=1 t X i d ˆM i (s) so it is equivalent to cumulating the martingale residuals (in certain way). Also, by Taylor series expansion, n 1/2 U( ˆβ, t) n 1/2 U(β, t) (n 1 I( ˆβ, t))n 1/2 ( ˆβ β) n 1/2 ( ˆβ β) (n 1 I( ˆβ, τ)) 1 n 1/2 U(β, τ) Hence n 1/2 U( ˆβ, { } t) n 1/2 U(β, t) I( ˆβ, t)(i( ˆβ, τ)) 1 U(β, τ) Slide 4/51
41 Cumulative martingale residuals n 1/2 U( ˆβ, t) is thus asymptotically equivalent to the process ( ) n 1/2 M 1 (t) I(t, ˆβ)I 1 (τ, ˆβ)M 1 (τ), (12) where M 1 (t) = n M 1i (t) = i=1 n t i=1 (X i (s) e(s, β ))dm i (s) with e(t, β ) = lim p E(t, β ). Note: n 1/2 U( ˆβ, τ) = Slide 41/51
42 Cumulative martingale residuals The distribution of the process n 1/2 M 1 (t) (t [, τ]) is asymptotically equivalent to n 1/2 n i=1 t (X i (s) E(s, ˆβ))dN i (s)g i where G 1,..., G n are independent standard normals. Since M i s are i.i.d. with variance E(N i ). Alternatively to this resampling approach one may also ([?]), show n 1/2 n ˆM 1i (t)g i, with ˆM 1i (t) = i=1 t (X i (s) E(s, ˆβ))d ˆM i (s). Therefore it can be shown that n 1/2 U( ˆβ, t) is asymptotically equivalent to n 1/2 n i=1 ( ˆM1i (t) I(t, ˆβ)I 1 (τ, ˆβ) ˆM 1i (τ)) G i, for almost any sequence of the counting processes that determines ˆM and I. Slide 42/51
43 Cumulative martingale residuals To summarize the cumulative score process plots one may look at test statistics like sup U j ( ˆβ, t), j = 1..., p, t [,τ] Slide 43/51
44 PBC-data > library(timereg) > fit=cox.aalen(surv(days,status==1)~ prop(factor(sex))+prop(lthick)+ prop(ulc), data=melanoma) > summary(fit) Proportional Cox terms : Coef. SE Robust SE D2log(L)^-1 z P-val prop(factor(sex)) prop(lthick) prop(ulc) Test for Proportionality sup hat U(t) p-value H_ prop(factor(sex)) prop(lthick) prop(ulc) > plot(fit,score=t) Slide 44/51
45 PBC-data prop(factor(sex))1 prop(lthick) Cumulative coefficients Cumulative coefficients Time Time prop(ulc) Cumulative coefficients Time Slide 45/51
46 Cumulative Residuals LWY also suggested : M c(t, z) = t K T z (s)d ˆM(s) where K z(t) is an n 1 matrix with elements I(X i1 (t) z) for i = 1,..., n focusing here on the first continuous covariate X 1, say. Thus cumulating residuals versus both time and the covariate values. A 1-dimensional process : M c(z) = τ which can be plotted against z. If X 1 (t) = X 1 then τ K T z (t)d ˆM(t) = i process in z. K T z (t)d ˆM(t), (13) ˆM i (τ)i(x i1 z) which is a This will reveal if the functional form is appropriate. These processes can also be resampled. Slide 46/51
47 Cumulative residuals > fit.gof=cox.aalen(surv(days,status==1)~ prop(lthick),melanoma, residuals=1) > > resids=cum.residuals(fit.gof,melanoma,cum.resid=1) > summary(resids) Test for cumulative MG-residuals Grouped cumulative residuals not computed, you must provide modelmatrix to get these (see help) Residual versus covariates consistent with model sup hat B(t) p-value H_: B(t)= Call: cum.residuals(fit.gof, melanoma, cum.resid = 1) Slide 47/51
48 Cumulative martingale residuals Cumulative residuals Slide 48/51
49 Cumulative residuals > fit.gof=cox.aalen(surv(days,status==1)~ prop(thick),melanoma, residuals=1) > resids=cum.residuals(fit.gof,melanoma,cum.resid=1) > summary(resids) Test for cumulative MG-residuals Grouped cumulative residuals not computed, you must provide modelmatrix to get these (see help) Residual versus covariates consistent with model sup hat B(t) p-value H_: B(t)= Call: cum.residuals(fit.gof, melanoma, cum.resid = 1) Slide 49/51
50 Cumulative martingale residuals Cumulative residuals Slide 5/51
51 Summary Cox s proportional hazards model. Nice model with easily interpretable parameters - relative risks! Is is used heavily in Biostatistcs. Are the relative risks really not depending on time? Check model carefully. Resampling techniques useful for evaluating asymptotics distributions. Useful goodness-of-fit based on cumulative residuals with p-values. Slide 51/51
Goodness-Of-Fit for Cox s Regression Model. Extensions of Cox s Regression Model. Survival Analysis Fall 2004, Copenhagen
Outline Cox s proportional hazards model. Goodness-of-fit tools More flexible models R-package timereg Forthcoming book, Martinussen and Scheike. 2/38 University of Copenhagen http://www.biostat.ku.dk
More informationSurvival analysis in R
Survival analysis in R Niels Richard Hansen This note describes a few elementary aspects of practical analysis of survival data in R. For further information we refer to the book Introductory Statistics
More informationSurvival analysis in R
Survival analysis in R Niels Richard Hansen This note describes a few elementary aspects of practical analysis of survival data in R. For further information we refer to the book Introductory Statistics
More informationOutline. Cox's regression model Goodness-of-t methods. Cox's proportional hazards model: Survival analysis
Outline Cox's regression model Goodness-of-t methods Giuliana Cortese gico@biostat.ku.dk Advanced Survival Analysis 21 Copenhagen Cox's proportional hazards model: assumptions Goodness-of-t for Cox's model:
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 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 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 information( t) Cox regression part 2. Outline: Recapitulation. Estimation of cumulative hazards and survival probabilites. Ørnulf Borgan
Outline: Cox regression part 2 Ørnulf Borgan Department of Mathematics University of Oslo Recapitulation Estimation of cumulative hazards and survival probabilites Assumptions for Cox regression and check
More informationOutline. Frailty modelling of Multivariate Survival Data. Clustered survival data. Clustered survival data
Outline Frailty modelling of Multivariate Survival Data Thomas Scheike ts@biostat.ku.dk Department of Biostatistics University of Copenhagen Marginal versus Frailty models. Two-stage frailty models: copula
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 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 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 informationADVANCED STATISTICAL ANALYSIS OF EPIDEMIOLOGICAL STUDIES. Cox s regression analysis Time dependent explanatory variables
ADVANCED STATISTICAL ANALYSIS OF EPIDEMIOLOGICAL STUDIES Cox s regression analysis Time dependent explanatory variables Henrik Ravn Bandim Health Project, Statens Serum Institut 4 November 2011 1 / 53
More informationβ j = coefficient of x j in the model; β = ( β1, β2,
Regression Modeling of Survival Time Data Why regression models? Groups similar except for the treatment under study use the nonparametric methods discussed earlier. Groups differ in variables (covariates)
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 informationLecture 12. Multivariate Survival Data Statistics Survival Analysis. Presented March 8, 2016
Statistics 255 - Survival Analysis Presented March 8, 2016 Dan Gillen Department of Statistics University of California, Irvine 12.1 Examples Clustered or correlated survival times Disease onset in family
More informationMore on Cox-regression
More on Cox-regression p. 1/45 More on Cox-regression STK4080 H16 1. Repetition 2. Left truncation 3. Time-dependent covariates 4. Stratified Cox-regression 5. Residuals - Model check 6. How to handle
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 informationPart [1.0] Measures of Classification Accuracy for the Prediction of Survival Times
Part [1.0] Measures of Classification Accuracy for the Prediction of Survival Times Patrick J. Heagerty PhD Department of Biostatistics University of Washington 1 Biomarkers Review: Cox Regression Model
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 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 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 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 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 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 informationLog-linearity for Cox s regression model. Thesis for the Degree Master of Science
Log-linearity for Cox s regression model Thesis for the Degree Master of Science Zaki Amini Master s Thesis, Spring 2015 i Abstract Cox s regression model is one of the most applied methods in medical
More informationLecture 7 Time-dependent Covariates in Cox Regression
Lecture 7 Time-dependent Covariates in Cox Regression So far, we ve been considering the following Cox PH model: λ(t Z) = λ 0 (t) exp(β Z) = λ 0 (t) exp( β j Z j ) where β j is the parameter for the the
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 informationProportional hazards regression
Proportional hazards regression Patrick Breheny October 8 Patrick Breheny Survival Data Analysis (BIOS 7210) 1/28 Introduction The model Solving for the MLE Inference Today we will begin discussing regression
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 informationTHESIS for the degree of MASTER OF SCIENCE. Modelling and Data Analysis
PROPERTIES OF ESTIMATORS FOR RELATIVE RISKS FROM NESTED CASE-CONTROL STUDIES WITH MULTIPLE OUTCOMES (COMPETING RISKS) by NATHALIE C. STØER THESIS for the degree of MASTER OF SCIENCE Modelling and Data
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 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 informationOutline. Frailty modelling of Multivariate Survival Data. Clustered survival data. Clustered survival data
Outline Frailty modelling of Multivariate Survival Data Thomas Scheike ts@biostat.ku.dk Department of Biostatistics University of Copenhagen Marginal versus Frailty models. Two-stage frailty models: copula
More informationLecture 10. Diagnostics. Statistics Survival Analysis. Presented March 1, 2016
Statistics 255 - Survival Analysis Presented March 1, 2016 Dan Gillen Department of Statistics University of California, Irvine 10.1 Are model assumptions correct? Is the proportional hazards assumption
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 informationRelative-risk regression and model diagnostics. 16 November, 2015
Relative-risk regression and model diagnostics 16 November, 2015 Relative risk regression More general multiplicative intensity model: Intensity for individual i at time t is i(t) =Y i (t)r(x i, ; t) 0
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 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 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 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 informationAccelerated Failure Time Models
Accelerated Failure Time Models Patrick Breheny October 12 Patrick Breheny University of Iowa Survival Data Analysis (BIOS 7210) 1 / 29 The AFT model framework Last time, we introduced the Weibull distribution
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 informationExtensions of Cox Model for Non-Proportional Hazards Purpose
PhUSE Annual Conference 2013 Paper SP07 Extensions of Cox Model for Non-Proportional Hazards Purpose Author: Jadwiga Borucka PAREXEL, Warsaw, Poland Brussels 13 th - 16 th October 2013 Presentation Plan
More informationLecture 8 Stat D. Gillen
Statistics 255 - Survival Analysis Presented February 23, 2016 Dan Gillen Department of Statistics University of California, Irvine 8.1 Example of two ways to stratify Suppose a confounder C has 3 levels
More informationChapter 4 Regression Models
23.August 2010 Chapter 4 Regression Models The target variable T denotes failure time We let x = (x (1),..., x (m) ) represent a vector of available covariates. Also called regression variables, regressors,
More informationThe coxvc_1-1-1 package
Appendix A The coxvc_1-1-1 package A.1 Introduction The coxvc_1-1-1 package is a set of functions for survival analysis that run under R2.1.1 [81]. This package contains a set of routines to fit Cox models
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 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 informationMAS3301 / MAS8311 Biostatistics Part II: Survival
MAS330 / MAS83 Biostatistics Part II: Survival M. Farrow School of Mathematics and Statistics Newcastle University Semester 2, 2009-0 8 Parametric models 8. Introduction In the last few sections (the KM
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 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 informatione 4β e 4β + e β ˆβ =0.765
SIMPLE EXAMPLE COX-REGRESSION i Y i x i δ i 1 5 12 0 2 10 10 1 3 40 3 0 4 80 5 0 5 120 3 1 6 400 4 1 7 600 1 0 Model: z(t x) =z 0 (t) exp{βx} Partial likelihood: L(β) = e 10β e 10β + e 3β + e 5β + e 3β
More informationTime-dependent coefficients
Time-dependent coefficients Patrick Breheny December 1 Patrick Breheny Survival Data Analysis (BIOS 7210) 1/20 Introduction As we discussed previously, stratification allows one to handle variables that
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 informationMultistate models and recurrent event models
Multistate models Multistate models and recurrent event models Patrick Breheny December 10 Patrick Breheny Survival Data Analysis (BIOS 7210) 1/22 Introduction Multistate models In this final lecture,
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 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 informationBeyond GLM and likelihood
Stat 6620: Applied Linear Models Department of Statistics Western Michigan University Statistics curriculum Core knowledge (modeling and estimation) Math stat 1 (probability, distributions, convergence
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 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 informationSemiparametric Regression
Semiparametric Regression Patrick Breheny October 22 Patrick Breheny Survival Data Analysis (BIOS 7210) 1/23 Introduction Over the past few weeks, we ve introduced a variety of regression models under
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 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 informationMultistate models and recurrent event models
and recurrent event models Patrick Breheny December 6 Patrick Breheny University of Iowa Survival Data Analysis (BIOS:7210) 1 / 22 Introduction In this final lecture, we will briefly look at two other
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 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 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 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 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 informationApproximation of Survival Function by Taylor Series for General Partly Interval Censored Data
Malaysian Journal of Mathematical Sciences 11(3): 33 315 (217) MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES Journal homepage: http://einspem.upm.edu.my/journal Approximation of Survival Function by Taylor
More informationTied survival times; estimation of survival probabilities
Tied survival times; estimation of survival probabilities Patrick Breheny November 5 Patrick Breheny Survival Data Analysis (BIOS 7210) 1/22 Introduction Tied survival times Introduction Breslow approximation
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 informationAnalysis of Time-to-Event Data: Chapter 4 - Parametric regression models
Analysis of Time-to-Event Data: Chapter 4 - Parametric regression models Steffen Unkel Department of Medical Statistics University Medical Center Göttingen, Germany Winter term 2018/19 1/25 Right censored
More informationChapter 7: Hypothesis testing
Chapter 7: Hypothesis testing Hypothesis testing is typically done based on the cumulative hazard function. Here we ll use the Nelson-Aalen estimate of the cumulative hazard. The survival function is used
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 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 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 informationNonparametric two-sample tests of longitudinal data in the presence of a terminal event
Nonparametric two-sample tests of longitudinal data in the presence of a terminal event Jinheum Kim 1, Yang-Jin Kim, 2 & Chung Mo Nam 3 1 Department of Applied Statistics, University of Suwon, 2 Department
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 informationST745: Survival Analysis: Cox-PH!
ST745: Survival Analysis: Cox-PH! Eric B. Laber Department of Statistics, North Carolina State University April 20, 2015 Rien n est plus dangereux qu une idee, quand on n a qu une idee. (Nothing is more
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 informationST495: Survival Analysis: Maximum likelihood
ST495: Survival Analysis: Maximum likelihood Eric B. Laber Department of Statistics, North Carolina State University February 11, 2014 Everything is deception: seeking the minimum of illusion, keeping
More informationP n. This is called the law of large numbers but it comes in two forms: Strong and Weak.
Large Sample Theory Large Sample Theory is a name given to the search for approximations to the behaviour of statistical procedures which are derived by computing limits as the sample size, n, tends to
More informationTime-dependent covariates
Time-dependent covariates Rasmus Waagepetersen November 5, 2018 1 / 10 Time-dependent covariates Our excursion into the realm of counting process and martingales showed that it poses no problems to introduce
More informationMulti-state Models: An Overview
Multi-state Models: An Overview Andrew Titman Lancaster University 14 April 2016 Overview Introduction to multi-state modelling Examples of applications Continuously observed processes Intermittently observed
More informationLecture 9. Statistics Survival Analysis. Presented February 23, Dan Gillen Department of Statistics University of California, Irvine
Statistics 255 - Survival Analysis Presented February 23, 2016 Dan Gillen Department of Statistics University of California, Irvine 9.1 Survival analysis involves subjects moving through time Hazard may
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 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 informationEmpirical Likelihood Confidence Band
University of Kentucky UKnowledge Theses and Dissertations--Statistics Statistics 2015 Empirical Likelihood Confidence Band Shihong Zhu University of Kentucky, shihong2015@gmail.com Click here to let us
More informationStatistics in medicine
Statistics in medicine Lecture 4: and multivariable regression Fatma Shebl, MD, MS, MPH, PhD Assistant Professor Chronic Disease Epidemiology Department Yale School of Public Health Fatma.shebl@yale.edu
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 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 informationLogistic regression model for survival time analysis using time-varying coefficients
Logistic regression model for survival time analysis using time-varying coefficients Accepted in American Journal of Mathematical and Management Sciences, 2016 Kenichi SATOH ksatoh@hiroshima-u.ac.jp Research
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 informationA Regression Model For Recurrent Events With Distribution Free Correlation Structure
A Regression Model For Recurrent Events With Distribution Free Correlation Structure J. Pénichoux(1), A. Latouche(2), T. Moreau(1) (1) INSERM U780 (2) Université de Versailles, EA2506 ISCB - 2009 - Prague
More informationOn a connection between the Bradley-Terry model and the Cox proportional hazards model
On a connection between the Bradley-Terry model and the Cox proportional hazards model Yuhua Su and Mai Zhou Department of Statistics University of Kentucky Lexington, KY 40506-0027, U.S.A. SUMMARY This
More informationCase-control studies
Matched and nested case-control studies Bendix Carstensen Steno Diabetes Center, Gentofte, Denmark b@bxc.dk http://bendixcarstensen.com Department of Biostatistics, University of Copenhagen, 8 November
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 information