Estimation of parameters for the truncated exponential distribution

We consider the right truncated exponential distribution where the truncation point is unknown and show that the ML equation has a unique solution over an extended parameter space. In the case of the estimation of the truncation point T we show that the asymptotic distribution of the MLE is not centered at T. A modified MLE is introduced which outperforms all other considered estimators including the minimum variance unbiased estimator. Asymptotic as well as small sample properties of different estimators are investigated and compared. The truncated exponential distribution has an increasing failure rate, ideally suited for use as a survival distribution for biological and industrial data.     

Simultaneous marginal survival estimators when doubly censored data is present

A doubly censoring scheme occurs when the lifetimes T being measured, from a well-known time origin, are exactly observed within a window [L, R] of observational time and are otherwise censored either from above (right-censored observations) or below (left-censored observations). Sample data consists on the pairs (U, δ) where U = min{R, max{T, L}} and δ indicates whether T is exactly observed (δ = 0), right-censored (δ = 1) or left-censored (δ = −1). We are interested in the estimation of the marginal behaviour of the three random variables T, L and R based on the observed pairs (U, δ). We propose new nonparametric simultaneous marginal estimators [^(S)]_T, [^(S)]_L{\hat S_{T}, \hat S_{L}} and [^(S)]_R{\hat S_{R}} for the survival functions of T, L and R, respectively, by means of an inverse-probability-of-censoring approach. The proposed estimators [^(S)]_T, [^(S)]_L{\hat S_{T}, \hat S_{L}} and [^(S)]_R{\hat S_{R}} are not computationally intensive, generalize the empirical survival estimator and reduce to the Kaplan-Meier estimator in the absence of left-censored data. Furthermore, [^(S)]_T{\hat S_{T}} is equivalent to a self-consistent estimator, is uniformly strongly consistent and asymptotically normal. The method is illustrated with data from a cohort of drug users recruited in a detoxification program in Badalona (Spain). For these data we estimate the survival function for the elapsed time from starting IV-drugs to AIDS diagnosis, as well as the potential follow-up time. A simulation study is discussed to assess the performance of the three survival estimators for moderate sample sizes and different censoring levels.     

Multivariate normal distribution approaches for dependently truncated data

Many statistical methods for truncated data rely on the independence assumption regarding the truncation variable. In many application studies, however, the dependence between a variable X of interest and its truncation variable L plays a fundamental role in modeling data structure. For truncated data, typical interest is in estimating the marginal distributions of (L, X) and often in examining the degree of the dependence between X and L. To relax the independence assumption, we present a method of fitting a parametric model on (L, X), which can easily incorporate the dependence structure on the truncation mechanisms. Focusing on a specific example for the bivariate normal distribution, the score equations and Fisher information matrix are provided. A robust procedure based on the bivariate t-distribution is also considered. Simulations are performed to examine finite-sample performances of the proposed method. Extension of the proposed method to doubly truncated data is briefly discussed.     

A note on inference for P(X < Y) for right truncated exponentially distributed data

In this paper, a likelihood based analysis is developed and applied to obtain confidence intervals and p values for the stress-strength reliability R  =  P(X  <  Y) with right truncated exponentially distributed data. The proposed method is based on theory given in Fraser et al. (Biometrika 86:249–264, 1999) which involves implicit but appropriate conditioning and marginalization. Monte Carlo simulations are used to illustrate the accuracy of the proposed method.     

On the explosion rate of maximum-bias functions

The maxbias function B_T() contains much information about the robustness properties of the estimate T. This function satisfies B_T(0)=0 and B_T()<for all 0<<whereis the breakdown point of T. Hampel (1974)pioneered the study of the limiting behaviour of B_T(?) where ? → 0. He computed and optimized the rate γ at which B_T(?) approaches 0 when ? → 0. This rate is now called the contamination sensitivity of T, and constitutes one of the cornerstones of the theory of robustness. We show that much can also be learned from the study of the limiting behaviour of B_T(?) when ? → ?*. A new robustness measure, called the relative explosion rate, can be obtained by studying the limiting relative maxbias behaviour of two extimates when approaches their common breakdown point ?*. Like the contamination sensitivity, the relative explosion rate can be readily derived from the estimate's score function. General formulae are given for M-estimates of scale and S-, MM- and τ-estimates of regression. We also show that the maxbias behaviour for large ? is largely determined by the curvature of the estimate's score function near zero. This motivates our definition and study of the local order of a score function.     

On the mean L1-error in the heteroscedastic deconvolution problem with compactly supported noises

We study the heteroscedastic deconvolution problem when random noises have compactly supported densities. In this context, the Fourier transforms of the densities can vanish on the real line. We propose a truncated type of estimator for target density and derive the convergence rate of the mean L¹-error uniformly over a class of target densities. A lower bound for the mean L¹-error is also established. Some simulations will be given to illustrate the performance of the proposed estimator.     

Using a Truncated C p Statistic for Variable Selection in Multiple Linear Regression

In multiple linear regression analysis each lower-dimensional subspace L of a known linear subspace M of ? ⁿ corresponds to a non empty subset of the columns of the regressor matrix. For a fixed subspace L, the C _p statistic is an unbiased estimator of the mean square error if the projection of the response vector onto L is used to estimate the expected response. In this article, we consider two truncated versions of the C _p statistic that can also be used to estimate this mean square error. The C _p statistic and its truncated versions are compared in two example data sets, illustrating that use of the truncated versions may result in models different from those selected by standard C _p .     

On generalized exponential transformations for proportions

ABSTRACT

Transformation of the response is a popular method to meet the usual assumptions of statistical methods based on linear models such as ANOVA and t-test. In this paper, we introduce new families of transformations for proportions or percentage data. Most of the transformations for proportions require 0 < x < 1 (where x denotes the proportion), which is often not the case in real data. The proposed families of transformations allow x = 0 and x = 1. We study the properties of the proposed transformations, as well as the performance in achieving normality and homoscedasticity. We analyze three real data sets to empirically show how the new transformation performs in meeting the usual assumptions. A simulation study is also performed to study the behavior of new families of transformations.     

Nonparametric analysis of interval censored and doubly truncated data

Doubly truncated data appear in a number of applications, including astronomy and survival analysis. For double-truncated data, the lifetime T is observable only when U≤T≤V, where U and V are the left-truncated and right-truncated time, respectively. In some situations, the lifetime T also suffers interval censoring. Using the EM algorithm of Turnbull [The empirical distribution function with arbitrarily grouped censored and truncated data, J. R. Stat. Soc. Ser. B 38 (1976), pp. 290–295] and iterative convex minorant algorithm [P. Groeneboom and J.A. Wellner, Information Bounds and Nonparametric Maximum Likelihood Estimation, Birkhäuser, Basel, 1992], we study the performance of the nonparametric maximum-likelihood estimates (NPMLEs) of the distribution function of T. Simulation results indicate that the NPMLE performs adequately for the finite sample.     

Robustness of weighted L p–depth and L p–median

Summary: L ^p –norm weighted depth functions are introduced and the local and global robustness of these weighted L ^p –depth functions and their induced multivariate medians are investigated via influence function and finite sample breakdown point. To study the global robustness of depth functions, a notion of finite sample breakdown point is introduced. The weighted L ^p –depth functions turn out to have the same low breakdown point as some other popular depth functions. Their influence functions are also unbounded. On the other hand, the weighted L ^p –depth induced medians are globally robust with the highest possible breakdown point for any reasonable estimator. The weighted L ^p –medians are also locally robust with bounded influence functions for suitable weight functions. Unlike other existing depth functions and multivariate medians, the weighted L ^p depth and medians are easy to calculate in high dimensions. The price for this advantage is the lack of affine invariance and equivariance of the weighted L ^p depth and medians, respectively.*The author thanks the referees for their very insightful and constructive comments and suggestions which led to corrections and substantial improvements. Supported in part by NSF Grants DMS-0071976 and DMS-0134628.     

Detecting linear sequences and subsequences

Greenwood [1946. The statistical study of infectious diseases (with discussion). J. Roy. Statist. Soc. 109, 85–110], using an _L2$L_2$ distance, and others have addressed the question of detecting a too-linear fit of the occurrence times _T0<_T1<?<_Tn$T_0<T_1<?<T_n$ of a sequence of random events. Two convenient distances are introduced here, then applied to the more challenging problem of detecting too-linear subsequences, where the multiple subsequence effect must be taken into account. Two interpretations of “linear subsequence” are considered.     

Shelf-life and its estimation in drug stability studies

One important property of any drug product is its stability over time. Drug stability studies are routinely carried out in the pharmaceutical industry in order to measure the degradation of an active pharmaceutical ingredient of a drug product. One important study objective is to estimate the shelf-life of the drug; the estimated shelf-life is required by the US Food and Drug Administration to be printed on the package label of the drug. This involves a suitable definition of the true shelf-life and the construction of an appropriate estimate of the true shelf-life. In this paper, the true shelf-life T_β is defined as the time point at which 100β% of all the individual dosage units (e.g. tablets) of the drug have the active ingredient content no less than the lowest acceptable limit L, where β and L are prespecified constants. The value of T_β depends on the parameters of the assumed degradation model of the active ingredient content and so is unknown. A lower confidence bound T?_β for T_β is then provided and used as the estimated shelf-life of the drug.     

Flexible modelling of simultaneously interval censored and truncated time‐to‐event data

This paper deals with the analysis of data from a HET‐CAM^VT experiment. From a statistical perspective, such data yield many challenges. First of all, the data are typically time‐to‐event like data, which are at the same time interval censored and right truncated. In addition, one has to cope with overdispersion as well as clustering. Traditional analysis approaches ignore overdispersion and clustering and summarize the data into a continuous score that can be analysed using simple linear models. In this paper, a novel combined frailty model is developed that simultaneously captures all of the aforementioned statistical challenges posed by the data. Copyright © 2015 John Wiley & Sons, Ltd.     

Lp (1 < p ⩽ 2) solutions of one-dimensional BSDEs whose generator is weakly monotonic in y and non-Lipschitz in z

ABSTRACT

In this paper, we start with establishing the existence of a minimal (maximal) L^p (1 < p ? 2) solution to a one-dimensional backward stochastic differential equation (BSDE), where the generator g satisfies a p-order weak monotonicity condition together with a general growth condition in y and a linear growth condition in z. Then, we propose and prove a comparison theorem of L^p (1 < p ? 2) solutions to one-dimensional BSDEs with q-order (1 ? q < p) weak monotonicity and uniform continuity generators. As a consequence, an existence and uniqueness result of L^p (1 < p ? 2) solutions is also given for BSDEs whose generator g is q-order (1 ? q < p) weakly monotonic with a general growth in y and uniformly continuous in z.     

Some asymptotic properties of point estimation with n-dimensional fuzzy data

Some asymptotic properties of point estimation with n-dimensional fuzzy data with respect to a special L ₂-metric ρ are investigated in this article. It is shown that the collection of all n-dimensional fuzzy data endowed with the ρ-metric is a complete and separable space. Some criterions for point estimation in such fuzzy environments are proposed, and the sample mean, variance and covariance with n-dimensional fuzzy data under these criteria are further studied.     

Jackknifing L-estimates

Linear functions of order statistics (“L-estimates”) of the form T_n =under jackknifing are investigated. This paper proves that with suitable conditions on the function J, the jackknifed version T_n of the L-estimate T_n has the same limit distribution as T_n. It is also shown that the jackknife estimate of the asymptotic variance of n^1/2 is consistent. Furthermore, the Berry-Esséen rate associated with asymptotic normality, and a law of the iterated logarithm of a class of jackknife L-estimates, are characterized.     

Multi - sample test of equal gamma distribution scale parameters in presence of unknown common shape parameter

Shiue and Bain proposed an approximate F statistic for testing equality of two gamma distribution scale parameters in presence of a common and unknown shape parameter. By generalizing Shiue and Bain's statistic we develop a new statistic for testing equality of L >= 2 gamma distribution scale parameters. We derive the distribution of the new statistic ESP for L = 2 and equal sample size situation. For other situations distribution of ESP is not known and test based on the ESP statistic has to be performed by using simulated critical values. We also derive a C(α) statistic CML and develop a likelihood ratio statistic, LR, two modified likelihood ratio statistics M and MLB and a quadratic statistic Q. The distribution of each of the statistics CML, LR, M, MLB and Q is asymptotically chi-square with L - 1 degrees of freedom. We then conducted a monte-carlo simulation study to compare the perfor- mance of the statistics ESP, LR, M, MLB, CML and Q in terms of size and power. The statistics LR, M, MLB and Q are in general liberal and do not show power advantage over other statistics. The statistic CML, based on its asymptotic chi-square distribution, in general, holds nominal level well. It is most powerful or nearly most powerful in most situations and is simple to use. Hence, we recommend the statistic CML for use in general. For better power the statistic ESP, based on its empirical distribution, is recommended for the special situation for which there is evidence in the data that λ₁ < … < λ_L and n₁ < … < n_L, where λ₁ …, λ_L are the scale parameters and n₁,…, n_L are the sample sizes.     

Detecting a Change in the Correlation Coefficient in a Sequence of Bivariate Normal Variables

We consider the problem of testing the hypothesis that the correlation coefficient is stable in a sequence of n observations of independent, bivariate normal random variables against the alternative that the correlation coefficient changes after an unknown point t(t < n). We propose an estimate of the changepoint t and report on power comparisons between the commonly used test for this problem and our proposed test. Some applications to finance are discussed.     

Double filter instrumental variable estimation of panel data models with weakly exogenous variables

In this article, we propose instrumental variables (IV) and generalized method of moments (GMM) estimators for panel data models with weakly exogenous variables. The model is allowed to include heterogeneous time trends besides the standard fixed effects (FE). The proposed IV and GMM estimators are obtained by applying a forward filter to the model and a backward filter to the instruments in order to remove FE, thereby called the double filter IV and GMM estimators. We derive the asymptotic properties of the proposed estimators under fixed T and large N, and large T and large N asymptotics where N and T denote the dimensions of cross section and time series, respectively. It is shown that the proposed IV estimator has the same asymptotic distribution as the bias corrected FE estimator when both N and T are large. Monte Carlo simulation results reveal that the proposed estimator performs well in finite samples and outperforms the conventional IV/GMM estimators using instruments in levels in many cases.