Estimating population size with truncated sampling

Let X₁, X₂,…,X_n be independent, indentically distributed random variables with density f(x,θ) with respect to a σ-finite measure μ. Let R be a measurable set in the sample space X. The value of X is observable if X ? (X?R) and not otherwise. The number J of observable X’s is binomial, N, Q, Q = 1?P(X ? R). On the basis of J observations, it is desired to estimate N and θ. Estimators considered are conditional and unconditional maximum likelihood and modified maximum likelihood using a prior weight function to modify the likelihood before maximizing. Asymptotic expansions are developed for the [Ncirc]’s of the form [Ncirc] = N + α√N + β + o_p(1), where α and β are random variables. All estimators have the same α, which has mean 0, variance σ² (a function of θ) and is asymptotically normal. Hence all are asymptotically equivalent by the usual limit distributional theory. The β’s differ and Eβ can be considered an “asymptotic bias”. Formulas are developed to compare the asymptotic biases of the various estimators. For a scale parameter family of absolutely continuous distributions with X = (0,∞) and R = (T,∞), special formuli are developed and a best estimator is found.     

A BETA-BINOMIAL MODEL FOR ESTIMATING THE SIZE OF A HETEROGENEOUS POPULATION

This paper compares the properties of various estimators for a beta‐binomial model for estimating the size of a heterogeneous population. It is found that maximum likelihood and conditional maximum likelihood estimators perform well for a large population with a large capture proportion. The jackknife and the sample coverage estimators are biased for low capture probabilities. The performance of the martingale estimator is satisfactory, but it requires full capture histories. The Gibbs sampler and Metropolis‐Hastings algorithm provide reasonable posterior estimates for informative priors.     

Parameter estimation of three-parameter Weibull distribution based on progressively Type-II censored samples

In this paper, the estimation of parameters for a three-parameter Weibull distribution based on progressively Type-II right censored sample is studied. Different estimation procedures for complete sample are generalized to the case with progressively censored data. These methods include the maximum likelihood estimators (MLEs), corrected MLEs, weighted MLEs, maximum product spacing estimators and least squares estimators. We also proposed the use of a censored estimation method with one-step bias-correction to obtain reliable initial estimates for iterative procedures. These methods are compared via a Monte Carlo simulation study in terms of their biases, root mean squared errors and their rates of obtaining reliable estimates. Recommendations are made from the simulation results and a numerical example is presented to illustrate all of the methods of inference developed here.     

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     

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.     

Estimation in a truncated bivariate poisson distribution

Moment estimators for parameters in a truncated bivariate Poisson distribution are derived in Hamdan (1972) for the special case of λ₁ = λ₂, Where λ₁, λ₂ are the marginal means. Here we derive the maximum likelihood estimators for this special case. The information matrix is also obtained which provides asymptotic covariance matrix of the maximum likelihood estimators. The asymptotic covariance matrix of moment estimators is also derived. The asymptotic efficiency of moment estimators is computed and found to be very low.     

Testing uniformity based on new entropy estimators

In this paper, we first introduce new entropy estimators for distributions with known and bounded supports. Our estimators are obtained by using constrained maximum likelihood estimation of cumulative distribution function for absolutely continuous distributions with known and bounded supports. We prove the consistency of our estimators. Then, we propose uniformity tests based on the proposed entropy estimators and compare their powers with the powers of other tests of uniformity. Our simulation results show that the proposed entropy estimators perform well in estimating entropy and testing uniformity.     

Odds ratio inference from stratified samples

One of the common problems encountered in applied statistics is that of comparing two proportions from stratified samples. One approach to this problem is via inference on the corresponding odds ratio. In this paper, the various point and interval estimators of and hypothesis testing procedures for a common odds ratio from multiple 2 ×2 tables are reviewed. Based On research to date, the conditional maximum likelihood and Mantel-Haenszel estimators are recommended as the point estimators of choice. Neither confidence intervals nor hypothesis testing metthods have been studied as well as the point estimators, but there is a confidence interval method associated with the Mantel-Haenszel estimator that is a good choice.     

Estimation of tail parameters under type i censoring

This paper is a continuation of previous work concerning the estimation of tail-parameters under Type II censoring (Weissman 1978). The same estimation problem is considered here, this truip under Type I censoring. A sample of size n is censored below aE a given level x₀it is assumed that che underlying distriibution .function (df)belogs to the domain of attraction of a known extreme-value distribution and that K - K(x_o) , the number of observed values, remains finite as on - ∞ . We offer here estimators, which are asymptotically maximum likelihood estimators (MLE's), for quantiles associated with the tail of F such as location and scale parameters, quantiles and F(x) itself (for x in the tail). The results are applied to two illustrative examples.     

Efficient estimators with categorical ranked set samples: estimation procedures for osteoporosis

Ranked set sampling (RSS) design as a cost-effective sampling is a powerful tool in situations where measuring the variable of interest is costly and time-consuming; however, ranking information about sampling units can be obtained easily through inexpensive and easy to measure characteristics at little or no cost. In this paper, we study RSS data for analysis of an ordinal population. First, we compare the problem of non-representative extreme samples under RSS and commonly-used simple random sampling. Using RSS data with tie information, we propose non-parametric and maximum likelihood estimators for population parameters. Through extensive numerical studies, we investigate the effect of various factors including ranking ability, tie generating mechanisms, the number of categories and population setting on the performance of the estimators. Finally, we apply the proposed methods to the bone disorder data to estimate the proportions of patients with osteopenia and osteoporosis status.     

Estimation of the parameters of the gamma distribution by means of a maximal invariant

The use of a scale invariance criterion allows estimation of the shape parameter of the two parameter gamma distribution without estimating the scale parameter. Simulation experiments are used to show that the resulting estimators of both parameters are better than the usual maximum likelihood estimators in terms of both bias and mean square error. Approximately unbiased versions of the maximal invariant based estimators are derived and are shown to be as good as approximately unbiased versions of the usual maximum likelihood estimators     

Statistical inference for some Volterra population processes in a random environment

We study the problem of maximum-likelihood estimation in some random-environment population models, defined through nonlinear ltǒ stochastic differential equations. It is shown that a criticality parameter can be estimated consistently. The asymptotic behavior of the estimators is analyzed, and a goodness-of-fit test is proposed.     

Computational algorithms for the factor model

Algorithms for computing the maximum likelihood estimators and the estimated covariance matrix of the estimators of the factor model are derived. The algorithms are particularly suitable for large matrices and for samples that give zero estimates of some error variances. A method of constructing estimators for reduced models is presented. The algorithms can also be used for the multivariate errors-in-variables model with known error covariance matrix.     

The combined dynamically weighted modified maximum likelihood estimators of the location and scale parameters

In this study, two new types of estimators of the location and scale parameters are proposed having high efficiency and robustness; the dynamically weighted modified maximum likelihood (DWMML) and the combined dynamically weighted modified maximum likelihood (CDWMML) estimators. Three pairs of the DWMML and two pairs of the CDWMML estimators of the location and scale parameters are produced, namely, the DWMML1, the DWMML2 and the DWMML3, and the CDWMML1 and the CDWMML2 estimators, respectively. Based on the simulation results, the DWMML1 estimators of the location and scale parameters are almost fully efficient (under normality) and robust at the same time. The DWMML3 estimators are asymptotically fully efficient and more robust than the M-estimators. The DWMML2 estimators are a compromise between efficiency and robustness. The CDWMML1 and CDWMML2 estimators are jointly very efficient and robust. Particularly, the CDWMML1 and CDWMML2 estimators of the scale parameter are superior compared to the other estimators of the scale parameter.     

Regression analysis with two mixed forms of heteroscedasticity

Regression models are here considered in which disturbances are related to both the expectation of the dependent variable and a linear conbination of certain auxiliary variables. The maximum likelihood and weighted least squares estimators are compared in estimating the form of heteroscedasticity and regression coefficients. Also a test for heteroscedasticity is discussed. Finally an example is worked out for the purpose of illustration.     

Inflated modified power series distributions with applications

In certain applications involving discrete data, it is sometimes found that X = 0 is observed with a frequency significantly higher than predicted by the assumed model. Zero inflated Poisson, binomial and negative binomial models have been employed in some clinical trials and in some regression analysis problems.

In this paper, we study the zero inflated modified power series distributions (IMPSD) which include among others the generalized Poisson and the generalized negative binomial distributions and hence the Poisson, binomial and negative binomial distributions. The structural properties along with the distribution of the sum of independent IMPSD variables are studied. The maximum likelihood estimation of the parameters of the model is examined and the variance-covariance matrix of the estimators is obtained. Finally, examples are presented for the generalized Poisson distribution to illustrate the results.     

Semiparametric Estimation of Stochastic Production Frontier Models

This article extends the linear stochastic frontier model proposed by Aigner, Lovell, and Schmidt to a semiparametric frontier model in which the functional form of the production frontier is unspecified and the distributions of the composite error terms are of known form. Pseudolikelihood estimators of the parameters characterizing the two error terms of the model are constructed based on kernel estimation of the conditional mean function. The Monte Carlo results show that the proposed estimators perform well in finite samples. An empirical application is presented. Extensions to a partially linear frontier function and to more flexible one-sided error distributions than the half-normal are discussed     

Estimation in Maxwell distribution with randomly censored data

In many practical situations, complete data are not available in lifetime studies. Many of the available observations are right censored giving survival information up to a noted time and not the exact failure times. This constitutes randomly censored data. In this paper, we consider Maxwell distribution as a survival time model. The censoring time is also assumed to follow a Maxwell distribution with a different parameter. Maximum likelihood estimators and confidence intervals for the parameters are derived with randomly censored data. Bayes estimators are also developed with inverted gamma priors and generalized entropy loss function. A Monte Carlo simulation study is performed to compare the developed estimation procedures. A real data example is given at the end of the study.     

A comparison of some robust,adaptive, and partially adaptive estimators of regression models

Numerous estimation techniques for regression models have been proposed. These procedures differ in how sample information is used in the estimation procedure. The efficiency of least squares (OLS) estimators implicity assumes normally distributed residuals and is very sensitive to departures from normality, particularly to "outliers" and thick-tailed distributions. Lead absolute deviation (LAD) estimators are less sensitive to outliers and are optimal for laplace random disturbances, but not for normal errors. This paper reports monte carlo comparisons of OLS,LAD, two robust estimators discussed by huber, three partially adaptiveestimators, newey's generalized method of moments estimator, and an adaptive maximum likelihood estimator based on a normal kernal studied by manski. This paper is the first to compare the relative performance of some adaptive robust estimators (partially adaptive and adaptive procedures) with some common nonadaptive robust estimators. The partially adaptive estimators are based on three flxible parametric distributions for the errors. These include the power exponential (Box-Tiao) and generalized t distributions, as well as a distribution for the errors, which is not necessarily symmetric. The adaptive procedures are "fully iterative" rather than one step estimators. The adaptive estimators have desirable large sample properties, but these properties do not necessarily carry over to the small sample case.

The monte carlo comparisons of the alternative estimators are based on four different specifications for the error distribution: a normal, a mixture of normals (or variance-contaminated normal), a bimodal mixture of normals, and a lognormal. Five hundred samples of 50 are used. The adaptive and partially adaptive estimators perform very well relative to the other estimation procedures considered, and preliminary results suggest that in some important cases they can perform much better than OLS with 50 to 80% reductions in standard errors.

     

Bias of maximum likelihood estimators of the response standard deviation in regression

The bias of maximum likelihood estimators of the standard deviation of the response in location/scale regression models is considered. Results are obtained for a very wide family of densities for the response variable. These are used to propose point estimators with improved mean square error properties and to demonstrate the importance of bias correction in statistical inference when samples are moderately small.