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.     

Estimation of a population size through capture-mark-recapture method: a comparison of various point and interval estimators

This article deals with the estimation of a fixed population size through capture-mark-recapture method that gives rise to hypergeometric distribution. There are a few well-known and popular point estimators available in the literature, but no good comprehensive comparison is available about their merits. Apart from the available estimators, an empirical Bayes (EB) estimator of the population size is proposed. We compare all the point estimators in terms of relative bias and relative mean squared error. Next, two new interval estimators – (a) an EB highest posterior distribution interval and (b) a frequentist interval estimator based on a parametric bootstrap method, are proposed. The comparison is then carried among the two proposed interval estimators and interval estimators derived from the currently available estimators in terms of coverage probability and average length (AL). Based on comprehensive numerical results, we rank and recommend the point estimators as well as interval estimators for practical use. Finally, a real-life data set for a green treefrog population is used as a demonstration for all the methods discussed.     

Between stability and higher-order asymptotics

We discuss the effects of model misspecifications on higher-order asymptotic approximations of the distribution of estimators and test statistics. In particular we show that small deviations from the model can wipe out the nominal improvements of the accuracy obtained at the model by second-order approximations of the distribution of classical statistics. Although there is no guarantee that the first-order robustness properties of robust estimators and tests will carry over to second-order in a neighbourhood of the model, the behaviour of robust procedures in terms of second-order accuracy is generally more stable and reliable than that of their classical counterparts. Finally, we discuss some related work on robust adjustments of the profile likelihood and outline the role of computer algebra in this type of research.     

Estimation of the shape and scale parameters of Pareto distribution using ranked set sampling

In the current paper, the estimation of the shape and location parameters α and c, respectively, of the Pareto distribution will be considered in cases when c is known and when both are unknown. Simple random sampling (SRS) and ranked set sampling (RSS) will be used, and several traditional and ad hoc estimators will be considered. In addition, the estimators of α, when c is known using an RSS version based on the order statistic that maximizes the Fisher information for a fixed set size, will be considered. These estimators will be compared in terms of their biases and mean square errors. The estimators based on RSS can be real competitors against those based on SRS.     

Uniqueness of permuted treatment estimator distributions in the proportional hazards model

We discuss findings regarding the permutation distributions of treatment effect estimators in the proportional hazards model. For fixed sample size n, we will prove that all uncensored and untied event times yield the same permutation distribution of treatment effect estimators in the proportional hazards model. In other words this distribution is irrelevant with respect to the actual event times. We will show several uniqueness properties under different conditions. These properties are useful for small sample permutation tests and also helpful to large sample cases.     

Factorial moment estimation for the bivariate generalized Waring distribution

In the context of accident theory, the bivariate generalized Waring distribution (Xekalaki, 1984) is known to offer the possibility of obtaining distinguishable estimates of the “contribution” of chance, risk exposure and proneness to an accident situation. In this paper an estimation procedure based on the first and second order factorial moments is discussed for fitting the distribution to data. Expressions for the asymptotic standard errors of the estimators of the distribution parameters as well as of the resulting estimators of the variance components that represent the roles of the above mentioned factors are given.     

A new absolute continuous bivariate generalized exponential distribution

The generalized exponential is the most commonly used distribution for analyzing lifetime data. This distribution has several desirable properties and it can be used quite effectively to analyse several skewed life time data. The main aim of this paper is to introduce absolutely continuous bivariate generalized exponential distribution using the method of Block and Basu (1974). In fact, the Block and Basu exponential distribution will be extended to the generalized exponential distribution. We call the new proposed model as the Block and Basu bivariate generalized exponential distribution, then, discuss its different properties. In this case the joint probability distribution function and the joint cumulative distribution function can be expressed in compact forms. The model has four unknown parameters and the maximum likelihood estimators cannot be obtained in explicit form. To compute the maximum likelihood estimators directly, one needs to solve a four dimensional optimization problem. The EM algorithm has been proposed to compute the maximum likelihood estimations of the unknown parameters. One data analysis is provided for illustrative purposes. Finally, we propose some generalizations of the proposed model and compare their models with each other.     

The microeconometric estimation of treatment effects—An overview

Summary The need to evaluate the performance of active labour market policies is not questioned any longer. Even though OECD countries spend significant shares of national resources on these measures, unemployment rates remain high or even increase. We focus on microeconometric evaluation which has to solve the fundamental evaluation problem and overcome the possible occurrence of selection bias. When using non-experimental data, different evaluation approaches can be thought of. The aim of this paper is to review the most relevant estimators, discuss their identifying assumptions and their (dis-)advantages. Thereby we will present estimators based on some form of exogeneity (selection on observables) as well as estimators where selection might also occur on unobservable characteristics. Since the possible occurrence of effect heterogeneity has become a major topic in evaluation research in recent years, we will also assess the ability of each estimator to deal with it. Additionally, we will also discuss some recent extensions of the static evaluation framework to allow for dynamic treatment evaluation. The authors thank Stephan L. Thomsen, Christopher Zeiss and one anonymous referee for valuable comments. The usual disclaimer applies.     

Gamma mixture of generalized error distribution

Abstract

A new symmetric heavy-tailed distribution, namely gamma mixture of generalized error distribution is defined by scaling generalized error distribution with gamma distribution, its probability density function, k-moment, skewness and kurtosis are derived. After tedious calculation, we also give the Fisher information matrix, moment estimators and maximum likelihood estimators for the parameters of gamma mixture of generalized error distribution. In order to evaluate the effectiveness of the point estimators and the stability of Fisher information matrix, extensive simulation experiments are carried out in three groups of parameters. Additionally, the new distribution is applied to Apple Inc. stock (AAPL) data and compared with normal distribution, F-S skewed standardized t distribution and generalized error distribution. It is found that the new distribution has better fitting effect on the data under the Akaike information criterion (AIC). To a certain extent, our results enrich the probability distribution theory and develop the scale mixture distribution, which will provide help and reference for financial data analysis.     

Goodness-of-fit testing for Weibull-type behavior

In the process of analyzing data, testing the fit of a model under consideration is a prerequisite for performing inference about the model parameters. In this paper we examine the goodness-of-fit testing problem for assessing whether a sample is consistent with the Weibull-type model. Inspired by the Jackson and the Lewis test statistics, originally proposed as goodness-of-fit tests for the exponential distribution, we introduce two new statistics for testing Weibull-type behavior, and study their asymptotic properties. Moreover, given that the statistics are ratios of estimators for the Weibull-tail coefficient, we obtain new estimators for the latter, and establish their consistency and asymptotic normality. The small sample behavior of our statistics and estimators is evaluated on the basis of a simulation study.     

Nonparametric Shrinkage Estimation for Aalen's Additive Hazards Model

Aalen's nonparametric additive model in which the regression coefficients are assumed to be unspecified functions of time is a flexible alternative to Cox's proportional hazards model when the proportionality assumption is in doubt. In this paper, we incorporate a general linear hypothesis into the estimation of the time‐varying regression coefficients. We combine unrestricted least squares estimators and estimators that are restricted by the linear hypothesis and produce James‐Stein‐type shrinkage estimators of the regression coefficients. We develop the asymptotic joint distribution of such restricted and unrestricted estimators and use this to study the relative performance of the proposed estimators via their integrated asymptotic distributional risks. We conduct Monte Carlo simulations to examine the relative performance of the estimators in terms of their integrated mean square errors. We also compare the performance of the proposed estimators with a recently devised LASSO estimator as well as with ridge‐type estimators both via simulations and data on the survival of primary billiary cirhosis patients.     

Bernstein conditional density estimation with application to conditional distribution and regression functions

In this paper we propose a smooth nonparametric estimation for the conditional probability density function based on a Bernstein polynomial representation. Our estimator can be written as a finite mixture of beta densities with data-driven weights. Using the Bernstein estimator of the conditional density function, we derive new estimators for the distribution function and conditional mean. We establish the asymptotic properties of the proposed estimators, by proving their asymptotic normality and by providing their asymptotic bias and variance. Simulation results suggest that the proposed estimators can outperform the Nadaraya–Watson estimator and, in some specific setups, the local linear kernel estimators. Finally, we use our estimators for modeling the income in Italy, conditional on year from 1951 to 1998, and have another look at the well known Old Faithful Geyser data.     

Shape Constrained Non‐parametric Estimators of the Baseline Distribution in Cox Proportional Hazards Model

Abstract. We investigate non‐parametric estimation of a monotone baseline hazard and a decreasing baseline density within the Cox model. Two estimators of a non‐decreasing baseline hazard function are proposed. We derive the non‐parametric maximum likelihood estimator and consider a Grenander type estimator, defined as the left‐hand slope of the greatest convex minorant of the Breslow estimator. We demonstrate that the two estimators are strongly consistent and asymptotically equivalent and derive their common limit distribution at a fixed point. Both estimators of a non‐increasing baseline hazard and their asymptotic properties are obtained in a similar manner. Furthermore, we introduce a Grenander type estimator for a non‐increasing baseline density, defined as the left‐hand slope of the least concave majorant of an estimator of the baseline cumulative distribution function, derived from the Breslow estimator. We show that this estimator is strongly consistent and derive its asymptotic distribution at a fixed point.     

Bayes Estimator of Inverse Gaussian Parameters Under General Entropy Loss Function Using Lindley's Approximation

In this article, we obtained Bayes estimators of parameters of Inverse Gaussian distributions under asymmetric loss function using Lindley's Approximation (L-Approximation). The proposed estimators have been compared with the corresponding estimators obtained under symmetric loss function and MLE for their risks. This comparison is illustrated using Monte-Carlo study of 2,000 simulated sample from the Inverse Gaussian distribution.     

Estimation of parameters of a right truncated exponential distribution

The maximum likelihood, moment and mixture of the estimators are for samples from the right truncated exponential distribution. The estimators are compared empirically when all the parameters are unknown; their bias and mean square error are investigated with the help of numerical technique. We have shown that these estimators are asymptotically unbiased. At the end, we conclude that mixture estimators are better than the maximum likelihood and moment estimators.     

Near universal consistency of the maximum pseudolikelihood estimator for discrete models

Maximum pseudolikelihood (MPL) estimators are useful alternatives to maximum likelihood (ML) estimators when likelihood functions are more difficult to manipulate than their marginal and conditional components. Furthermore, MPL estimators subsume a large number of estimation techniques including ML estimators, maximum composite marginal likelihood estimators, and maximum pairwise likelihood estimators. When considering only the estimation of discrete models (on a possibly countably infinite support), we show that a simple finiteness assumption on an entropy-based measure is sufficient for assessing the consistency of the MPL estimator. As a consequence, we demonstrate that the MPL estimator of any discrete model on a bounded support will be consistent. Our result is valid in parametric, semiparametric, and nonparametric settings.     

Estimating with percentile grouped and truncated date from a scale parameter distribution

Data which is grouped and truncated is considered. We are given numbers n₁<…<n_k=n and we observe X_ni ),i=1,…k, and the tottal number of observations available (N> n_k is unknown. If the underlying distribution has one unknown parameter θ which enters as a scale parameter, we examine the form of the equations for both conditional, unconditional and modified maximum likelihood estimators of θ and N and examine when these estimators will be finite, and unique. We also develop expressions for asymptotic bias and search for modified estimators which minimize the maximum asymptotic bias. These results are specialized tG the zxponential distribution. Methods of computing the solutions to the likelihood equatims are also discussed.     

An extension of Franklin's randomized response technique for multi-proportions

ABSTRACT

In this article, we introduce six estimators, three based on row averages and the remaining three on column averages of population proportions for trichotomous population when randomized response sampling with a normal randomizing distribution is used. The estimators have been obtained using the method of moments. All the proposed estimators are shown to be unbiased and their variances have been worked out. The percent relative efficiencies of the column total based estimators with respect to row total based estimators are investigated through empirical study.     

Bayes estimators of the reliability function and parameter of inverted exponential distribution using informative and non-informative priors

In this paper, we propose Bayes estimators of the parameter and reliability function of inverted exponential distribution under the general entropy loss function for complete, type I and type II censored samples. The proposed estimators have been compared with the corresponding maximum-likelihood estimators for their simulated risks (average loss over sample space).     

Statistical Properties of a Modified Accelerated Design for Phase I Cancer Clinical Trials

We consider the semiparametric profile likelihood inference for the distribution function under doubly censored data. For further developments of the statistical inference based on the profile likelihood ratio and alternative tools such as the score or Wald-type inference, we discuss the structures of the profile likelihood estimators and their derivatives included in the score function and the Fisher function of the profile likelihood, establishing the consistencies of their estimators.