Bias reduction for the maximum likelihood estimator of the doubly-truncated Poisson distribution

ABSTRACT

We derive an analytic expression for the bias of the maximum likelihood estimator of the parameter in a doubly-truncated Poisson distribution, which proves highly effective as a means of bias correction. For smaller sample sizes, our method outperforms the alternative of bias correction via the parametric bootstrap. Bias is of little concern in the positive Poisson distribution, the most common form of truncation in the applied literature. Bias appears to be the most severe in the doubly-truncated Poisson distribution, when the mean of the distribution is close to the right (upper) truncation.     

Estimation of the covariance matrix with two-step monotone missing data

Abstract

We suggest shrinkage based technique for estimating covariance matrix in the high-dimensional normal model with missing data. Our approach is based on the monotone missing scheme assumption, meaning that missing values patterns occur completely at random. Our asymptotic framework allows the dimensionality p grow to infinity together with the sample size, N, and extends the methodology of Ledoit and Wolf (2004) Ledoit, O., Wolf, M. (2004). A well-conditioned estimator for large dimensional covariance matrices. J. Multivariate Anal. 88:365–411.[Crossref], [Web of Science ®] , [Google Scholar] to the case of two-step monotone missing data. Two new shrinkage-type estimators are derived and their dominance properties over the Ledoit and Wolf (2004) Ledoit, O., Wolf, M. (2004). A well-conditioned estimator for large dimensional covariance matrices. J. Multivariate Anal. 88:365–411.[Crossref], [Web of Science ®] , [Google Scholar] estimator are shown under the expected quadratic loss. We perform a simulation study and conclude that the proposed estimators are successful for a range of missing data scenarios.     

Gamma shared frailty model based on reversed hazard rate

Abstract

Frailty models are used in survival analysis to account for unobserved heterogeneity in individual risks to disease and death. To analyze bivariate data on related survival times (e.g., matched pairs experiments, twin, or family data), shared frailty models were suggested. Shared frailty models are frequently used to model heterogeneity in survival analysis. The most common shared frailty model is a model in which hazard function is a product of random factor(frailty) and baseline hazard function which is common to all individuals. There are certain assumptions about the baseline distribution and distribution of frailty. In this paper, we introduce shared gamma frailty models with reversed hazard rate. We introduce Bayesian estimation procedure using Markov Chain Monte Carlo (MCMC) technique to estimate the parameters involved in the model. We present a simulation study to compare the true values of the parameters with the estimated values. Also, we apply the proposed model to the Australian twin data set.     

ML and Bayes estimation in a two-phase tandem queue with a second optional service and random feedback

ABSTRACT

In this article, we consider a two-phase tandem queueing model with a second optional service and random feedback. The first phase of service is essential for all customers and after the completion of the first phase of service, any customer receives the second phase of service with probability α, feedback to the tail of the first queue with probability β if the service is not successful and leaves the system with probability 1 ? α ? β. In this model, our main purpose is to estimate the parameters of the model, traffic intensity, and mean system size, in the steady state, via maximum likelihood and Bayesian methods. Furthermore, we find asymptotic confidence intervals for mean system size. Finally, by a simulation study, we compute the confidence levels and mean length for asymptotic confidence intervals of mean system size with a nominal level 0.95.     

Semiparametric inference on partially linear single-index model

ABSTRACT

We consider semiparametric inference on the partially linearsingle-index model (PLSIM). The generalized likelihood ratio (GLR) test is proposed to examine whether or not a family of new semiparametric models fits adequately our given data in the PLSIM. A new GLR statistic is established to deal with the testing of the index parameter α₀ in the PLSIM. The newly proposed statistic is shown to asymptotically follow a χ²-distribution with the scale constant and the degrees of freedom being independent of the nuisance parameters or function. Some finite sample simulations and a real example are used to illustrate our proposed methodology.     

A family of estimators of population variance in two-occasion rotation patterns

Abstract

In this article, we have considered the problem of estimation of population variance on current (second) occasion in two occasion successive (rotation) sampling. A class of estimators of population variance has been proposed and its asymptotic properties have been discussed. The proposed class of estimators is compared with the sample variance estimator when there is no matching from the previous occasion and the Singh et al. (2013) estimator. Optimum replacement policy is discussed. It has been shown that the suggested estimator is more efficient than the Singh et al. (2013) estimator and a usual unbiased estimator when there is no matching. An empirical study is carried out in support of the present study.     

Inference in log-alpha-power and log-skew-normal multivariate models

ABSTRACT

Random vectors with positive components are common in many applied fields, for example, in meteorology, when daily precipitation is measured through a region Marchenko and Genton (2010 Marchenko, Y., Genton, M. (2010). Multivariate log-skew-elliptical distributions with applications to precipitation data. Environmetrics 21:318–340.[Crossref], [Web of Science ®] , [Google Scholar]). Frequently, the log-normal multivariate distribution is used for modeling this type of data. This modeling approach is not appropriate for data with high asymmetry or kurtosis. Consequently, more flexible multivariate distributions than the log-normal multivariate are required. As an alternative to this distribution, we propose the log-alpha-power multivariate and log-skew-normal multivariate models. The first model is an extension for positive data of the fractional order statistics model Durrans (1992 Durrans, S. (1992). Distributions of fractional order statistics in hydrology. Water Resour. Res. 28:1649–1655.[Crossref], [Web of Science ®] , [Google Scholar]). The second one is an extension of the log-skew-normal model studied by Mateu-Figueras and Pawlowsky-Glahn (2007 Mateu-Figueras, G., Pawlowsky-Glahn, V. (2007). The skew-normal distribution on the simplex. Commun. Stat.-Theory Methods 36:1787–1802.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). We study parameter estimation for these models by means of pseudo-likelihood and maximum likelihood methods. We illustrate the proposal analyzing a real dataset.     

Inferential analysis for the reliability parameter based on the three-parameter Lindley distribution

ABSTRACT

In this article, we consider the estimation of R = P(Y < X), when Y and X are two independent three-parameter Lindley (LI) random variables. On the basis of two independent samples, the modified maximum likelihood estimator along its asymptotic behavior and conditional likelihood-based estimator are used to estimate R. We also propose sample-based estimate of R and the associated credible interval based on importance sampling procedure. A real life data set involving the times to breakdown of an insulating fluid is presented and analyzed for illustrative purposes.     

Regression Analysis of Left-truncated and Case I Interval-censored Data with the Additive Hazards Model

In recent years the analysis of interval-censored failure time data has attracted a great deal of attention and such data arise in many fields including demographical studies, economic and financial studies, epidemiological studies, social sciences, and tumorigenicity experiments. This is especially the case in medical studies such as clinical trials. In this article, we discuss regression analysis of one type of such data, Case I interval-censored data, in the presence of left-truncation. For the problem, the additive hazards model is employed and the maximum likelihood method is applied for estimations of unknown parameters. In particular, we adopt the sieve estimation approach that approximates the baseline cumulative hazard function by linear functions. The resulting estimates of regression parameters are shown to be consistent and efficient and have an asymptotic normal distribution. An illustrative example is provided.