An augmented inverse probability weighted survival function estimator
|
|
- Penelope Alberta Quinn
- 5 years ago
- Views:
Transcription
1 An augmented inverse probability weighted survival function estimator Sundarraman Subramanian & Dipankar Bandyopadhyay Abstract We analyze an augmented inverse probability of non-missingness weighted estimator of a survival function for a missing censoring indicator model, in the absence and presence of left truncation. The estimator improves upon its precursor but is still not the best in terms of achieving minimal asymptotic variance. KEY WORDS: Asymptotic variance; Dikta semiparametric estimator; Empirical estimator; Gaussian process; Gradient vector; Missing at random. 1 Introduction In this article, we focus on a missing censoring indicator MCI) model, with observed data {X i, ξ i, σ i ) 1 i n }, where X i = mint i, C i ), T is a lifetime of interest, C is an independent censoring variable, ξ i = 1 when δ i = IT i C i ) is observed and is otherwise, and σ i = ξ i δ i. The censoring indicators δ i are assumed to be missing at random MAR), which implies that P ξ = 1 X = x, δ = d) = P ξ = 1 X = x) = πx) van der Laan and McKeague, 1998). Writing px) = P δ = 1 X = x), MAR means that ξ and δ are conditionally independent given X: P σ = 1 X = x) = πx)px). We investigate semiparametric estimation of St), the survival function of T, using an augmented inverse probability weighted AIPW) scheme. Center for Applied Mathematics and Statistics, Department of Mathematical Sciences, New Jersey Institute of Technology, Newark, New Jersey 712, United States Department of Biostatistics, Bioinformatics and Epidemiology, Medical University of South Carolina, Charleston, SC 29425, United States 1
2 Inverse probability weighting IPW) has gained much currency of late. Horvitz and Thompson 1952) introduced IPW in the context of sample survey, while Koul, Susarla, and van Ryzin 1981) first employed the weighting scheme for censored regression. More recently, Robins and Rotnitzky 1992), Satten and Datta 21, 24), and Rotnitzky and Robins 25), among others, have also introduced IPW-based analysis in survival and related settings. However, the IPW estimators are not consistent under misspecification of the weighting probability, leading to the investigation of AIPW estimators Scharfstein, Rotnitzky and Robins 1999). The key element of such AIPW estimators is that they are doubly robust in the sense that they remain consistent if any one of px) and πx) is specified correctly. Robins, Rotnitzky and van der Laan 2), among others, have studied double robust estimators. See also Tsiatis 26) and van der Laan and Robins 23). 2 Estimation We first introduce some preliminary notation and present a brief review of two semiparametric estimators of St) when there are MCIs. These are the Dikta estimator Subramanian, 24a) and the IPW estimator. The latter was investigated recently for a left truncated version of the MCI model Subramanian and Bandyopadhyay, 28). 2.1 Preliminaries Let Ht) denote the distribution function of X and Ĥt) denote its empirical estimator. We focus on estimation over [, τ], where Hτ) < 1. Let Qt) = P X t, δ = 1) denote the subdistribution function corresponding to uncensored failures. Its empirical estimator, however, is inapplicable when there are MCIs. To address this, specify a model for pt). 2
3 That is, take pt) = pt, θ), where θ is a k dimensional unknown Euclidean parameter. The Dikta semiparametric estimator of Qt) is given by ˆQ D t) = ˆps)dĤs), 1) where ˆps) = ps, ˆθ), and ˆθ is the maximum likelihood estimate MLE) of θ. The MLE ˆθ can be obtained by maximizing a likelihood based on the complete cases ξ i = 1, see, for example, Subramanian 24a). Let θ denote the true value of θ and define p s) = ps, θ ). For r = 1,..., k, write p r s, θ ) = ps, θ)/ θ r θ=θ and let P s) = [p 1 s, θ ),..., p k s, θ )] T denote the gradient vector evaluated at θ. Let I = E[πX)P X)P T X)/p X)1 p X)))]. Also, let αs, t) = P T s)i 1 P t), and define Bt) 1 n B i t) = 1 n ξ i δ i p X i )) p X i )1 p X i )) αt, X i). 2) Let M be such that Ms)dHs) <. From Lemma 3.5 of Dikta 1998), we have the following representation, assuming that the model for pt) is correctly specified: n 1/2 {ˆpt) p t)} = n 1/2 Bt) + Rn t), 3) where, uniformly for t τ, the term n 1/2 R n t) is bounded above by Mt) ˆθ r θ,r )ˆθ s θ,s ) = Mt) O p n 1/2 ) + o p n 1/2 ) ) 2 = Mt)Op n 1 ). 1 r,s k Under the regularity conditions given by Dikta 1998), which we denote henceforth by R, it can be shown that Subramanian, 24a) n 1/2 ˆQ D t) Qt)) is asymptotically linear with influence function given by At) + Bt), where At) = p X)IX t) Qt), 4) Bt) = Bs)dHs) = ξ δ p X)) p X)1 p X)) αs, X)dHs). 5) 3
4 It follows that the process n 1/2 ˆQ D t) Qt)) converges weakly to a zero-mean Gaussian process whose variance at t [the second moment of At)+Bt)], the benchmark for comparing competing estimators of Qt) in this article, is given by V D t) = p s)dqs) Q 2 t) + αu, v)dhu)dhv). 6) Alternatively, specify πt) = πt, γ) parametrically, where γ is an unknown m dimensional parameter. Set ˆπs) = πs, ˆγ) and ˆγ is the MLE of γ. Let Ŵ t) denote the empirical estimator of W t) = P X t, ξ = 1, σ = 1). From an IPW representation for Qt) e.g., Subramanian, 26), the IPW estimator takes the form ˆQ I t) = 1 dŵ s). 7) ˆπs) Denote the gradient vector for this scenario by P s, θ), and by P s) when it is evaluated at the true value γ. Write βs, t) = P T s)j 1 E[ P X) P T X)/π X)1 π X)))]. Define P t), where the information matrix J = V t) 1 n V i t) = 1 n ξ i π X i )) π X i )1 π X i ) βt, X i). 8) Analogous to Eq. 3), assuming that the model for πt) is correctly specified, we have that n 1/2 ˆπt) π t)) = n 1/2 V t) + Rn t). 9) Assuming that πt) is bounded away from, uniformly for t [, τ], Eq. 8) and Eq. 9) and some basic analysis imply that n 1/2 ˆQ I t) Qt)) is asymptotically linear with influence function given by Ăt) + Bt), where IX t, ξ = 1, σ = 1) Ăt) = Qt), 1) π X) V s) Bt) = π s) dqs) = ξ π X)) t βs, X) dqs). 11) π X)1 π X)) π s) 4
5 It follows that the process n 1/2 ˆQ I t) Qt)) converges weakly to a zero-mean Gaussian process whose variance at t [the second moment of Ăt) + Bt)] is given by V I t) = 1 π s) dqs) Q2 t) βu, v) dqu)dqv). 12) π u)π v) Furthermore, it can be shown that V I t) V D t), see, for example, Subramanian and Bandyopadhyay 28), where an analogous result is derived for a left truncated version. Both estimators ˆQ D t) and ˆQ I t) work under the tacit assumption that the parametric models for pt) and πt) respectively are correctly specified, and are inconsistent when the putative models are ill-specified. The AIPW estimator of Qt) that we analyze in the next subsection is a double robust estimator Scharfstein et. al., 1999, Tsiatis, 26), in that it retains consistency when either or both of pt) and πt) is specified correctly. 2.2 AIPW subdistribution function estimator The AIPW estimator of Qt) takes the form e.g., Tsiatis, 26) ˆQ AI t) = 1 n { ξi δ i ˆπX i ) + 1 ξ ) } i ˆpX i ) IX i t). 13) ˆπX i ) We analyze ˆQ AI t) for three different scenarios. Scenario 1 represents the case only pt) correctly specified, Scenario 2 the case only πt) correctly specified, and Scenario 3 represents both pt) and πt) correctly specified. Clearly, ˆQAI t) = ˆQ D t) + T n t), where T n t) = 1 n { } ξi ˆπX i ) δ i ˆpX i )) IX i t), 14) and ˆQ D t) is defined by Eq. 1). We analyze this version of ˆQ AI t) for Scenario 1. On the other hand, we can also write ˆQ AI t) = ˆQ I t) + T n t), where ˆQ I t) is defined by Eq. 7) and T n t) = 1 n 1 ξ ) i ˆpX i )IX i t), 15) ˆπX i ) 5
6 which is the version that we will analyze for Scenario 2. We now state some basic assumptions, introduced by Kouassi and Singh 1997). We denote them henceforth by RP. When the model for pt) is misspecified Scenario 2), we assume that there exists θ IR k such that, with p t) = pt, θ ), the following hold: E ˆθ θ ) = o n 1/2), 16) ˆθ θ = O p n 1/2 ), 17) p t) p t) <. 18) Eq. 16) and Eq. 17) are introduced to ensure that the MLEs ˆθ and ˆpt) converge at the standard n 1/2 rates even when the parametric model is misspecified. Eq. 18) requires that the bias of ˆpt) is bounded when pt, θ) does not hold. Analogous conditions for the MLEs ˆγ, ˆπt), and π t) = πt, γ ) are also assumed denoted by RPI henceforth). Also, it will be assumed that ˆpt) P p t) and ˆπt) P π t) uniformly for t [, τ]. The last conditions of uniformity which are also included in RP or RPI) can be enforced through appropriate boundedness assumptions, see Eq. 3) and the equation following it. For Scenario 1, ˆpt) P p t) and ˆQ D t) P Qt), uniformly for t [, τ], see Eq. 3) or Dikta 1998). We can write [cf. Eq. 14)] T n t) = T n,1 t) + T n,2 t) + T n,3 t) where { ξi T n,1 t) = 1 } n π X i ) δ i p X i )) IX i t) 19) T n,2 t) = 1 { } ξi n π X i ) ˆpX i) p X i )) IX i t) 2) T n,3 t) = 1 { } ξi ˆπX i ) π X i )) δ n ˆπX i )π i ˆpX i )) IX i t). 21) X i ) By the strong law large numbers T n,1 t) a.s.. Next, T n,2 t) = O p sup t τ ˆpt) p t) ), 6
7 and therefore T n,2 t) P. Also, T n,3 t) = O p sup t τ ˆπt) π t) ) P. Therefore, when only pt) is specified correctly, it is clear that ˆQ AI t) = ˆQ D t)+t n t) retains consistency. ) ) Note, however, that n 1/2 ˆQAI t) Qt) and n 1/2 ˆQD t) Qt) do not have the same limiting distribution, since n 1/2 T n t) is not o p 1), see Theorem 1. For Scenario 2, we can show that ˆπt) P π t) and ˆQ I t) P Qt), uniformly for t [, τ]. We can write [cf. Eq. 15)] T n t) = T n,1 t) + T n,2 t) + T n,3 t) where T n,1 t) = 1 n T n,2 t) = 1 n T n,3 t) = 1 n 1 ξ i π X i ) { ξi ˆπX i ) π X i )) ˆπX i )π X i ) 1 ξ i π X i ) ) p X i )IX i t) 22) } ˆpX i ) IX i t) 23) ) ˆpX i ) p X i )) IX i t). 24) As before, we can show that T n t). P Therefore, for this case also, ˆQ AI t) = ˆQ I t) + T n t) ) ) retains consistency. Again, it may be noted that n 1/2 ˆQAI t) Qt) and n 1/2 ˆQI t) Qt) do not have the same limiting distribution, since n 1/2 Tn t) is not o p 1), see Theorem 2. Let πt) = 1 π t)/π t), which equals when πt) is correctly specified. Along with Eq. 4) and Eq. 5), we also define two other quantities Ct) = { } ξ π X) δ p X)) IX t), 25) Dt) = ξ δ p X)) p X)1 p X)) as well as the following variance function V 1 t): π s) αs, X)dHs), 26) π s) V 1 t) = p s)dqs) Q 2 t) + + πu)αu, v)dhu)dhv) π u) π u)) 2 1 p u))dqu). 27) The following theorem describes the asymptotic distribution of ˆQ AI t) for Scenario 1. 7
8 Theorem 1 Suppose that the regularity conditions R and RPI hold and that π t) is bounded ) away from, uniformly for t [, τ]. The process n 1/2 ˆQAI t) Qt) converges weakly in D[, τ] to a zero-mean Gaussian process Z 1 with variance function V 1 t) defined by Eq. 27). Proof We know that n 1/2 ˆQ D t) Qt)) is asymptotically linear with influence function At) + Bt). From Eq. 19), n 1/2 T n,1 t) is asymptotically linear with influence function Ct). Let ˆ W t) denote the empirical estimator of W t) = P X t, ξ = 1). For Eq. 2), analogous to calculations in Subramanian and Bean 28), we can show that, uniformly for t [, τ], T n,2 t) = ˆps) p s) d ˆ W s) = π s) ˆps) p s) d π s) W s) + o p n 1/2 ). It readily follows, from Eq. 3) and then Eq. 2), that n 1/2 T n,2 t) is also asymptotically linear with influence function Dt). Next, note that dw s) = p s)d W s). It can be shown that, uniformly for t [, τ], T n,3 t) defined by Eq. 21) is asymptotically negligible: T n,3 t) = = ˆπs) π s) ˆπs)π s) {dŵ s) ˆps)d ˆ W s) } ˆπs) π s) { dw s) p π s)) 2 s)d W s) } + o p n 1/2 ) = o p n 1/2 ). It follows that n 1/2 ˆQ AI t) Qt)) is asymptotically linear with influence function At) + Bt)+Ct)+Dt) [cf. Eqs. 4), 5), 25), and 26)]. After some routine variance calculations, it can be shown that the asymptotic variance is given by V 1 t) defined by Eq. 27). We next focus on Scenario 2. Let d p t) = p t) p t), which is when pt) is correctly specified. Along with Eq. 1) and Eq. 11), we also define two other quantities: Ct) = Dt) = 1 ξ ) p X)IX t) 28) π X) ξ π X) π X)1 π X)) 8 p s) βs, X)dHs), 29) π s)
9 as well as the following variance function V 2 t): V 2 t) = 1 π u) dqu) Q2 t) 1 π u) p u)dhu) π u) d p u) d p v) βu, v)dhv)dhu). 3) π u) π v) The following theorem describes the asymptotic distribution of ˆQ AI t) for Scenario 2. Theorem 2 Suppose that the regularity conditions R and RP hold and that π t) is bounded ) away from, uniformly for t [, τ]. The process n 1/2 ˆQAI t) Qt) converges weakly in D[, τ] to a zero-mean Gaussian process Z 2 with variance function V 2 t) defined by Eq. 3). Proof We know that n 1/2 ˆQ I t) Qt)) is asymptotically linear with influence function Ăt) + Bt). From Eq. 22), n 1/2 Tn,1 t) is asymptotically linear with influence function Ct). For T n,2 t) defined by Eq. 23), we can show that, uniformly for t [, τ], T n,2 t) = ˆπs) π s)) ˆps)d ˆ W s) = ˆπs)π s) ˆπs) π s)) p s)d πs) W s) + o 2 p n 1/2 ). It follows, from Eq. 9) and then Eq. 8), that n 1/2 Tn,2 t) is also asymptotically linear with influence function Dt). Furthermore, it can be shown that T n,3 t), defined by Eq. 24) is o p n 1/2 ). It follows that n 1/2 ˆQ AI t) Qt)) is asymptotically linear with influence function Ăt)+ Bt)+ Ct)+ Dt) [cf. Eqs. 1), 11), 28), and 29)]. After some variance calculations, it can be shown that the asymptotic variance is given by V 2 t) defined by Eq. 3). It is difficult to compare the variances V 1 t) and V 2 t) with V D t) [Eq. 6)], or for that matter V I t) [Eq. 12)]. However, for Scenario 3, a definitive conclusion can be made. Define V t) = Qt)1 Qt)) + 1 p u)) 1 π u) π u) dqu). 31) The following theorem describes the asymptotic distribution of ˆQ AI t) for Scenario 3. 9
10 Theorem 3 Suppose that the regularity condition R holds and that π t) is bounded away ) from, uniformly for t [, τ]. The process n 1/2 ˆQAI t) Qt) converges weakly in D[, τ] to a zero-mean Gaussian process Z with variance function V t) defined by Eq. 31). Proof Follow proof of Th. 1 with π t) = π t), or follow proof of Th. 2 with p t) = p t). In each case, the asymptotic variance reduces to V t), since πt) = and d p t) =. Note that when π and p are both correctly specified Scenario 3), the variance of the AIPW estimator equals the information bound for estimating Qt) Subramanian 24a), which was in fact derived from the efficient influence function for estimating St) van der Laan and McKeague, 1998). However, V t) V D t), see Subramanian 24a). In particular, we have the order relation V D t) V t) V I t). Thus, although the AIPW estimator ˆQ AI t) improves upon its IPW precursor, it is still inferior to the Dikta estimator ˆQ D t). 2.3 AIPW survival function estimator We now briefly describe the AIPW estimator of St). This estimator results from employing a series of compactly differentiable mappings of the basic estimators ˆQ AI t) and 1 Ĥt) Gill and Johansen, 199). Specifically, the estimator takes the form Ŝ AI t) = s t π 1 d ˆQs) ) 1 Ĥs ). 32) The weak convergence of ŜAIt) can be derived in a straightforward way through the functional delta method, see p of Gill and Johansen 199). Considering only Scenario 3, n 1/2 { ˆQ AI t) Qt), 1 Ĥt ) 1 Ht)} D Z, Z H ), in D[, τ] D [, τ], ), as n, where Z is the limiting Gaussian process defined in the statement of Th. 3. Here, is the supremum norm, and is the max supremum norm. We then have by 1
11 the functional delta method that n 1/2 {ŜAIt) St)} converges weakly in D[, τ], ) to W t) = St) dw s)/1 Hs)) as n, where W t) = Zt) Z Hs)dΛs). We can then show that the variance function V t) of W t) is given by V t) = St) 2 [ dλu) 1 Hu) + 1 p u)) 1 π u) π u) dλu) 1 Hu) ], 33) which is the asymptotic variance of the reduced data NPMLE van der Laan and McKeague, 1998), and of two competing estimators Subramanian 24b, 26) as well. However, V t) VD t) Subramanian, 24a), where V Dt) = St) 2 [ p s) 1 Hs) dλs) + ] αu, v) 1 Hu))1 Hv)) dhu)dhv) 34) is the asymptotic variance of the Dikta semiparametric estimator ŜDt). Note that V D t) can also be derived from Eq. 6) and the functional delta method described above. 2.4 Extension to a left truncated version The three scenarios investigated above can be analyzed in an analogous way for a left truncated version of the MCI model. For this version, we observe X, Z, ξ, σ) whenever X Z, where Z is a left truncation variable, ξ = 1 when δ is observed and otherwise, and σ = ξδ Subramanian and Bandyopadhyay, 28). When X < Z, however, nothing is observed. The observed data are {X i, Z i, ξ i, σ i } 1 i n, an i.i.d. sample from the distribution of X, Z, ξ, σ) which is observed i.e., when Z i X i ). Let W t) = P Z t X Z X) and W 1 t) = P X t, δ = 1 Z X). The cumulative hazard takes the form Λt) = a W dw 1 u)/ W u), where W denotes the d.f. of X and a W = inf{x : W x) > }. The 11
12 asymptotic variance of the AIPW estimator of St) for Scenario 3 is [compare with Eq. 33)] V T t) = St) 2 [ a W dλu) W u) + 1 p u)) 1 π u) a W π u) ] dλu), W u) where p t) = P δ = 1 X = t, Z X), and π t) = P ξ = 1 X = t, Z X). The Dikta semiparametric estimator for this version of the MCI model has asymptotic variance V T t) = a W p u) W u) dλu) + a W αu, v) a W p u)p v) dλu)dλv), see Subramanian and Bandyopadhyay 28). Moreover, we can show that V T t) V T t). References Dikta, G., On semiparametric random censorship models. J. Statist. Plann. Inference 66, Gill, R. D., Johansen, S., 199. A survey of product-integration with a view toward application in survival analysis. Ann. Statist. 18, Horvitz, D.G., Thompson, D.J., A generalization of sampling without replacement from a finite universe, J. Ameri. Statist. Assoc. 47, Kouassi, D. A., Singh, J., A semiparametric approach to hazard estimation with randomly censored observations. J. Ameri. Statist. Assoc. 92, Koul, H., Susarla, V., Van Ryzin, J., Regression analysis of randomly right censored data. Ann. Statist. 9, van der Laan, M. J., McKeague, I. W., Efficient estimation from right-censored data when failure indicators are missing at random. Ann. Statist. 26,
13 van der Laan, M.J., Robins, J.M., 23. Unified methods for censored longitudinal data and causality. Springer, New York. Robins, J.M., Rotnitzky, A., Recovery of information and adjustment for dependent censoring using surrogate markers. AIDS Epidemiology-Methodological Issues, Eds. N.Jewell, K.Dietz, V.Farewell, Boston: Birkhauser, Robins, J.M., Rotnitzky, A., and van der Laan, M.J., 2. Discussion of On profile likelihood by Murphy and van der Vaart. J. Ameri. Statist. Assoc. 95, Rotnitzky, A., Robins, J.M., 25. Inverse probability weighted estimation in survival analysis, Encyclopedia of Biostatistics, Eds. P. Armitage, T. Colton, New York: Wiley. Satten, G.A., Datta, S., 21. The Kaplan-Meier estimator as an inverse-probability-ofcensoring weighted average, Am. Stat. 55, Satten, G.A., Datta, S., 24. Marginal analysis of multistage data. Handbook of Statistics 23: Advances in Survival Analysis. Eds. N. Balakrishnan and C.R.Rao Scharfstein, D.O., Rotnitzky, A., Robins, J.M., Adjusting for nonignorable drop-out using semiparametric nonresponse models with discussion). J. Ameri. Statist. Assoc. 94, Subramanian, S., 24a. The missing censoring-indicator model of random censorship. Handbook of Statistics 23: Advances in Survival Analysis. Eds. N. Balakrishnan and C.R.Rao Subramanian, S., 24b. Asymptotically efficient estimation of a survival function in the 13
14 missing censoring indicator model. J. Nonparametr. Statist. 16, Subramanian, S., 26. Survival analysis for the missing censoring indicator model using kernel density estimation techniques. Stat. Methodol. 3, Subramanian, S., Bean, D., 28a. Hazard function estimation from homogeneous right censored data with missing censoring indicators. Stat. Methodol. in press). Subramanian, S., Bandyopadhyay, D., 28b. Semiparametric left truncation and right censorship models with missing censoring indicators. Statist. Probab. Lett. 78, Tsiatis, A. A. 26. Semiparametric theory and missing data. Springer, New York. 14
Median regression using inverse censoring weights
Median regression using inverse censoring weights Sundarraman Subramanian 1 Department of Mathematical Sciences New Jersey Institute of Technology Newark, New Jersey Gerhard Dikta Department of Applied
More informationModel checks via bootstrap when there are missing binary data.
1 INTRODUCTION Model checks via bootstrap when there are missing binary data. Gerhard Dikta a, Sundarraman Subramanian b, Thorsten Winkler a a Department of Medizintechnik und Technomathematik, Fachhochschule
More informationEstimation of the Bivariate and Marginal Distributions with Censored Data
Estimation of the Bivariate and Marginal Distributions with Censored Data Michael Akritas and Ingrid Van Keilegom Penn State University and Eindhoven University of Technology May 22, 2 Abstract Two new
More informationWeighting Methods. Harvard University STAT186/GOV2002 CAUSAL INFERENCE. Fall Kosuke Imai
Weighting Methods Kosuke Imai Harvard University STAT186/GOV2002 CAUSAL INFERENCE Fall 2018 Kosuke Imai (Harvard) Weighting Methods Stat186/Gov2002 Fall 2018 1 / 13 Motivation Matching methods for improving
More informationEfficient Semiparametric Estimators via Modified Profile Likelihood in Frailty & Accelerated-Failure Models
NIH Talk, September 03 Efficient Semiparametric Estimators via Modified Profile Likelihood in Frailty & Accelerated-Failure Models Eric Slud, Math Dept, Univ of Maryland Ongoing joint project with Ilia
More informationStatistical 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 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 informationNonparametric Likelihood and Doubly Robust Estimating Equations for Marginal and Nested Structural Models. Zhiqiang Tan 1
Nonparametric Likelihood and Doubly Robust Estimating Equations for Marginal and Nested Structural Models Zhiqiang Tan 1 Abstract. Drawing inferences about treatment effects is of interest in many fields.
More informationModification and Improvement of Empirical Likelihood for Missing Response Problem
UW Biostatistics Working Paper Series 12-30-2010 Modification and Improvement of Empirical Likelihood for Missing Response Problem Kwun Chuen Gary Chan University of Washington - Seattle Campus, kcgchan@u.washington.edu
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 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 informationCovariate Balancing Propensity Score for General Treatment Regimes
Covariate Balancing Propensity Score for General Treatment Regimes Kosuke Imai Princeton University October 14, 2014 Talk at the Department of Psychiatry, Columbia University Joint work with Christian
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 informationApproximate Self Consistency for Middle-Censored Data
Approximate Self Consistency for Middle-Censored Data by S. Rao Jammalamadaka Department of Statistics & Applied Probability, University of California, Santa Barbara, CA 93106, USA. and Srikanth K. Iyer
More informationInconsistency of the MLE for the joint distribution. of interval censored survival times. and continuous marks
Inconsistency of the MLE for the joint distribution of interval censored survival times and continuous marks By M.H. Maathuis and J.A. Wellner Department of Statistics, University of Washington, Seattle,
More informationSTAT Sample Problem: General Asymptotic Results
STAT331 1-Sample Problem: General Asymptotic Results In this unit we will consider the 1-sample problem and prove the consistency and asymptotic normality of the Nelson-Aalen estimator of the cumulative
More 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 informationHarvard University. Harvard University Biostatistics Working Paper Series
Harvard University Harvard University Biostatistics Working Paper Series Year 2015 Paper 197 On Varieties of Doubly Robust Estimators Under Missing Not at Random With an Ancillary Variable Wang Miao Eric
More informationLecture 22 Survival Analysis: An Introduction
University of Illinois Department of Economics Spring 2017 Econ 574 Roger Koenker Lecture 22 Survival Analysis: An Introduction There is considerable interest among economists in models of durations, which
More informationUniversity of California, Berkeley
University of California, Berkeley U.C. Berkeley Division of Biostatistics Working Paper Series Year 24 Paper 153 A Note on Empirical Likelihood Inference of Residual Life Regression Ying Qing Chen Yichuan
More 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 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 informationUniversity of California, Berkeley
University of California, Berkeley U.C. Berkeley Division of Biostatistics Working Paper Series Year 2009 Paper 251 Nonparametric population average models: deriving the form of approximate population
More informationEmpirical Processes & Survival Analysis. The Functional Delta Method
STAT/BMI 741 University of Wisconsin-Madison Empirical Processes & Survival Analysis Lecture 3 The Functional Delta Method Lu Mao lmao@biostat.wisc.edu 3-1 Objectives By the end of this lecture, you will
More informationUniversity of California, Berkeley
University of California, Berkeley U.C. Berkeley Division of Biostatistics Working Paper Series Year 2010 Paper 260 Collaborative Targeted Maximum Likelihood For Time To Event Data Ori M. Stitelman Mark
More informationConditional Empirical Likelihood Approach to Statistical Analysis with Missing Data
Conditional Empirical Likelihood Approach to Statistical Analysis with Missing Data by Peisong Han A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy
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 informationINVERSE PROBABILITY WEIGHTED ESTIMATION FOR GENERAL MISSING DATA PROBLEMS
IVERSE PROBABILITY WEIGHTED ESTIMATIO FOR GEERAL MISSIG DATA PROBLEMS Jeffrey M. Wooldridge Department of Economics Michigan State University East Lansing, MI 48824-1038 (517) 353-5972 wooldri1@msu.edu
More informationSemiparametric maximum likelihood estimation in normal transformation models for bivariate survival data
Biometrika (28), 95, 4,pp. 947 96 C 28 Biometrika Trust Printed in Great Britain doi: 1.193/biomet/asn49 Semiparametric maximum likelihood estimation in normal transformation models for bivariate survival
More informationOptimal Treatment Regimes for Survival Endpoints from a Classification Perspective. Anastasios (Butch) Tsiatis and Xiaofei Bai
Optimal Treatment Regimes for Survival Endpoints from a Classification Perspective Anastasios (Butch) Tsiatis and Xiaofei Bai Department of Statistics North Carolina State University 1/35 Optimal Treatment
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 informationGoodness-of-fit tests for the cure rate in a mixture cure model
Biometrika (217), 13, 1, pp. 1 7 Printed in Great Britain Advance Access publication on 31 July 216 Goodness-of-fit tests for the cure rate in a mixture cure model BY U.U. MÜLLER Department of Statistics,
More informationThis is the submitted version of the following book chapter: stat08068: Double robustness, which will be
This is the submitted version of the following book chapter: stat08068: Double robustness, which will be published in its final form in Wiley StatsRef: Statistics Reference Online (http://onlinelibrary.wiley.com/book/10.1002/9781118445112)
More informationTopics and Papers for Spring 14 RIT
Eric Slud Feb. 3, 204 Topics and Papers for Spring 4 RIT The general topic of the RIT is inference for parameters of interest, such as population means or nonlinearregression coefficients, in the presence
More informationNonparametric rank based estimation of bivariate densities given censored data conditional on marginal probabilities
Hutson Journal of Statistical Distributions and Applications (26 3:9 DOI.86/s4488-6-47-y RESEARCH Open Access Nonparametric rank based estimation of bivariate densities given censored data conditional
More informationPackage Rsurrogate. October 20, 2016
Type Package Package Rsurrogate October 20, 2016 Title Robust Estimation of the Proportion of Treatment Effect Explained by Surrogate Marker Information Version 2.0 Date 2016-10-19 Author Layla Parast
More information5 Methods Based on Inverse Probability Weighting Under MAR
5 Methods Based on Inverse Probability Weighting Under MAR The likelihood-based and multiple imputation methods we considered for inference under MAR in Chapters 3 and 4 are based, either directly or indirectly,
More informationSize and Shape of Confidence Regions from Extended Empirical Likelihood Tests
Biometrika (2014),,, pp. 1 13 C 2014 Biometrika Trust Printed in Great Britain Size and Shape of Confidence Regions from Extended Empirical Likelihood Tests BY M. ZHOU Department of Statistics, University
More informationRobustness to Parametric Assumptions in Missing Data Models
Robustness to Parametric Assumptions in Missing Data Models Bryan Graham NYU Keisuke Hirano University of Arizona April 2011 Motivation Motivation We consider the classic missing data problem. In practice
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 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 informationCalibration Estimation for Semiparametric Copula Models under Missing Data
Calibration Estimation for Semiparametric Copula Models under Missing Data Shigeyuki Hamori 1 Kaiji Motegi 1 Zheng Zhang 2 1 Kobe University 2 Renmin University of China Economics and Economic Growth Centre
More informationNonparametric estimation of linear functionals of a multivariate distribution under multivariate censoring with applications.
Chapter 1 Nonparametric estimation of linear functionals of a multivariate distribution under multivariate censoring with applications. 1.1. Introduction The study of dependent durations arises in many
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 informationLikelihood Construction, Inference for Parametric Survival Distributions
Week 1 Likelihood Construction, Inference for Parametric Survival Distributions In this section we obtain the likelihood function for noninformatively rightcensored survival data and indicate how to make
More informationSurvival Analysis: Weeks 2-3. Lu Tian and Richard Olshen Stanford University
Survival Analysis: Weeks 2-3 Lu Tian and Richard Olshen Stanford University 2 Kaplan-Meier(KM) Estimator Nonparametric estimation of the survival function S(t) = pr(t > t) The nonparametric estimation
More informationUniversity of California, Berkeley
University of California, Berkeley U.C. Berkeley Division of Biostatistics Working Paper Series Year 2009 Paper 248 Application of Time-to-Event Methods in the Assessment of Safety in Clinical Trials Kelly
More informationApplication of Time-to-Event Methods in the Assessment of Safety in Clinical Trials
Application of Time-to-Event Methods in the Assessment of Safety in Clinical Trials Progress, Updates, Problems William Jen Hoe Koh May 9, 2013 Overview Marginal vs Conditional What is TMLE? Key Estimation
More informationEstimation for two-phase designs: semiparametric models and Z theorems
Estimation for two-phase designs:semiparametric models and Z theorems p. 1/27 Estimation for two-phase designs: semiparametric models and Z theorems Jon A. Wellner University of Washington Estimation for
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 informationarxiv: v2 [stat.me] 17 Jan 2017
Semiparametric Estimation with Data Missing Not at Random Using an Instrumental Variable arxiv:1607.03197v2 [stat.me] 17 Jan 2017 BaoLuo Sun 1, Lan Liu 1, Wang Miao 1,4, Kathleen Wirth 2,3, James Robins
More informationANALYSIS OF COMPETING RISKS DATA WITH MISSING CAUSE OF FAILURE UNDER ADDITIVE HAZARDS MODEL
Statistica Sinica 18(28, 219-234 ANALYSIS OF COMPETING RISKS DATA WITH MISSING CAUSE OF FAILURE UNDER ADDITIVE HAZARDS MODEL Wenbin Lu and Yu Liang North Carolina State University and SAS Institute Inc.
More informationBounded, Efficient, and Doubly Robust Estimation with Inverse Weighting
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 Biometrika (2008), 94, 2, pp. 1 22 C 2008 Biometrika Trust Printed
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 informationSession 9: Introduction to Sieve Analysis of Pathogen Sequences, for Assessing How VE Depends on Pathogen Genomics Part I
Session 9: Introduction to Sieve Analysis of Pathogen Sequences, for Assessing How VE Depends on Pathogen Genomics Part I Peter B Gilbert Vaccine and Infectious Disease Division, Fred Hutchinson Cancer
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 informationAssessing the effect of a partly unobserved, exogenous, binary time-dependent covariate on -APPENDIX-
Assessing the effect of a partly unobserved, exogenous, binary time-dependent covariate on survival probabilities using generalised pseudo-values Ulrike Pötschger,2, Harald Heinzl 2, Maria Grazia Valsecchi
More informationA Sampling of IMPACT Research:
A Sampling of IMPACT Research: Methods for Analysis with Dropout and Identifying Optimal Treatment Regimes Marie Davidian Department of Statistics North Carolina State University http://www.stat.ncsu.edu/
More informationHarvard University. Harvard University Biostatistics Working Paper Series
Harvard University Harvard University Biostatistics Working Paper Series Year 2008 Paper 85 Semiparametric Maximum Likelihood Estimation in Normal Transformation Models for Bivariate Survival Data Yi Li
More informationAsymptotic Properties of Nonparametric Estimation Based on Partly Interval-Censored Data
Asymptotic Properties of Nonparametric Estimation Based on Partly Interval-Censored Data Jian Huang Department of Statistics and Actuarial Science University of Iowa, Iowa City, IA 52242 Email: jian@stat.uiowa.edu
More informationPart III Measures of Classification Accuracy for the Prediction of Survival Times
Part III Measures of Classification Accuracy for the Prediction of Survival Times Patrick J Heagerty PhD Department of Biostatistics University of Washington 102 ISCB 2010 Session Three Outline Examples
More informationCalibration Estimation of Semiparametric Copula Models with Data Missing at Random
Calibration Estimation of Semiparametric Copula Models with Data Missing at Random Shigeyuki Hamori 1 Kaiji Motegi 1 Zheng Zhang 2 1 Kobe University 2 Renmin University of China Institute of Statistics
More informationMaximum likelihood: counterexamples, examples, and open problems. Jon A. Wellner. University of Washington. Maximum likelihood: p.
Maximum likelihood: counterexamples, examples, and open problems Jon A. Wellner University of Washington Maximum likelihood: p. 1/7 Talk at University of Idaho, Department of Mathematics, September 15,
More informationComment: Understanding OR, PS and DR
Statistical Science 2007, Vol. 22, No. 4, 560 568 DOI: 10.1214/07-STS227A Main article DOI: 10.1214/07-STS227 c Institute of Mathematical Statistics, 2007 Comment: Understanding OR, PS and DR Zhiqiang
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 informationSmooth nonparametric estimation of a quantile function under right censoring using beta kernels
Smooth nonparametric estimation of a quantile function under right censoring using beta kernels Chanseok Park 1 Department of Mathematical Sciences, Clemson University, Clemson, SC 29634 Short Title: Smooth
More informationPropensity Score Weighting with Multilevel Data
Propensity Score Weighting with Multilevel Data Fan Li Department of Statistical Science Duke University October 25, 2012 Joint work with Alan Zaslavsky and Mary Beth Landrum Introduction In comparative
More informationEmpirical Likelihood in Survival Analysis
Empirical Likelihood in Survival Analysis Gang Li 1, Runze Li 2, and Mai Zhou 3 1 Department of Biostatistics, University of California, Los Angeles, CA 90095 vli@ucla.edu 2 Department of Statistics, The
More informationST745: Survival Analysis: Nonparametric methods
ST745: Survival Analysis: Nonparametric methods Eric B. Laber Department of Statistics, North Carolina State University February 5, 2015 The KM estimator is used ubiquitously in medical studies to estimate
More informationENHANCED PRECISION IN THE ANALYSIS OF RANDOMIZED TRIALS WITH ORDINAL OUTCOMES
Johns Hopkins University, Dept. of Biostatistics Working Papers 10-22-2014 ENHANCED PRECISION IN THE ANALYSIS OF RANDOMIZED TRIALS WITH ORDINAL OUTCOMES Iván Díaz Johns Hopkins University, Johns Hopkins
More informationSimple design-efficient calibration estimators for rejective and high-entropy sampling
Biometrika (202), 99,, pp. 6 C 202 Biometrika Trust Printed in Great Britain Advance Access publication on 3 July 202 Simple design-efficient calibration estimators for rejective and high-entropy sampling
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 informationDiscussion of Missing Data Methods in Longitudinal Studies: A Review by Ibrahim and Molenberghs
Discussion of Missing Data Methods in Longitudinal Studies: A Review by Ibrahim and Molenberghs Michael J. Daniels and Chenguang Wang Jan. 18, 2009 First, we would like to thank Joe and Geert for a carefully
More informationEmpirical likelihood ratio with arbitrarily censored/truncated data by EM algorithm
Empirical likelihood ratio with arbitrarily censored/truncated data by EM algorithm Mai Zhou 1 University of Kentucky, Lexington, KY 40506 USA Summary. Empirical likelihood ratio method (Thomas and Grunkmier
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 informationCalibration Estimation of Semiparametric Copula Models with Data Missing at Random
Calibration Estimation of Semiparametric Copula Models with Data Missing at Random Shigeyuki Hamori 1 Kaiji Motegi 1 Zheng Zhang 2 1 Kobe University 2 Renmin University of China Econometrics Workshop UNC
More informationLarge sample theory for merged data from multiple sources
Large sample theory for merged data from multiple sources Takumi Saegusa University of Maryland Division of Statistics August 22 2018 Section 1 Introduction Problem: Data Integration Massive data are collected
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 informationRobust estimates of state occupancy and transition probabilities for Non-Markov multi-state models
Robust estimates of state occupancy and transition probabilities for Non-Markov multi-state models 26 March 2014 Overview Continuously observed data Three-state illness-death General robust estimator Interval
More informationCensoring and Truncation - Highlighting the Differences
Censoring and Truncation - Highlighting the Differences Micha Mandel The Hebrew University of Jerusalem, Jerusalem, Israel, 91905 July 9, 2007 Micha Mandel is a Lecturer, Department of Statistics, The
More informationCensoring Unbiased Regression Trees and Ensembles
Johns Hopkins University, Dept. of Biostatistics Working Papers 1-31-216 Censoring Unbiased Regression Trees and Ensembles Jon Arni Steingrimsson Department of Biostatistics, Johns Hopkins Bloomberg School
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 informationNonparametric Bayes Estimator of Survival Function for Right-Censoring and Left-Truncation Data
Nonparametric Bayes Estimator of Survival Function for Right-Censoring and Left-Truncation Data Mai Zhou and Julia Luan Department of Statistics University of Kentucky Lexington, KY 40506-0027, U.S.A.
More informationWeighted Likelihood for Semiparametric Models and Two-phase Stratified Samples, with Application to Cox Regression
doi:./j.467-9469.26.523.x Board of the Foundation of the Scandinavian Journal of Statistics 26. Published by Blackwell Publishing Ltd, 96 Garsington Road, Oxford OX4 2DQ, UK and 35 Main Street, Malden,
More informationDouble Robustness. Bang and Robins (2005) Kang and Schafer (2007)
Double Robustness Bang and Robins (2005) Kang and Schafer (2007) Set-Up Assume throughout that treatment assignment is ignorable given covariates (similar to assumption that data are missing at random
More informationData Integration for Big Data Analysis for finite population inference
for Big Data Analysis for finite population inference Jae-kwang Kim ISU January 23, 2018 1 / 36 What is big data? 2 / 36 Data do not speak for themselves Knowledge Reproducibility Information Intepretation
More informationTargeted Maximum Likelihood Estimation for Adaptive Designs: Adaptive Randomization in Community Randomized Trial
Targeted Maximum Likelihood Estimation for Adaptive Designs: Adaptive Randomization in Community Randomized Trial Mark J. van der Laan 1 University of California, Berkeley School of Public Health laan@berkeley.edu
More informationExtending causal inferences from a randomized trial to a target population
Extending causal inferences from a randomized trial to a target population Issa Dahabreh Center for Evidence Synthesis in Health, Brown University issa dahabreh@brown.edu January 16, 2019 Issa Dahabreh
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 informationAsymptotic Distributions for the Nelson-Aalen and Kaplan-Meier estimators and for test statistics.
Asymptotic Distributions for the Nelson-Aalen and Kaplan-Meier estimators and for test statistics. Dragi Anevski Mathematical Sciences und University November 25, 21 1 Asymptotic distributions for statistical
More 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 informationNONPARAMETRIC BAYES ESTIMATOR OF SURVIVAL FUNCTIONS FOR DOUBLY/INTERVAL CENSORED DATA
Statistica Sinica 14(2004, 533-546 NONPARAMETRIC BAYES ESTIMATOR OF SURVIVAL FUNCTIONS FOR DOUBLY/INTERVAL CENSORED DATA Mai Zhou University of Kentucky Abstract: The non-parametric Bayes estimator with
More informationSimultaneous confidence bands using model assisted Cox regression
Simultaneous confidence bands using model assisted Cox regression Shoubhik Mondal and Sundarraman Subramanian Center for Applied Mathematics and Statistics Department of Mathematical Sciences New Jersey
More informationStep-Stress Models and Associated Inference
Department of Mathematics & Statistics Indian Institute of Technology Kanpur August 19, 2014 Outline Accelerated Life Test 1 Accelerated Life Test 2 3 4 5 6 7 Outline Accelerated Life Test 1 Accelerated
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 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 informationInvestigation of goodness-of-fit test statistic distributions by random censored samples
d samples Investigation of goodness-of-fit test statistic distributions by random censored samples Novosibirsk State Technical University November 22, 2010 d samples Outline 1 Nonparametric goodness-of-fit
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 informationEMPIRICAL LIKELIHOOD ANALYSIS FOR THE HETEROSCEDASTIC ACCELERATED FAILURE TIME MODEL
Statistica Sinica 22 (2012), 295-316 doi:http://dx.doi.org/10.5705/ss.2010.190 EMPIRICAL LIKELIHOOD ANALYSIS FOR THE HETEROSCEDASTIC ACCELERATED FAILURE TIME MODEL Mai Zhou 1, Mi-Ok Kim 2, and Arne C.
More informationUniversity of California, Berkeley
University of California, Berkeley U.C. Berkeley Division of Biostatistics Working Paper Series Year 2011 Paper 288 Targeted Maximum Likelihood Estimation of Natural Direct Effect Wenjing Zheng Mark J.
More information