A Convolution Estimator for the Density of Nonlinear Regression Observations

Abstract. The problem of estimating an unknown density function has been widely studied. In this article, we present a convolution estimator for the density of the responses in a nonlinear heterogenous regression model. The rate of convergence for the mean square error of the convolution estimator is of order n ^?1 under certain regularity conditions. This is faster than the rate for the kernel density method. We derive explicit expressions for the asymptotic variance and the bias of the new estimator, and further a data‐driven bandwidth selector is proposed. We conduct simulation experiments to check the finite sample properties, and the convolution estimator performs substantially better than the kernel density estimator for well‐behaved noise densities.     

Ranked-Weighted Utilities and Qualitative Convolution

For gambles—non-numerical consequences attached to uncertain chance events—analogues are proposed for the sum of independent random variables and their convolution. Joint receipt of gambles is the analogue of the sum of random variables. Because it has no unique expansion as a first-order gamble analogous to convolution, a definition of qualitative convolution is proposed. Assuming ranked, weighted-utility representations (RWU) over gains (and, separately, over losses, but not mixtures of both), conditions are given for the equivalence of joint receipt, qualitative convolution, and a utility expression like expected value. As background, some properties of RWU are developed.     

On univariate and bivariate generalized gamma convolutions

This paper has two parts. In the first part some results for generalized gamma convolutions (GGCs) are reviewed. A GGC is a limit distribution for sums of independent gamma variables. In the second part, bivariate gamma distributions and bivariate GGCs are considered. New bivariate gamma distributions are derived from shot-noise models. The remarkable property hyperbolic complete monotonicity (HCM) for a function is considered both in the univariate case and in the bivariate case.     

Ordering Properties of Convolutions of Exponential Random Variables

Convolutions of independent random variables are usually compared. In this paper, after a synthetic comparison with respect to hazard rate ordering between sums of independent exponential random variables, we focus on the special case where one sum is identically distributed. So, for a given sum of n independent exponential random variables, we deduce the "best" Erlang-n bounds, with respect to each of the usual orderings: mean ordering, stochastic ordering, hazard rate ordering and likelihood ratio ordering.     

Evaluation of Bayesian multiple stage estimation under spatial CAR model variants

In this study, an evaluation of Bayesian hierarchical models is made based on simulation scenarios to compare single-stage and multi-stage Bayesian estimations. Simulated datasets of lung cancer disease counts for men aged 65 and older across 44 wards in the London Health Authority were analysed using a range of spatially structured random effect components. The goals of this study are to determine which of these single-stage models perform best given a certain simulating model, how estimation methods (single- vs. multi-stage) compare in yielding posterior estimates of fixed effects in the presence of spatially structured random effects, and finally which of two spatial prior models – the Leroux or ICAR model, perform best in a multi-stage context under different assumptions concerning spatial correlation. Among the fitted single-stage models without covariates, we found that when there is low amount of variability in the distribution of disease counts, the BYM model is relatively robust to misspecification in terms of DIC, while the Leroux model is the least robust to misspecification. When these models were fit to data generated from models with covariates, we found that when there was one set of covariates – either spatially correlated or non-spatially correlated, changing the values of the fixed coefficients affected the ability of either the Leroux or ICAR model to fit the data well in terms of DIC. When there were multiple sets of spatially correlated covariates in the simulating model, however, we could not distinguish the goodness of fit to the data between these single-stage models. We found that the multi-stage modelling process via the Leroux and ICAR models generally reduced the variance of the posterior estimated fixed effects for data generated from models with covariates and a UH term compared to analogous single-stage models. Finally, we found the multi-stage Leroux model compares favourably to the multi-stage ICAR model in terms of DIC. We conclude that the mutli-stage Leroux model should be seriously considered in applications of Bayesian disease mapping when an investigator desires to fit a model with both fixed effects and spatially structured random effects to Poisson count data.     

On The Method of Discrete Probability Distributions in Risk and Reliability Calculations–Application to Seismic Risk Assessment

If the point of view is adopted that in calculations of real-world phenomena we almost invariably have significant uncertainty in the numerical values of our parameters, then, in these calculations, numerical quantities should be replaced by probability distributions and mathematical operations between these quantities should be replaced by analogous operations between probability distributions. Also, practical calculations one way or another always require discretization or truncation. Combining these two thoughts leads to a numerical approach to probabilistic calculations having great simplicity, power, and elegance. The philosophy and technique of this approach is described, some pitfalls are pointed out, and an application to seismic risk assessment is outlined.     

Modeling material stress using integrated Gaussian Markov random fields

The equations of a physical constitutive model for material stress within tantalum grains were solved numerically using a tetrahedrally meshed volume. The resulting output included a scalar vonMises stress for each of the more than 94,000 tetrahedra within the finite element discretization. In this paper, we define an intricate statistical model for the spatial field of vonMises stress which uses the given grain geometry in a fundamental way. Our model relates the three-dimensional field to integrals of latent stochastic processes defined on the vertices of the one- and two-dimensional grain boundaries. An intuitive neighborhood structure of the said boundary nodes suggested the use of a latent Gaussian Markov random field (GMRF). However, despite the potential for computational gains afforded by GMRFs, the integral nature of our model and the sheer number of data points pose substantial challenges for a full Bayesian analysis. To overcome these problems and encourage efficient exploration of the posterior distribution, a number of techniques are now combined: parallel computing, sparse matrix methods, and a modification of a block update strategy within the sampling routine. In addition, we use an auxiliary variables approach to accommodate the presence of outliers in the data.     

Some bivariate families of lagrangian probability distributions

The bivariate Lagrange expansion, given by Poincare (1986), has been explained and slightly modified which gives bivariate Lagrangian probability models. A generalized bivariate Lagrangian Poisson distribution with six parameters has been obtained and studied. Also, the bivariate Lagrangian binomial, bivariate Lagrangian negative binomial and bivariate Lagrangian logarithmic series distribution have been obtained.     

Robust hierarchical Bayes estimation of exchangeable means

Estimation of the mean of a multivariate normal distribution is considered. The components of the mean vector θ are assumed to be exchangeable; this is modelled in a hierarchical fashion with independent Cauchy distributions as the first-stage prior. The resulting generalized Bayes estimator is calculated and shown to be robust with respect to the presence of outlying means. Alternative estimators that have similar behaviour but are cheaper to compute are also derived.     

广义卷积代数与Hun半群及（，＊_α）的半群结构

利用Ｈｕｎ半群理论，证明了任一正则广义卷积代数，按广义卷积运算和弱收敛拓扑构成一个可度量化，稳定，可模的Ｈｕｎ半群，且无除么元以外的幂等元，并研究了（，＊α）的半群结构．