Estimation and Inference of Quantile Regression. for Survival Data under Biased Sampling
|
|
- Catherine Hicks
- 6 years ago
- Views:
Transcription
1 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 sampling scheme Here, we provide the justification of (11) which is introduced in Section 3. Observe that f T, (t, = 1 A, Z) = P T (t, t + dt), = 1 A, Z/dt = I(t A)f A,V (A, t A Z)S c (t A Z) = I(t A)f(t Z)S c (t A Z) 1 µ(z) = I(t A)P (T (t, t + dt), C > t Z)/dt 1 µ(z) 1 = I(t A)f T, (t, = 1 Z) µ(z), f T, (t, = A, Z) = P ( T (t, t + dt), = A, Z)/dt = I(t A) t f A,V (A, s A Z)ds g c (t A Z) 1 = I(t A) f(s Z)ds g c (t A Z) t µ(z) = I(t A)P T > t, C (t, t + dt) Z/dt 1 µ(z) 1 = I(t A)f T, (t, = Z) µ(z). 1
2 It follows that that v i (t) = f T, ( T i, i Z i ) f T, ( T i, i Z i ) f T, (t, 1 Z i ) f T, (t, 1 Z i) = I(A i t). S2 Proof of Theorems 1 and 2 For the biased samplings introduced in Sections 2.2 and 3, let f T (t Z) be the conditional density function of T, i.e., f T (t Z) = δ,1 f T, (t, δ Z), where f T, (t, δ Z) is defined in (8) of Section 2.2. For a vector a, let a 2 denote aa, and a denote the Euclidean norm of a. For b R p, define m(b) = E ZN(e Z b ), m(b) = E Zv(e Z b )Y (e Z b ), m n (b) = 1 n Z i N i (e Z i b ), n m n (b) = 1 n Z i v i (e Z i b )Y i (e Z i b ), n B(b) = E Z 2 f T, (e Z b, 1 Z) exp(z b), J(b) = E Z 2 v(e Z b )f T (e Z b Z) exp(z b). For d >, define the set B(d) as B(d) = b R p : inf m(b) m(β (τ)) d, τ (,τ u] where β (τ) is the true parameter value, τ (, τ u ]; and τ u satisfies Condition C4 below. 2
3 Furthermore, consider the setting in Section 3.2. For b R p, define m (b) = E Z ψ( T )N(e Z b ), m (b) = E Z ψ(e Z b )v(e Z b )Y (e Z b ), m n(b) = 1 n Z i ψ( n T i )N i (e Z i b ), m n(b) = 1 n Z i ψ(e Z i b )v i (e Z i b )Y i (e Z i b ), n B (b) = E [Z ψ(e Z b )]Z f T, (e Z b, 1 Z) exp(z b), J (b) = E [Z ψ(e Z b )]Z v(e Z b )f T (e Z b Z) exp(z b). We assume the following conditions for theoretical derivation. C1: Z is uniformly bounded, i.e., sup i Z i < and the matrix E(Z 2 ) is positive definite. C2: mβ(τ) is a Lipschitz continuous function of τ (, τ u ], and f T (t z) and f T, (t, 1 z) are bounded above and continuous uniformly in t and z. C3: There exists d > such that for b B(d ) and any Z, f T, (e Z b, 1 Z) > and J(b)B(b) 1 is uniformly bounded. C4: inf τ [τl,τ u] eigmin Bβ (τ) > for any τ l (, τ u ], where eigmin( ) denotes the minimum eigenvalue of a matrix. C5: The class of weight functions (Z, A, T, ) v(e Z b, A, T,, Z); b R d is Glivenko- Cantelli and uniformly bounded; and (Z, A, T, ) v(e Z b, A, T,, Z); b B(d ) belongs to a Donsker class. C6: The weight function ψ(t) is positive, bounded and differentiable; (Z, A, T, ) ψ(e Z b, A, T,, Z); b R d is Glivenko-Cantelli and (Z, A, T, ) ψ(e Z b, A, T,, Z); b B(d ) belongs to a Donsker class. Matrix W (β, τ) is nonsingular, inf t [τl,τ] eigmin B β (t) 3
4 W (β, τ) 1 B β (t) > for any τ l (, τ] and J (b)b (b) W (β, τ) 1 B (b) 1 is uniformly bounded for b B(d ). Conditions C1 C4 are mild assumptions concerning the covariates Z, the underlying regression quantile parameter process β(τ), and the density functions associated with the observed data ( T, ). The boundedness assumption of covariates in C1 are often assumed in survival models although it can be further relaxed with extra technical complexity. For Conditions C2 and C3, the boundedness and uniform continuity assumption of f T (t z) and f T, (t, 1 z) is reasonable in many biased sampling problems. For instance, it is satisfied for case-cohort study if the density functions of survival and censoring times are bounded and uniformly continuous. Similarly, for the length-biased sampling introduced in Section 3, this is satisfied if the density functions of T and C are bounded and uniformly continuous. Under these conditions, it can been verified that m(b) and m(b) are differentiable and B(b) = m(b)/ b and J(b) = m(b)/ b are well defined. As in Peng and Huang (28), Condition C4 is assumed to ensure the identifiability of β (τ). The regularity of the weight function assumed in Condition C5 is satisfied for all the sampling schemes introduced in Section 3. In particular, for the case-cohort type designs (Examples 3 and 4), the weight function v does not depend on β and satisfies C5. For length-biased sampling with weight function in the form (14), v(e Z b ), b B(d ) is a VC class, which implies C5 (Theorems and in van der Vaart and Wellner, 1996). Condition C6 is similar as conditions C3-C5 and is satisfied for many weight functions. It is assumed to ensure the asymptotic normality of the efficient estimator. Proof of Theorem 1 The proof follows closely that in Peng and Huang (28) and we only present the key steps. We should note that, however, their results cannot be directly applied due to the associated random weight functions in the estimating equation. Let α (τ) = mβ (τ), ˆα(τ) = mˆβ(τ), and A(d) = m(b) : b B(d). Following the argument in Peng and Huang (28), m is a one-to-one mapping from B(d ) to A(d ). There 4
5 exists an inverse function κ : B(d ) A(d ) such that κm(b) = b, for any b B(d ). Uniform Convergence Let ξ n,k = m n ˆβ(τ k ) k m n ˆβ(u)dH(u). According to the estimating equation taking the general form of (1), we have sup k ξ n,k = O(1) sup i Z i /n = O(n 1 ). From mβ (τ k ) k mβ (s)dh(s) =, we have the following decomposition: mˆβ(τ k ) mβ (τ k ) [ ] k 1 = m n ˆβ(τ k ) mˆβ(τ k ) + j= j+1 k [ ] + m n ˆβ(s) mˆβ(s) dh(s) + ξ n,k. τ j [ ] mˆβ(s) mβ (s) dh(s) (S1) Under the condition that Z is uniformly bounded, Z i N(e Z i b ); b R p and Z i Y i (e Z i b ); b R p are Glivenko-Cantelli. By Condition C5 and Theorem 3 in van der Vaart and Wellner (2), Z i v i (e Z i b )Y i (e Z i b ); b R p is also Glivenko-Cantelli. Then, we have the following results: sup b R p m(b) m n (b) and sup b R p m(b) m n (b) almost surely. This implies that the first and third terms in (S1) are ignorable. Denote c n, := sup k [m n ˆβ(τ k ) mˆβ(τ k )]+ k [ m nˆβ(s) mˆβ(s)]dh(s), and it follows that c n, = o p (1). Under condition C2, there exists c 1 such that mβ (τ) mβ (τ ) < c 1 τ τ for any τ, τ (, τ u ]. For = τ τ < τ 1, since mˆβ() =, we have sup τ τ<τ 1 mˆβ(τ) mβ (τ) = sup τ τ<τ 1 mβ (τ) b n, := c 1 S L(n). Then for large enough n, we know b n, < d and ˆβ n (τ) B(d ) for τ [τ, τ 1 ). Thus, we can write sup τ τ<τ 1 mˆβ(τ) mβ (τ) = sup τ τ<τ 1 m[κ ˆα(τ)] m[κα(τ)]. Under Condition C3, there exists α(τ) between α (τ) and ˆα(τ) such that sup mˆβ(τ) mβ (τ) = sup J[κ α(τ)]b[κ α(τ)] 1 α(τ) α (τ) c 2 b n,, τ τ<τ 1 τ τ<τ 1 5
6 where c 2 is some constant. Next consider τ 1 τ < τ 2. We have sup mˆβ(τ) mβ (τ) mˆβ(τ 1 ) mβ (τ 1 ) + sup mβ (τ 1 ) mβ (τ). τ 1 τ<τ 2 τ 1 τ<τ 2 We know sup τ1 τ<τ 2 mβ (τ 1 ) mβ (τ) c 1 S L(n). Further from the decomposition (S1) that gives an upper bound for mˆβ(τ 1 ) mβ (τ 1 ), we have: sup mˆβ(τ) mβ (τ) b n,1 := c n, + c 3 n 1 + c 2 b n, (1 τ u ) 1 S L(n) + c 1 S L(n), τ 1 τ<τ 2 where c 3 is some big constant. Note that b 1,n and it is smaller than d for large enough n. This implies ˆβ n (τ) B(d ) for τ [τ 1, τ 2 ). Then an induction argument similarly to the Proof of Theorem 1 in Peng and Huang (28) gives k 1 sup mˆβ(τ) mβ (τ) b n,k := c n, +c 3 n 1 +c 2 b n,i (1 τ u ) 1 S L(n) +c 1 S L(n). τ k τ τ k+1 This implies that sup τ τ τ u mˆβ(τ) mβ (τ) p and for τ τ τ u and large enough n, ˆβ(τ) B(d ). Thus sup τ ˆα(τ) α (τ) = sup τ mˆβ(τ) mβ (τ) p. Take Taylor s expansion of mˆβ(τ) around α (τ), and we have sup ˆβ(τ) β (τ) τ [τ l,τ u] sup Bβ (τ) 1 ˆα(τ) α (τ) + τ [τ l,τ u] sup ɛ n(τ), τ [τ l,τ u] where ɛ n(τ) is the remainder term of the Taylor expansion and sup τ [τl,τ u] ɛ n(τ). By Condition C4, we have sup τ [τl,τ u] Bβ (τ) 1 ˆα(τ) α (τ) = O ˆα(τ) α (τ) and this implies the uniform consistency. 6
7 Weak Convergence Following Lemma B.1 of Peng and Huang (28), we have sup n 1/2 τ (,τ u] n [ Z i N i (e Z ˆβ(τ) i ) N i (e Z i β (τ) ) n 1/2 mˆβ(τ) mβ (τ)] = o p (1) (S2) and sup n 1/2 τ (,τ u] n ( Z i I ˆβ(τ)) ( ) Ti e Z i I Ti e Z i β (τ) n 1/2 [ mˆβ(τ) mβ (τ)] = o p (1). (S3) In addition, for the grid size chosen in Theorem 1, S n (ˆβ, τ) = o p (1) uniformly in τ (, τ u ]. Then, the following equations hold uniformly in τ (, τ u ], ] S n (β, τ) =n [mˆβ(τ) 1/2 mβ (τ) [ ] n 1/2 mˆβ(u) mβ (u) dh(u) + o p (1) ] =n [mˆβ(τ) 1/2 mβ (τ) n [ 1/2 Jβ (u)bβ (u) 1 + o p (1) ] [ ] mˆβ(u) mβ (u) dh(u) + o p (1). A discrete version of the above equation ] [ ] S n (β, τ) =n [m 1/2 n ˆβ(τ) m n β (τ) n 1/2 m n ˆβ(u) m n β (u) dh(u) + o p (1) 7
8 leads to the following recursive formula: S n (β, τ k ) S n (β, τ k 1 ) ] [ = n [m 1/2 n ˆβ(τ k ) m n β (τ k ) n 1/2 n m 1/2 n (ˆβ(τ k 1 )) m n (β (τ k 1 )) = Bβ (τ k )n 1/2 ˆβ(τ k ) β (τ k ) ] m n ˆβ(τ k 1 ) m n β (τ k 1 ) H(τ k ) H(τ k 1 ) [Bβ (τ k 1 ) + Jβ (τ k 1 )H(τ k ) H(τ k 1 )]n 1/2 ˆβ(τ k 1 ) β (τ k 1 ) + o p (1). The above recursive equation gives the following approximation result Bβ (τ k )n 1/2 ˆβ(τ k ) β (τ k ) = S n (β, τ k ) S n (β, τ k 1 ) [I + Jβ (τ k 1 )B 1 β (τ k 1 )H(τ k ) H(τ k 1 )] S n (β, τ k 1 ) S n (β, τ k 2 ) k [I + Jβ (τ h 1 )B 1 β (τ h 1 )H(τ h ) H(τ h 1 )] S n (β, τ 1 ) S n (β, τ ) h=2 + o p (1). (S4) Using the product integration theory (Gill and Johansen, 199; Andersen et al., 1993 II.6), we can write the above decomposition as n 1/2 [mˆβ(τ) mβ (τ)] = φ S n (β, τ) + o p (1), (S5) where φ is functional defined as follows. For any g G = g : [, τ u ] R p, g is left-continuous with right limit, g() = and product integral I(s, t) = π u (s,t] [I p +Jβ (u)bβ (u) 1 ]dh(u) 8
9 with π the product integral notation, φ(g)(τ) is defined as φ(g)(τ) = I(s, τ)dg(s). (S6) Note that Z i N i (e Z iβ (τ) ), τ [τ l, τ u ] is a Donsker class and uniformly bounded. From Condition C5, the class Z i v i (e Z i β (s) )Y i (e Z i β (s) ), s [τ l, τ u ] is also Donsker (see Example in Section 2.1 of van der Vaart and Wellner, 1996). Since v i(e Z i β (s) )Y i (e Z i β (s) )dh(s) is Libschitz in τ, from the permanence properties of the Lipschitz transformation of the Donsker class (Theorem 2.1.6, van der Vaart and Wellner, 1996), we can show that Z i N i (e Z iβ (τ) ) Z i v i (e Z i β (s) )Y i (e Z i β (s) )dh(s), τ [τ l, τ u ] is a Donsker class. Then the Donsker theorem implies that S n (β, τ) converges weakly to a Gaussian process U(τ) living on τ [τ l, τ u ] with mean and covariance matrix Σ(s, t), where Σ(s, t) = Eu i (s)u i (t) with u i (τ) = Z i N i (e Z iβ (τ) ) Z i v i(e Z i β (s) )Y i (e Z i β (s) )dh(s). Furthermore, since φ is a linear operator, φ S n (β, τ) converges weakly to φ U(τ) and φ U(τ) is also a Gaussian process (Römisch, 25). From (S5), apply Taylor s expansion and we have Bβ (τ)n 1/2 ˆβ(τ) β (τ) = φ S n (β, τ) + o p (1). Then n 1/2 ˆβ(τ) β (τ) converges weakly to the Gaussian process Bβ (τ) 1 φ U(τ) with covariance matrix Bβ (τ) 1 Σ [Bβ (τ) 1 ], where Σ (τ) denotes the limiting covariance matrix of φ S n (β, τ). Proof of Theorem 2 Next we prove the weak convergence of nˆβ eff (τ) β (τ). From the above definitions, we can write η(β, τ) and η (β, τ k ) in (18) and (2) as η(β, τ) = n 1/2 m nβ(τ) n 1/2 m nβ(u)dh(u) 9
10 and k 1 η (β, τ k ) = n 1/2 m nβ(τ k ) n 1/2 j= m nˆβ eff (τ j )H(τ j+1 ) H(τ j ). By the proposed estimation method for ˆβ eff (τ), the sequential estimator ˆβ(τ k ), 1 k L, minimizes n 1 η (β, τ k ) W (ˆβ int, τ) 1 η (β, τ k ). Since ˆβ int is taken as a consistent estimator of β, it can be shown that W (ˆβ int, τ) 1 W (β, τ) 1 in probability. Further note that B (b) = m (b) b and J (b) = m (b) b. As in the proof of Theorem 1, m n(b) and m n(b) converge to m (b) and m (b) uniformly in probability; then the sequential estimators ˆβ eff (τ k ) satisfy n 1/2 B ˆβ eff (τ k ) W (β, τ) 1 η (ˆβ eff, τ k ) = o p (1). Then from a similar argument as in the proof of Theorem 1 for the consistency of ˆβ( ), we have that ˆβ eff (t) is uniformly consistent in t τ. Similarly as Lemma B.1 of Peng and Huang (28), we have for β(t) such that sup t β(t) β (t) in probability, sup [m n β(t) m nβ (t)] [m β(t) m β (t)] = o p ( n 1/2 ) t (,τ] and Then we can write sup [ m n β(t) m nβ (t)] [ m β(t) m β (t)] = o p (n 1/2 ). t (,τ] η ( β, τ k ) = n 1/2 B β (τ k ) β(τ k ) β (τ k ) + o p (1) k 1 +n 1/2 m nβ (τ k ) n 1/2 1 j= m nˆβ eff (τ j )H(τ j+1 ) H(τ j ).
11 This implies the sequential estimator ˆβ eff (τ k ) satisfies B β (τ k ) W (β, τ) 1 η (ˆβ eff, τ k ) = o p (1), uniformly in probability, and furthermore, uniformly in t, B β (t)w (β, τ) 1 η(ˆβ eff, t) = o p (1). This equation plays a similar role in the current proof as that of the unweighted estimating equation S n (ˆβ, t) = o p (1) in the proof of Theorem 1 for the weak convergence of ˆβ. The weak convergence of ˆβ eff then follows from a similar argument. S3 Resampling method An alternative Resampling method A commonly used resampling method is the perturbationbased approach proposed by Jin et al. (23) for the AFT model; its modified version for the quantile regression was adopted in Peng and Huang (28) under unbiased sampling. This approach can be easily extended to the biased sampling case. In particular, consider the following stochastic perturbation of the observed estimating equations S n (β, τ) = n 1/2 n ζ i Z i N i (e Z i β(τ) ) v i (e Z i β(s) )Y i (e Z i β(s) )dh(s) =. (S7) where ζ i s are i.i.d. from a known nonnegative distribution with mean 1 and variance 1, such as the exponential distribution with rate 1. The above estimation can be implemented in the same fashion as discussed in Section 2.3. Conditional on the observed data, we independently generate the variates (ζ 1,..., ζ n ) M b times, with M b being a big number, and obtain the corresponding M b estimates ˆβ (τ) by solving (S7). The resampling method described above is consistent. Following the proof of Theorem 3 below, conditional on Z i, i, T i n, we have S n (ˆβ, τ) has the same asymptotic distribution 11
12 as S n (β, τ) in probability; moreover, for τ uniformly in (, τ u ], we have n 1/2 β (τ) ˆβ(τ) = B(β (τ)) 1 φ S n (ˆβ, τ) + o p (1), (S8) where function φ is defined as in (S6). Therefore, conditional on the observed data, n 1/2 ˆβ (τ) ˆβ(τ) converges weakly to the same limiting process of n 1/2 ˆβ(τ) β (τ) for τ [τ l, τ u ], where τ l (, τ u ). Proof of Theorem 3 We first show that, conditional on the observed data, Sn (ˆβ, τ) has the same asymptotic distribution as S n (β, τ) in probability. Consider the weighed summation defined in (S7) S n (β, τ) = n 1/2 n ζ i Z i N i (e Z i β(τ) ) v i (e Z i β(s) )Y i (e Z i β(s) )dh(s). Since E(ζ i ) = 1 and V ar(ζ i ) = 1, a similar argument as in the proof of Theorem 1 gives that sup S n (ˆβ, τ) = o p (1), τ (,τ u] where ˆβ is the solution of S n (β, τ) = using the proposed sequential algorithm. Because sup τ (,τu] S n (ˆβ, τ) = o p (1) as shown in the proof of Theorem 1, for τ uniformly in (, τ u ], S n (ˆβ, τ) can be expressed as S n (ˆβ, τ) = S n (ˆβ, τ) S n (ˆβ, τ) + o p (1) = n 1/2 n (1 ζ i )û i (τ) + o p (1), where û i (τ) = Z i N(e Z i ˆβ(τ) ) Z i v i (e Z ˆβ(s) i, T i, i )Y i (e Z ˆβ(s) i )dh(s). 12
13 Then we have for any s, t [τ l, τ u ], [ n n E n 1/2 (1 ζ i )û i (s) n 1/2 Zi (1 ζ i )û i (t), i, T i n n = n 1 û i (s)û i (t) Σ(s, t), ] in probability as n, where Σ(, ) is the corresponding limiting covariance function of S n (β, ) as defined in the proof of Theorem 1. Therefore, by a similar argument as in Lin, Wei, and Ying (1993), conditional on Z i, i, T i n, Sn (ˆβ, τ) has the same limiting covariance matrix and asymptotic distribution as S n (β, τ) in probability. From the decomposition in (S4), then we can show that conditional on the data, φ n S n (ˆβ, τ) has the same asymptotic distribution as φ n ( S n (β, τ)) in probability. Furthermore, since n 1/2 [mˆβ(τ) mβ (τ)] = φ n S n (β, τ)+o p (1), we obtain that conditional on the data, φ n S n (ˆβ, τ) converges weakly to the limiting distribution of n 1/2 [mˆβ(τ) mβ (τ)] in probability. Following from Zeng and Lin (28), we use the perturbation estimators ˆB and Ĵ to estimate B and J. The consistency property can be established as follows. Following Lemma B.1 of Peng and Huang (28), for β(t) such that sup t β(t) β (t) in probability, the following result holds uniformly in τ n 1/2 [m n β(τ) m n β (τ)] n 1/2 [m β(τ) mβ (τ)] = o p (1). This implies that n 1/2 [m n β(τ) m n β (τ)] = Bβ (τ)n 1/2 β(τ) β (τ) + o p (1). Therefore, the above approximation holds for β(τ) = ˆβ(τ) + n 1/2 γ, where γ follows a p- dimensional multivariate normal distribution with zero mean and identity covariance matrix. 13
14 For M generated γ s, let γ be the M p matrix (γ 1,, γ M ) and m be the m-dimensional vector [n 1/2 m n ˆβ(τ) + n 1/2 γ i n 1/2 m n ˆβ(τ); i = 1,, M]. Then we have Bb (τ) = (γ γ) 1 γ m + o p (1). Note that (γ γ) 1 exist with probability 1 for M > p. Therefore, the slope matrix estimator given by regressing the perturbed values n 1/2 m n ˆβ(τ) + n 1/2 γ i on γ i, i.e., ˆBˆβ(τ) = (γ γ) 1 γ m, is consistent. Similarly we have the consistency of Ĵ. This implies that ˆBˆβ(τ) 1 φ n S n (ˆβ, τ) has the same asymptotic distribution as n 1/2 ˆβ(τ) β (τ). This validates the proposed resampling method. References Andersen, P. K., Borgan, Ø., Gill, R. D., and Keiding, N. (1993), Statistical Models Based on Counting Processes, Springer, New York. Gill, R. D. and Johansen, S. (199), A Survey of Product-Integration with a View Toward Application in Survival Analysis, The Annals of Statistics, 18, Jin, Z., Lin, D. Y., Wei, L. J., and Ying, Z. (23), Rank-Based Inference for the Accelerated Failure Time Model, Biometrika, 9, Lin, D. Y., Wei, L.-J., and Ying, Z. (1993), Checking the Cox Model with Cumulative Sums of Martingale-Based Residuals, Biometrika, 8, Peng, L. and Huang, Y. (28), Survival Analysis with Quantile Regression Models, Journal of the American Statistical Association, 13,
15 Römisch, W. (25), Delta Method, Infinite Dimensional, Encyclopedia of Statistical Sciences. van der Vaart, A. and Wellner, J. A. (2), Preservation theorems for Glivenko-Cantelli and uniform Glivenko-Cantelli classes, in High dimensional probability II, Springer, pp van der Vaart, A. W. and Wellner, J. A. (1996), Weak Convergence, Springer. 15
Efficiency 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 Doubly Censored Data with Applications to Cystic Fibrosis Studies. Technical Report July 1, 2010
Quantile Regression for Doubly Censored Data with Applications to Cystic Fibrosis Studies by Shuang Ji 1, Limin Peng 1, Yu Cheng 2, and HuiChuan Lai 3 Technical Report 1-3 July 1, 21 1 Department of Biostatistics
More informationIntroduction to Empirical Processes and Semiparametric Inference Lecture 02: Overview Continued
Introduction to Empirical Processes and Semiparametric Inference Lecture 02: Overview Continued Michael R. Kosorok, Ph.D. Professor and Chair of Biostatistics Professor of Statistics and Operations Research
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 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 informationQuantile Processes for Semi and Nonparametric Regression
Quantile Processes for Semi and Nonparametric Regression Shih-Kang Chao Department of Statistics Purdue University IMS-APRM 2016 A joint work with Stanislav Volgushev and Guang Cheng Quantile Response
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 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 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 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 informationSmoothed and Corrected Score Approach to Censored Quantile Regression With Measurement Errors
Smoothed and Corrected Score Approach to Censored Quantile Regression With Measurement Errors Yuanshan Wu, Yanyuan Ma, and Guosheng Yin Abstract Censored quantile regression is an important alternative
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 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 informationarxiv:submit/ [math.st] 6 May 2011
A Continuous Mapping Theorem for the Smallest Argmax Functional arxiv:submit/0243372 [math.st] 6 May 2011 Emilio Seijo and Bodhisattva Sen Columbia University Abstract This paper introduces a version of
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 informationSOME CONVERSE LIMIT THEOREMS FOR EXCHANGEABLE BOOTSTRAPS
SOME CONVERSE LIMIT THEOREMS OR EXCHANGEABLE BOOTSTRAPS Jon A. Wellner University of Washington The bootstrap Glivenko-Cantelli and bootstrap Donsker theorems of Giné and Zinn (990) contain both necessary
More informationSupplementary Appendix: Difference-in-Differences with. Multiple Time Periods and an Application on the Minimum. Wage and Employment
Supplementary Appendix: Difference-in-Differences with Multiple Time Periods and an Application on the Minimum Wage and mployment Brantly Callaway Pedro H. C. Sant Anna August 3, 208 This supplementary
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 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 informationasymptotic normality of nonparametric M-estimators with applications to hypothesis testing for panel count data
asymptotic normality of nonparametric M-estimators with applications to hypothesis testing for panel count data Xingqiu Zhao and Ying Zhang The Hong Kong Polytechnic University and Indiana University Abstract:
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 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 informationSimple and Fast Overidentified Rank Estimation for Right-Censored Length-Biased Data and Backward Recurrence Time
Biometrics 74, 77 85 March 2018 DOI: 10.1111/biom.12727 Simple and Fast Overidentified Rank Estimation for Right-Censored Length-Biased Data and Backward Recurrence Time Yifei Sun, 1, * Kwun Chuen Gary
More informationA Resampling Method on Pivotal Estimating Functions
A Resampling Method on Pivotal Estimating Functions Kun Nie Biostat 277,Winter 2004 March 17, 2004 Outline Introduction A General Resampling Method Examples - Quantile Regression -Rank Regression -Simulation
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 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 informationQuantile Regression for Dynamic Panel Data
Quantile Regression for Dynamic Panel Data Antonio Galvao 1 1 Department of Economics University of Illinois NASM Econometric Society 2008 June 22nd 2008 Panel Data Panel data allows the possibility of
More informationMaximum likelihood estimation in the proportional hazards cure model
AISM DOI 1.17/s1463-7-12-x Maximum likelihood estimation in the proportional hazards cure model Wenbin Lu Received: 12 January 26 / Revised: 16 November 26 The Institute of Statistical Mathematics, Tokyo
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 informationRegression Calibration in Semiparametric Accelerated Failure Time Models
Biometrics 66, 405 414 June 2010 DOI: 10.1111/j.1541-0420.2009.01295.x Regression Calibration in Semiparametric Accelerated Failure Time Models Menggang Yu 1, and Bin Nan 2 1 Department of Medicine, Division
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 informationQuantile Regression for Recurrent Gap Time Data
Biometrics 000, 1 21 DOI: 000 000 0000 Quantile Regression for Recurrent Gap Time Data Xianghua Luo 1,, Chiung-Yu Huang 2, and Lan Wang 3 1 Division of Biostatistics, School of Public Health, University
More informationQuantile Regression: Inference
Quantile Regression: Inference Roger Koenker University of Illinois, Urbana-Champaign 5th RMetrics Workshop, Meielisalp: 28 June 2011 Roger Koenker (UIUC) Introduction Meielisalp: 28.6.2011 1 / 29 Inference
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 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 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 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 informationSTAT 992 Paper Review: Sure Independence Screening in Generalized Linear Models with NP-Dimensionality J.Fan and R.Song
STAT 992 Paper Review: Sure Independence Screening in Generalized Linear Models with NP-Dimensionality J.Fan and R.Song Presenter: Jiwei Zhao Department of Statistics University of Wisconsin Madison April
More informationNETWORK-REGULARIZED HIGH-DIMENSIONAL COX REGRESSION FOR ANALYSIS OF GENOMIC DATA
Statistica Sinica 213): Supplement NETWORK-REGULARIZED HIGH-DIMENSIONAL COX REGRESSION FOR ANALYSIS OF GENOMIC DATA Hokeun Sun 1, Wei Lin 2, Rui Feng 2 and Hongzhe Li 2 1 Columbia University and 2 University
More informationImproving linear quantile regression for
Improving linear quantile regression for replicated data arxiv:1901.0369v1 [stat.ap] 16 Jan 2019 Kaushik Jana 1 and Debasis Sengupta 2 1 Imperial College London, UK 2 Indian Statistical Institute, Kolkata,
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 informationStability of optimization problems with stochastic dominance constraints
Stability of optimization problems with stochastic dominance constraints D. Dentcheva and W. Römisch Stevens Institute of Technology, Hoboken Humboldt-University Berlin www.math.hu-berlin.de/~romisch SIAM
More informationExpecting the Unexpected: Uniform Quantile Regression Bands with an application to Investor Sentiments
Expecting the Unexpected: Uniform Bands with an application to Investor Sentiments Boston University November 16, 2016 Econometric Analysis of Heterogeneity in Financial Markets Using s Chapter 1: Expecting
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 informationUniversity of California San Diego and Stanford University and
First International Workshop on Functional and Operatorial Statistics. Toulouse, June 19-21, 2008 K-sample Subsampling Dimitris N. olitis andjoseph.romano University of California San Diego and Stanford
More informationVerifying Regularity Conditions for Logit-Normal GLMM
Verifying Regularity Conditions for Logit-Normal GLMM Yun Ju Sung Charles J. Geyer January 10, 2006 In this note we verify the conditions of the theorems in Sung and Geyer (submitted) for the Logit-Normal
More informationQuantile Regression: Inference
Quantile Regression: Inference Roger Koenker University of Illinois, Urbana-Champaign Aarhus: 21 June 2010 Roger Koenker (UIUC) Introduction Aarhus: 21.6.2010 1 / 28 Inference for Quantile Regression Asymptotics
More informationAsymptotic normality of Powell s kernel estimator
Asymptotic normality of Powell s kernel estimator Kengo Kato This version: July 3, 2009 Abstract In this paper, we establish asymptotic normality of Powell s kernel estimator for the asymptotic covariance
More informationSEMIPARAMETRIC LIKELIHOOD RATIO INFERENCE. By S. A. Murphy 1 and A. W. van der Vaart Pennsylvania State University and Free University Amsterdam
The Annals of Statistics 1997, Vol. 25, No. 4, 1471 159 SEMIPARAMETRIC LIKELIHOOD RATIO INFERENCE By S. A. Murphy 1 and A. W. van der Vaart Pennsylvania State University and Free University Amsterdam Likelihood
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 informationOn Estimation of Partially Linear Transformation. Models
On Estimation of Partially Linear Transformation Models Wenbin Lu and Hao Helen Zhang Authors Footnote: Wenbin Lu is Associate Professor (E-mail: wlu4@stat.ncsu.edu) and Hao Helen Zhang is Associate Professor
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 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 informationHomework # , Spring Due 14 May Convergence of the empirical CDF, uniform samples
Homework #3 36-754, Spring 27 Due 14 May 27 1 Convergence of the empirical CDF, uniform samples In this problem and the next, X i are IID samples on the real line, with cumulative distribution function
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 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 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 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 informationAnalysis of least absolute deviation
Analysis of least absolute deviation By KANI CHEN Department of Mathematics, HKUST, Kowloon, Hong Kong makchen@ust.hk ZHILIANG YING Department of Statistics, Columbia University, NY, NY, 10027, U.S.A.
More informationTESTINGGOODNESSOFFITINTHECOX AALEN MODEL
ROBUST 24 c JČMF 24 TESTINGGOODNESSOFFITINTHECOX AALEN MODEL David Kraus Keywords: Counting process, Cox Aalen model, goodness-of-fit, martingale, residual, survival analysis. Abstract: The Cox Aalen regression
More 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 informationAsymptotics of minimax stochastic programs
Asymptotics of minimax stochastic programs Alexander Shapiro Abstract. We discuss in this paper asymptotics of the sample average approximation (SAA) of the optimal value of a minimax stochastic programming
More informationMedian 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 informationMultivariate Least Weighted Squares (MLWS)
() Stochastic Modelling in Economics and Finance 2 Supervisor : Prof. RNDr. Jan Ámos Víšek, CSc. Petr Jonáš 23 rd February 2015 Contents 1 2 3 4 5 6 7 1 1 Introduction 2 Algorithm 3 Proof of consistency
More informationThe Contraction Method on C([0, 1]) and Donsker s Theorem
The Contraction Method on C([0, 1]) and Donsker s Theorem Henning Sulzbach J. W. Goethe-Universität Frankfurt a. M. YEP VII Probability, random trees and algorithms Eindhoven, March 12, 2010 joint work
More informationLecture 11: Regression Methods I (Linear Regression)
Lecture 11: Regression Methods I (Linear Regression) Fall, 2017 1 / 40 Outline Linear Model Introduction 1 Regression: Supervised Learning with Continuous Responses 2 Linear Models and Multiple Linear
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 informationEfficient Estimation of Censored Linear Regression Model
2 3 4 5 6 7 8 9 2 3 4 5 6 7 8 9 2 2 22 23 24 25 26 27 28 29 3 3 32 33 34 35 36 37 38 39 4 4 42 43 44 45 46 47 48 Biometrika (2), xx, x, pp. 4 C 28 Biometrika Trust Printed in Great Britain Efficient Estimation
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 informationLecture 11: Regression Methods I (Linear Regression)
Lecture 11: Regression Methods I (Linear Regression) 1 / 43 Outline 1 Regression: Supervised Learning with Continuous Responses 2 Linear Models and Multiple Linear Regression Ordinary Least Squares Statistical
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 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 informationASYMPTOTIC PROPERTIES OF BRIDGE ESTIMATORS IN SPARSE HIGH-DIMENSIONAL REGRESSION MODELS
ASYMPTOTIC PROPERTIES OF BRIDGE ESTIMATORS IN SPARSE HIGH-DIMENSIONAL REGRESSION MODELS Jian Huang 1, Joel L. Horowitz 2, and Shuangge Ma 3 1 Department of Statistics and Actuarial Science, University
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 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 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 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 informationM- and Z- theorems; GMM and Empirical Likelihood Wellner; 5/13/98, 1/26/07, 5/08/09, 6/14/2010
M- and Z- theorems; GMM and Empirical Likelihood Wellner; 5/13/98, 1/26/07, 5/08/09, 6/14/2010 Z-theorems: Notation and Context Suppose that Θ R k, and that Ψ n : Θ R k, random maps Ψ : Θ R k, deterministic
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 informationMonotonicity and Aging Properties of Random Sums
Monotonicity and Aging Properties of Random Sums Jun Cai and Gordon E. Willmot Department of Statistics and Actuarial Science University of Waterloo Waterloo, Ontario Canada N2L 3G1 E-mail: jcai@uwaterloo.ca,
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 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 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 informationTesting linearity against threshold effects: uniform inference in quantile regression
Ann Inst Stat Math (2014) 66:413 439 DOI 10.1007/s10463-013-0418-9 Testing linearity against threshold effects: uniform inference in quantile regression Antonio F. Galvao Kengo Kato Gabriel Montes-Rojas
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 informationAsymmetric least squares estimation and testing
Asymmetric least squares estimation and testing Whitney Newey and James Powell Princeton University and University of Wisconsin-Madison January 27, 2012 Outline ALS estimators Large sample properties Asymptotic
More informationASYMPTOTIC PROPERTIES OF BRIDGE ESTIMATORS IN SPARSE HIGH-DIMENSIONAL REGRESSION MODELS
The Annals of Statistics 2008, Vol. 36, No. 2, 587 613 DOI: 10.1214/009053607000000875 Institute of Mathematical Statistics, 2008 ASYMPTOTIC PROPERTIES OF BRIDGE ESTIMATORS IN SPARSE HIGH-DIMENSIONAL REGRESSION
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 informationNonconcave Penalized Likelihood with A Diverging Number of Parameters
Nonconcave Penalized Likelihood with A Diverging Number of Parameters Jianqing Fan and Heng Peng Presenter: Jiale Xu March 12, 2010 Jianqing Fan and Heng Peng Presenter: JialeNonconcave Xu () Penalized
More information5.1 Approximation of convex functions
Chapter 5 Convexity Very rough draft 5.1 Approximation of convex functions Convexity::approximation bound Expand Convexity simplifies arguments from Chapter 3. Reasons: local minima are global minima and
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 informationlarge number of i.i.d. observations from P. For concreteness, suppose
1 Subsampling Suppose X i, i = 1,..., n is an i.i.d. sequence of random variables with distribution P. Let θ(p ) be some real-valued parameter of interest, and let ˆθ n = ˆθ n (X 1,..., X n ) be some estimate
More informationUNIVERSITÄT POTSDAM Institut für Mathematik
UNIVERSITÄT POTSDAM Institut für Mathematik Testing the Acceleration Function in Life Time Models Hannelore Liero Matthias Liero Mathematische Statistik und Wahrscheinlichkeitstheorie Universität Potsdam
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 informationStochastic Optimization with Inequality Constraints Using Simultaneous Perturbations and Penalty Functions
International Journal of Control Vol. 00, No. 00, January 2007, 1 10 Stochastic Optimization with Inequality Constraints Using Simultaneous Perturbations and Penalty Functions I-JENG WANG and JAMES C.
More informationClosest Moment Estimation under General Conditions
Closest Moment Estimation under General Conditions Chirok Han and Robert de Jong January 28, 2002 Abstract This paper considers Closest Moment (CM) estimation with a general distance function, and avoids
More informationDelta Theorem in the Age of High Dimensions
Delta Theorem in the Age of High Dimensions Mehmet Caner Department of Economics Ohio State University December 15, 2016 Abstract We provide a new version of delta theorem, that takes into account of high
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 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 informationCENSORED QUANTILE REGRESSION WITH COVARIATE MEASUREMENT ERRORS
Statistica Sinica 21 211, 949-971 CENSORED QUANTILE REGRESSION WITH COVARIATE MEASUREMENT ERRORS Yanyuan Ma and Guosheng Yin Texas A&M University and The University of Hong Kong Abstract: We study censored
More information