Estimation of parameters in logistic and log-logistic distribution with grouped data

In many situations, instead of a complete sample, data are available only in grouped form. For example, grouped failure time data occur in studies in which subjects are monitored periodically to determine whether failure has occurred in the predetermined intervals. Here the model under consideration is the log-logistic distribution. This paper demonstrates the existence and uniqueness of the MLEs of the parameters of the logistic distribution under mild conditions with grouped data. The times with the maximum failure rate and the mode of the p.d.f. of the log-logistic distribution are also estimated based on the MLEs. The methodology is further studied with simulations and exemplified with a data set with artificially introduced grouping from a locomotive life test study.     

Practical considerations when analyzing discrete survival times using the grouped relative risk model

The grouped relative risk model (GRRM) is a popular semi-parametric model for analyzing discrete survival time data. The maximum likelihood estimators (MLEs) of the regression coefficients in this model are often asymptotically efficient relative to those based on a more restrictive, parametric model. However, in settings with a small number of sampling units, the usual properties of the MLEs are not assured. In this paper, we discuss computational issues that can arise when fitting a GRRM to small samples, and describe conditions under which the MLEs can be ill-behaved. We find that, overall, estimators based on a penalized score function behave substantially better than the MLEs in this setting and, in particular, can be far more efficient. We also provide methods of assessing the fit of a GRRM to small samples.     

On maximum likelihood estimation in parametric regression with missing covariates

We consider parametric regression problems with some covariates missing at random. It is shown that the regression parameter remains identifiable under natural conditions. When the always observed covariates are discrete, we propose a semiparametric maximum likelihood method, which does not require parametric specification of the missing data mechanism or the covariate distribution. The global maximum likelihood estimator (MLE), which maximizes the likelihood over the whole parameter set, is shown to exist under simple conditions. For ease of computation, we also consider a restricted MLE which maximizes the likelihood over covariate distributions supported by the observed values. Under regularity conditions, the two MLEs are asymptotically equivalent and strongly consistent for a class of topologies on the parameter set.     

Estimating the Pareto parameters under progressive censoring data for constant-partially accelerated life tests

Accelerated life testing is widely used in product life testing experiments since it provides significant reduction in time and cost of testing. In this paper, assuming that the lifetime of items under use condition follow the two-parameter Pareto distribution of the second kind, partially accelerated life tests based on progressively Type-II censored samples are considered. The likelihood equations of the model parameters and the acceleration factor are reduced to a single nonlinear equation to be solved numerically to obtain the maximum-likelihood estimates (MLEs). Based on normal approximation to the asymptotic distribution of MLEs, the approximate confidence intervals (ACIs) for the parameters are derived. Two bootstrap CIs are also proposed. The classical Bayes estimates cannot be obtained in explicit form, so we propose to apply Markov chain Monte Carlo method to tackle this problem, which allows us to construct the credible interval of the involved parameters. Analysis of a simulated data set has also been presented for illustrative purposes. Finally, a Monte Carlo simulation study is carried out to investigate the precision of the Bayes estimates with MLEs and to compare the performance of different corresponding CIs considered.     

On the existence and uniqueness of the maximum likelihood estimates of the parameters of Birnbaum–Saunders distribution based on Type-I,Type-II and hybrid censored samples

The Birnbaum–Saunders (BS) distribution is a positively skewed distribution and is a common model for analysing lifetime data. In this paper, we discuss the existence and uniqueness of the maximum likelihood estimates (MLEs) of the parameters of BS distribution based on Type-I, Type-II and hybrid censored samples. The line of proof is based on the monotonicity property of the likelihood function. We then describe the numerical iterative procedure for determining the MLEs of the parameters, and point out briefly some recently developed simple methods of estimation in the case of Type-II censoring. Some graphical illustrations of the approach are given for three real data from the reliability literature. Finally, for illustrative purpose, we also present an example in which the MLEs do not exist.     

Exact inference for a simple step-stress model with Type-I hybrid censored data from the exponential distribution

In reliability and life-testing experiments, the researcher is often interested in the effects of extreme or varying stress factors such as temperature, voltage and load on the lifetimes of experimental units. Step-stress test, which is a special class of accelerated life-tests, allows the experimenter to increase the stress levels at fixed times during the experiment in order to obtain information on the parameters of the life distributions more quickly than under normal operating conditions. In this paper, we consider the simple step-stress model under the exponential distribution when the available data are Type-I hybrid censored. We derive the maximum likelihood estimators (MLEs) of the parameters assuming a cumulative exposure model with lifetimes being exponentially distributed. The exact distributions of the MLEs of parameters are obtained through the use of conditional moment generating functions. We also derive confidence intervals for the parameters using these exact distributions, asymptotic distributions of the MLEs and the parametric bootstrap methods, and assess their performance through a Monte Carlo simulation study. Finally, we present two examples to illustrate all the methods of inference discussed here.     

Location and Scale Parameter Family of Distributions Attaining the Bhattacharyya Bound

It is well known that, under appropriate regularity conditions, the variance of an unbiased estimator of a real-valued function of an unknown parameter can coincide with the Cramér–Rao lower bound only if the family of distributions is a one-parameter exponential family. But it seems that the necessary conditions about the probability distribution for which there exists an unbiased estimator whose variance coincides with the Bhattacharyya lower bound are not completely known. The purpose of this paper is to specify the location, scale, and location-scale parameter family of distributions attaining the general order Bhattacharyya bound in certain class.     

Estimating the parameters of a doubly truncated normal distribution

This paper deals with the maximum likelihood estimation of parameters for a doubly truncated normal distribution when the truncation points are known. We prove, in this case, that the MLEs are nonexistent (become infinite) with positive probability. For estimators that exist with probability one, the class of Bayes modal estimators or modified maximum likelihood estimators is explored. Another useful estimating procedure, called mixed estimation, is proposed. Simulations compare the behavior of the MLEs, the modified MLEs, and the mixed estimators which reveal that the MLE, in addition to being nonexistent with positive probability, behaves poorly near the upper boundary of the interval of its existence. The modified MLEs and the mixed estimators are seen to be remarkably better than the MLE     

Step partially accelerated life tests under finite mixture models

In this paper, step partially accelerated life tests are considered when the lifetime of an item under use condition follows a finite mixture of distributions. The analysis is performed when each of the components follows a general class of distributions, which includes, among others, the Weibull, compound Weibull (or three-parameter Burr type XII), power function, Gompertz and compound Gompertz distributions. Based on type-I censoring, the maximum likelihood estimates (MLEs) of the mixing proportions, scale parameters and acceleration factor are obtained. Special attention is paid to a mixture of two exponential components. Simulation results are obtained to study the precision of MLEs.     

Parametric estimation of location and scale parameters in ranked set sampling

This paper develops parametric inference for the parameters of location-scale family of distributions based on a ranked set sample. Likelihood function incorporates within-set ranking errors into the model through a missing data mechanism. The maximum likelihood estimators of the location-scale and missing data model parameters are constructed and an EM-algorithm is provided. It is shown that the proposed estimator is robust against imperfect ranking error and provides higher efficiency over its competitors.     

Regular mles for nonregular distributions

Distributions whose extremity values of the support depend on unknown pa¬rameters are usually known as nonregular distributions. In most cases, the MLEs for these parameters cannot be obtained by differentiation. Familiar examples are the uniform distribution on the interval (0,0) and the truncated exponential distribution with truncation parameter 0. However, there exist distributions whose extremity points of the support depend on unknown pa¬rameters, which nevertheless are regular in the sense that the MLEs can be obtained by differentiation. This note provides a method of constructing such nonregular distributions with regular MLEs.     

Exact inference on multiple exponential populations under a joint type-II progressive censoring scheme

ABSTRACT

Recently Mondal and Kundu [Mondal S, Kundu D. A new two sample type-II progressive censoring scheme. Commun Stat Theory Methods. 2018. doi:10.1080/03610926.2018.1472781] introduced a Type-II progressive censoring scheme for two populations. In this article, we extend the above scheme for more than two populations. The aim of this paper is to study the statistical inference under the multi-sample Type-II progressive censoring scheme, when the underlying distributions are exponential. We derive the maximum likelihood estimators (MLEs) of the unknown parameters when they exist and find out their exact distributions. The stochastic monotonicity of the MLEs has been established and this property can be used to construct exact confidence intervals of the parameters via pivoting the cumulative distribution functions of the MLEs. The distributional properties of the ordered failure times are also obtained. The Bayesian analysis of the unknown model parameters has been provided. The performances of the different methods have been examined by extensive Monte Carlo simulations. We analyse two data sets for illustrative purposes.     

Statistical analysis of middle censored competing risks data with exponential distribution

In this paper, we consider some problems of estimation and reconstruction based on middle censored competing risks data. It is assumed that the lifetime distributions of the latent failure times are independent and exponential distributed with different parameters and also that the censoring mechanism is independent. The maximum likelihood estimators (MLEs) of the unknown parameters are obtained. We then use the asymptotic distribution of the MLEs to construct approximate confidence intervals. Based on gamma priors, Lindley's approximation method is applied to obtain the Bayesian estimates of the unknown parameters under squared error loss function. Since it is not possible to construct the credible intervals, we propose and implement the Gibbs sampling technique to construct the credible intervals. Several point reconstructors for failure time of censored units are provided. Finally, a simulation study is given by Monte-Carlo simulations to evaluate the performances of the different methods and a data set is analysed to illustrate the proposed procedures.     

RANKED SET SAMPLING FROM LOCATION-SCALE FAMILIES OF SYMMETRIC DISTRIBUTIONS

Statistical inference based on ranked set sampling has primarily been motivated by nonparametric problems. However, the sampling procedure can provide an improved estimator of the population mean when the population is partially known. In this article, we consider estimation of the population mean and variance for the location-scale families of distributions. We derive and compare different unbiased estimators of these parameters based on rindependent replications of a ranked set sample of size n.Large sample properties, along with asymptotic relative efficiencies, help identify which estimators are best suited for different location-scale distributions.     

Consistent estimation of parameters and quantiles of the three-parameter gamma distribution based on Type-II right-censored data

In this paper, we propose a method of estimation of parameters and quantiles of the three-parameter gamma distribution based on Type-II right-censored data. In the proposed method, under mild conditions, the estimates always exist uniquely, and the estimators have consistency over the entire parameter space. Through Monte Carlo simulations, we further show that the proposed method performs well compared with another prominent method of estimation in terms of bias and root mean-squared error in small-sample situations. Finally, two real data sets are used for illustrating the proposed method.     

A comparative study of the K-means algorithm and the normal mixture model for clustering: Univariate case

This paper gives a comparative study of the K-means algorithm and the mixture model (MM) method for clustering normal data. The EM algorithm is used to compute the maximum likelihood estimators (MLEs) of the parameters of the MM model. These parameters include mixing proportions, which may be thought of as the prior probabilities of different clusters; the maximum posterior (Bayes) rule is used for clustering. Hence, asymptotically the MM method approaches the Bayes rule for known parameters, which is optimal in terms of minimizing the expected misclassification rate (EMCR).     

On the estimation of the extreme value and normal distribution parameters based on progressive type-II hybrid-censored data

A progressive hybrid censoring scheme is a mixture of type-I and type-II progressive censoring schemes. In this paper, we mainly consider the analysis of progressive type-II hybrid-censored data when the lifetime distribution of the individual item is the normal and extreme value distributions. Since the maximum likelihood estimators (MLEs) of these parameters cannot be obtained in the closed form, we propose to use the expectation and maximization (EM) algorithm to compute the MLEs. Also, the Newton–Raphson method is used to estimate the model parameters. The asymptotic variance–covariance matrix of the MLEs under EM framework is obtained by Fisher information matrix using the missing information and asymptotic confidence intervals for the parameters are then constructed. This study will end up with comparing the two methods of estimation and the asymptotic confidence intervals of coverage probabilities corresponding to the missing information principle and the observed information matrix through a simulation study, illustrated examples and real data analysis.     

Fisher’s observed information matrix for location and scale parameters under type i censoring

An expression for Fisher's observed information matrix is given under type I censoring for any location-scale distribution under mild requirements. It is illustrated on a data set which has been analyzed by several authors.     

On estimating parameters of a progressively censored lognormal distribution

We consider the problem of making statistical inference on unknown parameters of a lognormal distribution under the assumption that samples are progressively censored. The maximum likelihood estimates (MLEs) are obtained by using the expectation-maximization algorithm. The observed and expected Fisher information matrices are provided as well. Approximate MLEs of unknown parameters are also obtained. Bayes and generalized estimates are derived under squared error loss function. We compute these estimates using Lindley's method as well as importance sampling method. Highest posterior density interval and asymptotic interval estimates are constructed for unknown parameters. A simulation study is conducted to compare proposed estimates. Further, a data set is analysed for illustrative purposes. Finally, optimal progressive censoring plans are discussed under different optimality criteria and results are presented.     

Inference based on progressive Type I interval censored data from log-normal distribution

This article considers inference for the log-normal distribution based on progressive Type I interval censored data by both frequentist and Bayesian methods. First, the maximum likelihood estimates (MLEs) of the unknown model parameters are computed by expectation-maximization (EM) algorithm. The asymptotic standard errors (ASEs) of the MLEs are obtained by applying the missing information principle. Next, the Bayes’ estimates of the model parameters are obtained by Gibbs sampling method under both symmetric and asymmetric loss functions. The Gibbs sampling scheme is facilitated by adopting a similar data augmentation scheme as in EM algorithm. The performance of the MLEs and various Bayesian point estimates is judged via a simulation study. A real dataset is analyzed for the purpose of illustration.