A theory of statistical models for Monte Carlo integration

Summary. The task of estimating an integral by Monte Carlo methods is formulated as a statistical model using simulated observations as data. The difficulty in this exercise is that we ordinarily have at our disposal all of the information required to compute integrals exactly by calculus or numerical integration, but we choose to ignore some of the information for simplicity or computational feasibility. Our proposal is to use a semiparametric statistical model that makes explicit what information is ignored and what information is retained. The parameter space in this model is a set of measures on the sample space, which is ordinarily an infinite dimensional object. None-the-less, from simulated data the base-line measure can be estimated by maximum likelihood, and the required integrals computed by a simple formula previously derived by Vardi and by Lindsay in a closely related model for biased sampling. The same formula was also suggested by Geyer and by Meng and Wong using entirely different arguments. By contrast with Geyer's retrospective likelihood, a correct estimate of simulation error is available directly from the Fisher information. The principal advantage of the semiparametric model is that variance reduction techniques are associated with submodels in which the maximum likelihood estimator in the submodel may have substantially smaller variance than the traditional estimator. The method is applicable to Markov chain and more general Monte Carlo sampling schemes with multiple samplers.     

Pricing power options with a generalized jump diffusion

In this article, the valuation of power option is investigated when the dynamic of the stock price is governed by a generalized jump-diffusion Markov-modulated model. The systematic risk is characterized by the diffusion part, and the non systematic risk is characterized by the pure jump process. The jumps are described by a generalized renewal process with generalized jump amplitude. By introducing NASDAQ Index Model, their risk premium is identified respectively. A risk-neutral measure is identified by employing Esscher transform with two families of parameters, which represent the two parts risk premium. In this article, the non systematic risk premium is considered, based on which the price of power option is studied under the generalized jump-diffusion Markov-modulated model. In the case of a special renewal process with log double exponential jump amplitude, the accurate expressions for the Esscher parameters and the pricing formula are provided. By numerical simulation, the influence of the non systematic risk’s price and the index of the power options on the price of the option is depicted.     

Extreme values identification in regression using a peaks-over-threshold approach

The problem of heavy tail in regression models is studied. It is proposed that regression models are estimated by a standard procedure and a statistical check for heavy tail using residuals is conducted as a tool for regression diagnostic. Using the peaks-over-threshold approach, the generalized Pareto distribution quantifies the degree of heavy tail by the extreme value index. The number of excesses is determined by means of an innovative threshold model which partitions the random sample into extreme values and ordinary values. The overall decision on a significant heavy tail is justified by both a statistical test and a quantile–quantile plot. The usefulness of the approach includes justification of goodness of fit of the estimated regression model and quantification of the occurrence of extremal events. The proposed methodology is supplemented by surface ozone level in the city center of Leeds.     

Efficient inference for parameters of unobservable periodic autoregressive time series

This paper considers estimating the model coefficients when the observed periodic autoregressive time series is contaminated by a trend. The proposed Yule–Walker estimators are obtained by a two-step procedure. In the first step, the trend is estimated by a weighted local polynomial, and the residuals are obtained by subtracting the trend estimates from the observations; in the second step, the model coefficients are estimated by the well-known Yule–Walker method via the residuals. It is shown that under certain conditions such Yule–Walker estimators are oracally efficient, i.e., they are asymptotically equivalent to those obtained from periodic autoregressive time series without a trend. An easy-to-use implementation procedure is provided. The performance of the estimators is illustrated by simulation studies and real data analysis. In particular, the simulation studies show that the proposed estimator outperforms that obtained from the residuals when the trend is estimated by kernel smoothing without taking the heteroscedasticity into consideration.     

A simple and asymptotically optimal test for the equality of k(≥2) exponential distributions based on type II censored samples

A test for the equality of two or more two-parameter exponential distributions is suggested. It is developed on an intuitive basis and is obtained by combining two independent tests by the Fisher method (1950, pp. 99-101). The test is simple for application and is optimal asymptotically in the sense of Bahadur efficiency (1960). A numerical example is discussed to illustrate its application in a real-world situation. The Monte Carlo simulation is used for calculating its power which is compared with that of the test suggested by Singh and Narayan (1983). The suggested test is found oftener more powerful.     

A precise estimator for the log-normal mean

The log-normal distribution is frequently encountered in applications. The uniformly minimum variance unbiased (UMVU) estimator for the log-normal mean is given explicitly by a formula found by Finney in 1941. In contrast to this the most commonly used estimator for a log-normal mean is the sample mean. This is possibly due to the complexity of the formula given by Finney. A modified maximum likelihood estimator which approximates the UMVU estimator is derived here. It is sufficiently simple to be implemented in elementary spreadsheet applications. An elementary approximate formula for the root-mean-square error of the suggested estimator and the UMVU estimator is presented. The suggested estimator is compared with the sample mean, the maximum likelihood, and the UMVU estimators by Monte Carlo simulation in terms of root-mean-square error.     

A Bounded Derivative Model for Prior Ignorance about a Real-valued Parameter

A new method is proposed for drawing coherent statistical inferences about a real-valued parameter in problems where there is little or no prior information. Prior ignorance about the parameter is modelled by the set of all continuous probability density functions for which the derivative of the log-density is bounded by a positive constant. This set is translation-invariant, it contains density functions with a wide variety of shapes and tail behaviour, and it generates prior probabilities that are highly imprecise. Statistical inferences can be calculated by solving a simple type of optimal control problem whose general solution is characterized. Detailed results are given for the problems of calculating posterior upper and lower means, variances, distribution functions and probabilities of intervals. In general, posterior upper and lower expectations are achieved by prior density functions that are piecewise exponential. The results are illustrated by normal and binomial examples     

A penalized likelihood approach to image warping

A warping is a function that deforms images by mapping between image domains. The choice of function is formulated statistically as maximum penalized likelihood, where the likelihood measures the similarity between images after warping and the penalty is a measure of distortion of a warping. The paper addresses two issues simultaneously, of how to choose the warping function and how to assess the alignment. A new, Fourier–von Mises image model is identified, with phase differences between Fourier-transformed images having von Mises distributions. Also, new, null set distortion criteria are proposed, with each criterion uniquely minimized by a particular set of polynomial functions. A conjugate gradient algorithm is used to estimate the warping function, which is numerically approximated by a piecewise bilinear function. The method is motivated by, and used to solve, three applied problems: to register a remotely sensed image with a map, to align microscope images obtained by using different optics and to discriminate between species of fish from photographic images.     

Computational methods for local regression

Local regression is a nonparametric method in which the regression surface is estimated by fitting parametric functions locally in the space of the predictors using weighted least squares in a moving fashion similar to the way that a time series is smoothed by moving averages. Three computational methods for local regression are presented. First, fast surface fitting and evaluation is achieved by building ak-d tree in the space of the predictors, evaluating the surface at the corners of the tree, and then interpolating elsewhere by blending functions. Second, surfaces are made conditionally parametric in any proper subset of the predictors by a simple alteration of the weighting scheme. Third degree-of-freedom quantities that would be extremely expensive to compute exactly are approximated, not by numerical methods, but through a statistical model that predicts the quantities from the trace of the hat matrix, which can be computed easily.     

最低工资标准下的个税免征额区域化设计--以北京市为例

目前,我国个人所得税免征额的设定采用的是全国统一的标准,这种方式虽然比较简便,但是忽略了我国各地区经济发展的差异性。对于不同生活水平的人群按照统一标准进行征税,显然缺乏合理性。有鉴于此,本文结合我国的实际,借鉴国际上一些国家比较成熟的做法,提出在免征额设定时可以通过将各地区加以分类,然后以最低工资标准为基础设计免征额,并以北京市为例进行了说明。     

A new Bayesian wavelet thresholding estimator of nonparametric regression

The methods of estimation of nonparametric regression function are quite common in statistical application. In this paper, the new Bayesian wavelet thresholding estimation is considered. The new mixture prior distributions for the estimation of nonparametric regression function by applying wavelet transformation are investigated. The reversible jump algorithm to obtain the appropriate prior distributions and value of thresholding is used. The performance of the proposed estimator is assessed with simulated data from well-known test functions by comparing the convergence rate of the proposed estimator with respect to another by evaluating the average mean square error and standard deviations. Finally by applying the developed method, density function of galaxy data is estimated.     

货币政策中介目标M2的修正和新生

M2是传统的货币政策中介目标,然而,随着经济和金融环境的变化,它的有效性备受质疑。在没有更好的替代指标出现以前,探讨修正的M2是学术界的热点。迪维西亚法和现金等价法都是基于流动性总量的测度来修正简单加总法测度下的M2。本文通过更为简洁的模型推导,首次在构建原理上论证,在测度M2的流动性总量时,现金等价法的精确性明显优于迪维西亚法。本文还基于货币需求函数的稳定性和选择货币政策中介目标的可控性、相关性原则,对简单加总法、迪维西亚法和现金等价法测度下的M2的合理性进行实证比较,结果表明,现金等价法测度下的M2在货币需求函数的稳定性以及作为货币政策中介目标的可控性、相关性上都明显优于简单加总法和迪维西亚法测度下的M2。因此,本文认为,应该用现金等价法来修正传统M2,使其能够准确反映流动性总量的变化,提升其作为货币政策中介目标的效力。     

Wavelet Threshold Estimators for Data with Correlated Noise

Wavelet threshold estimators for data with stationary correlated noise are constructed by applying a level-dependent soft threshold to the coefficients in the wavelet transform. A variety of threshold choices is proposed, including one based on an unbiased estimate of mean-squared error. The practical performance of the method is demonstrated on examples, including data from a neurophysiological context. The theoretical properties of the estimators are investigated by comparing them with an ideal but unattainable `bench-mark', that can be considered in the wavelet context as the risk obtained by ideal spatial adaptivity, and more generally is obtained by the use of an `oracle' that provides information that is not actually available in the data. It is shown that the level-dependent threshold estimator performs well relative to the bench-mark risk, and that its minimax behaviour cannot be improved on in order of magnitude by any other estimator. The wavelet domain structure of both short- and long-range dependent noise is considered, and in both cases it is shown that the estimators have near optimal behaviour simultaneously in a wide range of function classes, adapting automatically to the regularity properties of the underlying model. The proofs of the main results are obtained by considering a more general multivariate normal decision theoretic problem.     

On the Asymptotic Distribution of the Likelihood Ratio Statistic

Some results on the estimation of a symmetric density function are given. For the case when the point of symmetry, θ, is known it is shown that a symmetrized kernel estimator is, as measured by MISE, approximately as good as a non-symmetrized one based on twice as many observations. This result remains ’true if the estimated density is a normal one and θ is estimated by the sample mean. Some Monte Carlo results for several densities and sample sizes are given for the case when θ is estimated by the sample median.     

An alpha-power extension for the Birnbaum–Saunders distribution

Data on fatigue to exceed a critical value (or to grow to a critical level at which failure is likely to occur) is typically adjusted using the Birnbaum–Saunders (BS) distribution [see Birnbaum ZW, Saunders SC. A new family of life distributions. J Appl Probab. 1969a;6:319–327]. Although this type of distribution is asymmetric, in some cases the degree of skewness and/or kurtosis are outside the distributional range allowed by the BS distribution. Thus, a more adequate distribution model for better adjusting such unexpected deviations is called for. With this in mind, the main object of this paper is to propose an extension of the BS distribution based on the asymmetric alpha-power family of distributions [see Pewsey A, Gómez HW, Bolfarine H. Likelihood-based inference for power distributions. Test. 2012;21(4):775–789]. This extension is called the exponentiated BS distribution. We expect that by replacing the normal distribution by such more general family, a more flexible BS family is obtained. Asymmetry in the alpha-power family is controlled by a shape parameter, which also presents a similar performance in the extended BS family. The paper presents the density function for the extended BS and derives closed-form expressions for moments. Estimation is dealt with by using maximum likelihood estimators. Large sample inference can be conducted by using the Fisher information matrix derived in the paper. Estimation performance is studied by using a small scale simulation study. Results of a real application illustrates the good performance of the proposed approach.     

Nonparametric estimation of multivariate extreme-value copulas

Extreme-value copulas arise in the asymptotic theory for componentwise maxima of independent random samples. An extreme-value copula is determined by its Pickands dependence function, which is a function on the unit simplex subject to certain shape constraints that arise from an integral transform of an underlying measure called spectral measure. Multivariate extensions are provided of certain rank-based nonparametric estimators of the Pickands dependence function. The shape constraint that the estimator should itself be a Pickands dependence function is enforced by replacing an initial estimator by its best least-squares approximation in the set of Pickands dependence functions having a discrete spectral measure supported on a sufficiently fine grid. Weak convergence of the standardized estimators is demonstrated and the finite-sample performance of the estimators is investigated by means of a simulation experiment.     

1978年以来中国铁路运输区域差异统计分析与评价

区域经济发展的不平衡必然使得区域铁路运输具有差异性。利用1978-2005年的铁路运输数据，采用重心模型、变异系数和加权变异系数，从客运量、客运周转量、货运量和货运周转量4个方面，对中国铁路运输的区域差异进行了统计分析和评价。结果发现：1．客运量重心在逐步向西南方向移动，其中向西移动了2经度，向南移动了2．6纬度；铁路客运量区域差异在逐步减小。2．客运周转量重心呈直线形逐步向西南方向移动，其中，重心向西移动了2．3经度，向南移动了2．9纬度；铁路客运周转量区域差异整体上呈减小趋势，但其中波动较大，可以划分为5个阶段。3．货运量重心波动较大；货运量区域差异在逐渐增大，而且区域差异的增长幅度明显加大，两极分化越来越严重。4．货运周转量重心变化轨迹略呈指数型；货运周转量区域差异呈先减小后增大的趋势，2000年以前区域差异在减小，2000年以后区域差异在增大。     

Markov random-field models for estimating local labour markets

Summary. This work is motivated by data on daily travel-to-work flows observed between pairs of elemental territorial units of an Italian region. The data were collected during the 1991 population census. The aim of the analysis is to partition the region into local labour markets. We present a new method for this which is inspired by the Bayesian texture segmentation approach. We introduce a novel Markov random-field model for the distribution of the variables that label the local labour markets for each territorial unit. Inference is performed by means of Markov chain Monte Carlo methods. The issue of model hyperparameter estimation is also addressed. We compare the results with those obtained by applying a classical method. The methodology can be applied with minor modifications to other data sets.