June 9th 2015
Faculty of Economic and Business Sciences | Salón de Grados
Generalized exponential tilting
2015/06/09 – 12:30 h | Paul Eilers – Erasmus University Medical Centre, Rotterdam, the Netherlands
Paul Eilers has a guest position at the Erasmus University Medical Centre, where he retired in 2013 as professor of Genetical Statistics. Besides in genetics, he has been working in diverse fields: mechanical vibrations, aquatic ecology, air and soil pollution, medical statistics and the social sciences. He was trained as an electronic engineer. He is interested in statistical computation, time series analysis, high-dimensional data analysis and chemometrics. His real love is smoothing of complex data.
Abstract
Suppose we are given a density and we ask for the density that is closest to it, but has a different expectation. If we express closeness using the Kullback-Leibler distance, we find that we have to add a linear function to the logarithm of the given density. This is equivalent to multiplication with an exponential function, leading to the name exponential tilting.
We can go further and specify expected values of a set of functions. Then we find that we have to add a linear combination of these functions to the logarithm of the given density. The proper values of the coefficients can be computed efficiently with Newton-Raphson iterations.
Suppose now that we are given a set of densities. We can ask the reverse question: can we find a “mother density” such that the observed densities are (approximately) the results of exponential tilting? The answer is affirmative, but it is of little value if we work with observed data. However, starting from histograms with narrow bins and applying penalized Poisson regression (to get a smooth estimate) we can obtain excellent results. We call this exploratory exponential tilting (EET).
In EET it is assumed that the tilting functions are known. A more ambitious goal is to estimate them from the data in a semi-parametric way. As will be shown, this goal is attainable. It can lead to a parsimonious but accurate description of large sets of densities. Also the patterns in the tilting functions can teach us something about the underlying processes. We call this generalized exponential tilting (GET).
GET can be extended with penalties. This offers a new approach to two-dimensional density estimation.
Applications to real data show the usefulness of the propose variants of exponential tilting.
This is joint work with Giancarlo Camarda (Institut Nationale d’Études Démographiques, Paris) and Jutta Gampe (Max Plank Institute for Demographic Research, Rostock).
Gallery