Convergence rates for empirical bayes estimators of parameters in linear regressin models

In this paper, the convergence rates of the EB estimators of the regression coefficients and the error variance in a linear model are obtained. The rates can approximate to O(n¹) arbitrarily. The convergency of the EB estimators of the regression coefiicients and the variance components in a variance component model is also investigated. The investigation makes use of the results concerning the convergence rates of the EB estimators of the parameters in multi-parameter exponential families.     

Empirical Bayes estimation in continuous one-parameter exponential families under associated samples

In this paper, we study the empirical Bayes (EB) estimation in continuous one-parameter exponential families under negatively associated (NA) samples and positively associated (PA) samples. Under certain regularity conditions, it is shown that the convergence rates of proposed EB estimators under NA or PA samples are the same as those of EB estimators under independent observations, which significantly improve the existing results in EB estimation under associated samples.     

On empirical Bayes estimation of variance components in random effects model

In this paper, Bayes estimators of variance components are derived for the one-way random effects model, and empirical Bayes (EB) estimators are constructed by the kernel estimation method of a multivariate density and its mixed partial derivatives. It is shown that the EB estimators are asymptotically optimal and convergence rates are established. Finally, an example concerning the main results is given.     

Analysis of Type-II hybrid censored competing risks data

Kundu and Gupta [Analysis of hybrid life-tests in presence of competing risks. Metrica. 2007;65:159–170] provided the analysis of Type-I hybrid censored competing risks data, when the lifetime distributions of the competing cause of failures follows exponential distribution. In this paper, we consider the analysis of Type-II hybrid censored competing risks data. It is assumed that latent lifetime distributions of the competing causes of failures follow independent exponential distributions with different scale parameters. It is observed that the maximum likelihood estimators of the unknown parameters do not always exist. We propose the modified estimators of the scale parameters, which coincide with the corresponding maximum likelihood estimators when they exist, and asymptotically they are equivalent. We obtain the exact distribution of the proposed estimators. Using the exact distributions of the proposed estimators, associated confidence intervals are obtained. The asymptotic and bootstrap confidence intervals of the unknown parameters are also provided. Further, Bayesian inference of some unknown parametric functions under a very flexible Beta-Gamma prior is considered. Bayes estimators and associated credible intervals of the unknown parameters are obtained using the Monte Carlo method. Extensive Monte Carlo simulations are performed to see the effectiveness of the proposed estimators and one real data set has been analysed for the illustrative purposes. It is observed that the proposed model and the method work quite well for this data set.     

Reducing component estimation for varying coefficient models

The estimation problem for varying coefficient models has been studied by many authors. We consider the problem in the case that the unknown functions admit different degrees of smoothness. In this paper we propose a reducing component local polynomial method to estimate the unknown functions. It is shown that all of our estimators achieve the optimal convergence rates. The asymptotic distributions of our estimators are also derived. The established asymptotic results and the simulation results show that our estimators outperform the the existing two-step estimators when the coefficient functions admit different degrees of smoothness. We also develop methods to speed up the estimation of the model and the selection of the bandwidths.     

Componentwise B-spline estimation for varying coefficient models with longitudinal data

A componentwise B-spline method is proposed for estimating the unknown functions in the varying-coefficient models with longitudinal data. Different amounts of smoothing are used for different individual coefficient functions and the estimators of different coefficient functions are obtained by different minimization operations. The local asymptotic bias and variance of the estimators are derived. It is shown that our estimators achieve the local and global optimal convergence rates even if the coefficient functions belong to different smoothness families. The asymptotic distributions of the estimators are also established and are used to construct approximate pointwise confidence intervals for coefficient functions. Finite sample properties of our procedures are studied through Monte Carlo simulations.     

Flexible empirical Bayes estimation for wavelets

Wavelet shrinkage estimation is an increasingly popular method for signal denoising and compression. Although Bayes estimators can provide excellent mean-squared error (MSE) properties, the selection of an effective prior is a difficult task. To address this problem, we propose empirical Bayes (EB) prior selection methods for various error distributions including the normal and the heavier-tailed Student t -distributions. Under such EB prior distributions, we obtain threshold shrinkage estimators based on model selection, and multiple-shrinkage estimators based on model averaging. These EB estimators are seen to be computationally competitive with standard classical thresholding methods, and to be robust to outliers in both the data and wavelet domains. Simulated and real examples are used to illustrate the flexibility and improved MSE performance of these methods in a wide variety of settings.     

BAYES EMPIRICAL BAYES APPROACH TO ESTIMATION OF THE FAILURE RATE IN EXPONENTIAL DISTRIBUTION

ABSTRACT

In the empirical Bayes (EB) decision problem consisting of squared error estimation of the failure rate in exponential distribution, a prior Λ is placed on the gamma family of prior distributions to produce Bayes EB estimators which are admissible. A subclass of such estimators is shown to be asymptotically optimal (a.o.). The results of a Monte Carlo study are presented to demonstrate the a.o. property of the Bayes EB estimators.     

Nonparametric smooth estimation of the expected inactivity time function

The expected inactivity time (EIT) function (also known as the mean past lifetime function) is a well known reliability function which has application in many disciplines such as survival analysis, actuarial studies and forensic science, to name but a few. In this paper, we use a fixed design local polynomial fitting technique to obtain estimators for the EIT function when the lifetime random variable has an unknown distribution. It will be shown that the proposed estimators are asymptotically unbiased, consistent and also, when standardized, has an asymptotic normal distribution. An optimal bandwidth, which minimizes the AMISE (asymptotic mean integrated squared error) of the estimator, is derived. Numerical examples based on simulated samples from various lifetime distributions common in reliability studies will be presented to evaluate the performances of these estimators. Finally, three real life applications will also be presented to further illustrate the wide applicability of these estimators.     

Inference Under Right Censoring in a Discrete Setup

In this article, we consider the right random censoring scheme in a discrete setup when the lifetime and censoring variables are independent and have geometric distributions with means 1/θ₁ and 1/θ₂, respectively. We first obtain the Maximum Likelihood and Method of Moment estimators of the unknown parameters. We also find the Bayes and Posterior Regret Gamma Minimax estimators of the parameters for the two cases when the prior distributions are dependent and independent, assuming a squared error loss function. We then discuss the Proportional Hazard model, and obtain Maximum Likelihood estimators of the unknown parameters and derive the Bayes estimators assuming squared error loss using Markov Chain Monte Carlo methods.     

Bivariate inverse Weibull distribution

ABSTRACT

Recently it is observed that the inverse Weibull (IW) distribution can be used quite effectively to analyse lifetime data in one dimension. The main aim of this paper is to define a bivariate inverse Weibull (BIW) distribution so that the marginals have IW distributions. It is observed that the joint probability density function and the joint cumulative distribution function can be expressed in compact forms. Several properties of this distribution such as marginals, conditional distributions and product moments have been discussed. We obtained the maximum likelihood estimates for the unknown parameters of this distribution and their approximate variance– covariance matrix. We perform some simulations to see the performances of the maximum likelihood estimators. One data set has been re-analysed and it is observed that the bivariate IW distribution provides a better fit than the bivariate exponential distribution.     

Wavelet‐based estimators for mixture regression

We consider a process that is observed as a mixture of two random distributions, where the mixing probability is an unknown function of time. The setup is built upon a wavelet‐based mixture regression. Two linear wavelet estimators are proposed. Furthermore, we consider three regularizing procedures for each of the two wavelet methods. We also discuss regularity conditions under which the consistency of the wavelet methods is attained and derive rates of convergence for the proposed estimators. A Monte Carlo simulation study is conducted to illustrate the performance of the estimators. Various scenarios for the mixing probability function are used in the simulations, in addition to a range of sample sizes and resolution levels. We apply the proposed methods to a data set consisting of array Comparative Genomic Hybridization from glioblastoma cancer studies.     

Inference on constant stress accelerated life tests for log-location-scale lifetime distributions with type-I hybrid censoring

In this paper, we consider a constant stress accelerated life test terminated by a hybrid Type-I censoring at the first stress level. The model is based on a general log-location-scale lifetime distribution with mean life being a linear function of stress and with constant scale. We obtain the maximum likelihood estimators (MLE) and the approximate maximum likelihood estimators (AMLE) of the model parameters. Approximate confidence intervals, likelihood ratio tests and two bootstrap methods are used to construct confidence intervals for the unknown parameters of the Weibull and lognormal distributions using the MLEs. Finally, a simulation study and two illustrative examples are provided to demonstrate the performance of the developed inferential methods.     

On length and area-biased Maxwell distributions

This article addresses the various properties and different methods of estimation of the unknown parameter of length and area-biased Maxwell distributions. Although, our main focus is on estimation from both frequentist and Bayesian point of view, yet, various mathematical and statistical properties of length and area-biased Maxwell distributions (such as moments, moment-generating function (mgf), hazard rate function, mean residual lifetime function, residual lifetime function, reversed residual life function, conditional moments and conditional mgf, stochastic ordering, and measures of uncertainty) are derived. We briefly describe different frequentist approaches, namely, maximum likelihood estimator, moments estimator, least-square and weighted least-square estimators, maximum product of spacings estimator and compare them using extensive numerical simulations. Next we consider Bayes estimation under different types of loss function (symmetric and asymmetric loss functions) using inverted gamma prior for the scale parameter. Furthermore, Bayes estimators and their respective posterior risks are computed and compared using Markov chain Monte Carlo (MCMC) algorithm. Also, bootstrap confidence intervals using frequentist approaches are provided to compare with Bayes credible intervals. Finally, a real dataset has been analyzed for illustrative purposes.     

On Estimating The Quantile Function For Distributions Constrained by Star-Shaped Ordering

Consider distributions F and G such that G ^-1 F is star-shaped. In the problem of estimating the quantile functions for lifetime distributions, the estimators developed by Rojo (1998) are compared with the commonly used empirical quantile function. Both the one-sample and the two-sample methods of estimation are considered for a wide class of lifetime distributions. In addition, the behavior of the estimators is examined for star-shaped ordered lifetime distributions of the important class of coherent k- out-of-n reliability systems. Results of a Monte Carlo study are presented which compare the behavior of the new estimators with that of the empirical quantile function interms of bias and mean-squared error. As the behavior of these estimators typically depends on the tail behavior of the underlying distributions, the examples presented here include distributions with short, medium and long tails. A formula for the inverse of the Kaplan-Meier estimator is provided and used to generate the simulations in the case of censored data.     

Berry-Esseen-type bounds for the kernel estimator of conditional distribution and conditional quantiles

Berry-Esseen bounds of order O(n^−1/2) have been obtained for several classes of statistics. In this paper, the rates of convergence in central limit theorem for conditional empirical functions and conditional sample quantiles based on kernel estimators are studied for both conditional and unconditional distributions.     

Large sample properties and confidence bands for component-wise varying-coefficient regression with longitudinal dependent variable

Longitudinal studies with repeatedly measured dependent variable (out-come) and time-invariant covariates are common in biomedical and epidemi-ological studies. A useful statistical tool to evaluate the effects of covariates on the outcome variable over time is the varying-coefficient regression, which considers a linear relationship between the covariates and the outcome at a specific time point but assumes the linear coefficients to be smooth curves over time. In order to provide adequate smoothing for each coefficient curve, Wu and Chiang ( 1999 ) proposed a class of component-wise kernel estimators and determined the large sample convergence rates and some of the constant terms of the mean squared errors of their estimators. In this paper we calcu¬late the explicit large sample mean squared errors, including the convergence rates and ail the constant terms, and the asymptotic distributions of the kernel estimators of Wu and Chiang ( 1999 ). These asymptotic distributions are used to construct point-wise confidence intervals and Bonferroni-type confidence bands for the coefficient curves. Through a Monte Carlo simulation, wre show that our confidence regions have adequate coverage probabilities. Applying our procedures to a NIH fetal growth study, we show that our procedures are useful to determine the effects of maternal height, cigarette smoking and al¬cohol consumption on the growth of fetal abdominal circumference over time during pregnancy.     

Analysis of left censored data from the generalized exponential distribution

The generalized exponential distribution proposed by Gupta and Kundu [Gupta, R.D and Kundu, D., 1999, Generalized exponential distributions. Australian and New Zealand Journal of Statistics, 41(2), 173–188.] is an important lifetime distribution in survival analysis. In this paper, we consider the maximum likelihood estimation procedure of the parameters of the generalized exponential distribution when the data are left censored. We obtain the maximum likelihood estimators of the unknown para-meters and the Fisher information matrix. Simulation studies are carried out to observe the performance of the estimators in small sample.     

On cusp location estimation for perturbed dynamical systems

We consider the problem of parameter estimation in the case of observation of the trajectory of the diffusion process. We suppose that the drift coefficient has a singularity of cusp type and that the unknown parameter corresponds to the position of the point of the cusp. The asymptotic properties of the maximum likelihood estimator and Bayesian estimators are described in the asymptotic of small noise, that is, as the diffusion coefficient tends to zero. The consistency, limit distributions, and the convergence of moments of these estimators are established.     

The odd log-logistic Lindley Poisson model for lifetime data

We define two new lifetime models called the odd log-logistic Lindley (OLL-L) and odd log-logistic Lindley Poisson (OLL-LP) distributions with various hazard rate shapes such as increasing, decreasing, upside-down bathtub, and bathtub. Various structural properties are derived. Certain characterizations of OLL-L distribution are presented. The maximum likelihood estimators of the unknown parameters are obtained. We propose a flexible cure rate survival model by assuming that the number of competing causes of the event of interest has a Poisson distribution and the time to event has an OLL-L distribution. The applicability of the new models is illustrated by means real datasets.