The exact solution to the mle of a decreasing density based on arbitrarily right censored data

Nonparametric maximum likelihood estimation of decreasing and unimodal density functions, based on observations subject to arbitrary right censorship, was considered by McNichols and Padgett(1982). In order to compute their estimators, however, nonlinear equations with linear constraints had to be maximized using numerical techniques. The exact solution to this problem can now be found. An example illustrates the simplicity of the method.     

A nonparametric estimation procedure for the Hawkes process: comparison with maximum likelihood estimation

In earlier work, Kirchner [An estimation procedure for the Hawkes process. Quant Financ. 2017;17(4):571–595], we introduced a nonparametric estimation method for the Hawkes point process. In this paper, we present a simulation study that compares this specific nonparametric method to maximum-likelihood estimation. We find that the standard deviations of both estimation methods decrease as power-laws in the sample size. Moreover, the standard deviations are proportional. For example, for a specific Hawkes model, the standard deviation of the branching coefficient estimate is roughly 20% larger than for MLE – over all sample sizes considered. This factor becomes smaller when the true underlying branching coefficient becomes larger. In terms of runtime, our method clearly outperforms MLE. The present bias of our method can be well explained and controlled. As an incidental finding, we see that also MLE estimates seem to be significantly biased when the underlying Hawkes model is near criticality. This asks for a more rigorous analysis of the Hawkes likelihood and its optimization.     

The Statistical Properties of the Equity Estimator

In this article we consider the Equity estimator proposed by Krishnamurthi and Rangaswamy. We show that this estimator is inconsistent and does not necessarily improve on the mean squared error (MSE) of the least squares (LS) estimator. We perform a Monte Carlo experiment based on the price-promotion model used in marketing research, with marketing data, comparing the MSE of the Equity estimator to that of two empirical Bayes estimators and the LS estimator. We find that the empirical Bayes estimators have substantially smaller MSE than the Equity estimator in almost every case.     

Pairwise-Difference Rank Estimation of the Transformation Model

This article considers pairwise-difference rank estimators of the coefficient vector in a transformation model. These estimators, like other existing rank estimators, require no subjective bandwidth choice. Monte Carlo simulations, numerical asymptotic efficiency comparisons, and two empirical applications suggest that the proposed estimators perform well in comparison with existing semiparametric estimators.     

A note on the greenwood variance estimator in cohort life tables: Isolating the effect of one source of bias

A modification of the Greenwood variance estimator is defined and shown to be free of bias whenever its constitu­ent interval estimators are conditionally unbiased, given the sample size at the start of the interval. Using the modified estimator as a standard of comparison, the original Greenwood estimator is seen to have an intrinsic positive bias.Under­estimation of variances through the use of Greenwood's formula must be due to bias in the constituent interval estimators and/or, with fixed interval bounds, due to disregarding the random character of the total number of life table intervals to exhaustion of ttje sample. Some easy to prove properties of the modified and the original Greenwood estimators are stated that apply in the absence of censoring. A suggest­ion is made for reducing the bias of the interval variance estimators.     

Semiparametric estimation of employment duration models

There are a variety of economic areas, such as studies of employment duration and of the durability of capital goods, in which data on important variables typically are censored. The standard techinques for estimating a model from censored data require the distributions of unobservable random components of the model to be specified a priori up to a finite set of parameters, and misspecification of these distributions usually leads to inconsistent parameter estimates. However, economic theory rarely gives guidance about distributions and the standard estimation techniques do not provide convenient methods for identifying distributions from censored data. Recently, several distribution-free or semiparametric methods for estimating censored regression models have been developed. This paper presents the results of using two such methods to estimate a model of employment duration. The paper reports the operating characteristics of the semiparametric estimators and compares the semiparametric estimates with those obtained from a standard parametric model.     

A new lifetime distribution for a series-parallel system: properties,applications and estimations under progressive type-II censoring

A compound class of zero truncated Poisson and lifetime distributions is introduced. A specialization is paved to a new three-parameter distribution, called doubly Poisson-exponential distribution, which may represent the lifetime of units connected in a series-parallel system. The new distribution can be obtained by compounding two zero truncated Poisson distributions with an exponential distribution. Among its motivations is that its hazard rate function can take different shapes such as decreasing, increasing and upside-down bathtub depending on the values of its parameters. Several properties of the new distribution are discussed. Based on progressive type-II censoring, six estimation methods [maximum likelihood, moments, least squares, weighted least squares and Bayes (under linear-exponential and general entropy loss functions) estimations] are used to estimate the involved parameters. The performance of these methods is investigated through a simulation study. The Bayes estimates are obtained using Markov chain Monte Carlo algorithm. In addition, confidence intervals, symmetric credible intervals and highest posterior density credible intervals of the parameters are obtained. Finally, an application to a real data set is used to compare the new distribution with other five distributions.     

Non-parametric Estimation of the Death Rate in Branching Diffusions

We consider finite systems of diffusing particles in with branching and immigration. Branching of particles occurs at position dependent rate. Under ergodicity assumptions, we estimate the position-dependent branching rate based on the observation of the particle process over a time interval [0, t ]. Asymptotics are taken as t  → ∞. We introduce a kernel-type procedure and discuss its asymptotic properties with the help of the local time for the particle configuration. We compute the minimax rate of convergence in squared-error loss over a range of Hölder classes and show that our estimator is asymptotically optimal.     

A fortran implementation of univariate fourier series density estimation

A set of Fortran-77 subroutines is described which compute a nonparametric density estimator expressed as a Fourier series. In addition, a subroutine is given for the estimation of a cumulative distribution. Performance measures are given based on samples from a Weibull distribution. Due to small size and modest space demands, these subroutines are easily implemented on most small computers.     

Consistency of the least median of squares estimator in nonlinear regression

Recent research by Sakata and White (1995) presents the consistency and asymptotic normality of S-estimators in nonlinear regression. It is well known from research in linear regression that it is important to use a consistent high breakdown estimator as an initial estimate when computing an S-estimate. This paper presents the proof of the weak consistency of the least median of squares estimator in a nonlinear regression setting, thus suggesting that it is a reasonable choice for the starting value for computing S-estimates in nonlinear regression.     

An algorithm for computing nonparametric estimators of change‐points

An efficient FORTRAN subroutine for computing three nonparametric point estimators of change-points is offered. Estimates computed from the subroutine are obtained for two classical data sets     

Generalized confidence intervals compatible with the Min test for simultaneous comparisons of one subpopulation to several other subpopulations

A problem where one subpopulation is compared with several other subpopulations in terms of means with the goal of estimating the smallest difference between the means commonly arises in biology, medicine, and many other scientific fields. A generalization of Strass-burger-Bretz-Hochberg approach for two comparisons is presented for cases with three and more comparisons. The method allows constructing an interval estimator for the smallest mean difference, which is compatible with the Min test. An application to a fluency-disorder study is illustrated. Simulations confirmed adequate probability coverage for normally distributed outcomes for a number of designs.     

Unbiased Estimation of Central Moments by using U-statistics

We obtain an estimator of the r th central moment of a distribution, which is unbiased for all distributions for which the first r moments exist. We do this by finding the kernel which allows the r th central moment to be written as a regular statistical functional. The U-statistic associated with this kernel is the unique symmetric unbiased estimator of the r th central moment, and, for each distribution, it has minimum variance among all estimators which are unbiased for all these distributions.     

On the Estimation of Parameters of Kumaraswamy-G Distributions

Cordeiro and de Castro proposed a new family of generalized distributions based on the Kumaraswamy distribution (denoted as Kw-G). Nadarajah et al. showed that the density function of the new family of distributions can be expressed as a linear combination of the density of exponentiated family of distributions. They derived some properties of Kw-G distributions and discussed estimation of parameters using the maximum likelihood (ML) method. Cheng and Amin and Ranneby introduced a new method of estimating parameters based on Kullback–Leibler divergence (the maximum spacing (MSP) method). In this article, the estimates of parameters of Kw-G distributions are obtained using the MSP method. For some special Kw-G distributions, the new estimators are compared with ML estimators. It is shown by simulations and a real data application that MSP estimators have better properties than ML estimators.     

Recursive estimation of the mode of a multivariate density

Let f be an unknown possibly multimodal density on R^d and let X₁, X₂, … be a sequence of independent random vectors with density f. Several recursive estimates of the mode of f are proposed, and sufficient conditions ensuring their weak and strong consistency are established.     

Tests of Linear Hypotheses and Lav Estimation: A Monte Carlo Comparison

A Monte Carlo simulation is used to study the performance of the Wald, likelihood ratio and Lagrange multiplier tests for regression coefficients when least absolute value regression is used. The simulation results provide support for use of the Lagrange multiplier test, especially when certain computational advantages are considered.     

Moving truncations barrier-function methos for estimation in three-parameter lognormal models

A method of centres algorithm for maximum likelihood estimation in the three-parameter lognormal model is presented and discussed, The algorithm is a member of the class of moving truncations algorithms for solving nonlinear programming problems and is able to move the numerical search out of the region of the infinite maximum of the conditional likelihood function, thereby permitting convergence to an interior relative maximum of this function. The algorithm also includes an optimality test to locate the primary relative maximum of the likelihood function.     

Estimation of Stigmatized Characteristics of a Hidden Gang in a Finite Population

This paper considers the problem of estimating the size and mean value of a stigmatized quantitative character of a hidden gang in a finite population. The proposed method may be applied to solve domestic problems in a particular country or across countries: for example, a government may be interested in estimating the average income of victims or perpetrators of domestic violence. The proposed method is based on the technique introduced by Warner (1965) to estimate the proportion of a sensitive attribute in a finite population without threatening the privacy of the respondents. Expressions for the bias and variance of the proposed estimators are given, to a first order of approximation. Circumstances in which the method can be applied are studied and illustrated using a numerical example.