The Uniform Asymptotics of the Overshoot of a Random Walk with Light-Tailed Increments

Let {S _n : n ≥ 0} be a random walk with light-tailed increments and negative drift, and let τ(x) be the first time when the random walk crosses a given level x ≥ 0. Tang (2007 Tang , Q. ( 2007 ). The overshoot of a random walk with negative drift . Statist. Probab. Lett. 77 : 158 – 165 .[Crossref], [Web of Science ®] , [Google Scholar]) obtained the asymptotics of P(S _τ(x) ? x > y, τ(x) < ∞) as x → ∞, which is uniform for y ≥ f(x) for any positive function f(x) → ∞ as x → ∞. In this article, the uniform asymptotics of P(S _τ(x) ? x > y, τ(x) < ∞) as x → ∞, for 0 ≤ y ≤ N for any positive number N will be given. Using the above two results, the uniform asymptotics of P(S _τ(x) ? x > y, τ(x) < ∞) as x → ∞, for y ≥ 0, is presented.     

Bivariate Aging Properties under Archimedean Dependence Structures

Let X  = (X, Y) be a pair of lifetimes whose dependence structure is described by an Archimedean survival copula, and let X _t  = [(X ? t, Y ? t) | X > t, Y > t] denotes the corresponding pair of residual lifetimes after time t ≥ 0. Multivariate aging notions, defined by means of stochastic comparisons between X and X _t , with t ≥ 0, were studied in Pellerey (2008 Pellerey , F. ( 2008 ). On univariate and bivariate aging for dependent lifetimes with Archimedean survival copulas . Kybernetika 44 : 795 – 806 .[Web of Science ®] , [Google Scholar]), who considered pairs of lifetimes having the same marginal distribution. Here, we present the generalizations of his results, considering both stochastic comparisons between X _t and X _t+s for all t, s ≥ 0 and the case of dependent lifetimes having different distributions. Comparisons between two different pairs of residual lifetimes, at any time t ≥ 0, are discussed as well.     

Reliability Properties of Reversed Residual Lifetime

Abstract

If the random variable X denotes the lifetime (X ≥ 0, with probability one) of a unit, then the random variable X _t  = (t ? X|X ≤ t), for a fixed t > 0, is known as `time since failure', which is analogous to the residual lifetime random variable used in reliability and survival analysis. The reversed hazard rate function, which is related to the random variable X _t , has received the attention of many researchers in the recent past [(cf. Shaked, M., Shanthikumar, J. G., 1994 Shaked, M. and Shanthikumar, J. G. 1994. Stochastic Orders and Their Applications New York: Academic Press.  [Google Scholar]). Stochastic Orders and Their Applications. New York: Academic Press]. In this paper, we define some new classes of distributions based on the random variable X _t and study their interrelations. We also define a new ordering based on the mean of the random variable X_t and establish its relationship with the reversed hazard rate ordering.     

On asymptotics of multivariate integrals with applications to records

Let { X _n ,n≥1} be a sequence of iid. Gaussian random vectors in R ^d , d≥2, with nonsingular distribution function F. In this paper the asymptotics for the sequence of integrals I _F,n (G _n )?n∫ _{R ^d} G _n ^n?1( X ) dF( X ) is considered with G _n some distribution function on R ^d . In the case G _n =F the integral I _F,n (F)/n is the probability that a record occurs in X ₁,…, X _n at index n. [1] Gnedin, A.V. 1998. Records from a Multivariate Normal Sample. Statist. Probab. Lett., 39: 11–15. [Crossref], [Web of Science ®] , [Google Scholar] obtained lower and upper asymptotic bounds for this case, whereas [2] Ledford, W.A. and Twan, A.J. 1998. On the Tail Concomitant Behaviour for Extremes. Adv. Appl. Probab., 30: 197–215. [Crossref], [Web of Science ®] , [Google Scholar] showed the rate of convergence if d=2. In this paper we derive the exact rate of convergence of I _F,n (G _n ) for d≥2 under some restrictions on the distribution function G _n . Some related results for multivariate Gaussian tails are discussed also.     

Simultaneous Canonization of Linear Models

Canonical form plays a similar role in linear models to spectral decomposition in matrix analysis. Let X = (X ₁,…, X _n )′ be a random vector with expectation Aβ and the variance–covariance matrix σV, where V is positive definite and let rank(A) = r. Then there exists a nonsingular linear transformation from X to T = (T ₁,…, T _n )′, such that ET _i  = η _i , for i = 1,…, r and zero for i > r, while cov(T _i , T _j ) = δ _ij σ. This canonical form, introduced by Ko?odziejczyk (1935 Ko?odziejczyk , S. ( 1935 ). On an important class of statistical hypotheses . Biometrika 27 : 161 – 190 .[Crossref] , [Google Scholar]), was used, among others, by Scheffé (1959 Scheffé , H. ( 1959 ). Analysis of Variance . New York : Wiley . [Google Scholar]) and by Lehmann (1959, 1986 Lehmann , E. L. (1959, 1986 ). Testing Statistical Hypotheses . New York : Wiley . [Google Scholar]). This technique is extended here for arbitrary (possibly singular) V and for simultaneous canonization of two models of this type.     

Smoothed Mann–Whitney–Wilcoxon Procedure for Two-Sample Location Problem

This study is mainly concerned with estimating a shift parameter in the two-sample location problem. The proposed Smoothed Mann–Whitney–Wilcoxon method smooths the empirical distribution functions of each sample by using convolution technique, and it replaces unknown distribution functions F(x) and G(x ? Δ₀) with the new smoothed distribution functions F _s (x) and G _s (x ? Δ₀), respectively. The unknown shift parameter Δ₀ is estimated by solving the gradient function S _n (Δ) with respect to an arbitrary variable Δ. The asymptotic properties of the new estimator are established under some conditions that are similar to the Generalized Wilcoxon procedure proposed by Anderson and Hettmansperger (1996 Anderson , G. F. , Hettmansperger , T. P. ( 1996 ). Generalized Wilcoxon methods for the one and two-sample location models . In: Brunner , E. , Denker , M. , eds. Research Developments in Probability and Statistics: Festschrift in Honor of Madan L. Puri on the Occasion of his 65th Birthday . Zeist, The Netherlands : VSP BV , pp. 303 – 317 . [Google Scholar]). Some of these properties are asymptotic normality, asymptotic level confidence interval, and hypothesis testing for Δ₀. Asymptotic relative efficiency of the proposed method with respect to the least squares, Generalized Wilcoxon and Hodges and Lehmann (1963 Hodges , J. L. , Lehmann , E. L. ( 1963 ). Estimates of location based on rank tests . Ann. Mathemat. Statist. 34 : 598 – 611 .[Crossref] , [Google Scholar]) procedures are also calculated under the contaminated normal model.     

New Generalizations of Cauchy Distribution

The generalized skew-normal distribution introduced by Balakrishnan (2002 Balakrishnan , N. ( 2002 ). Discussion on ‘Skew multivariate models related to hidden truncation and/or selective reporting’ by B. C. Arnold and R. J. Beaver . Test 11 : 37 – 39 .[Web of Science ®] , [Google Scholar]) is used to obtain new generalizations of univariate Cauchy distribution with two parameters, denoted by GC _{m, n} (a, b) with m and n non-negative integer numbers and a, b ∈ R. For cases (m, n) = (1, 2), (m, n) = (2, 1), (m, n) = (0, 3) and (m, n) = (3, 0) explicit forms of the density functions are derived and compared to previous generalizations of Cauchy and skew-Cauchy distributions.     

Linear Combinations,Products and Ratios of Simplicial or Spherical Variates

The density level sets of the two types of measures under consideration are l _{2, p} -circles with p = 1 and p = 2, respectively. The intersection-percentage function (ipf) of such a measure reflects the percentages which the level set corresponding to the p-radius r shares for each r > 0 with a set to be measured. The geometric measure representation formulae in Richter (2009 Richter , W.-D. (2009). Continuous l _{n, p} -symmetric distributions. Lithuanian Mathemat. J. 49:93–108.[Crossref], [Web of Science ®] , [Google Scholar]) is based upon these ipf's and will be used here for evaluating exact cdf's and pdf's for the linear combination, the product, and the ratio of the components of two-dimensional simplicial or spherically distributed random vectors.     

A self-normalized central limit theorem for a ρ-mixing stationary sequence

We investigate a self-normalized central limit theorem for a ρ-mixing stationary sequence {X_i, i ? 1} of random variables such that L(x) ? E(X²₁I{|X₁| ? x}) is a slowly varying function as x → ∞. The results obtained generalize the results of Gine, Gotze, and Mason (1997) and Mason (2005 Mason, D. M. 2005. The asymptotic distribution of self-normalized triangular arrays. Journal of Theoretical Probability 18 (4):853–70.[Crossref], [Web of Science ®] , [Google Scholar]) to ρ-mixing sequences.     

Natural Exponential Families and Generalized Hypergeometric Measures

Let ν be a positive Borel measure on ?ⁿ and _pF_q(a₁,…, a_p; b₁,…, b_q; s) be a generalized hypergeometric series. We define a generalized hypergeometric measure, μ_p,q := _pF_q(a₁,…, a_p; b₁,…, b_q;ν), as a series of convolution powers of the measure ν, and we investigate classes of probability distributions which are expressible as such a measure. We show that the Kemp (1968 Kemp , A. W. ( 1968 ). A wide class of discrete distributions and the associated differential equations . Sankhyā, Ser. A 30 : 401 – 410 . [Google Scholar]) family of distributions is an example of μ_p,q in which ν is a Dirac measure on ?. For the case in which ν is a Dirac measure on ?ⁿ, we relate μ_p,q to the diagonal natural exponential families classified by Bar-Lev et al. (1994 Bar-Lev , S. K. , Bshouty , D. , Enis , P. , Letac , G. , Lu , I. , Richards , D. ( 1994 ). The diagonal natural exponential families on ? ⁿ and their classification . J. Theoret. Probab. 7 : 883 – 929 .[Crossref] , [Google Scholar]). For p < q, we show that certain measures μ_p,q can be expressed as the convolution of a sequence of independent multi-dimensional Bernoulli trials. For p = q, q + 1, we show that the measures μ_p,q are mixture measures with the Dufresne and Poisson-stopped-sum probability distributions as their mixing measures.     

On Optimal Estimating Functions in the Presence of Nuisance Parameters

When there is only one interesting parameter θ₁ and one nuisance parameter θ₂, Godambe and Thompson (1974 Godambe , V. P. , Thompson , M. E. ( 1974 ). Estimating equations in the presence of a nuisance parameter . Ann. Statist. 2 : 568 – 571 .[Crossref], [Web of Science ®] , [Google Scholar]) showed that the optimal estimating function for θ₁ essentially is a linear function of the θ₁-score, the square of the θ₂-score, and the derivative of θ₂-score with respect to θ₂. Mukhopadhyay (2000b) generalized this result to m nuisance parameters. Mukhopadhyay (2000 Mukhopadhyay , P. ( 2000 ). On some lower bounds for the variance of an estimating function . Int. J. Math. Statist. Sci. 9 ( 2 ). [Google Scholar] 2002a Mukhopadhyay , P. ( 2002a ). On estimating functions in the presence of a nuisance parameter . Commun. Statist. Theor. Mem. 31 ( 1 ): 31 – 36 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar] b Mukhopadhyay , P. ( 2002b ): Some lower bounds on variance of estimating functions . J. Statist. Res. 36 ( 2 ): 189 – 197 . [Google Scholar]) obtained lower bounds to the variance of regular estimating functions in the presence of nuisance parameters. Taking cues from these results we propose a method of finding optimal estimating function for θ₁ by taking the multiple regression equation on θ₁ score and Bhattacharyya's (1946 Bhattacharyya , A. ( 1946 ). On some analogues of the amount of information and their use in statistical estimation . Sankhya A 8 : 1 – 14 . [Google Scholar]) scores with respect to θ₂. The result is extended to the case of m nuisance parameters.     

Quantile-based Chernoff distance for truncated random variables

Several probability distributions such as power-Pareto distribution (see Gilchrist 2000 Gilchrist, W. 2000. Statistical modelling with quantile functions. Boca Raton, FL: Chapman and Hall/CRC.[Crossref] , [Google Scholar] and Hankin and Lee 2006 Hankin, R. K. S., and A. Lee. 2006. A new family of non-negative distributions. Australian and New Zealand Journal of Statistics 48:67–78.[Crossref], [Web of Science ®] , [Google Scholar]), various forms of lambda distributions (see Ramberg and Schmeiser 1974 Ramberg, J. S., and B. W. Schmeiser. 1974. An appropriate method for generating asymmetric random variables. Communications of the ACM 17:78–82.[Crossref], [Web of Science ®] , [Google Scholar] and Freimer et al. 1988 Freimer, M., S. Mudholkar, G. Kollia, and C. T. Lin. 1988. A study of the generalized lambda family. Communications in Statistics - Theory and Methods 17:3547–67.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), Govindarajulu distribution (see Nair, Sankaran, and Vineshkumar 2012 Nair, U. N., P. G. Sankaran, and B. Vineshkumar. 2012. The Govindarajulu distribution: some properties and applications. Communications in Statistics—Theory and Methods 41:4391–406.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), etc., do not have manageable distribution functions, though they have tractable quantile functions. Hence, analytical study of the properties of Chernoff distance of two random variables associated with these distributions via traditional distribution function-based tool becomes difficult. To make this simple, in this paper, we introduce quantile-based Chernoff distance for (left or right) truncated random variables and study its various properties. Some useful bounds as well as characterization results are obtained.     

D-Optimal Chemical Balance Weighing Designs with N ≡ 0 (mod 4) and 3 Objects

In this article, the problem of estimation of the individual weights of three objects using a chemical balance weighing design is considered. We use the criterion of D-optimality. We assume that the covariance matrix of errors is the matrix of first-order autoregressive process. Such problems were discussed in Li and Yang (2005 Li , C. H. , Yang , S. Y. ( 2005 ). On a conjecture in D-optimal designs with n ≡ 0 (mod 4) . Lin. Alg. Applic. 400 : 279 – 290 .[Crossref], [Web of Science ®] , [Google Scholar]) and also in Yeh and Lo Huang (2005 Yeh , H. G. , Lo Huang , M. N. ( 2005 ). On exact D-optimal designs with 2 two-level factors and n autocorrelated observations . Metrika 61 : 261 – 275 .[Crossref], [Web of Science ®] , [Google Scholar]). We present some results of D-optimal designs in certain class of designs with the design matrix X  ∈ M _n×3(±1) such that each column of matrix X has at least one 1 and one ?1.     

Circular First- and Second-Order Balanced Repeated Measurements Designs

Sharma (1977 Sharma , V. K. ( 1977 ). Change-over designs with complete balance for first and second order residual effect . Canad. J. Statist. 5 : 121 – 132 .[Crossref] , [Google Scholar]) and Aggarwal et al. (2006 Aggarwal , M. L. , Deng , L.-Y. , Jha , M. K. ( 2006 ). Balanced residual treatment effects designs of first and second order . Statist. Probab. Lett. 76 : 597 – 600 .[Crossref], [Web of Science ®] , [Google Scholar]) considered non circular construction of first- and second-order balanced repeated measurements designs. Sharma et al. (2002 Sharma , V. K. , Varghese , C. , Jaggi , S. ( 2002 ). On optimality of change-over designs balanced for first and second order residual effects . Metron 60 : 153 – 162 . [Google Scholar]) constructed circular first- and second-order balanced repeated measurements designs only for a class with parameters (v, p = 3n, n = v ²) and also showed its universal optimality. In this article, we consider circular construction of first- and second-order balanced repeated measurements designs and strongly balanced repeated measurements designs by using the method of cyclic shifts. Some new circular designs with parameters (v, p, n) for cases p = v, p < v and p > v are given.     

An Index of Local Sensitivity to Nonignorability for a Pseudolikelihood Method

For regression analysis of data with non response, sensitivity analysis is usually recommended. An index of local sensitivity to non ignorability (ISNI) (Troxel et al., 2004 Troxel , A. B. , Ma , G. , Heitjan , D. F. (2004). An index of local sensitivity to nonignorability. Statistica Sinica 14:1221–1237.[Web of Science ®] , [Google Scholar]) was derived to detect the sensitivity of maximum likelihood estimates to small departures from ignorability. However, ISNI requires specification of a parametric model for the missing-data mechanism. In this article, a local sensitivity index for a pseudolikelihood (PL) method that does not require specification of the mechanism is proposed. For bivariate data (x, y), when the non response mechanism is an arbitrary function of x + λy, this new index is defined as the first derivative of the PL estimate with respect to λ at λ = 0. The closed form was derived for normal regression data when the density function of the predictor x approximated by a kernel estimator in the PL method. The utility of this new local sensitivity index was illustrated through application on one dataset.     

Peaks-Over-Threshold Modeling Under Random Censoring

Recently, the topic of extreme value under random censoring has been considered. Different estimators for the index have been proposed (see Beirlant et al., 2007 Beirlant , J. , Guillou , A. , Dierckx , G. , Fils-Villetard , A. ( 2007 ). Estimation of the extreme value index and extreme quantiles under random censoring . Extremes 10 : 151 – 174 .[Crossref] , [Google Scholar]). All of them are constructed as the classical estimators (without censoring) divided by the proportion of non censored observations above a certain threshold. Their asymptotic normality was established by Einmahl et al. (2008 Einmahl , J. H. J. , Fils-Villetard , A. , Guillou , A. ( 2008 ). Statistics of extremes under random censoring . Bernoulli 14 ( 1 ): 207 – 227 . [Google Scholar]). An alternative approach consists of using the Peaks-Over-Threshold method (Balkema and de Haan, 1974 Balkema , A. , de Haan , L. ( 1974 ). Residual life at great age . Ann. Probab. 2 : 792 – 804 .[Crossref], [Web of Science ®] , [Google Scholar]; Smith, 1987 Smith , R. L. ( 1987 ). Estimating tails of probability distributions . Ann. Statist. 15 : 1174 – 1207 .[Crossref], [Web of Science ®] , [Google Scholar]) and to adapt the likelihood to the context of censoring. This leads to ML-estimators whose asymptotic properties are still unknown. The aim of this article is to propose one-step approximations, based on the Newton-Raphson algorithm. Based on a small simulation study, the one-step estimators are shown to be close approximations to the ML-estimators. Also, the asymptotic normality of the one-step estimators has been established, whereas in case of the ML-estimators it is still an open problem. The proof of our result, whose approach is new in the Peaks-Over-Threshold context, is in the spirit of Lehmann's theory (1991 Lehmann , E. L. ( 1991 ). Theory of Point Estimation . Pacific Grove , CA : Wadsworth & Brooks/Cole Advanced Books & Software .[Crossref] , [Google Scholar]).     

Nonparametric Testing of an Exclusion Restriction in Quantile Regression

Using the framework proposed by Bickel et al. (2006 Bickel , P. J. , Ritov , Y. , Stoker , T. ( 2006 ). Tailor-made tests for goodness-of-fit to semiparametric hypotheses . Ann. Stat. 34 ( 2 ): 721 – 741 . [Google Scholar]), we provide a score-based testing method to check the exclusion restriction in quantile regression, i.e., H: ν_α(Y|U, V) = ν_α(Y|U) w.p.1, where ν_α denotes the αth (0 < α < 1) quantile. A subsampling method is suggested to acquire the critical values and justified. The tests are all found to be consistent against fixed alternatives and have discriminating power against local alternatives at root-n scale. We address this particular problem as a representative among a wide family of semiparametric model checking problems. The methodology can be carried over to other goodness-of-fit testing of semiparametric models, possibly involve non smooth functions.     

Preservation of Some Aging Properties and Stochastic Orders by Weighted Distributions

By assuming that a random variable X possesses an aging property, we provide conditions under which the corresponding weighted version X ₁, with weight function w ₁(·), would also possess this aging property. Similarly, by assuming that two random variables X and Y are ordered with respect to a stochastic order, we provide conditions under which the corresponding weighted versions X ₁ and Y ₂, with weight functions w ₁(·) and w ₂(·), respectively, preserve this stochastic ordering. We also point out fallacies in the similar results claimed by Jain et al. (1989 Jain , K. , Singh , H. , Bagai , I. ( 1989 ). Relations for reliability measures of weighted distributions . Commun. Statist. Theory Meth. 18 : 4393 – 4412 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) and Bartoszewicz and Skolimowska (2006 Bartoszewicz , J. , Skolimowska , M. ( 2006 ). Preservation of classes of life distribution and stochastic orders under weighting . Statist. Probab. Lett. 76 : 587 – 596 .[Crossref] , [Google Scholar]) and correct them.     

Regression Analysis with Covariates Missing at Random: A Piece-wise Nonparametric Model for Missing Covariates

Statistical analysis for the regression model f _β(y | x, z) with missing values in the covariate vector X requires modeling of the covariate distribution g(x | z). Likelihood methods, including Ibrahim (1990 Ibrahim , J. G. ( 1990 ). Incomplete data in generalized linear models . J. Amer. Statist. Assoc. 85 : 765 – 769 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), Chen (2004 Chen , H. Y. (2004). Nonparametric and semiparametric models for missing covariates in parametric regression. J. Amer. Statist. Assoc. 99:1176–1189.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), and Zhao (2005 Zhao , Y. ( 2005 ). Design and Efficient Estimation in Regression Analysis with Missing Data in Two-Phase Studies. Ph.D. thesis , University of Waterloo . [Google Scholar]), need either X or Z to be discrete. This article considers extending the likelihood methods to deal with cases where both X and Z may be continuous. We propose modeling the covariate distribution g(x | z) using a piece-wise nonparametric model, then a maximum likelihood estimate (MLE) of β can be computed following the maximum likelihood estimating procedure of Chen (2004 Chen , H. Y. (2004). Nonparametric and semiparametric models for missing covariates in parametric regression. J. Amer. Statist. Assoc. 99:1176–1189.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) or Zhao (2005 Zhao , Y. ( 2005 ). Design and Efficient Estimation in Regression Analysis with Missing Data in Two-Phase Studies. Ph.D. thesis , University of Waterloo . [Google Scholar]). The resulting estimation method is easy to implement and the asymptotic properties of the MLE follow under certain conditions. Extensive simulation studies for different models indicate that the proposed method is acceptable for practical implementation. A real data example is used to illustrate the method.