Linear life expectancy regression with censored data
|
|
- Chrystal Copeland
- 5 years ago
- Views:
Transcription
1 Linear life expectancy regression with censored data By Y. Q. CHEN Program in Biostatistics, Division of Public Health Sciences, Fred Hutchinson Cancer Research Center, Seattle, Washington 98109, U.S.A. and S. CHENG Department of Epidemiology & Biostatistics, University of California, San Francisco, California 94143, U.S.A. Summary In the statistical literature, the life expectancy is usually characterised by the mean residual life function. Regression models are thus needed to study the association between the mean residual life functions and their covariates. In this article, we consider a linear mean residual life model and further developed inference procedures in the presence of potential censoring. The new model and inference procedures are applied to the Stanford heart transplant data. Additional semiparametric efficiency calculations and information bound are also considered. Some key words: Additive model; Counting process; Estimating equation; Mean residual life; Semiparametric model 1
2 1 Introduction Suppose that failure time T is a nonnegative random variable on a probability space {Ω, F, P}. In the statistical literature, the life expectancy of T at time t > 0, e(t) say, is usually characterised and studied by way of its mean residual life function, m(t) = E(T t T > t) = e(t) t. When covariates are present, the proportional mean residual life model of Oakes & Dasu (1990) can be used to study the association between m(t) and the covariates: m(t Z) =m 0 (t) exp(β T Z), (1) where Z are the p-vector covariates and β are the associated parameters. To estimate β in (1), Maguluri & Zhang (1994) and Oakes & Dasu (2003) studied some estimation procedures when there is no censoring. Recently, Chen & Cheng (2005) developed quasi partial score estimating equations when censoring is present and m 0 ( ) is unspecified. As a result of the nature of e(t), the shape of m(t) has an intrinsic constraint, in that e(t) = m(t)+ t should be monotonically nondecreasing in t 0. In the proportional mean residual life model, however, m(t Z)+t = m 0 (t) exp(β T Z)+t may not always satisfy this constraint for certain β R p, unless m 0 ( ) itself is monotonically nondecreasing, as pointed out in Oakes & Dasu (1990). A monotonically nondecreasing m 0 ( ), although plausible mathematically, may not be compatible with the actual underlying process, in the case of human ageing, for example. To cope with this constraint, we consider instead a linear mean residual life model, m(t Z) =m 0 (t)+β T Z. (2) It is apparent the additive structure in model (2) complies with the intrinsic constraint. Moreover, model (2) is equivalent to e(t Z) =e 0 (t)+β T Z. Thus, in practice, the parameter β in model (2) can be interpreted as the average difference in life expectancies per unit change in Z. For example, when Z = 0/1 is the treatment indicator in a randomised clinical trial, β can be considered as the treatment effect measured in terms of life expectancy. In general, m 0 (t)+β T Z is required to be nonnegative. 2
3 To apply the new linear model in real applications, it is desirable that m 0 ( ) be unspecified without restrictive parametric assumptions, and the fact that the survival outcomes are often censored has to be accommodated. The rest of this article aims to develop and study some inference procedures for the new linear model. 2 Estimation and Inference 2.1 Estimation of baseline function Suppose that there are n subjects in the dataset, and let T i and C i be the failure and censoring times, respectively, for i = 1, 2,...,n. The observed dataset thus consists of n independent identically distributed copies of {(X i, i,z i ); i =1, 2,...,n}, where X i = min(t i,c i ), i = I(T i C i ) and Z i are the covariates, respectively. Given Z i, T i and C i are assumed independent. Without confusion in notation, subscripts may be occasionally suppressed. Let N(t) =I(X t, = 1) and Y (t) =I(X t). Since the baseline residual life function in model (2) is unspecified, we first develop an estimator for m 0 (t) as if the true β were known. Let the true values of β and m 0 ( ) beβ and m ( ), respectively; and consider the filtration defined by F t = σ{n i (u),y i (u),z i ;0 u t, i =1, 2,...,n}. Then E{dN(t) F t ; β,m } = Y (t)dλ(t Z; β,m ), where Λ( ) denotes the cumulative hazard function. Let dm i (t; β,m )=dn i (t) Y i (t)dλ(t Z i ; β,m ), i =1, 2,...,n, where {M i (t; β,m )} are martingales with respect to F t. Applying the inversion formula in Oakes & Dasu (1990) to the linear model, we know that the survival function of Z is S(t Z; β,m )= m { (0) + β T Z t } m (t)+β T Z exp du. 0 m (u)+β T Z As a result, {m (t)+β T Z}dΛ(t Z; β,m )=d{t+m (t)}. Thus, it is reasonable to estimate m 0 ( ) by the estimating equation, { m 0 (t)+β T Z i }dn i (t) 3 Y i (t)d{t + m 0 (t)} =0, (3)
4 which is in fact m 0 (t) i dn i(t)/ i Y i(t)+d m 0 (t) = i {βt Z i dn i (t) Y i (t)dt}/ i Y i(t). Then this first-order ordinary differential equation yields a closed-form solution of m 0 (t; β )=Ŝ(t) 1 i Ŝ(u) {Y i(u)du β T Z i dn i (u)} t i Y, i(u) where Ŝ(t) = exp{ t 0 i dn i(u)/ i Y i(u)}. Here τ = inf{t : pr(x > t)=0} < to avoid technical discussion about the tail behaviour of censored data. By the convention of Reid (1981) and James (1986), would be redefined as 1 for the last observation to ensure a meaningful m 0 ( ). Since i {m (t)+β T Z i}dn i (t) i Y i(t)d{t+m (t)} = i {m (t)+β T Z}dM i(t; β,m ), by subtracting it from (3), we obtain that i { m 0(t; β ) m (t)}dn i (t)+ i Y i(t)d{ m 0 (t; β ) m (t)} = i {m (t)+β T Z i }dm i (t; β,m ). Thus, this leads to m 0 (t; β ) m (t) = Ŝ(t) 1 i Ŝ(u) {m (u)+β T Z i }dm i (u; β,m ) t i Y, (4) i(u) and hence the consistency and asymptotic normality of m 0 (t; β ) hold according to the standard martingale theory for counting processes. In addition, if we let µ (t; β) = m 0 (t; β)/ β, it follows that µ (t; β ) = Ŝ(t) 1 τ Ŝ(u) t j Z jdn j (u)/ j Y j(u) = µ (t) +o p (1), where µ (t) =E {ZS (t Z)} /E {S (t Z)} and S (t Z) =pr{x t Z}, respectively. 2.2 Extended Buckley-James estimation The linear mean residual life model (2) implies that E(T Z) =m 0 (0)+β T Z. It is thus closely related to the standard linear model in Miller (1976) and Buckley & James (1979), which assumes that T γ 0 γ1 T Z has the identical zero-mean distribution with E(T Z) =γ 0 +γ1 T Z. Since the Buckley-James estimation procedure (Buckley & James, 1979) has been shown to be a reliable estimation procedure in the linear regression models with censored data (Miller & Halpern, 1982; Lin & Wei, 1992), we consider an extension of this procedure to the estimation of the linear mean residual life model. 4
5 Let α(t,z; β) =Z{T m 0 (0) β T Z}, and let that ξ(x,,z; β) = α(x, Z; β)+(1 )E{α(T,Z; β) T > X,Z}. Then E(ξ Z) =E(α Z) = 0, which leads to the unbiased estimating equations g 1 (β) =n 1/2 ξ(x i, i,z i ; β) =0, (5) as in the Buckley-James procedure for the standard linear regression model. However, both α(t,z; β) and E{α(T,Z; β) T > X,Z} depend on the unknown m 0 ( ), in addition to β. In the linear regression model, Buckley & James (1979) required that the residuals of T m 0 (0) β T Z share the same distribution, which was then estimated by the Kaplan-Meier product-limit estimator. However, this is usually not satisfied in the linear mean residual life model. It is thus less straightforward to use Efron s (1967) self-consistency representation of the Kaplan-Meier estimator to simplify g 1 (β). Instead, we extend the Buckley-James procedure to the model (2) using the proposed baseline estimator in 2.1. Consider the natural estimator of α(t,z; β) defined by α(t,z; β) =Z{T m 0 (0) β T Z}. In addition, by the inversion formula, S(t Z; β) is estimated by Ŝ(t Z; β) ={ m 0(0) + β T Z}{ m 0 (t) +β T Z} 1 exp[ t { m 0 0(u)+β T Z} 1 du]. Therefore, E{α(T,Z; β) T > X,Z} can be estimated by Ê{α(T,Z; β) T > X,Z} = Ŝ(X Z; β) 1 τ α(u, Z; β)dŝ(u Z; β). X By substituting α(t,z; β) and E{α(T,Z; β) T > X,Z} by these estimators in g 1 (β), we obtain the final estimating equations ĝ 1 (β) =n 1/2 [ i α(x i,z i ; β)+(1 i )Ê {α(t,z i; β) T>X i } ] =0. (6) Let D 1 (β) = lim n n 1/2 ĝ 1 (β)/ β. When β = β, D 1 (β) = E[ α(x, Z; β)/ β + (1 ) Ê{α(T,Z; β) T>X,Z}/ β]. In addition, by a Multivariate Central Limit Theorem, n 1/2 ĝ 1 (β ) approaches a zero-mean normal distribution with covariance matrix V 1 = E { α(x, Z; β ) 2 }+E[(1 )Ê {α(t,z; β ) T>X,Z} 2 ], asymptotically. Denote by β BJ the solution to ĝ 1 (β) = 0. Following an application of Taylor expansion, β BJ is asymptotically normal: n 1/2 ( β BJ β ) N(0,D 1 1 V 1D 1 1 ), 5
6 in distribution, as n, provided that D 1 is nonsingular. Here V 1 and D 1 can be estimated by their empirical estimators, V 1 = n 1 i [ i α(x i,z i ; β BJ ) 2 +(1 i )Ê{α(T,Z i ; β BJ ) T> X i,z i } 2 ] and D 1 = n 1 i [ i α(x i,z i ; β BJ )/ β+(1 i ) Ê{α(T,Z i; β BJ ) T>X i,z i }/ β]. β BJ Here, α( ) is an ad hoc choice for the estimating functions, and therefore the estimator is not necessarily efficient. In fact, in the original Buckley-James estimation, when ɛ is normal with known variance in the standard linear regression model of Y = γ 0 + γ T 1 Z + ɛ, α(t,z; γ) =Z(T γ 0 γ T 1 Z) is the score function of the regression parameter γ 1. Its efficient use for distributions other than normal is however unclear. Later, in James (1986), α(t,z; γ) was naturally generalised to be log f(t Z; γ)/ γ, where f( ) is the known density function of T. Similarly for the current linear mean residual life model (2), one possibly more efficient choice for α is [ log f(t Z; β) 1 T ] α(t,z; β) = = Z β m 0 (T )+β T Z {1+m 0(u)}du, (7) 0 {m 0 (u)+β T Z} 2 instead of α(t,z; β) =Z{T m 0 (0) β T Z}. In particular, when the baseline mean residual life function is constant of m 0 (t) α 0, then α(t,z; β) =Z(α 0 + βz) 2 (T α 0 βz) = Zm(t Z) 2 (T α 0 βz). More discussion of efficient estimation can be found in Quasi partial score estimation An alternative approach is to construct similar quasi partial score equations to those in Chen & Cheng (2005). Note that the estimating functions of g 2 (β ; m )=n 1/2 n β T Z i}dn i (t) Y i (t){1+m (t)dt}] are unbiased. It is therefore natural to replace m ( ) with 0 Z i[{m (t)+ its estimator, m 0 ( ; β), and use ĝ 2 {β; m 0 (β)} = 0 to solve for β to provide an estimator. Since d{t + m 0 (t; β)} = n { m 0(t; β)+β T Z i }dn i (t)/ i Y i(t), as in (3), ĝ 2 {β; m 0 (β)} =0 is simplified as { n 1/2 Zi Z(t) } { m 0 (t; β)+β T Z i }dn i (t) =0. (8) 0 Here, Z(t) = i Y i(t)z i / i Y i(t) has the same limit as µ (t). Therefore, n 1/2 ĝ 2 (β )/ β tends to D 2 = E[ {Z µ 0 (t)} 2 dn(t)] almost surely, as n. Consider a decomposition 6
7 of ĝ 2 {β ; m 0 (β )} of the form ĝ 2 {β ; m 0 (β )} =[ĝ 2 {β ; m 0 (β )} ĝ 2 (β; m )] + ĝ 2 (β ; m ), as in Tsiatis (1991). By the martingale representation of m(t; β ) m (t) in(4), the first term in this decomposition equals n 1/2 0 Z i Z(t) [ Ŝ(t) j=1 t Ŝ(u) ] i Y i(u) {m (u)+β T Z j }dm j (u) dn i (t). (9) Let Z(t) =Ŝ(t) t 0 Ŝ(u) 1 j {Z j Z(u)}dN j (u)/ j Y j(t). Then (9) can be further written as n 1/2 τ Z(t){m i 0 (t) +β TZ i}dm i (t). Since the second term of ĝ 2 (β ; m )in the decomposition equals n 1/2 τ i {Z 0 i Z(t)}{m (t) +β TZ i}dm i (t), we obtain that ĝ 2 (β ) = ĝ 2 {β ; m 0 (β )} = n 1/2 τ i {Z 0 i Z(t) Z(t)}{m (t) +β T Z i }dm i (t). As a result, ĝ 2 (β ) tends in distribution to a zero-mean normal distribution with the asymptotic variance V 2 (β )=E[ Y 0 1(t){Z µ (t) µ (t)} 2 {1+m (t)}{m (t)+β T Z}dt], where µ (t) = lim n Z(t), as n. Let βqp be the solution to ĝ 2 (β) = 0. By Taylor expansion, n 1/2 ( β QP β )={ n 1/2 ĝ 2 (β )/ β} 1 ĝ 2 (β )+o p (1). Therefore, n 1/2 ( β QP β ) N(0,D 1 2 V 2D 1 2 ), in distribution, where D 2 and V 2 can be consistently estimated by D 2 = n 1 n {Z 0 i Z(t)} 2 dn i (t), and V 2 = n 1 n Y 0 i(t){z i Z(t) Z(t)} 2 {1+ m 0(t)}{ m 0 (t)+ β QPZ T i }dt. In addition, user-defined weight functions can be also included in the proposed quasi partial estimating functions, as in n 1/2 τ i W (t){z 0 i Z(t)} 2 { m 0 (t; β)+β T Z i }dn i (t) =0, where the F t predictable W ( ) converges to a deterministic function of w( ). Then the asymptotic properties of its solution, β w, say, can be similarly derived. 2.4 Efficiency considerations Given the ad hoc nature of the estimating equations for the proposed β BJ and β QP, more efficient estimation procedures may be of interest. To gain some insight into the efficient estimation of the semiparametric linear mean residual life model, we follow the procedure in Lai & Ying (1992) and Lin & Ying (1994) in calculating a semiparametric information 7
8 bound for this model. Let β 0 R p be a parameter and letη( ) be a fixed function. Consider the parametric submodels of (2), defined by m(t Z) =m 0 (t)+β T s Z s (t), (10) where β s =(β0 T,β T ) T and Z s ( ) =(η( ) T,Z T ) T, respectively. Then the associated loglikelihood function of β s is l(β s )= i 0 log[{m 0(t) +θ T η (t)}/{m 0 (t) +βs T Z si (t)}]dn i (t) Y i (t)[{1 +m 0 (t) η (t)}/{m +θt 0 (t)+βs TZ si(t)}]dt. Therefore, the score functions at β s = (0,β T ) T are l β(β s ) = l(β s )/ β βs=β s = i Z 0 i/{m 0 (t) +β T Z i }[dn i (t) Y i (t){1 + m 0(t)}/{m 0 (t) +β T Z i }dt] and l β 0 (β s ) = l(β s )/ β 0 βs=β s = i 0 [η (t)/{1 +m 0(t)} η(t)/{m 0 (t)+β T Z i }][dn i (t) Y i (t){1+m 0(t)}/{m 0 (t)+β T Z i }dt], respectively. We therefore obtain that, in distribution, n 1/2 l β (β s ) l β 0 (β s ) N{0,I(η)}, where I(η) = I ββ(η) I β0 β(η) I ββ0 (η) I β0 β 0 (η) under proper regularity conditions similar to those in Lai & Ying (1992), where I ββ (η) = E{l ββ (β s )}, I ββ0 (η) =I β0 β(η) T = E{l ββ 0 (β s )} and I β0 β 0 (η) =E{l β 0 β 0 ) (β s )}. Actual calculation leads to that I ββ (η) = lim n n 1 I ββ0 (η) = lim n n 1 I β0 β 0 (η) = lim n n E [ Yi (t){1+m 0(t)}Z 2 ] i dt, {m 0 (t)+β T Z i } 3 [ Yi (t){1+m 0 E (t)}z { } i η T ] (t) {m 0 (t)+β T Z i } 2 1+m 0(t) η(t) dt, m 0 (t)+β T Z i [ { } ] Y i (t){1+m E 0(t)} η 2 (t) {m 0 (t)+β T Z i } 1+m 0(t) η(t) dt. m 0 (t)+β T Z i We shall say that matrix A is larger than matrix B if A B is nonnegative definite. Then, by an application of the Cauchy-Schwarz inequality, I 1 β (η) ={I ββ(η) I ββ0 (η)i 1 β 0 β 0 (η)i β0 β(η) T } 1 reaches its maximum at η = η 0. Here, η 0 ( ) satisfies the ordinary differential equation η 0(t)E[Y (t)/{1 +m 0(t)}] η 0 (t)e[y (t)/{m 0 (t)+β T Z}] =E[Y (t)z/{m 0 (t)+β T Z}], with the solution of η 0 (t) =P (t) 1 P (u)q(u)du. Here, P (t) = exp( t E[Y (u)/{m t 0 0(u) +, 8
9 β T Z}]/E[Y (u)/{1+m 0(u)}]du) and Q(t) =E[Y (t)z/{m 0 (t)+β T Z}]/E[Y (t)/{1+m 0(t)}], respectively. Therefore, the upper parametric information bound for β is I 1 β (η 0)= ( lim n n 1 0 [ Yi (t){z i µ 0 (t; β )} 2 ] 1 E dt), (11) {m 0 (t)+β T Z i } 2 where µ 0 (t; β) =lim n E[ i Y i(t)z i /{m 0 (t)+β T Z i }]/E[ i Y i(t)/{m 0 (t)+β T Z i }]. The sample-splitting technique in Lai & Ying and Lin & Ying (1994) can be used to construct a regular semiparametric estimator with its asymptotic variance reaching I 1 β (η 0). The entire data would be randomly divided into two evenly disjoint groups. One has the first n 0 =[n/2] samples and the other has the remaining n 1 = n n 0 samples. Let β j and m 0j ( ) be the ad hoc estimators of β and m 0 ( ) estimated from the (j + 1)st group for j =0, 1. For example, they can be obtained from either the extended Buckley-James estimation procedure or the quasi partial estimation procedure. Consider the following estimating functions for β: g(β) = n 1 j j=0 0 Z i Ẑ(t; β 1 j, m 0,1 j ) m 0,1 j (t)+ β T 1 j Z i [ dn i (t) Y ] i(t)d{t + m 0,1 j (t)}, m 0,1 j (t)+β T Z i where Ẑ(t; β,m 0)= i [Y i(t)z i /{m 0 (t) +β T Z i }]/ i [Y i(t)/{m 0 (t) +β T Z i }]. By standard counting process arguments, the asymptotic behaviour for one group is not affected by its estimates of β and m 0 from another group, but they asymptotically approach the true values. Thus the asymptotic variance of the semiparametric estimator β such that g( β) = 0 equals I 1 β (η 0). According to Bickel et al. (1993), the upper parametric information bound for β in the model (10) atβ s would lead to the semiparametric information bound for β at β in the model (2) as well. Therefore the semiparametric information bound of the linear mean residual life model is I 1 β (η 0); that is, by the Cramér-Rao inequality, the variance of any regular semiparametric estimator β in the linear mean residual life model is larger than (I ββ I βθ I 1 θθ I T βθ) 1 for any η, ifn 1/2 ( β β ) converges to a zero-mean normal. 9
10 3 Examples We consider the well-known Stanford heart transplantation data analysed in Miller & Halpern (1982). The time-to-event outcome of this dataset is the lifetime since first heart transplantation between October 1967 and February Two covariates were originally used in analysis, namely, age at the time of first transplant and the T5 mismatch score, which measures the degree of tissue incompatibility between the hearts of the initial donor and recipient with respect to hla antigens. To provide a comparison with their results, we apply the linear mean residual life model to the lifetimes with these two covariates. The results are tabulated in Table 1, along with those from the partial likelihood estimates of the Cox proportional hazards models and from the Buckley-James estimates of the linear regression models for the base-10 log-transformed lifetimes. [Table 1 about here] All the estimates in the linear mean residual life models suggest that both age and the T5 mismatch score are negatively associated with the patients life expectancy, although none of them is significant for the T5 mismatch scores. The covariate age is not only a significant predictor for the patients hazard and their average survival lifetimes, but also for their life expectancy throughout time. Among the different estimation procedures of the same linear mean residual life model, it is not surprising to see that the efficient estimation procedure yields the smallest variance. In addition, as demonstrated in Miller & Halpern (1982), a quadratic relationship to the covariate of age is plausible for both the Cox proportional hazards model and the linear regression model. We therefore fitted the linear mean residual life model with the T5 mismatch scores omitted because their nonsignificance as shown in Table 1. Table 2 shows that, both age and squared age are significant predictors of the life expectancy. However, their effects are of different signs: age itself is positively correlated with the life expectancy while the quadratic age has a negative influence on life expectancy. This demonstrates similar patterns to those shown by the Cox model and the linear regression 10
11 model. Again, the efficient estimation procedure yields smaller variances than do the other ad hoc procedures. [Table 2 about here] We can assess the fit of the linear mean residual life model using the goodness-of-fit test of Gill & Schumacher (1987). We consider two different weight functions, W 1 (t) and W 2 (t), say, in the weighted quasi partial estimating equations of g w1 (β) and g w2 (β), respectively. Their corresponding estimators are denoted by β w1 and β w2, respectively. Then a significant difference of β w1 β w2 should suggest lack of fit. To implement this, we adapt the reliable variation suggested in Wei et al. (1990), involving the goodness-of-fit test statistic defined by κ(w 1,W 2 ) = min β N ( βw1 ) {(g w1(β) T,g w2 (β) T )Σ 1 w ( β w1 )(g w1 (β) T,g w2 (β) T ) T }, where N ( β w1 ) is a neighbourhood of β w1. Here Σ w ( β w1 ) is the covariance matrix of (g w1 ( β w1),g T w2 ( β w1) T T ) T, which can be estimated by computer-intensive methods, such as the bootstrap. Under the null hypothesis that the linear mean residual life model is adequate, κ(w 1,W 2 ) is approximately distributed as χ 2 p. For the Stanford data, κ(w 1,W 2 )=4.81 (p =0.09) for the linear mean residual life model based on age and T5, which suggests a mild lack of fit. For the linear mean residual life model based on age and squared age, the test statistic is 1.17 (p >0.5) which does not reject the model. 4 Discussion Some simulations were conducted to assess the finite-sample properties of the proposed methods. We took the sample size n to be 100 or 200 with two covariates Z =(Z 1,Z 2 ) T for each of n subjects. The covariate Z 1 is a Ber(0.5) random variable and Z 2 Un[0, 1]. We chose m (t) =t + 1, which corresponds to the Pareto distribution with the baseline survival function being (1 + t) 2. Failure times were generated according to the linear mean residual life model. The true parameters were (0, 0) T or (1, 1) T. Independent censoring times were generated from the uniform distribution so as to allow about 30% of the observations to be 11
12 censored. The simulation results are summarised in Table 3, with each entry based on 1000 simulated datasets. The results show that the estimators of β are virtually unbiased and the nominal 95% confidence intervals have acceptable coverage probabilities. More extensive simulation studies are needed to evaluate the efficiency of the proposed methods. [Table 3 about here] There is, however, the critical challenge of the tail behaviour of underlying distributions of failure times. When the underlying failure times are heavily right-skewed and censored early, as in the case of long-term survivors on cancer treatment or subjects in HIV/AIDS prevention/vaccine trials, it is impossible to estimate the mean residual life function on the whole positive real line without extra assumptions, although some techniques such those as in Gill (1983) can be extended. In general, it is difficult to determine a reasonable upper limit τ without the robustness of the proposed methodologies being compromised. One possible approach is modify the fully unspecified m 0 ( ) by including a parametric component in the tail. For instance, let τ be a prespecified truncation time. Then it is assumed that m 0 (t) =m 0 (t)i(t < τ)+m r I(t > τ), where m r is some positive constant. This means that the baseline mean residual life function is unspecified up to the truncation time τ, while it becomes exponential after τ. Then it is straightforward to extend the proposed methodologies to the whole positive real line. Acknowledgement The authors would like to thank the associate editor and an anonymous reviewer for their constructive comments and suggestions that led to a significant improvement of this paper. They would also thank Professors Zhiliang Ying, Chen-Hsin Chen and Louis Chen for fruitful discussion during the first author s visit to the National University of Singapore in , which was partially funded by the Institute for Mathematical Sciences, National University of Singapore. 12
13 References Buckley, J. & James, I. (1979). Linear regression with censored data. Biometrika 66, Bickel, P. J., Klassen, C. A. J., Ritov, Y. & Wellner, J. A. (1993). Efficient and Adaptive Estimation for Semiparametric Models. Baltimore, MD: Johns Hopkins University Press. Chen, Y. Q. & Cheng, S. (2005). Semiparametric regression analysis of mean residual life with censored survival data. Biometrika 92, Efron, B. (1967). The two sample problem with censored data. Proc. 5th Berkeley Symp. Math. Stat. Prob. 4, Berkeley, CA: UC Berkeley Press. Gill, R. D. (1983). Large sample behaviour of the product-limit estimator on the whole line. Ann. Statist. 11, Gill, R. D.& Schumacher, M. (1987). assumption. Biometrika 74, A simple test of the proportional hazards James, I. R. (1986). On estimating equations with censored data. Biometrika 73, Lai, T. L. & Ying, Z. (1992). Asymptotic efficiency estimation in censored and truncated regression models. Statist. Sinica 2, Lin, D. Y. & Ying, Z. (1994). Biometrika 81, Semiparametric analysis of the additive risk model. Lin, J. S. & Wei, L.-J. (1992). Linear regression analysis based on Buckley-James estimating equations. Biometrics 48, Maguluri, G. & Zhang, C.-H. (1994). Estimation in the mean residual life regression model. J. R. Statist. Soc. B 56,
14 Miller, R. & Halpern, J. (1982). Regression with censored data. Biometrika 69, Oakes, D. & Dasu, T. (1990). A note on residual life. Biometrika 77, Oakes, D. & Dasu, T. (2003). Inference for the proportional mean residual life model. In Crossing Boundaries: Statistical Essays in Honor of Jack Hall, Institute of Mathematical Statistics Lecture Notes Monograph Series, Vol. 43, Ed. J. E. Kolassa and D. Oakes, pp Hayward, CA: Institute of Mathematical Statistics. Reid, N. (1981). Influence functions for censored data. Ann. Statist. 9, Tsiatis, A. A. (1981). A large sample study of Cox s regression model. Ann. Statist. 9, Wei, L. J., Ying, Z. & Lin, D. Y. (1990). Linear regression analysis of censored survival data based on rank tests. Biometrika 77,
15 Table 1: Regression estimates with standard errors for 157 Stanford heart transplant patients Coefficients Age T5 Score β 1 s.e. β2 s.e. Cox model Linear regression model Linear mean residual life model Buckley-James Quasi partial score Efficient estimation Cox model, estimates obtained in original time-scale; Linear regression model, Buckley- James estimates obtained in 10-base log-transformed scale; linear mean residual life model, estimates obtained in original time-scale; s.e., estimated standard error 15
16 Table 2: Regression estimates with standard errors for 152 Stanford heart transplant patients who survived at least 10 days Coefficients Age Age 2 β 1 s.e. β2 s.e. Cox model Linear regression model Linear mean residual life model Buckley-James Quasi partial score Efficient estimation Cox model, estimates obtained in original time-scale; Linear regression model, Buckley- James estimates obtained in 10-base log-transformed scale; linear mean residual life model, estimates obtained in original time-scale; s.e., estimated standard error 16
17 Table 3: Summary of simulation studies Z 1 Z 2 β n Estimation Cov. Prob. Bias Cov. Prob. Bias (0,0) 100 BJ (0,0) 100 QP (0,0) 200 BJ (0,0) 200 QP (1,1) 100 BJ (1,1) 100 QP (1,1) 200 BJ (1,1) 200 QP BJ, the Buckley-James approach; QP, the quasi partial score equation approach; Cov. Prob., coverage probabilities of the 95% nominal confidence intervals; Bias, mean difference between the estimate and the true value 17
University 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 T = β T Z + ɛ Zi Z(u; β) } dn i (ue βzi ) = 0,
Accelerated failure time model: log T = β T Z + ɛ β estimation: solve where S n ( β) = n i=1 { Zi Z(u; β) } dn i (ue βzi ) = 0, Z(u; β) = j Z j Y j (ue βz j) j Y j (ue βz j) How do we show the asymptotics
More informationAttributable Risk Function in the Proportional Hazards Model
UW Biostatistics Working Paper Series 5-31-2005 Attributable Risk Function in the Proportional Hazards Model Ying Qing Chen Fred Hutchinson Cancer Research Center, yqchen@u.washington.edu Chengcheng Hu
More 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 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 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 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 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 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 informationAnalysis of transformation models with censored data
Biometrika (1995), 82,4, pp. 835-45 Printed in Great Britain Analysis of transformation models with censored data BY S. C. CHENG Department of Biomathematics, M. D. Anderson Cancer Center, University of
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 informationA note on L convergence of Neumann series approximation in missing data problems
A note on L convergence of Neumann series approximation in missing data problems Hua Yun Chen Division of Epidemiology & Biostatistics School of Public Health University of Illinois at Chicago 1603 West
More informationA COMPARISON OF POISSON AND BINOMIAL EMPIRICAL LIKELIHOOD Mai Zhou and Hui Fang University of Kentucky
A COMPARISON OF POISSON AND BINOMIAL EMPIRICAL LIKELIHOOD Mai Zhou and Hui Fang University of Kentucky Empirical likelihood with right censored data were studied by Thomas and Grunkmier (1975), Li (1995),
More informationSTAT 331. Accelerated Failure Time Models. Previously, we have focused on multiplicative intensity models, where
STAT 331 Accelerated Failure Time Models Previously, we have focused on multiplicative intensity models, where h t z) = h 0 t) g z). These can also be expressed as H t z) = H 0 t) g z) or S t z) = e Ht
More informationPublished online: 10 Apr 2012.
This article was downloaded by: Columbia University] On: 23 March 215, At: 12:7 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 172954 Registered office: Mortimer
More informationLinear rank statistics
Linear rank statistics Comparison of two groups. Consider the failure time T ij of j-th subject in the i-th group for i = 1 or ; the first group is often called control, and the second treatment. Let n
More informationAccelerated Failure Time Models: A Review
International Journal of Performability Engineering, Vol. 10, No. 01, 2014, pp.23-29. RAMS Consultants Printed in India Accelerated Failure Time Models: A Review JEAN-FRANÇOIS DUPUY * IRMAR/INSA of Rennes,
More informationEfficiency of Profile/Partial Likelihood in the Cox Model
Efficiency of Profile/Partial Likelihood in the Cox Model Yuichi Hirose School of Mathematics, Statistics and Operations Research, Victoria University of Wellington, New Zealand Summary. This paper shows
More informationQuantile Regression for Residual Life and Empirical Likelihood
Quantile Regression for Residual Life and Empirical Likelihood Mai Zhou email: mai@ms.uky.edu Department of Statistics, University of Kentucky, Lexington, KY 40506-0027, USA Jong-Hyeon Jeong email: jeong@nsabp.pitt.edu
More informationEstimation of Conditional Kendall s Tau for Bivariate Interval Censored Data
Communications for Statistical Applications and Methods 2015, Vol. 22, No. 6, 599 604 DOI: http://dx.doi.org/10.5351/csam.2015.22.6.599 Print ISSN 2287-7843 / Online ISSN 2383-4757 Estimation of Conditional
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 informationImproving Efficiency of Inferences in Randomized Clinical Trials Using Auxiliary Covariates
Improving Efficiency of Inferences in Randomized Clinical Trials Using Auxiliary Covariates Anastasios (Butch) Tsiatis Department of Statistics North Carolina State University http://www.stat.ncsu.edu/
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 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 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 informationLeast Absolute Deviations Estimation for the Accelerated Failure Time Model. University of Iowa. *
Least Absolute Deviations Estimation for the Accelerated Failure Time Model Jian Huang 1,2, Shuangge Ma 3, and Huiliang Xie 1 1 Department of Statistics and Actuarial Science, and 2 Program in Public Health
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 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 informationEfficiency Comparison Between Mean and Log-rank Tests for. Recurrent Event Time Data
Efficiency Comparison Between Mean and Log-rank Tests for Recurrent Event Time Data Wenbin Lu Department of Statistics, North Carolina State University, Raleigh, NC 27695 Email: lu@stat.ncsu.edu Summary.
More informationAFT Models and Empirical Likelihood
AFT Models and Empirical Likelihood Mai Zhou Department of Statistics, University of Kentucky Collaborators: Gang Li (UCLA); A. Bathke; M. Kim (Kentucky) Accelerated Failure Time (AFT) models: Y = log(t
More informationFull likelihood inferences in the Cox model: an empirical likelihood approach
Ann Inst Stat Math 2011) 63:1005 1018 DOI 10.1007/s10463-010-0272-y Full likelihood inferences in the Cox model: an empirical likelihood approach Jian-Jian Ren Mai Zhou Received: 22 September 2008 / Revised:
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 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 Accelerated Failure Time Model Under Biased. Sampling
The Accelerated Failure Time Model Under Biased Sampling Micha Mandel and Ya akov Ritov Department of Statistics, The Hebrew University of Jerusalem, Israel July 13, 2009 Abstract Chen (2009, Biometrics)
More informationEfficient Semiparametric Estimators via Modified Profile Likelihood in Frailty & Accelerated-Failure Models
NIH Talk, September 03 Efficient Semiparametric Estimators via Modified Profile Likelihood in Frailty & Accelerated-Failure Models Eric Slud, Math Dept, Univ of Maryland Ongoing joint project with Ilia
More informationOn the Breslow estimator
Lifetime Data Anal (27) 13:471 48 DOI 1.17/s1985-7-948-y On the Breslow estimator D. Y. Lin Received: 5 April 27 / Accepted: 16 July 27 / Published online: 2 September 27 Springer Science+Business Media,
More informationMAXIMUM LIKELIHOOD METHOD FOR LINEAR TRANSFORMATION MODELS WITH COHORT SAMPLING DATA
Statistica Sinica 25 (215), 1231-1248 doi:http://dx.doi.org/1.575/ss.211.194 MAXIMUM LIKELIHOOD METHOD FOR LINEAR TRANSFORMATION MODELS WITH COHORT SAMPLING DATA Yuan Yao Hong Kong Baptist University Abstract:
More informationOn least-squares regression with censored data
Biometrika (2006), 93, 1, pp. 147 161 2006 Biometrika Trust Printed in Great Britain On least-squares regression with censored data BY ZHEZHEN JIN Department of Biostatistics, Columbia University, New
More informationFrom semi- to non-parametric inference in general time scale models
From semi- to non-parametric inference in general time scale models Thierry DUCHESNE duchesne@matulavalca Département de mathématiques et de statistique Université Laval Québec, Québec, Canada Research
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 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 informationProduct-limit estimators of the survival function with left or right censored data
Product-limit estimators of the survival function with left or right censored data 1 CREST-ENSAI Campus de Ker-Lann Rue Blaise Pascal - BP 37203 35172 Bruz cedex, France (e-mail: patilea@ensai.fr) 2 Institut
More informationLSS: An S-Plus/R program for the accelerated failure time model to right censored data based on least-squares principle
computer methods and programs in biomedicine 86 (2007) 45 50 journal homepage: www.intl.elsevierhealth.com/journals/cmpb LSS: An S-Plus/R program for the accelerated failure time model to right censored
More informationSTAT331. Combining Martingales, Stochastic Integrals, and Applications to Logrank Test & Cox s Model
STAT331 Combining Martingales, Stochastic Integrals, and Applications to Logrank Test & Cox s Model Because of Theorem 2.5.1 in Fleming and Harrington, see Unit 11: For counting process martingales with
More 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 informationFULL LIKELIHOOD INFERENCES IN THE COX MODEL: AN EMPIRICAL LIKELIHOOD APPROACH
FULL LIKELIHOOD INFERENCES IN THE COX MODEL: AN EMPIRICAL LIKELIHOOD APPROACH Jian-Jian Ren 1 and Mai Zhou 2 University of Central Florida and University of Kentucky Abstract: For the regression parameter
More 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 informationIntroduction to Empirical Processes and Semiparametric Inference Lecture 25: Semiparametric Models
Introduction to Empirical Processes and Semiparametric Inference Lecture 25: Semiparametric Models Michael R. Kosorok, Ph.D. Professor and Chair of Biostatistics Professor of Statistics and Operations
More informationSurvival Analysis 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 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 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 informationChapter 4 Fall Notations: t 1 < t 2 < < t D, D unique death times. d j = # deaths at t j = n. Y j = # at risk /alive at t j = n
Bios 323: Applied Survival Analysis Qingxia (Cindy) Chen Chapter 4 Fall 2012 4.2 Estimators of the survival and cumulative hazard functions for RC data Suppose X is a continuous random failure time with
More 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 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 informationLecture 3. Truncation, length-bias and prevalence sampling
Lecture 3. Truncation, length-bias and prevalence sampling 3.1 Prevalent sampling Statistical techniques for truncated data have been integrated into survival analysis in last two decades. Truncation in
More informationPairwise rank based likelihood for estimating the relationship between two homogeneous populations and their mixture proportion
Pairwise rank based likelihood for estimating the relationship between two homogeneous populations and their mixture proportion Glenn Heller and Jing Qin Department of Epidemiology and Biostatistics Memorial
More informationand Comparison with NPMLE
NONPARAMETRIC BAYES ESTIMATOR OF SURVIVAL FUNCTIONS FOR DOUBLY/INTERVAL CENSORED DATA and Comparison with NPMLE Mai Zhou Department of Statistics, University of Kentucky, Lexington, KY 40506 USA http://ms.uky.edu/
More informationEstimation and Inference of Quantile Regression. for Survival Data under Biased Sampling
Estimation and Inference of Quantile Regression for Survival Data under Biased Sampling Supplementary Materials: Proofs of the Main Results S1 Verification of the weight function v i (t) for the lengthbiased
More informationA Recursive Formula for the Kaplan-Meier Estimator with Mean Constraints and Its Application to Empirical Likelihood
Noname manuscript No. (will be inserted by the editor) A Recursive Formula for the Kaplan-Meier Estimator with Mean Constraints and Its Application to Empirical Likelihood Mai Zhou Yifan Yang Received:
More informationCompeting risks data analysis under the accelerated failure time model with missing cause of failure
Ann Inst Stat Math 2016 68:855 876 DOI 10.1007/s10463-015-0516-y Competing risks data analysis under the accelerated failure time model with missing cause of failure Ming Zheng Renxin Lin Wen Yu Received:
More informationMultivariate Survival Data With Censoring.
1 Multivariate Survival Data With Censoring. Shulamith Gross and Catherine Huber-Carol Baruch College of the City University of New York, Dept of Statistics and CIS, Box 11-220, 1 Baruch way, 10010 NY.
More information1 Local Asymptotic Normality of Ranks and Covariates in Transformation Models
Draft: February 17, 1998 1 Local Asymptotic Normality of Ranks and Covariates in Transformation Models P.J. Bickel 1 and Y. Ritov 2 1.1 Introduction Le Cam and Yang (1988) addressed broadly the following
More informationStatistical Analysis of Competing Risks With Missing Causes of Failure
Proceedings 59th ISI World Statistics Congress, 25-3 August 213, Hong Kong (Session STS9) p.1223 Statistical Analysis of Competing Risks With Missing Causes of Failure Isha Dewan 1,3 and Uttara V. Naik-Nimbalkar
More informationPart III. Hypothesis Testing. III.1. Log-rank Test for Right-censored Failure Time Data
1 Part III. Hypothesis Testing III.1. Log-rank Test for Right-censored Failure Time Data Consider a survival study consisting of n independent subjects from p different populations with survival functions
More informationPoint and Interval Estimation for Gaussian Distribution, Based on Progressively Type-II Censored Samples
90 IEEE TRANSACTIONS ON RELIABILITY, VOL. 52, NO. 1, MARCH 2003 Point and Interval Estimation for Gaussian Distribution, Based on Progressively Type-II Censored Samples N. Balakrishnan, N. Kannan, C. T.
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 informationLikelihood ratio confidence bands in nonparametric regression with censored data
Likelihood ratio confidence bands in nonparametric regression with censored data Gang Li University of California at Los Angeles Department of Biostatistics Ingrid Van Keilegom Eindhoven University of
More 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 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 informationRank Regression Analysis of Multivariate Failure Time Data Based on Marginal Linear Models
doi: 10.1111/j.1467-9469.2005.00487.x Published by Blacwell Publishing Ltd, 9600 Garsington Road, Oxford OX4 2DQ, UK and 350 Main Street, Malden, MA 02148, USA Vol 33: 1 23, 2006 Ran Regression Analysis
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 informationTESTS FOR LOCATION WITH K SAMPLES UNDER THE KOZIOL-GREEN MODEL OF RANDOM CENSORSHIP Key Words: Ke Wu Department of Mathematics University of Mississip
TESTS FOR LOCATION WITH K SAMPLES UNDER THE KOIOL-GREEN MODEL OF RANDOM CENSORSHIP Key Words: Ke Wu Department of Mathematics University of Mississippi University, MS38677 K-sample location test, Koziol-Green
More 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 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 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 informationSurvival Prediction Under Dependent Censoring: A Copula-based Approach
Survival Prediction Under Dependent Censoring: A Copula-based Approach Yi-Hau Chen Institute of Statistical Science, Academia Sinica 2013 AMMS, National Sun Yat-Sen University December 7 2013 Joint work
More informationLEAST ABSOLUTE DEVIATIONS ESTIMATION FOR THE ACCELERATED FAILURE TIME MODEL
Statistica Sinica 17(2007), 1533-1548 LEAST ABSOLUTE DEVIATIONS ESTIMATION FOR THE ACCELERATED FAILURE TIME MODEL Jian Huang 1, Shuangge Ma 2 and Huiliang Xie 1 1 University of Iowa and 2 Yale University
More informationFrailty Models and Copulas: Similarities and Differences
Frailty Models and Copulas: Similarities and Differences KLARA GOETHALS, PAUL JANSSEN & LUC DUCHATEAU Department of Physiology and Biometrics, Ghent University, Belgium; Center for Statistics, Hasselt
More 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 informationGOODNESS-OF-FIT TESTS FOR ARCHIMEDEAN COPULA MODELS
Statistica Sinica 20 (2010), 441-453 GOODNESS-OF-FIT TESTS FOR ARCHIMEDEAN COPULA MODELS Antai Wang Georgetown University Medical Center Abstract: In this paper, we propose two tests for parametric models
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 informationGoodness-of-Fit Tests With Right-Censored Data by Edsel A. Pe~na Department of Statistics University of South Carolina Colloquium Talk August 31, 2 Research supported by an NIH Grant 1 1. Practical Problem
More informationPanel Count Data Regression with Informative Observation Times
UW Biostatistics Working Paper Series 3-16-2010 Panel Count Data Regression with Informative Observation Times Petra Buzkova University of Washington, buzkova@u.washington.edu Suggested Citation Buzkova,
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 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 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 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 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 informationBayesian Nonparametric Inference Methods for Mean Residual Life Functions
Bayesian Nonparametric Inference Methods for Mean Residual Life Functions Valerie Poynor Department of Applied Mathematics and Statistics, University of California, Santa Cruz April 28, 212 1/3 Outline
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 informationSemiparametric Mixed Effects Models with Flexible Random Effects Distribution
Semiparametric Mixed Effects Models with Flexible Random Effects Distribution Marie Davidian North Carolina State University davidian@stat.ncsu.edu www.stat.ncsu.edu/ davidian Joint work with A. Tsiatis,
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 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 informationOn consistency of Kendall s tau under censoring
Biometria (28), 95, 4,pp. 997 11 C 28 Biometria Trust Printed in Great Britain doi: 1.193/biomet/asn37 Advance Access publication 17 September 28 On consistency of Kendall s tau under censoring BY DAVID
More informationIntroduction to Empirical Processes and Semiparametric Inference Lecture 01: Introduction and Overview
Introduction to Empirical Processes and Semiparametric Inference Lecture 01: Introduction and Overview Michael R. Kosorok, Ph.D. Professor and Chair of Biostatistics Professor of Statistics and Operations
More informationSome New Methods for Latent Variable Models and Survival Analysis. Latent-Model Robustness in Structural Measurement Error Models.
Some New Methods for Latent Variable Models and Survival Analysis Marie Davidian Department of Statistics North Carolina State University 1. Introduction Outline 3. Empirically checking latent-model robustness
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 informationSEMIPARAMETRIC METHODS FOR ESTIMATING CUMULATIVE TREATMENT EFFECTS IN THE PRESENCE OF NON-PROPORTIONAL HAZARDS AND DEPENDENT CENSORING
SEMIPARAMETRIC METHODS FOR ESTIMATING CUMULATIVE TREATMENT EFFECTS IN THE PRESENCE OF NON-PROPORTIONAL HAZARDS AND DEPENDENT CENSORING by Guanghui Wei A dissertation submitted in partial fulfillment of
More informationPOWER AND SAMPLE SIZE DETERMINATIONS IN DYNAMIC RISK PREDICTION. by Zhaowen Sun M.S., University of Pittsburgh, 2012
POWER AND SAMPLE SIZE DETERMINATIONS IN DYNAMIC RISK PREDICTION by Zhaowen Sun M.S., University of Pittsburgh, 2012 B.S.N., Wuhan University, China, 2010 Submitted to the Graduate Faculty of the Graduate
More informationA Recursive Formula for the Kaplan-Meier Estimator with Mean Constraints
Noname manuscript No. (will be inserted by the editor) A Recursive Formula for the Kaplan-Meier Estimator with Mean Constraints Mai Zhou Yifan Yang Received: date / Accepted: date Abstract In this note
More information