Non parametric confidence intervals for the extinction probability of a Galton–Watson branching process

It is demonstrated that the confidence intervals (CIs) for the probability of eventual extinction and other parameters of a Galton–Watson branching process based upon the maximum likelihood estimators can often have substantially lower coverage when compared to the desired nominal confidence coefficient, especially in small, more realistic sample sizes. The same conclusion holds for the traditional bootstrap CIs. We propose several adjustments to these CIs, which greatly improves coverage in most cases. We also make a correction in an asymptotic variance formula given in Stigler (1971 Stigler, S.M. (1971). The estimation of the probability of extinction and other parameters associated with branching processes. Biometrika 58(3):499–508.[Crossref], [Web of Science ®] , [Google Scholar]). The focus here is on implementation of the CIs which have good coverage, in a wide variety of cases. We also consider expected CI lengths. Some recommendations are made.     

On the genealogy and coalescence times of Bienaymé–Galton–Watson branching processes

Coalescence processes have received a lot of attention in the context of conditional branching processes with fixed population size and non-overlapping generations. Here, we focus on similar problems in the context of the standard unconditional Bienaymé–Galton–Watson branching processes, either (sub)-critical or supercritical. Using an analytical tool, we derive the structure of some counting aspects of the ancestral genealogy of such processes, including: the transition matrix of the ancestral count process and an integral representation of various coalescence times distributions, such as the time to most recent common ancestor of a random sample of arbitrary size, including full size. We illustrate our results on two important examples of branching mechanisms displaying either finite or infinite reproduction mean, their main interest being to offer a closed form expression for their probability generating functions at all times. Large time behaviors are investigated.     

Some notes on the fisher–bingham family on the sphere

The fisher–Bingham family is potentially useful class of spherical distributions, but its practical application has been hindered by various problems, including; the identification of the form of the distribution, given the parameter values; and finding an effective technique for calculating the normalising constant, a requirement for maximum likelihood estimation. It is explained how these difficulties can be resolved. and then parameter estimation and hypothesis testing are discussed. A practical example is given.     

Some properties of Lin–Wong divergence on the past lifetime data

Measures of statistical divergence are used to assess mutual similarities between distributions of multiple variables through a variety of methodologies including Shannon entropy and Csiszar divergence. Modified measures of statistical divergence are introduced throughout the present article. Those modified measures are related to the Lin–Wong (LW) divergence applied on the past lifetime data. Accordingly, the relationship between Fisher information and the LW divergence measure was explored when applied on the past lifetime data. Throughout this study, a number of relations are proposed between various assessment methods which implement the Jensen–Shannon, Jeffreys, and Hellinger divergence measures. Also, relations between the LW measure and the Kullback–Leibler (KL) measures for past lifetime data were examined. Furthermore, the present study discusses the relationship between the proposed ordering scheme and the distance interval between LW and KL measures under certain conditions.     

New results on the Ristić–Balakrishnan family of distributions

ABSTRACT

Adding new shape parameters to expand a model into a larger family of distributions to provide significantly skewed and heavy-tails plays a fundamental role in distribution theory. For any continuous baseline G distribution, Risti? and Balakrishnan (2012 Risti?, M.M., Balakrishnan, N. (2012). The gamma exponentiated exponential distribution. J. Stat. Comput. Simul. 82:1191–1206.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) proposed the gamma-generated family of distributions with an extra positive shape parameter. They presented some special models of their family but did not study its properties. This paper examines some general mathematical properties of this family which hold for any baseline model. Some distributions are studied and a number of existing results in the literature can be recovered as special cases. We estimate the model parameters by maximum likelihood and illustrate the importance of the family by means of an application to a real data set.     

On the asymptotic behaviour of the recursive Nadaraya–Watson estimator associated with the recursive sliced inverse regression method

We investigate the asymptotic behaviour of the recursive Nadaraya–Watson estimator for the estimation of the regression function in a semiparametric regression model. On the one hand, we make use of the recursive version of the sliced inverse regression method for the estimation of the unknown parameter of the model. On the other hand, we implement a recursive Nadaraya–Watson procedure for the estimation of the regression function which takes into account the previous estimation of the parameter of the semiparametric regression model. We establish the almost sure convergence as well as the asymptotic normality for our Nadaraya–Watson estimate. We also illustrate our semiparametric estimation procedure on simulated data.     

Some contributions on the multivariate Poisson–Skellam probability distribution

In this article, we introduce a new form of distribution whose components have the Poisson or Skellam marginal distributions. This new specification allows the incorporation of relevant information on the nature of the correlations between every component. In addition, we present some properties of this distribution. Unlike the multivariate Poisson distribution, it can handle variables with positive and negative correlations. It should be noted that we are only interested in modeling covariances of order 2, which means between all pairs of variables. Some simulations are presented to illustrate the estimation methods. Finally, an application of soccer teams data will highlight the relationship between number of points per season and the goal differential by some covariates.     

On the effect of item nonresponse on the estimation of a two–panel–waves wage equation

Summary: In this paper, we present results of the estimation of a two–panel–waves wage equation based on completely observed units and on a multiply imputed data set. In addition to the survey information, reliable income data is available from the register. These external data are used to assess the reliability of wage regressions that suffer from item nonresponse. The findings reveal marked differences between the complete case analyses and both versions of multiple imputation analyses. We argue that the results based on the multiply imputed data sets are more reliable than those based on the complete case analysis.* We would like to thank Statistics Finland for providing the data. We are also very grateful to Susanna Sandström and Marjo Pyy–Martikainen for their helpful advice using the Finnish data. Helpful comments from Joachim Winter and participants of the Workshop on Item Nonresponse and Data Quality in Large Social Surveys, Basel, October, 2003, on an earlier version of the paper are greatfully acknowledged. Further, we would like to thank three anonymous referees and the editor for helpful comments and suggestions.     

Learning the two parameters of the Poisson–Dirichlet distribution with a forensic application

In forensic science, the rare type match problem arises when the matching characteristic from the suspect and the crime scene is not in the reference database; hence, it is difficult to evaluate the likelihood ratio that compares the defense and prosecution hypotheses. A recent solution consists of modeling the ordered population probabilities according to the two-parameter Poisson–Dirichlet distribution, which is a well-known Bayesian nonparametric prior, and plugging the maximum likelihood estimates of the parameters into the likelihood ratio. We demonstrate that this approximation produces a systematic bias that fully Bayesian inference avoids. Motivated by this forensic application, we consider the need to learn the posterior distribution of the parameters that governs the two-parameter Poisson–Dirichlet using two sampling methods: Markov Chain Monte Carlo and approximate Bayesian computation. These methods are evaluated in terms of accuracy and efficiency. Finally, we compare the likelihood ratio that is obtained by our proposal with the existing solution using a database of Y-chromosome haplotypes.     

Projection Mean–Variance Bounds on Expectations of kth Record Values from Restricted Families

We present sharp mean–variance bounds for expectations of kth record values based on distributions coming from restricted families of distributions. These families are defined in terms of convex or star ordering with respect to generalized Pareto distribution. The bounds for expectations of kth record values from DD, DFR, DDA, and DFRA families are special cases of our results. The bounds are derived by application of the projection method.     

Comparisons of ten estimation methods for the parameters of Marshall–Olkin extended exponential distribution

The aim of this article is to compare via Monte Carlo simulations the finite sample properties of the parameter estimates of the Marshall–Olkin extended exponential distribution obtained by ten estimation methods: maximum likelihood, modified moments, L-moments, maximum product of spacings, ordinary least-squares, weighted least-squares, percentile, Crámer–von-Mises, Anderson–Darling, and Right-tail Anderson–Darling. The bias, root mean-squared error, absolute and maximum absolute difference between the true and estimated distribution functions are used as criterion of comparison. The simulation study reveals that the L-moments and maximum products of spacings methods are highly competitive with the maximum likelihood method in small as well as in large-sized samples.     

A simplified proof of a theorem on the difference of the moore–penrose inverses of two positive semi–difinte matrices

Simplified proofs are given of a standard result that establishes positive semi–definiteness of the difference of the inverses of two non–singular matrices, and of the extension of this result by Milliken and Akdeniz (1977) to the difference of the Moore–Penrose inverse of two singular matrices.     

Tukey’s M-estimator of the Poisson parameter with a special focus on small means

We treat robust M-estimators for independent and identically distributed Poisson data. We introduce modified Tukey M-estimators with bias correction and compare them to M-estimators based on the Huber function as well as to weighted likelihood and other estimators by simulation in case of clean data and data with outliers. In particular, we investigate the problem of combining robustness and high efficiencies at small Poisson means caused by the strong asymmetry of such Poisson distributions and propose a further estimator based on adaptive trimming. The advantages of the constructed estimators are illustrated by an application to smoothing count data with a time varying mean and level shifts.     

Forecasting the Penetration of a New Product—A Bayesian Approach

We adopt a Bayesian approach to forecast the penetration of a new product into a market. We incorporate prior information from an existing product and/or management judgments into the data analysis. The penetration curve is assumed to be a nondecreasing function of time and may be under shape constraints. Markov-chain Monte Carlo methods are proposed and used to compute the Bayesian forecasts. An example on forecasting the penetration of color television using the information from black-and-white television is provided. The models considered can also be used to address the general bioassay and reliability stress-testing problems.     

On the rate of convergence of the Robbins–Monro's algorithm in a linear stochastic ill-posed problem with α-mixing data

In this paper we consider a recursive method of Robbins–Monro type to estimate the solution of the linear problem Ax = u, in which the second member is measured with α-mixing errors. We also show the almost complete convergence (a.co) of this algorithm specifying its convergence rate.     

Effect of individual observations on the Box–Cox transformation

In this paper, we consider the influence of individual observations on inferences about the Box–Cox power transformation parameter from a Bayesian point of view. We compare Bayesian diagnostic measures with the ‘forward’ method of analysis due to Riani and Atkinson. In particular, we look at the effect of omitting observations on the inference by comparing particular choices of transformation using the conditional predictive ordinate and the k _d measure of Pettit and Young. We illustrate the methods using a designed experiment. We show that a group of masked outliers can be detected using these single deletion diagnostics. Also, we show that Bayesian diagnostic measures are simpler to use to investigate the effect of observations on transformations than the forward search method.     

On Transition and First Hitting Time Densities and Moments of the Ornstein–Uhlenbeck Process

This paper derives transition and first hitting time densities and moments for the Ornstein–Uhlenbeck Process (OUP) between exponential thresholds. The densities are obtained by simplifying the process via Doob’s representation into Brownian motion between affine thresholds. The densities in this paper also offer easy-to-use and fast small-time approximations for the densities of OUP between constant thresholds given that exponential thresholds are virtually constant for a small time. This is of interest for estimation with high-frequency data given that extant approaches for constant thresholds impose a large demand on computing power. The moments of the transition distribution up to order n are derived within a closed-form recursive formula that offers valuable information for management. Expressions for the moments of the first hitting time distribution are also obtained in closed form by simplifying integrals via series expansions.     

Integration of model-based recursive partitioning with bias reduction estimation: a case study assessing the impact of Oliver’s four factors on the probability of winning a basketball game

AStA Advances in Statistical Analysis - In this contribution, we investigate the importance of Oliver’s Four Factors, proposed in the literature to identify a basketball team’s...     

Estimation of the scale parameter of the Rayleigh distribution with multiply type–II censored sample

Based on a multiply type-II censored sample, the maximum likelihood estimator (MLE) and Bayes estimator for the scale parameter and the reliability function of the Rayleigh distribution are derived. However, since the MLE does not exist an explicit form, an approximate MLE which is the maximizer of an approximate likelihood function will be given. The comparisons among estimators are investigated through Monte Carlo simulations. An illustrative example with the real data concerning the 23 ball bearing in the life test is presented.