Smoothing the Nelson-Aalen Estimtor Biostat 277 presentation Chi-hong Tseng
|
|
- Arabella Nicholson
- 6 years ago
- Views:
Transcription
1 Smoothing the Nelson-Aalen Estimtor Biostat 277 presentation Chi-hong seng Reference: 1. Andersen, Borgan, Gill, and Keiding (1993). Statistical Model Based on Counting Processes, Springer-Verlag, p Muller and Wang (1994). Hazard Rate Estimation Under Random Cnsoring with Varying Kernels and Bandwidths, Biometrics 50, Survival Data and Nelson-Aalen Estimator 1. Survival time 2. Censoring time C 3. Oservational time X = min(, C) 4. Censoring indicator D = I( C) 5. Hazard function: α(t) = lim Pr(t < t + t t) t 6. Cumulative hazard function : A(t) = t 0 α(u)du 7. Nelson-Aalen estimator Â(t) = D i n i t i t 8. n i = numer of suject at risk at time i 9. at-risk process Y i (t) = I(X i t), Y (t) = Y i (t) 10. counting process N i (t) = I(X i t, D i = 1), N(t) = N i (t) 11. Nelson-Aalen estimator Â(t) = t 12. J(s) = I(Y (s) > 0) 2 Kernel Function Estimator 0 J(u) n Y (u) dn(u) = he kernel function estimator is given y, ˆα(t) = K( t s )dâ(s) = j i=1 K( t j )(Y ( j )) 1. K: Kernel function, 1 K(s)ds = 1 e.g. K(x) = 0.75(1 x 2 ), x : andwidth 3. at each t, {j, t j t + } contriute to the sum. 1 t 0 J(u) Y (u) dn i(u)
2 4. Eˆα(t) = Consider α (t) = α(t) K( t s )Pr(Y (s) > 0)dA(s) = = K( t s )da (s) K( t s )J(s)dA(s) (almost = Eˆα(t)) K( t s )da(s) we can see that ˆα(t) α (t) is a martingale: ˆα(t) α (t) = = K( t s )d(â A )(s) K( t s and its expected predictale variation process: ) J(s)dM(s). Y (s) τ 2 (t) = E(ˆα(t) α (t)) 2 = 2 K 2 ( t s )E( J(s) Y (s) )α(s)ds 1 = K 2 J(t u) (u)e( )α(t u)du Y (t u) therefore an estimator of variance of ˆα(t): ˆτ 2 (t) = 2 K 2 ( t s ) J(s) Y 2 (s) dn(s) 3 Mean Intergrated Squared Error and Optimal Bandwith o otain the optimal gloal andwidth, we choose the one minimizing the mean intergrated squared error(mise), which is defined as, MISE(ˆα(t)) = E t 2 (ˆα(t) α(t)) 2 dt = E t 2 [(ˆα(t) α(t)) + ( α(t) α(t))] 2 dt = t 2 E(ˆα(t) α(t)) 2 dt + t 2 ( α(t) α(t)) 2 dt +2 t 2 (Eˆα(t) α(t))( α(t) α(t))dt = ( variance term) + (squared ias term) + R 2 (t) 2
3 cumulative hazard BM data time andwidth=100 smoothed hazard time 3
4 smoothed hazard andwidth= time andwidth=150 smoothed hazard time 4
5 Let with Assume t2 t2 R 2 = 2 (Eˆα(t) α(t))( α(t) α(t))dt = 2 R 1 (t)( α(t) α(t))dt R 1 (t) = (Eˆα(t) α(t)) = K(t s )Pr(Y (s) = 0)dA(s) = 1 K(u)Pr(Y (t nu) = 0)α(t n u)du R 1 (t) C 1 sup t [t1 c,t 2 +c] Pr(Y (t) = 0) R 2 (t) C 2 sup t [t1 c,t 2 +c] Pr(Y (t) = 0) 1. 1 K(u)du = 1, 1 uk(u)du = 0, 1 u2 K(u)du = k 2 > 0 2. there exist a sequence {a n }, a n as n (in most case, a n = n) 3. there exists a continuous function y such that E( a2 nj(t) Y (t) ) 1 y(t) uniformly on [ c, t 2 + c], as n 4. sup t [t1 c,t 2 +c] Pr(Y (t) = 0) = o(a 2 n ) Since α(t) α(t) = 1 K(u)(α(t nu) α(t))du = n α (t) 1 uk(u)du nα (t) 1 u2 K(u)du + o( 2 n) = nα (t)k 2 + o( 2 n) the squared ias term = t 2 ( α(t) α(t)) 2 dt = nk2 2 t2 (α (t)) 2 dt + o( 4 n) For the variance term: E(ˆα(t) α(t)) 2 = E(ˆα(t) α (t)) 2 + E(α (t) α(t)) 2 +2E(ˆα(t) α (t))(α (t) α(t)) the first term on the right-hand side (see ˆτ 2 (t)) can e expressed y, (a 2 n n ) α(t) y(t) 1 and the second term on the right-hand side: K 2 (u)du + o((a 2 n n ) ) α (t) α(t) = 1 K(u)I(Y (t n u) = 0)α(t n u)du = o(a 2 n ) 5
6 y Cauchy -Schwarz inequality, the last term on the right-hand side is of the form o((a 2 n 1/2 n ) ). so the variance term can e written as, t2 All together, E(ˆα(t) α(t)) 2 dt = (a 2 n n) MISE(ˆα) = nk 2 2 t2 1 (α (t)) 2 dt+a 2 n n K 2 (t)dt 1 t2 t2 K 2 (t)dt α(t) y(t) dt + o((a2 n n) ) α(t) y(t) dt+o(4 n)+o((a 2 n n ) ) to minimize the sum of the first two terms, the optimal andwidth is [ 1 t2 ] 1/5 [ n = an 2/5 k 2/5 2 K 2 α(t) t2 (t)dt y(t) dt (α (t)) 2 dt wo quantities are remained to estimated: α (t) and t 2 α(t) y(t) dt ˆα (t) = 3 K ( t s )dâ(s) and, an estimator of t 2 α(t) y(t) dt is ] /5 Because the ias is a 2 n t2 dâ(t) Y (t) dt. E(ˆα(t)) α(t) α (t)k 2 the ias-corrected estimate is ˆα(t) nˆα (t)k Cross Validation Method We can expressed the MISE, MISE() = E t2 (ˆα(t)) 2 dt 2E t2 ˆα(t)α(t)dt + t2 (α(t)) 2 dt, it s a function of ; the first term is easy to caluclate, the third term deos not depend on, and a consistent estimator of the second term is the cross -validation estimator: 2 i j 1 K( i j ) N( i) N( j ) Y ( i ) Y ( j ) and plotting MISE() against will find the optimal 6
7 3.2 Boundary Correction and Local Bandwidth Now we allow oth andwidth = (t) and kernel K = K t depend on t: ˆα(t, (t)) = (t) j K t ( t j (t) )(Y ( j)) Assuming the support of the survival time distriution is on [0, R]. interior points I = t : (t) t R (t) left oundary region B L = t : 0 t < (t) right oundary region B R = t : R (t) < t R the kernels K t are polynomials K t (z) = { K + (1, z) t I, K + (t/(t), z) t B L, K ((R t)/(t), z) t B R, where K ± : [0, 1] [, 1] R are ounded kernel function with K ± (q, z)dz = 1, K ± (q, z)zdz = 0, K ± (q, z)z 2 dz 0, K (q, z) = K + (q, z) and some differentiaility conditions. e.g, K + (q, z) = 12 (1 + q) 4(z + 1)[z(1 2q) + (3q2 2q + 1)/2] and K + (1, z) = 3 4 (1 x2 ) is the Epanechnikov kernel For local andwidth, we consider to minimize local mean square error to otain andwidth (t) MSE(t, (t)) = ˆv(t, (t)) + ˆβ 2 (t, (t)) ˆv and ˆβ are variance and ias estimatrs, ˆv(t, (t)) = 1 K 2 ˆα(t (t)u) t (u)( n(t) ŷ(t (t)u) )du ˆβ(t, (t)) = ˆα(t (t)u)k t (u)du ˆα(t) he following is the recommended algorithm to otain the α(t) with (t). 7
8 1 choose K + (q, z) and initial andwidth 0 to otain ˆα(t) 1a alternatively, choose a parametric model and otain MLE ˆα(t) 2 choose equidistant grid of m l points t i, i = 1 m l, etween 0 and R. choose a grid of andwidth j, j = 1,, l, etween some 1 and 2. compute ˆv(t i ) and ˆβ(t i ) with all andwidth j s and otain minimizer (t i ) of MSE(t i ). 3 otain final andwidth using oundary -modified smoother with andwidth 0 or ˆ(t) = Kt ( t t i ) (t i )/ K t ( t t i ) otain final estimate ˆα(t) with (t) = ˆ(t) 4 Large Sample Properties Now we discuss the large sample properties of the kernel function estimator with gloal andwidth. heoremn IV.2.1 pointwise consistency Assume 1. t e an interior point of 2. α is continuous at t 3. andwith n 0 as n 4. there exist an ǫ > 0 such that inf s [t ǫ,t+ǫ] n Y (s) p as n hen ˆα p α as n Proof. ˆα(t) α(t) ˆα(t) α (t) + α (t) α(t) + α(t) α(t) we have to show ˆα(t) α (t) p 0 α (t) α(t) p 0 α(t) α(t) p 0 8
9 Pr( ˆα(t) α (t) > η) Pr( sup n t n s t+ n 1 δ η 2 + Pr( n y Lenglart s inequality K( t u ) J(u)dM(u) t n n Y (u) > η) s K 2 (u) J(t nu)α(t n u)du ) > δ) Y (t n u) since α and K are ounded, and n Y (t) p in a neighorhood of t, the last term on the right-hand side can e aritrary small. so ˆα(t) α (t) p 0 α (t) α(t) 1 α(t) α(t) K(u) {1 J(t n u)}α(t n u)du p 0 1 K(u) α(t n u) α(t) du p 0 heoremn IV.2.2 uniform conisitency Assume 1. t e an interior point of 2. 0 < < t 2 < t e fixed numers 3. α is continuous on [0, t] 4. the kernel K is of ounded variation 5. andwith n 0 as n 6. 2 n t 0 J(s)α(s)ds Y (s) p 0 7. t 0 (1 J(s))α(s)ds p 0 hen as n, Proof. sup ˆα(s) α(s) p 0 s [,t 2 ] ˆα(t) α (t) = K( t s )d(â A )(s) 2 V (K) sup Â(s) t [0,t] A (s) where V (K) is the toal variation of K(t). similar to the proof of the consistency of Nelson-Aalen estimator, with assumption (6) the last term converges to zero. 9
10 sup t [t1,t 2 ] α (t) α(t) p 0 with assumption (7) sup t [t1,t 2 ] α(t) α(t) p 0 follows y the oundedness of K and continuity of α. Remark Sufficient conditions for consistency uniform consistency of Nelson Aalen estimator : inf s [0,t] Y (s) p pointwise consistency of kernel function estimator : inf s [0,t] Y (s) p uniform consistency of kernel function estimator : inf s [0,t] 2 Y (s) p heorem IV.2.4 Asymptotic normality Assume 1. t e an interior point of 2. α is continuous at t 3. there exits positive constants {a n }, increasing to as n, 4. n 0, a 2 n n as n 5. there exists a function y, positive and continuous at t such that hen as n, for an ǫ sup a 2 n Y (s) y(s) p 0 s [t ǫ,t+ǫ] 1. a n 1/2 n (ˆα(t) α(t)) d N(0, τ 2 (t)), τ 2 (t) = α(t) y(t) 2. a 2 n nˆτ 2 (t) p τ for t 2, ˆα( ) and ˆα(t 2 ) are asymptotic independent. Proof. A. a n n 1/2 (ˆα(t) α (t)) = with H(s) = a n n /2 K( t s ) J(s) n Y (s). Using Reolledo s martingale CL, H(s)dM(s), K 2 (u)du inf H 2 (s)y (s)α(s)ds = a2 n n K 2 ( t s n ) J(s)α(s)ds Y (s) = 1 a2 nk 2 (u) J(t nu)α(t n u)du Y (t n u) p α(t) y(t) 1 K2 (u)du, 10 as n
11 With any ε > 0, I( H(s) > ε) = I( K( t s ) J(s)a2 n n Y (s) > εa nn 1/2 ) converges to zero uniformly ecause a n 1 n /2, therefore Also inf H2 (s)y (s)α(s)i( H(s) > ε)ds 0 n a n 1/2 n (α (t) α(t)) 0 B. since a n n 1/2 (ˆα(t i ) α (t i )) = with H i (s) = a n n /2 K( t i s ) J(s) Y (s) therefore, H i (s)dm(s), H 1(s)H 2 (s)y (s)α(s)ds = a2 n n K( s n )K( t 2 s n ) J(s)α(s)ds Y (s) 0 heorem IV.2.5 Assume all conditions in heorem IV.2.4, and 1. α is twice continuously differentiale in a neighourhood of t 2. 1 K(u)du = 1, 1 uk(u)du = 0, 1 u2 K(u)du = k 2 > 0 3. limsup n a 2/5 n n < hen a n 1/2 n (ˆα(t) α(t) 2 2 nα (t)k 2 ) d N(0, τ 2 (t)) Proof. With aylor expansion (t is etween t n u and t) a n n 1/2 ( α(t) [ α(t) 2 2 nα (t)k 2 ) ] = a n n 1/2 1 K(u)α(t nu)du α(t) 2 2 n α (t)k 2 = 1 2 a nn 2/5 [ ] u 2 K(u)α(t )du α (t)k
STAT Sample Problem: General Asymptotic Results
STAT331 1-Sample Problem: General Asymptotic Results In this unit we will consider the 1-sample problem and prove the consistency and asymptotic normality of the Nelson-Aalen estimator of the cumulative
More 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 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 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 informationMAXIMUM LOCAL PARTIAL LIKELIHOOD ESTIMATORS FOR THE COUNTING PROCESS INTENSITY FUNCTION AND ITS DERIVATIVES
Statistica Sinica 21 (211), 17-128 MAXIMUM LOCAL PARTIAL LIKELIHOOD ESTIMATORS FOR THE COUNTING PROCESS INTENSITY FUNCTION AND ITS DERIVATIVES Feng Chen University of New South Wales Abstract: We propose
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 informationA GENERALIZED ADDITIVE REGRESSION MODEL FOR SURVIVAL TIMES 1. By Thomas H. Scheike University of Copenhagen
The Annals of Statistics 21, Vol. 29, No. 5, 1344 136 A GENERALIZED ADDITIVE REGRESSION MODEL FOR SURVIVAL TIMES 1 By Thomas H. Scheike University of Copenhagen We present a non-parametric survival model
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 informationEstimation of a quadratic regression functional using the sinc kernel
Estimation of a quadratic regression functional using the sinc kernel Nicolai Bissantz Hajo Holzmann Institute for Mathematical Stochastics, Georg-August-University Göttingen, Maschmühlenweg 8 10, D-37073
More informationAMS 147 Computational Methods and Applications Lecture 13 Copyright by Hongyun Wang, UCSC
Lecture 13 Copyright y Hongyun Wang, UCSC Recap: Fitting to exact data *) Data: ( x j, y j ), j = 1,,, N y j = f x j *) Polynomial fitting Gis phenomenon *) Cuic spline Convergence of cuic spline *) Application
More informationEVALUATIONS OF EXPECTED GENERALIZED ORDER STATISTICS IN VARIOUS SCALE UNITS
APPLICATIONES MATHEMATICAE 9,3 (), pp. 85 95 Erhard Cramer (Oldenurg) Udo Kamps (Oldenurg) Tomasz Rychlik (Toruń) EVALUATIONS OF EXPECTED GENERALIZED ORDER STATISTICS IN VARIOUS SCALE UNITS Astract. We
More informationNonparametric two-sample tests of longitudinal data in the presence of a terminal event
Nonparametric two-sample tests of longitudinal data in the presence of a terminal event Jinheum Kim 1, Yang-Jin Kim, 2 & Chung Mo Nam 3 1 Department of Applied Statistics, University of Suwon, 2 Department
More informationA General Kernel Functional Estimator with Generalized Bandwidth Strong Consistency and Applications
A General Kernel Functional Estimator with Generalized Bandwidth Strong Consistency and Applications Rafael Weißbach Fakultät Statistik, Universität Dortmund e-mail: Rafael.Weissbach@uni-dortmund.de March
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 informationAdaptive Kernel Estimation of The Hazard Rate Function
Adaptive Kernel Estimation of The Hazard Rate Function Raid Salha Department of Mathematics, Islamic University of Gaza, Palestine, e-mail: rbsalha@mail.iugaza.edu Abstract In this paper, we generalized
More informationSupplemental Material for KERNEL-BASED INFERENCE IN TIME-VARYING COEFFICIENT COINTEGRATING REGRESSION. September 2017
Supplemental Material for KERNEL-BASED INFERENCE IN TIME-VARYING COEFFICIENT COINTEGRATING REGRESSION By Degui Li, Peter C. B. Phillips, and Jiti Gao September 017 COWLES FOUNDATION DISCUSSION PAPER NO.
More informationSTAT 331. Martingale Central Limit Theorem and Related Results
STAT 331 Martingale Central Limit Theorem and Related Results In this unit we discuss a version of the martingale central limit theorem, which states that under certain conditions, a sum of orthogonal
More informationThe covariate order method for nonparametric exponential regression and some applications in other lifetime models.
The covariate order method for nonparametric exponential regression and some applications in other lifetime models. Jan Terje Kvaløy and Bo Henry Lindqvist Department of Mathematics and Natural Science
More informationKernel Density Estimation
Kernel Density Estimation Univariate Density Estimation Suppose tat we ave a random sample of data X 1,..., X n from an unknown continuous distribution wit probability density function (pdf) f(x) and cumulative
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 informationKernel density estimation in R
Kernel density estimation in R Kernel density estimation can be done in R using the density() function in R. The default is a Guassian kernel, but others are possible also. It uses it s own algorithm to
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 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 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 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 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 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 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 informationORTHOGONAL SERIES REGRESSION ESTIMATORS FOR AN IRREGULARLY SPACED DESIGN
APPLICATIONES MATHEMATICAE 7,3(000), pp. 309 318 W.POPIŃSKI(Warszawa) ORTHOGONAL SERIES REGRESSION ESTIMATORS FOR AN IRREGULARLY SPACED DESIGN Abstract. Nonparametric orthogonal series regression function
More informationCUSUM TEST FOR PARAMETER CHANGE IN TIME SERIES MODELS. Sangyeol Lee
CUSUM TEST FOR PARAMETER CHANGE IN TIME SERIES MODELS Sangyeol Lee 1 Contents 1. Introduction of the CUSUM test 2. Test for variance change in AR(p) model 3. Test for Parameter Change in Regression Models
More informationNonparametric Density Estimation
Nonparametric Density Estimation Econ 690 Purdue University Justin L. Tobias (Purdue) Nonparametric Density Estimation 1 / 29 Density Estimation Suppose that you had some data, say on wages, and you wanted
More informationNONPARAMETRIC CONFIDENCE INTERVALS FOR MONOTONE FUNCTIONS. By Piet Groeneboom and Geurt Jongbloed Delft University of Technology
NONPARAMETRIC CONFIDENCE INTERVALS FOR MONOTONE FUNCTIONS By Piet Groeneboom and Geurt Jongbloed Delft University of Technology We study nonparametric isotonic confidence intervals for monotone functions.
More informationSTAT331 Lebesgue-Stieltjes Integrals, Martingales, Counting Processes
STAT331 Lebesgue-Stieltjes Integrals, Martingales, Counting Processes This section introduces Lebesgue-Stieltjes integrals, and defines two important stochastic processes: a martingale process and a counting
More informationCALCULATION METHOD FOR NONLINEAR DYNAMIC LEAST-ABSOLUTE DEVIATIONS ESTIMATOR
J. Japan Statist. Soc. Vol. 3 No. 200 39 5 CALCULAION MEHOD FOR NONLINEAR DYNAMIC LEAS-ABSOLUE DEVIAIONS ESIMAOR Kohtaro Hitomi * and Masato Kagihara ** In a nonlinear dynamic model, the consistency and
More informationAsymptotic distributions of nonparametric regression estimators for longitudinal or functional data
Journal of Multivariate Analysis 98 (2007) 40 56 www.elsevier.com/locate/jmva Asymptotic distriutions of nonparametric regression estimators for longitudinal or functional data Fang Yao Department of Statistics,
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 informationSurvival analysis in R
Survival analysis in R Niels Richard Hansen This note describes a few elementary aspects of practical analysis of survival data in R. For further information we refer to the book Introductory Statistics
More informationVariance Function Estimation in Multivariate Nonparametric Regression
Variance Function Estimation in Multivariate Nonparametric Regression T. Tony Cai 1, Michael Levine Lie Wang 1 Abstract Variance function estimation in multivariate nonparametric regression is considered
More informationNonparametric Estimation of Smooth Conditional Distributions
Nonparametric Estimation of Smooth Conditional Distriutions Bruce E. Hansen University of Wisconsin www.ssc.wisc.edu/~hansen May 24 Astract This paper considers nonparametric estimation of smooth conditional
More informationPhilosophy and Features of the mstate package
Introduction Mathematical theory Practice Discussion Philosophy and Features of the mstate package Liesbeth de Wreede, Hein Putter Department of Medical Statistics and Bioinformatics Leiden University
More informationIntroduction to Regression
Introduction to Regression p. 1/97 Introduction to Regression Chad Schafer cschafer@stat.cmu.edu Carnegie Mellon University Introduction to Regression p. 1/97 Acknowledgement Larry Wasserman, All of Nonparametric
More informationResearch Article Almost Periodic Viscosity Solutions of Nonlinear Parabolic Equations
Hindawi Publishing Corporation Boundary Value Problems Volume 29, Article ID 873526, 15 pages doi:1.1155/29/873526 Research Article Almost Periodic Viscosity Solutions of Nonlinear Parabolic Equations
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 Analysis: Weeks 2-3. Lu Tian and Richard Olshen Stanford University
Survival Analysis: Weeks 2-3 Lu Tian and Richard Olshen Stanford University 2 Kaplan-Meier(KM) Estimator Nonparametric estimation of the survival function S(t) = pr(t > t) The nonparametric estimation
More informationProduct-limit estimators of the survival function with left or right censored data
Product-limit estimators of the survival function with left or right censored data 1 CREST-ENSAI Campus de Ker-Lann Rue Blaise Pascal - BP 37203 35172 Bruz cedex, France (e-mail: patilea@ensai.fr) 2 Institut
More 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 informationAn idea how to solve some of the problems. diverges the same must hold for the original series. T 1 p T 1 p + 1 p 1 = 1. dt = lim
An idea how to solve some of the problems 5.2-2. (a) Does not converge: By multiplying across we get Hence 2k 2k 2 /2 k 2k2 k 2 /2 k 2 /2 2k 2k 2 /2 k. As the series diverges the same must hold for the
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 informationGeneralized Seasonal Tapered Block Bootstrap
Generalized Seasonal Tapered Block Bootstrap Anna E. Dudek 1, Efstathios Paparoditis 2, Dimitris N. Politis 3 Astract In this paper a new lock ootstrap method for periodic times series called Generalized
More informationOn the Multi-Dimensional Controller and Stopper Games
On the Multi-Dimensional Controller and Stopper Games Joint work with Yu-Jui Huang University of Michigan, Ann Arbor June 7, 2012 Outline Introduction 1 Introduction 2 3 4 5 Consider a zero-sum controller-and-stopper
More informationModel Specification Testing in Nonparametric and Semiparametric Time Series Econometrics. Jiti Gao
Model Specification Testing in Nonparametric and Semiparametric Time Series Econometrics Jiti Gao Department of Statistics School of Mathematics and Statistics The University of Western Australia Crawley
More informationIntroduction to Regression
Introduction to Regression Chad M. Schafer May 20, 2015 Outline General Concepts of Regression, Bias-Variance Tradeoff Linear Regression Nonparametric Procedures Cross Validation Local Polynomial Regression
More informationNonparametric Regression. Changliang Zou
Nonparametric Regression Institute of Statistics, Nankai University Email: nk.chlzou@gmail.com Smoothing parameter selection An overall measure of how well m h (x) performs in estimating m(x) over x (0,
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 informationProperties of an infinite dimensional EDS system : the Muller s ratchet
Properties of an infinite dimensional EDS system : the Muller s ratchet LATP June 5, 2011 A ratchet source : wikipedia Plan 1 Introduction : The model of Haigh 2 3 Hypothesis (Biological) : The population
More informationDESIGN-ADAPTIVE MINIMAX LOCAL LINEAR REGRESSION FOR LONGITUDINAL/CLUSTERED DATA
Statistica Sinica 18(2008), 515-534 DESIGN-ADAPTIVE MINIMAX LOCAL LINEAR REGRESSION FOR LONGITUDINAL/CLUSTERED DATA Kani Chen 1, Jianqing Fan 2 and Zhezhen Jin 3 1 Hong Kong University of Science and Technology,
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 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 informationIntroduction to Regression
Introduction to Regression Chad M. Schafer cschafer@stat.cmu.edu Carnegie Mellon University Introduction to Regression p. 1/100 Outline General Concepts of Regression, Bias-Variance Tradeoff Linear Regression
More informationAppendix. Proof of Theorem 1. Define. [ ˆΛ 0(D) ˆΛ 0(t) ˆΛ (t) ˆΛ. (0) t. X 0 n(t) = D t. and. 0(t) ˆΛ 0(0) g(t(d t)), 0 < t < D, t.
Appendix Proof of Theorem. Define [ ˆΛ (D) X n (t) = ˆΛ (t) D t ˆΛ (t) ˆΛ () g(t(d t)), t < t < D X n(t) = [ ˆΛ (D) ˆΛ (t) D t ˆΛ (t) ˆΛ () g(t(d t)), < t < D, t where ˆΛ (t) = log[exp( ˆΛ(t)) + ˆp/ˆp,
More informationSpiking problem in monotone regression : penalized residual sum of squares
Spiking prolem in monotone regression : penalized residual sum of squares Jayanta Kumar Pal 12 SAMSI, NC 27606, U.S.A. Astract We consider the estimation of a monotone regression at its end-point, where
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 informationLecture 6 PREDICTING SURVIVAL UNDER THE PH MODEL
Lecture 6 PREDICTING SURVIVAL UNDER THE PH MODEL The Cox PH model: λ(t Z) = λ 0 (t) exp(β Z). How do we estimate the survival probability, S z (t) = S(t Z) = P (T > t Z), for an individual with covariates
More informationNonparametric Model Construction
Nonparametric Model Construction Chapters 4 and 12 Stat 477 - Loss Models Chapters 4 and 12 (Stat 477) Nonparametric Model Construction Brian Hartman - BYU 1 / 28 Types of data Types of data For non-life
More information( t) Cox regression part 2. Outline: Recapitulation. Estimation of cumulative hazards and survival probabilites. Ørnulf Borgan
Outline: Cox regression part 2 Ørnulf Borgan Department of Mathematics University of Oslo Recapitulation Estimation of cumulative hazards and survival probabilites Assumptions for Cox regression and check
More 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 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 informationPractice Exam 1. (A) (B) (C) (D) (E) You are given the following data on loss sizes:
Practice Exam 1 1. Losses for an insurance coverage have the following cumulative distribution function: F(0) = 0 F(1,000) = 0.2 F(5,000) = 0.4 F(10,000) = 0.9 F(100,000) = 1 with linear interpolation
More 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 informationFixed-b Inference for Testing Structural Change in a Time Series Regression
Fixed- Inference for esting Structural Change in a ime Series Regression Cheol-Keun Cho Michigan State University imothy J. Vogelsang Michigan State University August 29, 24 Astract his paper addresses
More informationJournal of Inequalities in Pure and Applied Mathematics
Journal of Inequalities in Pure and Applied Mathematics ON A HYBRID FAMILY OF SUMMATION INTEGRAL TYPE OPERATORS VIJAY GUPTA AND ESRA ERKUŞ School of Applied Sciences Netaji Subhas Institute of Technology
More informationAsymptotic statistics using the Functional Delta Method
Quantiles, Order Statistics and L-Statsitics TU Kaiserslautern 15. Februar 2015 Motivation Functional The delta method introduced in chapter 3 is an useful technique to turn the weak convergence of random
More informationExample: Bipolar NRZ (non-return-to-zero) signaling
Baseand Data Transmission Data are sent without using a carrier signal Example: Bipolar NRZ (non-return-to-zero signaling is represented y is represented y T A -A T : it duration is represented y BT. Passand
More informationEMPIRICAL ENVELOPE MLE AND LR TESTS. Mai Zhou University of Kentucky
EMPIRICAL ENVELOPE MLE AND LR TESTS Mai Zhou University of Kentucky Summary We study in this paper some nonparametric inference problems where the nonparametric maximum likelihood estimator (NPMLE) are
More 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 informationMaximum Likelihood Estimation for the Proportional Odds Model with Random Effects
Maximum Likelihood Estimation for the Proportional Odds Model with Random Effects DONGLIN ZENG, D. Y. LIN, and GUOSHENG YIN In this article, we study the semiparametric proportional odds model with random
More informationSurvival Analysis for Case-Cohort Studies
Survival Analysis for ase-ohort Studies Petr Klášterecký Dept. of Probability and Mathematical Statistics, Faculty of Mathematics and Physics, harles University, Prague, zech Republic e-mail: petr.klasterecky@matfyz.cz
More informationLecture 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 informationA Comparison of Different Approaches to Nonparametric Inference for Subdistributions
A Comparison of Different Approaches to Nonparametric Inference for Subdistributions Johannes Mertsching Johannes.Mertsching@gmail.com Master Thesis Supervision: dr. Ronald B. Geskus prof. Chris A.J. Klaassen
More informationFrontier estimation based on extreme risk measures
Frontier estimation based on extreme risk measures by Jonathan EL METHNI in collaboration with Ste phane GIRARD & Laurent GARDES CMStatistics 2016 University of Seville December 2016 1 Risk measures 2
More informationSemiparametric Trending Panel Data Models with Cross-Sectional Dependence
he University of Adelaide School of Economics Research Paper No. 200-0 May 200 Semiparametric rending Panel Data Models with Cross-Sectional Dependence Jia Chen, Jiti Gao, and Degui Li Semiparametric rending
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 informationA note on the decomposition of number of life years lost according to causes of death
A note on the decomposition of number of life years lost according to causes of death Per Kragh Andersen Research Report 12/2 Department of Biostatistics University of Copenhagen A note on the decomposition
More informationUnderstanding product integration. A talk about teaching survival analysis.
Understanding product integration. A talk about teaching survival analysis. Jan Beyersmann, Arthur Allignol, Martin Schumacher. Freiburg, Germany DFG Research Unit FOR 534 jan@fdm.uni-freiburg.de It is
More informationSurvival Analysis: Counting Process and Martingale. Lu Tian and Richard Olshen Stanford University
Survival Analysis: Counting Process and Martingale Lu Tian and Richard Olshen Stanford University 1 Lebesgue-Stieltjes Integrals G( ) is a right-continuous step function having jumps at x 1, x 2,.. b f(x)dg(x)
More informationThe Bernstein and Nikolsky inequalities for trigonometric polynomials
he ernstein and Nikolsky ineualities for trigonometric polynomials Jordan ell jordanbell@gmailcom Department of Mathematics, University of oronto January 28, 2015 1 Introduction Let = R/2πZ For a function
More informationLAN property for sde s with additive fractional noise and continuous time observation
LAN property for sde s with additive fractional noise and continuous time observation Eulalia Nualart (Universitat Pompeu Fabra, Barcelona) joint work with Samy Tindel (Purdue University) Vlad s 6th birthday,
More informationStatistical Properties of Numerical Derivatives
Statistical Properties of Numerical Derivatives Han Hong, Aprajit Mahajan, and Denis Nekipelov Stanford University and UC Berkeley November 2010 1 / 63 Motivation Introduction Many models have objective
More information7 Day 3: Time Varying Parameter Models
7 Day 3: Time Varying Parameter Models References: 1. Durbin, J. and S.-J. Koopman (2001). Time Series Analysis by State Space Methods. Oxford University Press, Oxford 2. Koopman, S.-J., N. Shephard, and
More informationONLINE SUPPLEMENT FOR EVALUATING FACTOR PRICING MODELS USING HIGH FREQUENCY PANELS
ONLINE SUPPLEMEN FOR EVALUAING FACOR PRICING MODELS USING HIGH FREQUENCY PANELS YOOSOON CHANG, YONGOK CHOI, HWAGYUN KIM AND JOON Y. PARK hi online upplement contain ome ueful lemma and their proof and
More informationI N S T I T U T D E S T A T I S T I Q U E B I O S T A T I S T I Q U E E T S C I E N C E S A C T U A R I E L L E S (I S B A)
I N S T I T U T D E S T A T I S T I Q U E B I O S T A T I S T I Q U E E T S C I E N C E S A C T U A R I E L L E S (I S B A) UNIVERSITÉ CATHOLIQUE DE LOUVAIN D I S C U S S I O N P A P E R 2011/2 ESTIMATION
More informationNatural boundary and Zero distribution of random polynomials in smooth domains arxiv: v1 [math.pr] 2 Oct 2017
Natural boundary and Zero distribution of random polynomials in smooth domains arxiv:1710.00937v1 [math.pr] 2 Oct 2017 Igor Pritsker and Koushik Ramachandran Abstract We consider the zero distribution
More informationReflected Brownian Motion
Chapter 6 Reflected Brownian Motion Often we encounter Diffusions in regions with boundary. If the process can reach the boundary from the interior in finite time with positive probability we need to decide
More informationAdditive functionals of infinite-variance moving averages. Wei Biao Wu The University of Chicago TECHNICAL REPORT NO. 535
Additive functionals of infinite-variance moving averages Wei Biao Wu The University of Chicago TECHNICAL REPORT NO. 535 Departments of Statistics The University of Chicago Chicago, Illinois 60637 June
More informationMinimax Rate of Convergence for an Estimator of the Functional Component in a Semiparametric Multivariate Partially Linear Model.
Minimax Rate of Convergence for an Estimator of the Functional Component in a Semiparametric Multivariate Partially Linear Model By Michael Levine Purdue University Technical Report #14-03 Department of
More informationInfinitely Imbalanced Logistic Regression
p. 1/1 Infinitely Imbalanced Logistic Regression Art B. Owen Journal of Machine Learning Research, April 2007 Presenter: Ivo D. Shterev p. 2/1 Outline Motivation Introduction Numerical Examples Notation
More informationNonparametric Trending Regression with Cross-Sectional Dependence
Nonparametric Trending Regression with Cross-Sectional Dependence Peter M. Robinson London School of Economics January 5, 00 Abstract Panel data, whose series length T is large but whose cross-section
More informationdy dx dx = 7 1 x dx dy = 7 1 x dx e u du = 1 C = 0
1. = 6x = 6x = 6 x = 6 x x 2 y = 6 2 + C = 3x2 + C General solution: y = 3x 2 + C 3. = 7 x = 7 1 x = 7 1 x General solution: y = 7 ln x + C. = e.2x = e.2x = e.2x (u =.2x, du =.2) y = e u 1.2 du = 1 e u
More informationSupport Vector Hazard Regression (SVHR) for Predicting Survival Outcomes. Donglin Zeng, Department of Biostatistics, University of North Carolina
Support Vector Hazard Regression (SVHR) for Predicting Survival Outcomes Introduction Method Theoretical Results Simulation Studies Application Conclusions Introduction Introduction For survival data,
More informationGraphical Diagnostic for Mortality Data Modeling
Preprints www.preprints.org NOT PEER-REVIEWED Posted: 18 July 2016 Article Graphical Diagnostic for Mortality Data Modeling M.L. Gamiz,, M.D. Martinez-Miranda and R. Raya-Miranda Department of Statistics
More information