Histogram-kernel error and its application for bin width selection in histograms

Histogram and kernel estimators are usually regarded as the two main classical data-based non- parametric tools to estimate the underlying density functions for some given data sets. In this paper we will integrate them and define a histogram-kernel error based on the integrated square error between histogram and binned kernel density estimator, and then exploit its asymptotic properties. Just as indicated in this paper, the histogram-kernel error only depends on the choice of bin width and the data for the given prior kernel densities. The asymptotic optimal bin width is derived by minimizing the mean histogram-kernel error. By comparing with Scott’s optimal bin width formula for a histogram, a new method is proposed to construct the data-based histogram without knowledge of the underlying density function. Monte Carlo study is used to verify the usefulness of our method for different kinds of density functions and sample sizes.     

A posteriori finite element error estimators for indefinite elliptic boundary value problems ∗

A general construction technique is presented for a posteriori error estimators of finite element solutions of elliptic boundary value problems that satisfy a Gång inequality. The estimators are obtained by an element–by–element solution of ‘weak residual’ with or without considering element boundary residuals. There is no order restriction on the finite element spaces used for the approximate solution or the error estimation; that is, the design of the estimators is applicable in connection with either one of the h–p–, or hp– formulations of the finite element method. Under suitable assumptions it is shown that the estimators are bounded by constant multiples of the true error in a suitable norm. Some numerical results are given to demonstrate the effectiveness and efficiency of the approach.     

A discontinuous Galerkin method for the Rosenau equation

A priori error estimates for the Rosenau equation, which is a K-dV like Rosenau equation modelled to describe the dynamics of dense discrete systems, have been studied by one of the authors. But since a priori error bounds contain the unknown solution and its derivatives, it is not effective to control error bounds with only a given step size. Thus we need to estimate a posteriori errors in order to control accuracy of approximate solutions using variable step sizes. A posteriori error estimates of the Rosenau equation are obtained by a discontinuous Galerkin method and the stability analysis is discussed for the dual problem. Numerical results on a posteriori error and wave propagation are given, which are obtained by using various spatial and temporal meshes controlled automatically by a posteriori error.     

Second-order nonlinear least squares estimation

The ordinary least squares estimation is based on minimization of the squared distance of the response variable to its conditional mean given the predictor variable. We extend this method by including in the criterion function the distance of the squared response variable to its second conditional moment. It is shown that this “second-order” least squares estimator is asymptotically more efficient than the ordinary least squares estimator if the third moment of the random error is nonzero, and both estimators have the same asymptotic covariance matrix if the error distribution is symmetric. Simulation studies show that the variance reduction of the new estimator can be as high as 50% for sample sizes lower than 100. As a by-product, the joint asymptotic covariance matrix of the ordinary least squares estimators for the regression parameter and for the random error variance is also derived, which is only available in the literature for very special cases, e.g. that random error has a normal distribution. The results apply to both linear and nonlinear regression models, where the random error distributions are not necessarily known.     

Guaranteed a posteriori error estimates for nonconforming finite element approximations to a singularly perturbed reaction–diffusion problem

We derive guaranteed a posteriori error estimates for nonconforming finite element approximations to a singularly perturbed reaction–diffusion problem. First, an abstract a posteriori error bound is derived under a special equilibration condition. Based on conservative flux reconstruction, two error estimators are proposed and provide actual upper error bounds in the usual energy norm without unknown constants, one of which can be directly constructed without solving local Neumann problems and provide practical computable error bounds. The error estimators also provide local lower bounds but with the multiplicative constants dependent on the diffusion coefficient and mesh size, where the constants can be bounded for enough small mesh size comparable with the square root of the diffusion coefficient. By adding edge jumps with weights to the energy norm, two modified error estimators with additional edge tangential jumps are shown to be robust with respect to the diffusion coefficient and provide guaranteed upper bounds on the error in the modified norm. Finally, the performance of the estimators are illustrated by the numerical results.     

Probability density estimation with surrogate data and validation sample

The probability density estimation problem with surrogate data and validation sample is considered. A regression calibration kernel density estimator is defined to incorporate the information contained in both surrogate variates and validation sample. Also, we define two weighted estimators which have less asymptotic variances but have bigger biases than the regression calibration kernel density estimator. All the proposed estimators are proved to be asymptotically normal. And the asymptotic representations for the mean squared error and mean integrated square error of the proposed estimators are established, respectively. A simulation study is conducted to compare the finite sample behaviors of the proposed estimators.     

Gradient recovery type a posteriori error estimates of virtual element method for an elliptic variational inequality of the second kind

In this paper, gradient recovery type a posteriori error estimators of virtual element discretization are derived for a simplified friction problem, which is a typical elliptic variational inequality of the second kind. Both the reliability and the efficiency of the error estimators are proved. In addition, one numerical example is presented to show the efficiency of the adaptive VEM based on the derived error estimators.     

Residual-based a posteriori error estimates of nonconforming finite element method for elliptic problems with Dirac delta source terms

Two residual-based a posteriori error estimators of the nonconforming Crouzeix-Raviart element are derived for elliptic problems with Dirac delta source terms.One estimator is shown to be reliable and efficient,which yields global upper and lower bounds for the error in piecewise W1,p seminorm.The other one is proved to give a global upper bound of the error in Lp-norm.By taking the two estimators as refinement indicators,adaptive algorithms are suggested,which are experimentally shown to attain optimal convergence orders.     

A posteriori FE error control for p-Laplacian by gradient recovery in quasi-norm

A posteriori error estimators based on quasi-norm gradient recovery are established for the finite element approximation of the p-Laplacian on unstructured meshes. The new a posteriori error estimators provide both upper and lower bounds in the quasi-norm for the discretization error. The main tools for the proofs of reliability are approximation error estimates for a local approximation operator in the quasi-norm.

     

Kernel density estimation for linear processes: Asymptotic normality and optimal bandwidth derivation

The problem of estimating the marginal density of a linear process by kernel methods is considered. Under general conditions, kernel density estimators are shown to be asymptotically normal. Their limiting covariance matrix is computed. We also find the optimal bandwidth in the sense that it asymptotically minimizes the mean square error of the estimators. The assumptions involved are easily verifiable.Supported in part by NSF grant DMS-9403718.     

鞅差误差序列下半参数EV回归模型的近邻估计

本文研究了误差为鞅差序列的条件下的一维半参数EV回归模型.利用两步估计的方法构造了参数分量和非参数分量的近邻估计,并且分别证明了估计量的L2相合性和强相合性,从而推广了在普通半参数回归模型已有的相关结论.     

Local a posteriori estimates for pointwise gradient errors in finite element methods for elliptic problems

We prove local a posteriori error estimates for pointwise gradient errors in finite element methods for a second-order linear elliptic model problem. First we split the local gradient error into a computable local residual term and a weaker global norm of the finite element error (the ``pollution term'). Using a mesh-dependent weight, the residual term is bounded in a sharply localized fashion. In specific situations the pollution term may also be bounded by computable residual estimators. On nonconvex polygonal and polyhedral domains in two and three space dimensions, we may choose estimators for the pollution term which do not employ specific knowledge of corner singularities and which are valid on domains with cracks. The finite element mesh is only required to be simplicial and shape-regular, so that highly graded and unstructured meshes are allowed.

     

对流扩散方程迎风有限元的自适应方法

本文对二维发展型对流扩散方程的迎风有限元格式给出了显式后验误差估计,证明了真实误差被后验误差估计器上下界定;并通过误差估计器建立了相应的自适应算法,数值例子表明了方法的有效性.     

When is the inverse regression estimator MSE-superior to the standard regression estimator in multivariate controlled calibration situations?

We assume as model a standard multivariate regression of y on x, fitted to a controlled calibration sample and used to estimate unknown x′s from observed y-values. The standard weighted least squares estimator (‘classical’, regress y on x and ‘solve’ for x) and the biased inverse regression estimator (regress x on y) are compared with respect to mean squared error. The regions are derived where the inverse regression estimator yields the smaller MSE. For any particular component of x this region is likely to contain ‘most’ future values in usual practice. For simultaneous estimation this needs not be true, however.     

A posteriori error estimates for discontinuous Galerkin methods of obstacle problems

We present a posteriori error analysis of discontinuous Galerkin methods for solving the obstacle problem, which is a representative elliptic variational inequality of the first kind. We derive reliable error estimators of the residual type. Efficiency of the estimators is theoretically explored and numerically confirmed.     

AN ADAPTIVE VERSION OF GLIMM＇S SCHEME

This article describes a local error estimator for Glimm＇s scheme for hyperbolic systems of conservation laws and uses it to replace the usual random choice in Glimm＇s scheme by an optimal choice. As a by-product of the local error estimator, the procedure provides a global error estimator that is shown numerically to be a very accurate estimate of the error in L1 （R） for all times. Although there is partial mathematical evidence for the error estimator proposed, at this stage the error estimator must be considered ad- hoc. Nonetheless, the error estimator is simple to compute, relatively inexpensive, without adjustable parameters and at least as accurate as other existing error estimators. Numerical experiments in 1-D for Burgers＇ equation and for Euler＇s system are performed to measure the asymptotic accuracy of the resulting scheme and of the error estimator.     

特征值问题协调/非协调有限元后验误差分析

本文研究对称椭圆特征值问题的有限元后验误差估计,包括协调元和非协调元,具有下列特色:(1)对协调/非协调元建立了有限元特征函数uh的误差与相应的边值问题有限元解的误差在局部能量模意义下的恒等关系式,该边值问题的右端为有限元特征值λh与uh的乘积,有限元解恰好为uh.从而边值问题有限元解在能量模意义下的局部后验误差指示子,包括残差型和重构型后验误差指示子,成为有限元特征函数在能量模意义下的局部后验误差指示子.(2)讨论了协调有限元特征函数的基于插值后处理的梯度重构型后验误差估计,对有限元特征函数的导数得到了最大模意义下的渐近准确局部后验误差指示子.     

AN ADAPTIVE VERSION OF GLIMM'S SCHEME

This article describes a local error estimator for Glimm's scheme for hyperbolic systems of conservation laws and uses it to replace the usual random choice in Glimm's scheme by an optimal choice.As a by-product of the local error estimator,the procedure provides a global error estimator that is shown numerically to be a very accurate estimate of the error in L1(R) for all times.Although there is partial mathematical evidence for the error estimator proposed,at this stage the error estimator must be considered adhoc.Nonetheless,the error estimator is simple to compute,relatively inexpensive,without adjustable parameters and at least as accurate as other existing error estimators.Numerical experiments in 1-D for Burgers' equation and for Euler's system are performed to measure the asymptotic accuracy of the resulting scheme and of the error estimator.     

Optimal Control of Population Dynamics

This paper studies the role of projection algorithms in conditional set membership estimation. These algorithms are known to be suboptimal in terms of the worst-case estimation error. A tight upper bound on the error of central projection estimators and interpolatory projection estimators is computed as a function of the conditional radius of information. Since the radius of information represents the minimum achievable error, the derived bound provides a measure of the reliability level of the suboptimal algorithms. The results are derived in a general deterministic setting, which allows the consideration of linearly parametrized approximations of a compact set of feasible problem elements.     

��Ԫ���ṹBerkson�������ģ�͵ķ���

Berkson测量误差模型在工业、农业、流行病学、经济学 等领域有广泛的应用. 本文考虑了一类多元超结构Berkson测量误差模型, 给出了该模型中参数的相合估计, 推导了估计的渐近分布, 并把该方法应用到一元超结构Berkson测量误差模型中, 最后给出了模拟结果和实例分析来说明本文所提出的估计方法的表现.