Control variate method for stationary processes

The sample mean is one of the most natural estimators of the population mean based on independent identically distributed sample. However, if some control variate is available, it is known that the control variate method reduces the variance of the sample mean. The control variate method often assumes that the variable of interest and the control variable are i.i.d. Here we assume that these variables are stationary processes with spectral density matrices, i.e. dependent. Then we propose an estimator of the mean of the stationary process of interest by using control variate method based on nonparametric spectral estimator. It is shown that this estimator improves the sample mean in the sense of mean square error. Also this analysis is extended to the case when the mean dynamics is of the form of regression. Then we propose a control variate estimator for the regression coefficients which improves the least squares estimator (LSE). Numerical studies will be given to see how our estimator improves the LSE.     

Predictive estimation of finite population mean using product estimator

If the predictive approach advocated byBasu [1971] is adopted for estimating the mean of a finite population, it is observed that the use of mean per unit estimator, regression estimator and ratio estimator as a predictor for the mean of unobserved units in the population result in the corresponding customary estimators of the mean of the whole population. Whereas if the product estimator is used as a predictor for the mean of unobserved units in the population, the resulting estimator of the mean of the whole population is different from the customary product estimator. The new estimator so obtained is compared with the customary product estimator.     

Conditional implicit mean and the law of iterated integrals

This paper presents a new framework which generalizes the concept of conditional expectation to mean values which are implicitly defined as unique solutions to some functional equation. We call such a mean value an implicit mean. The implicit mean and its very special example, the quasi-linear mean, have been extensively applied to economics and decision theory. This paper provides a procedure of defining the conditional implicit mean and then analyzes its properties. In particular, we show that the conditional implicit mean is in general “biased” in the sense that an analogue of the law of iterated expectations does not hold and we characterize the quasi-linear mean as the only implicit mean which is “unbiased”.     

Goodness-of-fit tests for long memory moving average marginal density

This paper addresses the problem of fitting a known density to the marginal error density of a stationary long memory moving average process when its mean is known and unknown. In the case of unknown mean, when mean is estimated by the sample mean, the first order difference between the residual empirical and null distribution functions is known to be asymptotically degenerate at zero, and hence can not be used to fit a distribution up to an unknown mean. In this paper we show that by using a suitable class of estimators of the mean, this first order degeneracy does not occur. We also investigate the large sample behavior of tests based on an integrated square difference between kernel type error density estimators and the expected value of the error density estimator based on errors. The asymptotic null distributions of suitably standardized test statistics are shown to be chi-square with one degree of freedom in both cases of the known and unknown mean. In addition, we discuss the consistency and asymptotic power against local alternatives of the density estimator based test in the case of known mean. A finite sample simulation study of the test based on residual empirical process is also included.     

Estimation of bias and standard error of an improved estimator of normal mean

This paper considers an improved estimator of normal mean which is obtained by considering a feasible version of minimum mean squared error estimator. The exact expression for the bias and the mean squared error are fairly complicated and do not provide any guidelines as how to estimate the standard error of improved estimator. As is well known that any estimator without a formula for standard error has little practical utility. We therefore derive unbiased estimators for the bias and mean squared error of the improved estimator. Incidently, they turn out to be minimum variance unbiased estimators. Further, this exercise yields a simple formula for estimating the standard error. Based on the criterion of estimated standard error, the efficiency of the improved estimator with respect to the traditional unbiased estimator (i.e., sample mean) is examined numerically. The relationship with asymptotic standard error is also studied.     

Estimation of a proportion in survey sampling using the logratio approach

The estimation of a mean of a proportion is a frequent task in statistical survey analysis, and often such ratios are estimated from compositions such as income components, wage components, tax components, etc. In practice, the weighted arithmetic mean is regularly used to estimate the center of the data. However, this estimator is not appropriate if the ratios are estimated from compositions, because the sample space of compositional data is the simplex and not the usual Euclidean space. We demonstrate that the weighted geometric mean is useful for this purpose. Even for different sampling designs, the weighted geometric mean shows excellent behavior.     

Time-consistent investment policies in Markovian markets: A case of mean–variance analysis

The optimal investment policy for a standard multi-period mean–variance model is not time-consistent because the variance operator is not separable in the sense of the dynamic programming principle. With a nested conditional expectation mapping, we develop an investment model with time consistency in Markovian markets. Furthermore, we examine the differences of the investment policies with a riskless asset from those without a riskless asset. Analytical solutions for time-consistent optimal investment policies and the resulting mean–variance efficient frontier are obtained. Finally, using numerical examples, we show that the optimal investment policy derived from our model is more efficient than that of the standard mean–variance model in which the trade-off is determined between the mean and variance of the terminal wealth.     

A note on the estimation of mean in normal population

An estimator for the mean of a Normal population is presented and its properties are analyzed. The efficiency with respect to the conventional sample mean is also examined.     

Quantile-locating functions and the distance between the mean and quantiles

Given a random variable X with finite mean, for each 0 < p < 1, a new sharp bound is found on the distance between a p -quantile of X and its mean in terms of the central absolute first moment of X . The new bounds strengthen the fact that the mean of X is within one standard deviation of any of its medians, as well as a recent quantile-generalization of this fact by O'Cinneide.     

Cusums for tracking arbitrary functionals

Cusum charts are widely used for detecting deviations of a process about a target value and also for finding evidence of change in the mean of a process. The testing theory approximates the process by a Wiener process or a Brownian bridge, respectively. For quality control, it is important that other aspects are monitored in addition to or instead of the mean. Here, we show that cusum theory is easily adapted when the target is not the mean but some other aspect of the distribution.     

Sequential estimation of normal mean under asymmetric loss function with a shrinkage stopping rule

The problem of estimating a normal mean with unknown variance is considered under an asymmetric loss function such that the associated risk is bounded from above by a known quantity. In the absence of a fixed sample size rule, a sequential stopping rule and two sequential estimators of the mean are proposed and second-order asymptotic expansions of their risk functions are derived. It is demonstrated that the sample mean becomes asymptotically inadmissible, being dominated by a shrinkage-type estimator. Also a shrinkage factor is incorporated in the stopping rule and similar inadmissibility results are established. Received September 1997     

Sequential point estimation of normal mean under LINEX loss function

A sequential point estimation of the mean of a normal distribution is considered under LINEX loss function. The regret of sequential procedures are obtained. Furthermore, it is shown that a sequential procedure with the sample mean as an estimate is asymptotically inadmissible. An accerelated stopping time is also considered. Received: December 1999     

Valuation under imperfect information: Bayesian learning from the performance of the firm and the market

This paper develops a valuation model for a project or firm in the presence of uncertainty about the mean of the probability distribution of the cash flows generated by the project. Its major point is that in the presence of parameter uncertainty the value of the project is smaller than in the case where the mean cash flows is perfectly known. The second point is that when there is a known covariance between project cash flows and aggregate market cash flows investors can learn about the unknown mean cash flows by observing the market. This is referred to as ‘learning from the market’.     

Bias in the estimation of the mean reversion parameter in continuous time models

It is well known that for continuous time models with a linear drift standard estimation methods yield biased estimators for the mean reversion parameter both in finite discrete samples and in large in-fill samples. In this paper, we obtain two expressions to approximate the bias of the least squares/maximum likelihood estimator of the mean reversion parameter in the Ornstein–Uhlenbeck process with a known long run mean when discretely sampled data are available. The first expression mimics the bias formula of Marriott and Pope (1954) for the discrete time model. Simulations show that this expression does not work satisfactorily when the speed of mean reversion is slow. Slow mean reversion corresponds to the near unit root situation and is empirically realistic for financial time series. An improvement is made in the second expression where a nonlinear correction term is included into the bias formula. It is shown that the nonlinear term is important in the near unit root situation. Simulations indicate that the second expression captures the magnitude, the curvature and the non-monotonicity of the actual bias better than the first expression.     

Estimation in Two-Stage Models with Heteroscedasticity

A surprising number of important problems can be cast in the framework of estimating a mean and variance using data arising from a two-stage structure. The first stage is a random sampling of "units" with some quantity of interest associated with the unit. The second stage produces an estimate of that quantity and usually, but not always, an estimated standard error, which may change considerably across units. Heteroscedasticity in the estimates over different units can arise for a number of reasons, including variation associated with the unit and changing sampling effort over units. This paper presents a broad discussion of the problem of making inferences for the population mean and variance associated with the unobserved true values at the first stage of sampling. A careful discussion of the causes of heteroscedasticity is given, followed by an examination of ways in which inferences can be carried out in a manner that is robust to the nature of the within unit heteroscedasticity. Among the conclusions are that under any type of heteroscedasticity, an unbiased estimate of the mean and the variance of the estimated mean can be obtained by using the estimates as if they were true unobserved values from the first stage. The issue of using the mean versus a weighted average which tries to account for the heteroscedasticity is also discussed. An unbiased estimate of the population variance is given and the variance of this estimate and its covariance with the estimated mean is provided under various types of heteroscedasticity. The two-stage setting arises in many contexts including the one-way random effects models with replication, meta-analysis, multi-stage sampling from finite populations and random coefficients models. We will motivate and illustrate the problem with data arising from these various contexts with the goal of providing a unified framework for addressing such problems.     

On parameter estimation of partly observed bilinear discrete-time stochastic systems

The parameter estimation problem of a partly observed nonlinear discrete-time stochastic system is considered. The unobserved component of the system is a q-dimensional stable autoregressive process of the pth order with random parameters, observed in the presence of multiplicative and additive noises. The distributions of all the noises of the system are supposed to be unknown. The problem is to estimate the mean of the drifting parameters of the object and variances of the additive noises of the system. Asymptotic correlation estimators of all these parameters are investigated and sequential estimators with given mean square accuracy of the mean of the drifting autoregressive parameters are obtained.     

Mean estimation bias in least squares estimation of autoregressive processes

Estimation of the parameters of an autoregressive process with a mean that is a function of time is considered. Approximate expressions for the bias of the least squares estimator of the autoregressive parameters that is due to estimating the unknown mean function are derived. For the case of a mean function that is a polynomial in time, a reparameterization that isolates the bias is given. Using the approximate expressions, a method of modifying the least squares estimator is proposed. A Monte Carlo study of the second-order autoregressive process is presented. The Monte Carlo results agree well with the approximate theory and, generally speaking, the modified least squares estimators performed better than the least squares estimator. For the second-order process we also considered the empirical properties of the estimated generalized least squares estimator of the mean function and the error made in predicting the process one, two and three periods in the future.     

Mean squared error of the empirical best linear unbiased predictor in an orthogonal finite discrete spectrum linear regression model

The mean squared error (MSE) of the empirical best linear unbiased predictor in an orthogonal finite discrete spectrum linear regression model is derived and a comparison with the MSE of the best linear unbiased predictor in this model is made. It is shown that under weak conditions these two mean square errors are asymptotically the same.     

A product Pareto distribution

The Pareto distributions are becoming increasing prominent in several applied areas. In this note, a new Pareto distribution is introduced. It takes the form of the product of two Pareto probability density functions. Various structural properties of this distribution are derived, including its cumulative distribution function, moments, mean deviation about the mean, mean deviation about the median, entropy, asymptotic distribution of the extreme order statistics, method of moments estimates, maximum likelihood estimates and the Fisher information matrix. The calculations involve the use of several special functions.     

Mean, Median and Mode in Binomial Distributions

Summary  While studying the median of the binomial distribution we discovered that the mean median-mode inequality, recently discussed in. Statistics Neerlandica by R unnen - burg 141 and V an Z wet [7] for continuous distributions, does not hold for the binomial distribution. If the mean is an integer, then mean = median = mode. In theorem 1 a sufficient condition is given for mode = median = rounded mean. If median and mode differ, the mean lies in between.