Estimation of the Parameters of the Birnbaum–Saunders Distribution

Some alternative estimators to the maximum likelihood estimators of the two parameters of the Birnbaum–Saunders distribution are proposed. Most have high efficiencies as measured by root mean square error and are robust to departure from the model as well as to outliers. In addition, the proposed estimators are easy to compute. Both complete and right-censored data are discussed. Simulation studies are provided to compare the performance of the estimators.     

A modification of the Whittaker–Henderson method of graduation

This paper presents a modified Whittaker–Henderson (WH) Method of Graduation. After giving a closed-form solution, we show that it is of practical use because it provides not only a smoothed series identical to that of the WH graduation, but also an extrapolation beyond the sample limit of current data. In addition, we introduce two other penalized least squares problems and show that they provide the same results as those of the modified WH graduation.     

An improved method of estimation for the parameters of the Birnbaum–Saunders distribution

The Birnbaum–Saunders (BS) distribution is a positively skewed distribution, frequently used for analysing lifetime data. In this paper, we propose a simple method of estimation for the parameters of the two-parameter BS distribution by making use of some key properties of the distribution. Compared with the maximum likelihood estimators and the modified moment estimators, the proposed method has smaller bias, but having the same mean square errors as these two estimators. We also discuss some methods of construction of confidence intervals. The performance of the estimators is then assessed by means of Monte Carlo simulations. Finally, an example is used to illustrate the method of estimation developed here.     

Explicit formulas for the smoother weights of the Whittaker–Henderson graduation of order 1

The Whittaker–Henderson (WH) graduation is a widely applied smoothing method. This paper contributes to the literature by providing explicit formulas for the smoother weights of the WH graduation of order 1 along with some related results, which leads to a richer understanding of the filter.     

Default Bayesian analysis of the Behrens–Fisher problem

In the Bayesian approach, the Behrens–Fisher problem has been posed as one of estimation for the difference of two means. No Bayesian solution to the Behrens–Fisher testing problem has yet been given due, perhaps, to the fact that the conventional priors used are improper. While default Bayesian analysis can be carried out for estimation purposes, it poses difficulties for testing problems. This paper generates sensible intrinsic and fractional prior distributions for the Behrens–Fisher testing problem from the improper priors commonly used for estimation. It allows us to compute the Bayes factor to compare the null and the alternative hypotheses. This default procedure of model selection is compared with a frequentist test and the Bayesian information criterion. We find discrepancy in the sense that frequentist and Bayesian information criterion reject the null hypothesis for data, that the Bayes factor for intrinsic or fractional priors do not.     

Distribution of the mean reversion estimator in the Ornstein–Uhlenbeck process

ABSTRACT

We derive the exact distribution of the maximum likelihood estimator of the mean reversion parameter (κ) in the Ornstein–Uhlenbeck process using numerical integration through analytical evaluation of a joint characteristic function. Different scenarios are considered: known or unknown drift term, fixed or random start-up value, and zero or positive κ. Monte Carlo results demonstrate the remarkably reliable performance of our exact approach across all the scenarios. In comparison, misleading results may arise under the asymptotic distributions, including the advocated infill asymptotic distribution, which performs poorly in the tails when there is no intercept in the regression and the starting value of the process is nonzero.     

A note on the equality of the OLSE and the BLUE of the parametric function in the general Gauss–Markov model

In this note we consider the equality of the ordinary least squares estimator (OLSE) and the best linear unbiased estimator (BLUE) of the estimable parametric function in the general Gauss–Markov model. Especially we consider the structures of the covariance matrix V for which the OLSE equals the BLUE. Our results are based on the properties of a particular reparametrized version of the original Gauss–Markov model.      

A Bayesian analysis of the Conway–Maxwell–Poisson cure rate model

The purpose of this paper is to develop a Bayesian analysis for the right-censored survival data when immune or cured individuals may be present in the population from which the data is taken. In our approach the number of competing causes of the event of interest follows the Conway–Maxwell–Poisson distribution which generalizes the Poisson distribution. Markov chain Monte Carlo (MCMC) methods are used to develop a Bayesian procedure for the proposed model. Also, some discussions on the model selection and an illustration with a real data set are considered.     

Rate of convergence for asymptotic variance of the Horvitz–Thompson estimator

Drawing distinct units without replacement and with unequal probabilities from a population is a problem often considered in the literature (e.g. Hanif and Brewer, 1980, Int. Statist. Rev. 48, 317–355). In such a case, the sample mean is a biased estimator of the population mean. For this reason, we use the unbiased Horvitz–Thompson estimator (1951). In this work, we focus our interest on the variance of this estimator. The variance is cumbersome to compute because it requires the calculation of a large number of second-order inclusion probabilities. It would be helpful to use an approximation that does not need heavy calculations. The Hájek (1964) variance approximation provides this advantage as it is free of second-order inclusion probabilities. Hájek (1964) proved that this approximation is valid under restrictive conditions that are usually not fulfilled in practice. In this paper, we give more general conditions and we show that this approximation remains acceptable for most practical problems.     

On the Behrens–Fisher problem from the spatial median point of view

We present a non-parametric affine-invariant test for the multivariate Behrens–Fisher problem. The proposed method based on the spatial medians is asymptotic and does not require normality of the data. To improve its finite sample performance, we apply a correction of the type which was already used in a similar test based on trimmed means, however, our simulations show that in the case of heavy-tailed distributions our method performs better. Also in a simulation comparison with a recently published rank-based test our test yields satisfactory results.     

The Volume of Transactions and the Circulation of Money in the United States, 1950–1979

In the identity of exchange I distinguish between currency and bank payments on one side and several types of transactions and the transfer of idle money on the other. An attempt is made to measure these variables, with varying success. On the payments side I argue that currency velocity is constant (and low) and that the vast rise of bank money velocity is largely due to increased short-term investment of idle funds. The results suggest an upward shift in the level of transactions in 1968–1972, which I attribute to changes in the international role of the dollar.     

Adaptive optimal scaling of Metropolis–Hastings algorithms using the Robbins–Monro process

ABSTRACT

We present an adaptive method for the automatic scaling of random-walk Metropolis–Hastings algorithms, which quickly and robustly identifies the scaling factor that yields a specified overall sampler acceptance probability. Our method relies on the use of the Robbins–Monro search process, whose performance is determined by an unknown steplength constant. Based on theoretical considerations we give a simple estimator of this constant for Gaussian proposal distributions. The effectiveness of our method is demonstrated with both simulated and real data examples.     

Wilks’ factorization of the matrix-variate Dirichlet–Riesz distributions

The Riesz distributions on positive definite symmetric matrices are used to introduce a class of Dirichlet–Riesz distributions. In addition, several distributional properties are stated. Essentially, we show the relationship between the Dirichlet–Riesz distributions of the first kind and the second kind, respectively. We derive Wilks’ factorization of the matrix-variate Dirichlet–Riesz. Further, several results on the product of Riesz and beta–Riesz matrices with a set of Dirichlet–Riesz matrices of the first kind have been derived.     

A decomposition of the Bradley–Blackwood paired-samples omnibus test

We show that the Bradley–Blackwood simultaneous test for equal means and equal variances in paired-samples additively decomposes into separate tests of these hypotheses. The test of equal variances in the decomposition is the standard Pitman–Morgan procedure. The test of equal means in the decomposition is based on a t-ratio with (n ? 2) degrees of freedom and has the additional restriction that the variances are equal.     

Estimation of the quantile function using Bernstein–Durrmeyer polynomials

This paper studies quantile estimation using Bernstein–Durrmeyer polynomials in terms of its mean squared error and integrated mean squared error including rates of convergence as well as its asymptotic distribution. Whereas the rates of convergence are achieved for i.i.d. samples, we also show that the consistency more or less directly follows from the consistency of the sample quantiles, such that our proposal can also be applied for risk measurement in finance and insurance. Furthermore, an improved estimator based on an error-correction approach is proposed for which a general consistency result is established. A crucial issue is how to select the degree of Bernstein–Durrmeyer polynomials. We propose a novel data-adaptive approach that controls the number of modes of the corresponding density estimator. Its consistency including an uniform error bound as well as its limiting distribution in the sense of a general invariance principle are established. The finite sample properties are investigated by a Monte Carlo study. Finally, the results are illustrated by an application to photovoltaic energy research.