More on Cox-regression
|
|
- Erick Goodman
- 5 years ago
- Views:
Transcription
1 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 departures from prop.haz. assumption
2 More on Cox-regression p. 2/45 Repetition Cox-regression With covariate x i the survival time T i has hazard α(t x i ) = exp(β x i )α 0 (t) where α 0 (t) is called basis (underlying) hazard and whereβ is a regressionsparameter. With D i indicator for death for individual i, T i right censoring time and R(t) = {i : T i t} = the risk set (right before) time t we estimate β by maximizing the Cox partial likelihood [ ] n exp(β Di x i ) n [ ] exp(β L(β) = Di x i ) = k R( T i ) exp(β x k ) S (0) (β,t) i=1 where S (0) (β,t) = k R( T i ) exp(β x k ) i=1
3 if the l -model hasq more parameters than theˆl-model and that ˆl-model is true (H0 -model). More on Cox-regression p. 3/45 Rep. Cox-regr., contd. In counting process notation we saw that we may express the score U(β) = log(l(β)) β = n i=1 [ ] x i S(1) (β,t) dn S (0) i (t). (β,t) Also L(β) may be treated as a regular likelihood. Thus ˆβ N(β,I(ˆβ) 1 ). where I(β) = 2 log(l(β)) and for two nested models with β 2 maximum log-partial likelihood respectively l andˆl we have LRT = 2(l ˆl) χ 2 q
4 More on Cox-regression p. 4/45 Estimation of cumulative hazarda 0 (t) = t 0 α 0(s)ds The common estimator fora 0 (t) is the Breslow-estimator  0 (t) = Ti t D i k R( T i ) exp(ˆβ x k ) = t 0 dn (s) S (0) (ˆβ,s) Note the similarity with the Nelson-Aalen estimator. Given Â0(t) it is simple to estimate cumulative hazard for an individual with hazardα(t x i ) = exp(β x i )α 0 (t) as Â(t x i ) = exp(ˆβ x i )Â0(t) The survival function S(t x i ) = P(T i > t x i ) can be estimated by Ŝ(t x i ) = exp( Â(t x i)) = exp( Â0(t)exp(ˆβ x i )) = Ŝ0(t) exp(ˆβ x i ).
5 No arguments have used the property of non-decreasing Y i (t) (might be a good exercise to check). More on Cox-regression p. 5/45 Cox-regression for left truncated and rightcens. data We may use the same likelihood expression in counting process notation. We may estimate β by maximizing the partial likelihood, or solving U(β) = log(l(β)) β = n i=1 [ ] x i S(1) (β,t) dn S (0) i (t) = 0 (β,t) where S (0) (β,t) = n i=1 Y i (t)exp(β x i ). We need only remember that Y i (t) need not be non-decreasing, but Y i (t) = I(L i < t T i ).
6 More on Cox-regression p. 6/45 R: Cox-reg. with l.trunc. and r.cens. The program need information on left truncation time L i in addition to ( T i,d i,x i ). Example: Time until death psychiatric patients. Data: x i = 1 or 2 for men / women L i = age of first time admitted to psychiatric ward (year) T i = age at death/censoring (year) D i = indicator of death. > coxph(surv(ageonset,agedeath,death) sex) coef exp(coef) se(coef) z p sex Likelihood ratio test=0.43 on 1 df, p=0.514 n= 26
7 More on Cox-regression p. 7/45 Time dependent covariates x i (t) Risk factors may depend on time: x i (t) = smoker at time t (yes/no) x i (t) = cum. no. cigarettes (pack years) smoked at age t x i (t) = no. years since quitting smoking It may furthermore be that the proportional hazards model does not fit the data well with the included covariates, but that a valid model is given by. α(t x i ) = exp(β 1 x i +β 2 x i t)α 0 (t) so that the effect of x i becomes larger (smaller)according to β 2 > 0(< 0). May code this by introducing new covariates, for instance x i2 (t) = tx i.
8 Cox-regression allows for time dependentx i (t)! Model: α i (t) = exp(β x i (t))α 0 (t). Will simply solve U(β) = log(l(β)) n [ ] = x i S(1) (β,t) dn β S (0) i (t) = 0 (β,t) i=1 (as before) with the difference that S (0) (β,t) = n i=1 Y i (t)exp(β x i (t)). depend on time dependent x i (t). Similarly to left truncation all theory goes through exactly the same way as before. For the sake of repetition we check E[U(β)] = 0. However: We need to assume that x i (t) is predictable! More on Cox-regression p. 8/45
9 More on Cox-regression p. 9/45 U(β) martingale with expectation 0. As before S (0) (β,t) β score = n i=1 x i(t)y i (t)exp(β x i (t)) = S (1) (β,t). Thus the U(β) = n i=1 [x i (t) S(1) (β,t) S (0) (β,t) ]dn i(t) With predictable covariatesx i (t) the model can be written dn i (s) = Y i (t)exp(β x i (t))α 0 (t)dt+dm i (t) (with standard interpretations) which gives U(β) = n i=1 [xi (t) S(1) (β,t) ]Y S (0) (β,t) i(t)exp(β x i (t))α 0 (t)dt + n i=1 [xi (t) S(1) (β,t) ]dm S (0) (β,t) i(t) [xi (t) S(1) (β,t) ]dm S (0) (β,t) i(t) = n i=1
10 More on Cox-regression p. 10/45 because n i=1 [x i (t) S(1) (β,t) S (0) (β,t) ]Y i(t)exp(β x i (t))α 0 (t)dt = 0. This since n i=1 [x i(t) S(1) (β,t) S (0) (β,t) ]Y i(t)exp(β x i (t)) and so = S (1) (β,t) S(1) (β,t) S (0) (β,t) S(0) (β,t) = 0 U(β) = n i=1 [x i (t) S(1) (β,t) S (0) (β,t) ]dm i(t) is a sum of integrals wrt. martingales, and itself a martingale with E[U(β)] = 0. just as for time constant and right censored data!
11 Time dependent covariates in R R only allows for time dependent covariates constant on intervals, i.e. step functions. Assume thatx i (t) = x j on interval < L ij,u ij ] for j = 1,2,...,J i. Need to represent this individual J i times in the data fil as left truncated data with L ij as left trunkcation time U ij as right censoring time D ij = D i I( event in interval j) as indicator x j as covariate value Note: Did not assume L i,j+1 = U ij, and individuals may well disappear and reenter later. More on Cox-regression p. 11/45
12 More on Cox-regression p. 12/45 Stratified Cox-regression Assume that the population is divided into s = 1, 2,..., S strata so that person i = 1,2,...,n s within stratumswith covariatex is has hazard α(t x is ) = exp(β x is )α 0s (t) where α 0s (t) is the baseline in stratum s (typically α 0s (t)/α 0s (t) varies with t.) Effects of covariatesβ are the same in all strata We get a partial likelihood from each stratum: [ n s exp(β x is ) L s (β) = k R s (T is ) exp(β x ks ) i=1 ] Dis with D is = indicator for individual is in stratum s etc.
13 More on Cox-regression p. 13/45 Stratified Cox-regression, contd. All strata give information on β. If we assume that the strata are independent (weak assumption) we may then combine this information by maximizing the stratified partial likelihood L(β) = S s=1 L s (β) The corresponding score function becomes U(β) = log(l(β)) β = S s=1 log(l s (β)) β = S s=1 U s (β) where U s (β) is the score function from stratum s. These all have expectation 0, i.e. E[U(β)] = 0.
14 More on Cox-regression p. 14/45 Stratified Cox-regression, III We find information I s from stratumsandi(β) total information becomes, due to independence, S Var[U(β)] = E[I s (β)] = E[I(β)] s=1 Thus may use the stratified partial likelihood as a regular likelihood. The stratified partial likelihood is useful when the proportional model does not hold for a categorical variable. Stratify on this variable and keep regression model for other covariates In particular this is useful when the stratification variable is a confounder and the main interest is on the other variables.
15 More on Cox-regression p. 15/45 Stratified Cox-regression i R Uses the Melanoma-data. Stratifies on grouped tumor thickness by command strata: coxph(surv(lifetime,dead) ulcer+sex+age+strata(grthick),data=mel) coef exp(coef) se(coef) z p ulcer sex age Likelihood ratio test=13.2 on 3 df, p= n= 205
16 More on Cox-regression p. 16/45 What did we gain from this? Compares with the Cox-regression where thickness is a categorical variable > coxph(surv(lifetime,dead) ulcer+sex+age+factor(grthick),data=mel) coef exp(coef) se(coef) z p ulcer sex age factor(grthick) factor(grthick) Likelihood ratio test=45.3 on 5 df, p=1.27e-08 n= 205 The results are only marginally different, but we have modeled in a much more flexible way. There is hardly a change in the standard errors, and the added flexibility did not lead to loss in efficiency.
17 More on Cox-regression p. 17/45 Baseline in each stratum in R. > nycox<-coxph(surv(lifetime,dead) ulcer+sex+age+strata(grthick),data=me > plot(survfit(nycox),fun="cum.haz") The proportional hazards assumption might be checked from these plots (or rather the log-cumulative hazard plots). Better methods are presented later in the lectures.
18 More on Cox-regression p. 18/45 How can the Cox-model fail? α(t x i ) = exp(β x i )α 0 (t)? The Cox-model is flexible wrt the baseline α 0 (t), but otherwise strict with respect to how the hazard depend on covariates, for instance We may have specified covariatex ik wrong, correct alternative may be f.ex. x ik = log(x ik) or x ik = x ik. We do not have a proportional model. The effect may vary with time, f.ex. α(t x i ) = exp(β(t) x i )α 0 (t) whereβ(t) is a function of time.
19 More on Cox-regression p. 19/45 Martingale residuals With a specified model for the hazard, say a proportional hazards modelα i (t) = exp(β x i )α 0 (t) we that M i (t) = N i (t) t 0 Y i (s)exp(β x i )α 0 (s)ds is a martingale we expectation zero. Inserting the Cox-estimator ˆβ forβ, the Breslow estimatordâ0(s) = dn (s) S (0) (ˆβ,s) forα 0(s)ds and the maximal right censored survival time τ we get the so called martingale residuals ˆM i = N i (τ) τ 0 Y i (s)exp(ˆβ x i )dâ0(s) which have the structure "Observed"-"Expected".
20 Example: Melanoma data We saw thatlog(tumorsize) was a better covariate than tumorsize directly. We will check if this can be discovered from martingale residuals. coxfit<-coxph(surv(lifetime,status==1) sex+ulcer+thickn, data=melanoma) martres<-coxfit$residuals plot(melanoma$thickn,martres) martres melanoma$thickn More on Cox-regression p. 20/45
21 Example: Melanoma data, contd. As often it is difficult to read off residual plots directly. We calculate mean mart. resid. for groups of tumor thickness > grthickn<-melanoma$grthick > lm(martres factor(grthickn)-1) Coefficients: factor(grthickn)1 factor(grthickn)2 factor(grthickn) > summary(lm(martres factor(grthickn))) Estimate Std. Error t value Pr(> t ) (Intercept) factor(grthickn) * factor(grthickn) F-statistic: on 2 and 202 DF, p-value: Note: This analysis is not strictly correct, from ABG one may figure out a correct test. But it shows that the martingale residuals are largest for the second group. More on Cox-regression p. 21/45
22 More on Cox-regression p. 22/45 GAM: Generalized Additive Models (Hastie & Tibshirani (1990): By smoothing techniques we may fit Linear model: Y = α+f(x)+ε Logistic model: log( p 1 p ) = α+f(x) where f(x) is some smooth function. This may be done in R with the library gam that may need to be downloaded from CRAN. Similar models may be fitted for survival data under specifications α(t x) = exp(f(x))α 0 (t)
23 More on Cox-regression p. 23/45 Example: Melanoma data, GAM-plot We add a gam-curve with 95% CI to the scatter plot of martingale residuals vs. thickness. thick<-melanoma$thickness) library(gam) plot(gam(martres s(thick)),se=t,ylim=c(min(martres),max(martres))) points(thick,martres) s(thick) thick
24 More on Cox-regression p. 24/45 GAM for Cox-regression in R Regression-spline: R-syntax coxph(surv(time,status) ns(z,df=4)) Cubic smoothing-spline (penalized partial likelihood): R-syntax: coxph(surv(time,status) pspline(z,df=4))
25 Penalized log-partial-likelihood: Maximize, for given smoothing parameterλ, log(l(f)) λ (f (x)) 2 dx where L(f) = Partial likelihood λ = "penalty"-term for curvature of f(x) In particular: λ = : No curvature,f(x) = βx straight line λ = 0 : No smoothing Maximation problem actually has a simple numerical solution. The penalty term λ corresponds to a certain degree of freedom and may be interpreted as df = no. covariates in mod., though df may be any real number > 0. More on Cox-regression p. 25/45
26 More on Cox-regression p. 26/45 Example: Penalized log-partial-likelihood: Melanoma-data, thickness, R-commands: newcoxfit<-coxph(surv(lifetime,status==1) pspline(thickn,df=4),data=mela newmel<-melanoma thckspacing<-min(melanoma$thickn) +(1:205)*(max(melanoma$thickn)-min(melanoma$thickn))/205 newmel$thickn<-thckspacing newpred<-predict(newcoxfit,newmel,type="terms",term=1,se=t) mi<-min(newpred$fit-1.96*newpred$se.fit) ma<-max(newpred$fit+1.96*newpred$se.fit) plot(thckspacing,newpred$fit,type="l",ylim=c(mi,ma),xlab="thickness") lines(thckspacing,newpred$fit+1.96*newpred$se.fit,lty=2) lines(thckspacing,newpred$fit-1.96*newpred$se.fit,lty=2)
27 Example: Penalized log-partial-likelihood: Melanoma-data, thickness, Plot of ˆf(x) against x = thickness: newpred$fit thickness The plot shows what we also found with grouping of thickness, it is the smallest tumors that have smaller risk. More on Cox-regression p. 27/45
28 May smooth with more covariates α(t x) = exp(β 1 x 1 +β 2 x 2 +β 3 x 3 +f 4 (x 4 ))α 0 (t) Ex: Melanoma,x 1 = sex,x 2 = ulc., x 3 = age,x 4 = tumorth. newcoxfit<-coxph(surv(lifetime,status==1) sex+ulcer+age +pspline(thickn,df=4),data=mel newpred<-predict(newcoxfit,newmel,type="terms",term=4,se=t) newpred$fit spline log thickness More on Cox-regression p. 28/45
29 More on Cox-regression p. 29/45 Non-proportional hazards The old-fashioned way of checking departure from proportionality is based on the following: With α(t x) = exp(βx)α 0 (t) we have log(a(t x)) = β x+log(a 0 (t)) i.e. log(a(t x)) for differentxand log(a 0 (t)) should be parallel lines. Thus if x is the only covariate and is categorical we may plot log of Nelson-Aalen for every level ofx If the lines are parallel then proportionality OK
30 More on Cox-regression p. 30/45 Ex. Melanoma: Checks tumor-thickness and ulceration Ulcer Tumortykkelse log(h) log(h) tid tid
31 More on Cox-regression p. 31/45 Multivariable models If we have decided thatx 1,x 2,...,x p should be included in the model and want to check if the categorical covariatex p+1 satisfies proportionality we should rather Fit a stratified Cox-model with x p+1 -levels as strata and x 1,x 2,...,x p as covariates. Plotlog(Âs(t, x)) against t for different levelssof x p+1 Are the lines parallel? Again, it might be an advantage to look at log(âs(t, x)) log(â1(t, x)).
32 More on Cox-regression p. 32/45 Schoenfeld residuals The Schoenfeld-residuals are given by, for T j such thatd j = 1 x jk n i=1 x iky i ( T j )exp(ˆβ x i ) n i=1 Y i( T j )exp(ˆβ x i ) There is thus one residual for event time and for each component of the covariates. Since Y i ( T j )exp(ˆβ x i ) ni=1 Y i ( T j )exp(ˆβ x i ) sums to one they may be thought of as point mass probabilities for some distribution. The interpretation of the distribution is that of the covariatesx j given that individual j experienced the event. And so n i=1 x iky i ( T j )exp(ˆβ x i ) n j=1 Y = x k ( T j ) i( T j )exp(ˆβ x i )
33 More on Cox-regression p. 33/45 More Schoenfeld becomes the expectation in this distribution. The Schoenfeldt-residuals at time T j may also be written x j x( T j ). Furthermore, with U(β) = the scorefunction, 0 = U(ˆβ) = [x j x( T j )] T j :D j =1 a sensible property for a residual. Sometimes component k for x jk x k ( T j ) shows a clear tendency for positive values over intervals (and negative over others). In the "positive" intervals there is the greater risk connected to component k than in the "negative".
34 More on Cox-regression p. 34/45 Schoenfeld, contd. We could have done a local Cox-regression within the "positive" interval. In such case the x(t) would tend to be larger. Thus the Schoenfeld-residuals give information as to whether the prop.haz. assumption holds. They may even be used to estimate how the hazard ratio varies over time. Let V(t) = S(2) (ˆβ,t) S (0) (ˆβ,t) [ S (1) (ˆβ,t) S (0) (ˆβ,t) The observed information is given as V(t)dN(t) and in particular V(t) can be interpreted as the variance of x given event att. ] 2
35 More on Cox-regression p. 35/45 Scaled Schoenfeld-residual If the true model equalsα(t x) = exp(β(t)x)α 0 (t) for some function β(t) we have as a 1.order approximation β( T j ) ˆβ +V( T j ) 1 (x j x( T j )) i.e. varying around the scaled Schoenfeld-residual V( T j ) 1 (x j x( T j )) that measures the departure from mean risk ˆβ. Plot of this quantity - smoothed over time - may give a picture of how the risk varies.
36 More on Cox-regression p. 36/45 Ex: Melanoma data Calculates and smooths scaled Schoenfeld-residuals for all 4 covariates in the melanom data. coxfit<-coxph(surv(lifetime,dead) sex+age+ulcer+logthick,data=mel) plot(cox.zph(coxfit)) Beta(t) for sex Beta(t) for age Time Time Beta(t) for ulcer Beta(t) for logthick Time Time
37 More on Cox-regression p. 37/45 Test for prop. assumption can be based on Schoenfeld-residuals: The tests are directed to departures from the model given by α(t x) = exp((β 0 +θg(t))x)α 0 (t) for specified functions g(t). Actually they are score-tests for a new covariate x (t) = g(t)x. If the model is extended by a new time-dependent covariate x p+1 (t) = x p g(t) the score for the coefficientβ p+1 is given U p+1 = j g( T j )[x p x p ( T j )] as a weighted sum of the Schoenfeld-residuals.
38 More on Cox-regression p. 38/45 Test for proportionality Example: Melanoma data: KM transform. Sex Ulceration Age χ 2 1 = 0.5 p=0.46 χ 2 1 = 0.7 p=0.40 χ 2 1 = 2.8 p=0.09 log(tumor thickness) χ 2 1 = 4.1 p=0.04 Indication for departure for tumor thickness R-syntax: coxfit<-coxph(surv(lifetime,dead) sex+age+ulcer+logthick,data=mel) cox.zph(coxfit) rho chisq p sex age ulcer logthick GLOBAL NA
39 More on Cox-regression p. 39/45 Strategies when proportional hazard fails Stratified Cox-regression Separate analyzes on disjoint time intervals Time-dependent covariates Alternative regression models Accelerated failure time models Additive models
40 More on Cox-regression p. 40/45 Ex. Stratified Cox-regression Weak departure wrt. thickness. Stratifies on grthick: > coxstrat<-coxph(surv(lifetime,dead) sex+age+ulcer+strata(grthick),data > coxstrat coef exp(coef) se(coef) z p sex age ulcer Likelihood ratio test=13.2 on 3 df, p= n= 205 > cox.zph(coxstrat) rho chisq p sex age ulcer GLOBAL NA But possibly the stratification changed other estimates somewhat?
41 Separate intervals We may split the time interval i 2 and make separate Cox-regressions within each interval: Ex: Melanoma data. Half of death before τ = 3 Analysis on [0,3 >: Uses events only if lifetime < 3 coxph(surv(lifetime,dead*(lifetime<3)) sex+age+ulcer+logthick,data=mel) coef exp(coef) se(coef) z p sex age ulcer logthick Analysis on [3, >: Uses only events with lifetime> 3 coxph(surv(lifetime,dead*(lifetime>3)) sex+age+ulcer+logthick,data=mel) coef exp(coef) se(coef) z p sex age ulcer logthick More on Cox-regression p. 41/45
42 More on Cox-regression p. 42/45 Time dependent covariats If similar parameter estimates for age, sex and ulceration on these intervals we may fit a common model log( α(t x) α 0 (t) ) = β 1x 1 +β 2 x 2 +β 3 x 3 +β 4 x 4 I(t < 3)+β 5 x 4 I(t 3) i.e. Cox-regression with time dependent covariates x 4 I(t < 3) x 4 I(t 3)
43 Cox-regression with time dependent covariates Need to set up data-frame with time dependent covariates: mel2 <- data.frame(intime = c(rep(0,205),rep(3,167))) mel2$outtime <- c(pmin(mel$lifetime,3),mel$lifetime[mel$lifetime>3]) mel2$indi <- c(mel$dead*(mel$lifetime<3),mel$dead[mel$lifetime>3]) mel2$sex mel2$ulcer mel2$age <- c(mel$sex,mel$sex[mel$lifetime>3]) <- c(mel$ulcer,mel$ulcer[mel$lifetime>3]) <- c(mel$age,mel$age[mel$lifetime>3]) mel2$logtha <- c(mel$logthick,rep(0,167)) mel2$logthb <- c(rep(0,205),mel$logthick[mel$lifetime>3]) coxph(surv(intime,outtime,indi) sex+ulcer+age+logtha+logthb,data=mel2) coef exp(coef) se(coef) z p sex ulcer age logtha logthb Likelihood ratio test=49 on 5 df, p=2.17e-09 n= 372 More on Cox-regression p. 43/45
44 More on Cox-regression p. 44/45 Advantages/disadvantages with strategies 1. Stratification Easy More difficult to show effect of stratification variable Allows for only a few problem covariates 2. Separate intervals Relatively easy Choice of interval difficult/arbitrary Looses power for covariates where the assumption is OK Many parameter estimates
45 More on Cox-regression p. 45/45 Advantages/disadvantages with strategies 3. Time dependent covariates Somewhat awkward to arrange (in R(?)) Difficult choice of interval Only helpful when prop.haz. OK for most covar. Consequences of departure from proportionality biased estimates of coefficients both for covariates where the assumption hold and fail biased survival estimates
( 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 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 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 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 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 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 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 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 informationFaculty of Health Sciences. Cox regression. Torben Martinussen. Department of Biostatistics University of Copenhagen. 20. september 2012 Slide 1/51
Faculty of Health Sciences Cox regression Torben Martinussen Department of Biostatistics University of Copenhagen 2. september 212 Slide 1/51 Survival analysis Standard setup for right-censored survival
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 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 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 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 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 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 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 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 informationUSING MARTINGALE RESIDUALS TO ASSESS GOODNESS-OF-FIT FOR SAMPLED RISK SET DATA
USING MARTINGALE RESIDUALS TO ASSESS GOODNESS-OF-FIT FOR SAMPLED RISK SET DATA Ørnulf Borgan Bryan Langholz Abstract Standard use of Cox s regression model and other relative risk regression models for
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 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 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 informationReduced-rank hazard regression
Chapter 2 Reduced-rank hazard regression Abstract The Cox proportional hazards model is the most common method to analyze survival data. However, the proportional hazards assumption might not hold. The
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 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 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 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 informationFaculty of Health Sciences. Regression models. Counts, Poisson regression, Lene Theil Skovgaard. Dept. of Biostatistics
Faculty of Health Sciences Regression models Counts, Poisson regression, 27-5-2013 Lene Theil Skovgaard Dept. of Biostatistics 1 / 36 Count outcome PKA & LTS, Sect. 7.2 Poisson regression The Binomial
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 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 informationGoodness-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 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 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 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 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 informationCorrelation and regression
1 Correlation and regression Yongjua Laosiritaworn Introductory on Field Epidemiology 6 July 2015, Thailand Data 2 Illustrative data (Doll, 1955) 3 Scatter plot 4 Doll, 1955 5 6 Correlation coefficient,
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 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 informationDynamic Prediction of Disease Progression Using Longitudinal Biomarker Data
Dynamic Prediction of Disease Progression Using Longitudinal Biomarker Data Xuelin Huang Department of Biostatistics M. D. Anderson Cancer Center The University of Texas Joint Work with Jing Ning, Sangbum
More informationLecture 2: Poisson and logistic regression
Dankmar Böhning Southampton Statistical Sciences Research Institute University of Southampton, UK S 3 RI, 11-12 December 2014 introduction to Poisson regression application to the BELCAP study introduction
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 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 informationLecture 11. Interval Censored and. Discrete-Time Data. Statistics Survival Analysis. Presented March 3, 2016
Statistics 255 - Survival Analysis Presented March 3, 2016 Motivating Dan Gillen Department of Statistics University of California, Irvine 11.1 First question: Are the data truly discrete? : Number of
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 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 informationLecture 5: Poisson and logistic regression
Dankmar Böhning Southampton Statistical Sciences Research Institute University of Southampton, UK S 3 RI, 3-5 March 2014 introduction to Poisson regression application to the BELCAP study introduction
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 12: Effect modification, and confounding in logistic regression
Lecture 12: Effect modification, and confounding in logistic regression Ani Manichaikul amanicha@jhsph.edu 4 May 2007 Today Categorical predictor create dummy variables just like for linear regression
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 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 informationStat 642, Lecture notes for 04/12/05 96
Stat 642, Lecture notes for 04/12/05 96 Hosmer-Lemeshow Statistic The Hosmer-Lemeshow Statistic is another measure of lack of fit. Hosmer and Lemeshow recommend partitioning the observations into 10 equal
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 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 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 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 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 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 information9. Estimating Survival Distribution for a PH Model
9. Estimating Survival Distribution for a PH Model Objective: Another Goal of the COX model Estimating the survival distribution for individuals with a certain combination of covariates. PH model assumption:
More 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 information11 November 2011 Department of Biostatistics, University of Copengen. 9:15 10:00 Recap of case-control studies. Frequency-matched studies.
Matched and nested case-control studies Bendix Carstensen Steno Diabetes Center, Gentofte, Denmark http://staff.pubhealth.ku.dk/~bxc/ Department of Biostatistics, University of Copengen 11 November 2011
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 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 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 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 informationMüller: Goodness-of-fit criteria for survival data
Müller: Goodness-of-fit criteria for survival data Sonderforschungsbereich 386, Paper 382 (2004) Online unter: http://epub.ub.uni-muenchen.de/ Projektpartner Goodness of fit criteria for survival data
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 informationCausal Hazard Ratio Estimation By Instrumental Variables or Principal Stratification. Todd MacKenzie, PhD
Causal Hazard Ratio Estimation By Instrumental Variables or Principal Stratification Todd MacKenzie, PhD Collaborators A. James O Malley Tor Tosteson Therese Stukel 2 Overview 1. Instrumental variable
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 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 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 I (CHL5209H)
Survival Analysis Dalla Lana School of Public Health University of Toronto olli.saarela@utoronto.ca February 3, 2015 21-1 Time matching/risk set sampling/incidence density sampling/nested design 21-2 21-3
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 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 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 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 informationCIMAT Taller de Modelos de Capture y Recaptura Known Fate Survival Analysis
CIMAT Taller de Modelos de Capture y Recaptura 2010 Known Fate urvival Analysis B D BALANCE MODEL implest population model N = λ t+ 1 N t Deeper understanding of dynamics can be gained by identifying variation
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 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 informationExtensions of Cox Model for Non-Proportional Hazards Purpose
PhUSE 2013 Paper SP07 Extensions of Cox Model for Non-Proportional Hazards Purpose Jadwiga Borucka, PAREXEL, Warsaw, Poland ABSTRACT Cox proportional hazard model is one of the most common methods used
More informationBuilding a Prognostic Biomarker
Building a Prognostic Biomarker Noah Simon and Richard Simon July 2016 1 / 44 Prognostic Biomarker for a Continuous Measure On each of n patients measure y i - single continuous outcome (eg. blood pressure,
More informationSSUI: Presentation Hints 2 My Perspective Software Examples Reliability Areas that need work
SSUI: Presentation Hints 1 Comparing Marginal and Random Eects (Frailty) Models Terry M. Therneau Mayo Clinic April 1998 SSUI: Presentation Hints 2 My Perspective Software Examples Reliability Areas that
More informationLecture 14: Introduction to Poisson Regression
Lecture 14: Introduction to Poisson Regression Ani Manichaikul amanicha@jhsph.edu 8 May 2007 1 / 52 Overview Modelling counts Contingency tables Poisson regression models 2 / 52 Modelling counts I Why
More informationModelling counts. Lecture 14: Introduction to Poisson Regression. Overview
Modelling counts I Lecture 14: Introduction to Poisson Regression Ani Manichaikul amanicha@jhsph.edu Why count data? Number of traffic accidents per day Mortality counts in a given neighborhood, per week
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 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 informationA fast routine for fitting Cox models with time varying effects
Chapter 3 A fast routine for fitting Cox models with time varying effects Abstract The S-plus and R statistical packages have implemented a counting process setup to estimate Cox models with time varying
More informationMultivariate Survival Analysis
Multivariate Survival Analysis Previously we have assumed that either (X i, δ i ) or (X i, δ i, Z i ), i = 1,..., n, are i.i.d.. This may not always be the case. Multivariate survival data can arise in
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 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 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 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 informationChapter 17. Failure-Time Regression Analysis. William Q. Meeker and Luis A. Escobar Iowa State University and Louisiana State University
Chapter 17 Failure-Time Regression Analysis William Q. Meeker and Luis A. Escobar Iowa State University and Louisiana State University Copyright 1998-2008 W. Q. Meeker and L. A. Escobar. Based on the authors
More informationMultivariable Fractional Polynomials
Multivariable Fractional Polynomials Axel Benner May 17, 2007 Contents 1 Introduction 1 2 Inventory of functions 1 3 Usage in R 2 3.1 Model selection........................................ 3 4 Example
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 informationGeneralized Linear Models
Generalized Linear Models Lecture 7. Models with binary response II GLM (Spring, 2018) Lecture 7 1 / 13 Existence of estimates Lemma (Claudia Czado, München, 2004) The log-likelihood ln L(β) in logistic
More informationBiostatistics for physicists fall Correlation Linear regression Analysis of variance
Biostatistics for physicists fall 2015 Correlation Linear regression Analysis of variance Correlation Example: Antibody level on 38 newborns and their mothers There is a positive correlation in antibody
More information8 Analysis of Covariance
8 Analysis of Covariance Let us recall our previous one-way ANOVA problem, where we compared the mean birth weight (weight) for children in three groups defined by the mother s smoking habits. The three
More informationJoint Modeling of Longitudinal Item Response Data and Survival
Joint Modeling of Longitudinal Item Response Data and Survival Jean-Paul Fox University of Twente Department of Research Methodology, Measurement and Data Analysis Faculty of Behavioural Sciences Enschede,
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 information