A bayesian framework for the ratio of two Poisson rates in the context of vaccine efficacy trials

Stéphane Laurent; Catherine Legrand

ESAIM: Probability and Statistics (2012)

  • Volume: 16, page 375-398
  • ISSN: 1292-8100

Abstract

top
In many applications, we assume that two random observations x and yare generated according to independent Poisson distributions ( λ S ) x1d4ab;(λS) and ( μ T ) x1d4ab;(μT) and we are interested in performing statistical inference on the ratio φ = λ / μ of the two incidence rates. In vaccine efficacy trials, x and y are typically the numbers of cases in the vaccine and the control groups respectively, φ is called the relative risk and the statistical model is called ‘partial immunity model’. In this paper we start by defining a natural semi-conjugate family of prior distributions for this model, allowing straightforward computation of the posterior inference. Following theory on reference priors, we define the reference prior for the partial immunity model when φ is the parameter of interest. We also define a family of reference priors with partial information on μ while remaining uninformative about φ. We notice that these priors belong to the semi-conjugate family. We then demonstrate using numerical examples that Bayesian credible intervals for φ enjoy attractive frequentist properties when using reference priors, a typical property of reference priors.

How to cite

top

Laurent, Stéphane, and Legrand, Catherine. "A bayesian framework for the ratio of two Poisson rates in the context of vaccine efficacy trials." ESAIM: Probability and Statistics 16 (2012): 375-398. <http://eudml.org/doc/273609>.

@article{Laurent2012,
abstract = {In many applications, we assume that two random observations x and yare generated according to independent Poisson distributions $(\lambda S)$x1d4ab;(λS) and $(\mu T)$x1d4ab;(μT) and we are interested in performing statistical inference on the ratio φ = λ / μ of the two incidence rates. In vaccine efficacy trials, x and y are typically the numbers of cases in the vaccine and the control groups respectively, φ is called the relative risk and the statistical model is called ‘partial immunity model’. In this paper we start by defining a natural semi-conjugate family of prior distributions for this model, allowing straightforward computation of the posterior inference. Following theory on reference priors, we define the reference prior for the partial immunity model when φ is the parameter of interest. We also define a family of reference priors with partial information on μ while remaining uninformative about φ. We notice that these priors belong to the semi-conjugate family. We then demonstrate using numerical examples that Bayesian credible intervals for φ enjoy attractive frequentist properties when using reference priors, a typical property of reference priors.},
author = {Laurent, Stéphane, Legrand, Catherine},
journal = {ESAIM: Probability and Statistics},
keywords = {Poisson rates; relative risk; vaccine efficacy; partial immunity model; semi-conjugate family; reference prior; Jeffreys’ prior; frequentist coverage; beta prime distribution; beta-negative binomial distribution; Jeffreys prior},
language = {eng},
pages = {375-398},
publisher = {EDP-Sciences},
title = {A bayesian framework for the ratio of two Poisson rates in the context of vaccine efficacy trials},
url = {http://eudml.org/doc/273609},
volume = {16},
year = {2012},
}

TY - JOUR
AU - Laurent, Stéphane
AU - Legrand, Catherine
TI - A bayesian framework for the ratio of two Poisson rates in the context of vaccine efficacy trials
JO - ESAIM: Probability and Statistics
PY - 2012
PB - EDP-Sciences
VL - 16
SP - 375
EP - 398
AB - In many applications, we assume that two random observations x and yare generated according to independent Poisson distributions $(\lambda S)$x1d4ab;(λS) and $(\mu T)$x1d4ab;(μT) and we are interested in performing statistical inference on the ratio φ = λ / μ of the two incidence rates. In vaccine efficacy trials, x and y are typically the numbers of cases in the vaccine and the control groups respectively, φ is called the relative risk and the statistical model is called ‘partial immunity model’. In this paper we start by defining a natural semi-conjugate family of prior distributions for this model, allowing straightforward computation of the posterior inference. Following theory on reference priors, we define the reference prior for the partial immunity model when φ is the parameter of interest. We also define a family of reference priors with partial information on μ while remaining uninformative about φ. We notice that these priors belong to the semi-conjugate family. We then demonstrate using numerical examples that Bayesian credible intervals for φ enjoy attractive frequentist properties when using reference priors, a typical property of reference priors.
LA - eng
KW - Poisson rates; relative risk; vaccine efficacy; partial immunity model; semi-conjugate family; reference prior; Jeffreys’ prior; frequentist coverage; beta prime distribution; beta-negative binomial distribution; Jeffreys prior
UR - http://eudml.org/doc/273609
ER -

References

top
  1. [1] N. Balakrishnan, N.L. Johnson and S. Kotz, Continuous Univariate Distributions, 2nd edition. John Wiley, New York 1 (1995). Zbl0821.62001
  2. [2] M.J. Bayarri and J. Berger, The interplay of Bayesian and frequentist analysis. Stat. Sci.19 (2004) 58–80. Zbl1062.62001MR2082147
  3. [3] J.O. Berger and J.M. Bernardo, Ordered Group Reference Priors With Applications to Multinomial and Variance Component Problems. Technical Report Dept. of Statistics, Purdue University (1989). 
  4. [4] J.O. Berger and J.M. Bernardo, Estimating a product of means : Bayesian analysis with reference priors. J. Amer. Statist. Assoc.84 (1989) 200–207. Zbl0682.62018MR999679
  5. [5] J.O. Berger and J.M. Bernardo, Ordered group reference priors, with applications to multinomial problems. Biometrika79 (1992) 25–37. Zbl0763.62014MR1158515
  6. [6] J.O. Berger and J.M. Bernardo, On the development of reference priors, edited by J.M. Bernardo, J.O. Berger, A.P. Dawid and A.F.M. Smith, Bayesian Statistics. University Press, Oxford (with discussion) 4 (1992) 35–60. Zbl1194.62019MR1380269
  7. [7] J.O. Berger and D. Sun, Reference priors with partial information. Biometrika85 (1998) 55–71. Zbl1067.62521MR1627242
  8. [8] J.O. Berger and R. Yang, A catalog of noninformative priors. ISDS Discussion Paper, Duke Univ. (1997) 97–42. 
  9. [9] J.O. Berger, J.M. Bernardo and D. Sun, The formal definition of reference priors. Ann. Stat. 37 (2009). Zbl1162.62013MR2502655
  10. [10] J.M. Bernardo, Reference posterior distributions for Bayesian inference (with discussion). J. R. Stat. Soc. B41 (1979) 113–148. Zbl0428.62004MR547240
  11. [11] J.M. Bernardo, Noninformative priors do not exist : a discussion. (with discussion) J. Stat. Plann. Inference 65 (1997) 159–189. MR1619672
  12. [12] J.M. Bernardo, Reference Analysis, edited by D.K. Dey and C.R. Rao. Handbook of Stat.25 (2005) 17–90. MR2490522
  13. [13] J.M. Bernardo, Intrinsic credible regions : an objective Bayesian approach to interval estimation (with discussion). Test14 (2005) 317–384. Zbl1087.62036MR2211385
  14. [14] J.M. Bernardo and J.M. Ramon, An introduction to Bayesian reference analysis : inference on the ratio of multinomial parameters. J. R. Stat. Soc. D47 (1998) 101–135. 
  15. [15] J.M. Bernardo and A.F.M. Smith, Bayesian Theory. Wiley, Chichester (1994). Zbl0943.62009MR1274699
  16. [16] D.A. Berry, M.C. Wolff and D. Sack, Decision making during a phase III randomized controlled trial. Control. Clin. Trials15 (1994) 360–378. 
  17. [17] L.D. Brown, T.T. Cai and A. DasGupta, Interval estimation for a binomial proportion (with discussion). Stat. Sci.16 (2001) 101–133. Zbl1059.62533MR1861069
  18. [18] L.D. Brown, T.T. Cai and A. DasGupta, Confidence intervals for a binomial proportion and edgeworth expansions. Ann. Stat.30 (2002) 160–201. Zbl1012.62026MR1892660
  19. [19] H. Chu and M.E. Halloran, Bayesian estimation of vaccine efficacy. Clin. Trials1 (2004) 306–314. 
  20. [20] R.D. Cousins, Improved central confidence intervals for the ratio of Poisson means. Nucl. Instrum. Methods Phys. Res. A417 (1998) 391–399. 
  21. [21] G.S. Datta and R. Mukerjee, Probability Matching Priors : Higher Order Asymptotics. Springer, New-York (2004). Zbl1044.62031MR2053794
  22. [22] M. Ewell, Comparing methods for calculating confidence intervals for vaccine efficacy. Stat. Med.15 (1996) 2379–2392. 
  23. [23] M.E. Halloran, I.M. Jr. Longini and C.J. Struchiner, Design and interpretation of vaccine field studies. Epidemiol. Rev.21 (1999) 73–88. 
  24. [24] N.L. Johnson, A.W. Kemp and S. Kotz, Univariate Discrete Distributions, 3rd edition. John Wiley, New York (2005). Zbl0773.62007MR2163227
  25. [25] R.E. Kass and L. Wasserman, The selection of prior distributions by formal rules. J. Am. Statist. Assoc.91 (1996) 1343–1370. Zbl0884.62007
  26. [26] C. Kleiber and S. Kotz, Statistical Size Distributions in Economics and Actuarial Sciences, Wiley (2003). Zbl1044.62014MR1994050
  27. [27] K. Krishnamoorthy and M. Lee, Inference for functions of parameters in discrete distributions based on fiducial approach : Binomial and Poisson cases. J. Statist. Plann. Inference140 (2009) 1182–1192. Zbl1181.62028MR2581121
  28. [28] K. Krishnamoorthy and J. Thomson, A more powerful test for comparing two Poisson means. J. Statist. Plann. Inference119 (2004) 23–35. Zbl1031.62013MR2018448
  29. [29] B. Lecoutre, And if you were a Bayesian without knowing it? Bayesian inference and maximum entropy methods in science and engineering. AIP Conf. Proc. 872 (2006) 15–22. Zbl1014.00013
  30. [30] E.L. Lehmann and J.P. Romano, Testing Statistical Hypotheses, 3rd edition. Springer, New York (2005). Zbl1076.62018MR2135927
  31. [31] B. Liseo, Elimination of Nuisance Parameters with Reference Noninformative Priors. Technical Report #90-58C, Purdue University, Department of Statistics (1990). 
  32. [32] R.M. Price and D.G. Bonett, Estimating the ratio of two Poisson rates. Comput. Stat. Data Anal.34 (2000) 345–356. Zbl1061.62523
  33. [33] C. Robert, The Bayesian Choice : From Decision-Theoretic Foundations to Computational Implementation, 2nd edition. Springer Texts in Statistics (2001). Zbl1129.62003MR1835885
  34. [34] J. Robins and L. Wasserman, Conditioning, likelihood and coherence : A review of some foundational concepts. J. Amer. Statist. Assoc.95 (2000) 1340–1346. Zbl1072.62507MR1825290
  35. [35] H. Sahai and A. Khurshid, Confidence intervals for the ratio of two Poisson means. Math. Sci.18 (1993) 43–50. Zbl0770.62022MR1227240
  36. [36] J.D. Stamey, D.M. Young, T.L. Bratcher, Bayesian sample-size determination for one and two Poisson rate parameters with applications to quality control. J. Appl. Stat.33 (2006) 583–594. Zbl1118.62315MR2240928

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.