Utilisation d’une fonction de perte quadratique pondérée dans les approches bayésiennes : application à la cartographie d’indicateur de santé

Léa Fortunato; Adeline Liverneaux; Denis Hémon; Chantal Guihenneuc-Jouyaux

Journal de la société française de statistique (2007)

  • Volume: 148, Issue: 3, page 71-86
  • ISSN: 1962-5197

Abstract

top
Background In the context of ecological model, we compare different Bayes estimates of relative risks by disease mapping model. Methods Bayesian estimates are based on the choice of loss function. A usual choice is the quadratic loss function and its minimization leads to the posterior mean. In this study, a discussion is done on the advantages of this choice, in particular by considering different objectives as focussing on extreme relative risks in disease mapping. These objectives determine the choice of the loss function. Three different loss functions are considered: the quadratic loss, the absolute loss and a weigthed quadratic loss function with weights having an exponential shape. The last one corresponds to an unusual choice with important weights on extreme relative risks. Results A study based on simulations is done with three different cases of relative risks in France: two first cases correspond to discontinuous risk surfaces while the last case presents continuous risk decrease following a geographic gradient. Estimates performances are evaluated by mean quadratic error and bias as well as ability of units detection with high relative risks. Results show that the different bayesian estimates are coherent with good performances. Nevertheless, often significative improvements of estimates based on exponential loss function are found in the case of discontinuous risk surface. Conclusion The choice of loss function and the study of the associated estimate are rarely done and nevertheless interesting in bayesian approach. This choice is based on a specific objective as, in geographical analyse, a particular attention on extreme relative risks. In this study, slight differences between the three considered loss functions are found showing in this sense their robustness but with often significant improvements of weighted estimate. Different configurations of relative risks and/or expected number of cases could give more important differences.

How to cite

top

Fortunato, Léa, et al. "Utilisation d’une fonction de perte quadratique pondérée dans les approches bayésiennes : application à la cartographie d’indicateur de santé." Journal de la société française de statistique 148.3 (2007): 71-86. <http://eudml.org/doc/93466>.

@article{Fortunato2007,
abstract = {Position du problème Dans le cadre des modèles écologiques, nous comparons différents estimateurs bayésiens des risques relatifs en utilisant un modèle de cartographie hiérarchique bayésien. Méthodes Les estimateurs ponctuels bayésiens dépendent du choix de la fonction de perte. La plus « classique » est la fonction de perte quadratique et sa minimisation amène à l’estimation par moyenne a posteriori. Une discussion est menée sur les avantages de cette approche, notamment en envisageant différents objectifs comme l’intérêt particulier porté sur les valeurs extrêmes des risques relatifs lors de cartographie. Ces objectifs guident le choix de la fonction de perte. Trois fonctions sont envisagées : la fonction de perte quadratique, la fonction de perte en valeur absolue et une fonction de perte quadratique pondérée avec des poids de type exponentiel. La dernière fonction correspond à un choix moins fréquent avec une pondération importante sur les risques relatifs extrêmes. Résultats Une étude par simulations a été réalisée en considérant trois modèles de risques relatifs. Les deux premiers correspondent à des surfaces de risques discontinues alors que le dernier cas présente une décroissance du risque sur le territoire selon un gradient géographique. Les performances des estimateurs sont évaluées par l’erreur quadratique et le biais relatif ainsi que la capacité à détecter les zones à risques élevés. Les résultats montrent que de manière générale, les différents estimateurs des risques relatifs sont convergents avec de bonnes performances. Néanmoins, cette étude permet de montrer que dans le cas de risques relatifs discontinus, le choix d’une fonction de perte de type exponentiel permet d’obtenir une amélioration souvent significative de l’estimation des risques relatifs. Conclusion Le choix d’une fonction de perte et l’étude de l’estimateur qui en découle constituent une démarche rarement utilisée et pourtant prometteuse dans les approches bayésiennes. Ce choix peut être guidé par l’objectif recherché comme ici, la mise en valeur des risques relatifs extrêmes. Cette étude a montré de légères différences entre les estimateurs, prouvant dans ce sens leur robustesse mais avec une amélioration souvent significative de l’estimateur pondéré. D’autres configurations de risques relatifs et/ou de nombre de cas attendus pourraient peut-être mettre à jour des différences plus importantes.},
author = {Fortunato, Léa, Liverneaux, Adeline, Hémon, Denis, Guihenneuc-Jouyaux, Chantal},
journal = {Journal de la société française de statistique},
keywords = {hierarchical bayesian model; loss functions; relative risks; ecological Poisson model},
language = {fre},
number = {3},
pages = {71-86},
publisher = {Société française de statistique},
title = {Utilisation d’une fonction de perte quadratique pondérée dans les approches bayésiennes : application à la cartographie d’indicateur de santé},
url = {http://eudml.org/doc/93466},
volume = {148},
year = {2007},
}

TY - JOUR
AU - Fortunato, Léa
AU - Liverneaux, Adeline
AU - Hémon, Denis
AU - Guihenneuc-Jouyaux, Chantal
TI - Utilisation d’une fonction de perte quadratique pondérée dans les approches bayésiennes : application à la cartographie d’indicateur de santé
JO - Journal de la société française de statistique
PY - 2007
PB - Société française de statistique
VL - 148
IS - 3
SP - 71
EP - 86
AB - Position du problème Dans le cadre des modèles écologiques, nous comparons différents estimateurs bayésiens des risques relatifs en utilisant un modèle de cartographie hiérarchique bayésien. Méthodes Les estimateurs ponctuels bayésiens dépendent du choix de la fonction de perte. La plus « classique » est la fonction de perte quadratique et sa minimisation amène à l’estimation par moyenne a posteriori. Une discussion est menée sur les avantages de cette approche, notamment en envisageant différents objectifs comme l’intérêt particulier porté sur les valeurs extrêmes des risques relatifs lors de cartographie. Ces objectifs guident le choix de la fonction de perte. Trois fonctions sont envisagées : la fonction de perte quadratique, la fonction de perte en valeur absolue et une fonction de perte quadratique pondérée avec des poids de type exponentiel. La dernière fonction correspond à un choix moins fréquent avec une pondération importante sur les risques relatifs extrêmes. Résultats Une étude par simulations a été réalisée en considérant trois modèles de risques relatifs. Les deux premiers correspondent à des surfaces de risques discontinues alors que le dernier cas présente une décroissance du risque sur le territoire selon un gradient géographique. Les performances des estimateurs sont évaluées par l’erreur quadratique et le biais relatif ainsi que la capacité à détecter les zones à risques élevés. Les résultats montrent que de manière générale, les différents estimateurs des risques relatifs sont convergents avec de bonnes performances. Néanmoins, cette étude permet de montrer que dans le cas de risques relatifs discontinus, le choix d’une fonction de perte de type exponentiel permet d’obtenir une amélioration souvent significative de l’estimation des risques relatifs. Conclusion Le choix d’une fonction de perte et l’étude de l’estimateur qui en découle constituent une démarche rarement utilisée et pourtant prometteuse dans les approches bayésiennes. Ce choix peut être guidé par l’objectif recherché comme ici, la mise en valeur des risques relatifs extrêmes. Cette étude a montré de légères différences entre les estimateurs, prouvant dans ce sens leur robustesse mais avec une amélioration souvent significative de l’estimateur pondéré. D’autres configurations de risques relatifs et/ou de nombre de cas attendus pourraient peut-être mettre à jour des différences plus importantes.
LA - fre
KW - hierarchical bayesian model; loss functions; relative risks; ecological Poisson model
UR - http://eudml.org/doc/93466
ER -

References

top
  1. [1] Besag J., York J. & Mollie A. (1991). Bayesian image restoration, with two applications in spatial statistics. Ann. I. Stat. Math. 1991 ; 43 : 1-59. Zbl0760.62029MR1105822
  2. [2] Clavel J., Goubin A., Auclerc M.F., Auvrignon A., Waterkeyn C., Patte C. et al. (2004). Incidence of childhood leukemia and non-Hodgkin’s lymphoma in France – national registry of childhood leukemia and lymphoma, 1990-1999, Eur J Cancer Prev ; 13 : 97-103. 
  3. [3] Elliott P., Cuzick J., English D. and Stern R. (1992). Geographical and Environmental Epidemiology : Methods for Small-Area Studies, Oxford University Press, Oxford. 
  4. [4] Elliott P., Wakefield J., Best N.G. & Briggs D. (2000). Spatial Epidemiology, methods and applications. Oxford : Oxford University Press. 
  5. [5] Evrard A.S., Hemon D., Billon S., Laurier D., Jougla E., Tirmarche M. & Clavel J. (2005). Ecological association between indoor Radon concentration and childhood leukemia incidence in France, 1990-1998. Eur. J. Cancer Prev. ; 14 : 147-57. 
  6. [6] Ghosh M. (1992). Constrained Bayes estimation with applications. J. Am. Stat. Assoc. ; 87 : 533-540. Zbl0781.62040MR1173817
  7. [7] Green P.J., Richardson S. (2002). Hidden Markov models and disease mapping. J. Am. Stat. Assoc. ; 97 : 1055-1070. Zbl1046.62117MR1951259
  8. [8] Kelsall J.E. & Wakefield J. (1999). Discussion of « Bayesian models for spatially correlated disease and exposure data » by Best et al. (1999). In Bayesian Statistics 6, eds J.M. Bernardo, J.O. Berger, A.P. Dawid, A.F.M. Smith. Oxford : Oxford University Press. 
  9. [9] Lawson A., Bigerri A., Bohning D., Lesaffre E., Viel J. F., Bertollini R. (1999). Disease Mapping and risk assessment for public health. John Wiley & Sons, UK. Zbl0942.00010
  10. [10] Louis T. (1984). Estimating a population of parameter values using Bayes and empirical Bayes methods. J. Am. Stat. Assoc ; 79 : 393-398. MR755093
  11. [11] Marshall E.C. and Spiegelhalter D.J. (2003). Approximate cross-validatory predictive checks in disease mapping models. Stat. Med. ; 22 : 1649-1660. 
  12. [12] Mollié A. (1996). Bayesian mapping of disease. In : Markov Chain Monte Carlo in Practice, Gilks W.R., Richardson S., Spiegelhalter D.J. Eds., Chapman & Hall ; Chapter 20, 359-380. Zbl0849.62060MR1397966
  13. [13] Spiegelhalter D.J., Thomas A. & Best N.G. (2001). WinBUGS Version 1.4, User Manual.. Medical Research Council Biostatistics Unit, Cambridge, and Imperial College School of Medicine, London, UK. 
  14. [14] Stern H. S. et N. CressieN. (2000). Posterior predictive model checks for disease mapping models. Stat. Med. ; 19 : 2377-2397. 
  15. [15] Wright D. L., Stern H. S., and Cressie N. (2003). Loss functions for estimation of extrema with an application to disease mapping. Can. J. Stat. ; 31 : 251-266. Zbl1042.62029MR2030123

NotesEmbed ?

top

You must be logged in to post comments.

To embed these notes on your page include the following JavaScript code on your page where you want the notes to appear.

Only the controls for the widget will be shown in your chosen language. Notes will be shown in their authored language.

Tells the widget how many notes to show per page. You can cycle through additional notes using the next and previous controls.

    
                

                
Note: Best practice suggests putting the JavaScript code just before the closing </body> tag.