Bayesian estimation based on ranked set sample from Morgenstern type bivariate exponential distribution when ranking is imperfect

In this paper we consider Bayes estimation based on ranked set sample when ranking is imperfect, in which units are ranked based on measurements made on an easily and exactly measurable auxiliary variable X which is correlated with the study variable Y. Bayes estimators under squared error loss function and LINEX loss function for the mean of the study variate Y, when (X, Y) follows a Morgenstern type bivariate exponential distribution, are obtained based on both usual ranked set sample and extreme ranked set sample. Estimation procedures developed in this paper are illustrated using simulation studies and a real data.     

The sample estimation

Statistical research, as a rule, is based on the sample. Therefore, it is important to evaluate the quality of the sample studied. Our sample evaluation is based on a new interpretation of such concepts as a random factor, disturbance factor, homogeneity, representativeness, and population.     

Berry-Esséen rate in asymptotic normality for perturbed sample quantiles

In estimating quantiles with a sample of sizeN obtained from a distributionF, the perturbed sample quantiles based on a kernel functionk have been investigated by many authors. It is well known that their behaviour depends on the choices of “window-width”, sayw _N. Under suitable and reasonably mild assumptions onF andk, Ralescu and Sun (1993) have recently proven that lim _N→∞ N ^1/4w_N=0 is the necessary and sufficient condition for the asymptotic normality of the perturbed sample quantiles. In this paper, their rate of convergence is investigated. It turns out that the optimal Berry-Esséen rate ofO(N^?1/2) can be achieved by choosing the window-width suitably, sayw _N=O(N^?1/2). The obtained results, in addition to being explicit enough to verify the sufficient condition for the asymptotic normality, improve Ralescu's (1992) result of which the rate is of order (logN)N ^?1/2.     

Bayesian inference in a sample selection model

This paper develops methods of Bayesian inference in a sample selection model. The main feature of this model is that the outcome variable is only partially observed. We first present a Gibbs sampling algorithm for a model in which the selection and outcome errors are normally distributed. The algorithm is then extended to analyze models that are characterized by nonnormality. Specifically, we use a Dirichlet process prior and model the distribution of the unobservables as a mixture of normal distributions with a random number of components. The posterior distribution in this model can simultaneously detect the presence of selection effects and departures from normality. Our methods are illustrated using some simulated data and an abstract from the RAND health insurance experiment.     

Bayesian prediction with cointegrated vector autoregressions

A complete procedure for calculating the joint predictive distribution of future observations based on the cointegrated vector autoregression is presented. The large degree of uncertainty in the choice of cointegration vectors is incorporated into the analysis via the prior distribution. This prior has the effect of weighing the predictive distributions based on the models with different cointegration vectors into an overall predictive distribution. The ideas of Litterman [Mimeo, Massachusetts Institute of Technology, 1980] are adopted for the prior on the short run dynamics of the process resulting in a prior which only depends on a few hyperparameters. A straightforward numerical evaluation of the predictive distribution based on Gibbs sampling is proposed. The prediction procedure is applied to a seven-variable system with a focus on forecasting Swedish inflation.     

Bayesian averaging,prediction and nonnested model selection

This paper studies the asymptotic relationship between Bayesian model averaging and post-selection frequentist predictors in both nested and nonnested models. We derive conditions under which their difference is of a smaller order of magnitude than the inverse of the square root of the sample size in large samples. This result depends crucially on the relation between posterior odds and frequentist model selection criteria. Weak conditions are given under which consistent model selection is feasible, regardless of whether models are nested or nonnested and regardless of whether models are correctly specified or not, in the sense that they select the best model with the least number of parameters with probability converging to 1. Under these conditions, Bayesian posterior odds and BICs are consistent for selecting among nested models, but are not consistent for selecting among nonnested models and possibly overlapping models. These findings have important bearing for applied researchers who are frequent users of model selection tools for empirical investigation of model predictions.     

Improved estimation of extreme quantiles in the multivariate Lomax (Pareto II) distribution

Estimation of a quantile of the common marginal distribution in a multivariate Lomax (Pareto II) distribution with unknown location and scale parameters is considered. For quadratic loss and specified extreme quantiles, it is established that the best affine equivariant procedure is inadmissible by constructing a better estimator.     

Bayesian prediction in doubly stochastic Poisson process

A stochastic marked point process model based on doubly stochastic Poisson process is considered in the problem of prediction for the total size of future marks in a given period, given the history of the process. The underlying marked point process \((T_{i},Y_{i})_{i\ge 1}\) , where \(T_{i}\) is the time of occurrence of the \(i\) th event and the mark \(Y_{i}\) is its characteristic (size), is supposed to be a non-homogeneous Poisson process on \(\mathbb {R}_{+}^{2}\) with intensity measure \(P\times \varTheta \) , where \(P\) is known, whereas \(\varTheta \) is treated as an unknown measure of the total size of future marks in a given period. In the problem of prediction considered, a Bayesian approach is used assuming that \(\varTheta \) is random with prior distribution presented by a gamma process. The best predictor with respect to this prior distribution is constructed under a precautionary loss function. A simulation study for comparing the behavior of the predictors under various criteria is provided.     

Bayesian point estimation of the cointegration space

A neglected aspect of the otherwise fairly well developed Bayesian analysis of cointegration is point estimation of the cointegration space. It is pointed out here that, due to the well known non-identification of the cointegration vectors, the parameter space is not Euclidean and the loss functions underlying the conventional Bayes estimators are therefore questionable. We present a Bayes estimator of the cointegration space which takes the curved geometry of the parameter space into account. This estimate has the interpretation of being the posterior mean cointegration space and is invariant to the order of the time series, a property not shared with many of the Bayes estimators in the cointegration literature. An overall measure of cointegration space uncertainty is also proposed. Australian interest rate data are used for illustration. A small simulation study shows that the new Bayes estimator compares favorably to the maximum likelihood estimator.     

Bayesian estimation of a random coefficient model

This paper carries out a Bayesian analysis of the Hildreth-Houck (1968) random coefficient model and applies it to some cross-section production function data. Posterior distributions for mean coefficients, actual coefficients, variances and variance ratios are derived. The variance ratio posteriors are largely uninformative but they do lead to relatively informative densities on the variances, and the problem of negative variance estimates, obtained with previous techniques, is overcome. Posterior densities for the mean coefficients are not extremely sensitive to the variance ratios.     

An ergodic theorem for proportions of observations that fall into random sets determined by sample quantiles

Assume that a sequence of observations \((X_n; n\ge 1)\) forms a strictly stationary process with an arbitrary univariate cumulative distribution function. We investigate almost sure asymptotic behavior of proportions of observations in the sample that fall into a random region determined by a given Borel set and a sample quantile. We provide sufficient conditions under which these proportions converge almost surly and describe the law of the limiting random variable.     

Bayesian U-type design for nonparametric response surface prediction

This paper deals with Bayesian design over U-type designs of n runs and s factors with q levels for nonparametric response surface prediction. The criterion is developed in terms of the asymptotic approach of Mitchell et al. (Ann Statist 22: 634–651, 1994) for a specific covariance kernel. An optimal design is given in approximate design theory over the all level combinations. A connection with orthogonality and aberration is established. A lower bound for the criterion is provided, and numerical results show that this lower bound is tight.     

Bayesian inference and prediction in the single server Markovian queue

The use of stochastic models and performance measures for the analysis of real life queuing scenarios are based on the fundamental premise that parameters values are known. This is a rarity since more often than not, parameters are usually unknown and require to be estimated. This paper presents techniques for the same from Bayesian perspective. The queue we intend to deal with is the M/M/1 queuing model. Several closed form expressions on posterior inference and prediction are presented which can be readily implemented using standard spreadsheet tools. Previous work in this direction resulted in non-existence of posterior moments. A way out is suggested. Interval estimates and tests of hypothesis on performance measures are also presented.     

A Bayesian estimation of macro and micro CES production functions

The paper presents a Bayesian estimation of CES production functions. The proposed estimation method is easier to use than methods so far developed and it allows a direct comparison with the maximum likelihood estimator. A Bayesian highest posterior density interval inference is made to examine the validity of the Cobb–Douglas representation. The method is applied to Japanese macro data and to micro data on two Japanese manufacturing plants. The results indicate that the current practice of employing the Cobb–Douglas form a priori ought to be corrected. In addition, the micro data present results paradoxical to the commonly held view of capital intensity.     

Small sample considerations in estimation from panel data

Recent interest in statistical inference for panel data has focused on the problem of unobservable, individual-specific, random effects and the inconsistencies they introduce in estimation when they are correlated with other exogenous variables. Analysis of this problem has always assumed the variance components to be known. In this paper, we re-examine some of these questions in finite samples when the variance components must be estimated. In particular, when the effects are uncorrelated with other explanatory variables, we show that (i) the feasible Gauss-Markov estimator is more efficient than the within groups estimator for all but the fewest degrees of freedom and its variance is never more than 17% above the Cramer-Rao bound, (ii) the asymptotic approximation to the variance of the feasible Gauss-Markov estimator is similarly within 17% of the true variance but remains significantly smaller for moderately large samples sizes, and (iii) more efficient estimators for the variance components do not necessarily yield more efficient feasible Gauss-Markov estimators.     

Large sample estimation and testing procedures for dynamic equation systems

In 1980, Palm and Zellner presented a number of joint or system estimation and testing procedures for the final equations, transfer functions and structural equations of a multivariate dynamic model with normally distributed vector moving average errors. In this paper, we show that there is a flaw in their argument which renders most of the procedures invalid. Possible alternative procedures are described.