×

Bayesian sequential design for copula models. (English) Zbl 1447.62091

Summary: Bayesian design requires determining the value of controllable variables in an experiment to maximise the information that will be obtained for subsequently collected data, with the majority of research in this field being focused on experiments that yield a univariate response. In this paper, a robust and computationally efficient Bayesian design approach is proposed to derive designs for experiments which yield bivariate discrete and mixed responses. To construct the joint distribution of responses, Copula models are considered, and a sequential Monte Carlo algorithm is adopted to reduce the computational effort required in deriving sequential designs. The total entropy utility function is considered to derive designs for the dual experimental goals of parameter estimation and model discrimination for Copula models. The results show that designs constructed within our framework are able to precisely estimate model parameters and that it is possible to discriminate between different competing Copula models. However, for experiments which yield binary and continuous data, it appears as though discriminating between Copula models can require a large number of data points, which may limit the general applicability of our methods and/or the range of experimental objectives that can be considered in experiments that yield multiple responses.

MSC:

62K05 Optimal statistical designs
62L05 Sequential statistical design
62L12 Sequential estimation
62H05 Characterization and structure theory for multivariate probability distributions; copulas

Software:

acebayes; CDVine
PDFBibTeX XMLCite
Full Text: DOI Link

References:

[1] Atkinson, AC, Robust optimum designs for transformation of the responses in a multivariate chemical kinetic model, Technometrics, 47, 4, 478-487 (2005)
[2] Borth, DM, A total entropy criterion for the dual problem of model discrimination and parameter estimation, J R Stat Soc Ser B (Methodol), 37, 1, 77-87 (1975) · Zbl 0297.62061
[3] Brechmann, EC; Schepsmeier, U., Modeling dependence with C-and D-vine Copulas: the R-package CDVine, J Stat Softw, 52, 3, 1-27 (2013)
[4] Carlin, BP; Kadane, JB; Gelfand, AE, Approaches for optimal sequential decision analysis in clinical trials, Biometrics, 54, 3, 964-975 (1998) · Zbl 1058.62548
[5] Chaloner, K.; Verdinelli, I., Bayesian experimental design: a review, Stat Sci, 10, 3, 273-304 (1995) · Zbl 0955.62617
[6] Clyde, M.; Chaloner, K., The equivalence of constrained and weighted designs in multiple objective design problems, J Am Stat Assoc, 91, 435, 1236-1244 (1996) · Zbl 0883.62079
[7] Craiu, RV; Rosenthal, JS, Bayesian computation via Markov chain Monte Carlo, Ann Rev Stat Appl, 1, 179-201 (2014)
[8] Del Moral, P.; Doucet, A.; Jasra, A., Sequential Monte Carlo samplers, J R Stat Soc Ser B (Stat Methodol), 68, 3, 411-436 (2006) · Zbl 1105.62034
[9] Deldossi L, Osmetti SA, Tommasi C (2018) Optimal design to discriminate between rival copula models for a bivariate binary response. TEST · Zbl 1420.62327
[10] Denman, NG; McGree, JM; Eccleston, JA; Duffull, SB, Design of experiments for bivariate binary responses modelled by Copula functions, Comput Stat Data Anal, 55, 4, 1509-1520 (2011) · Zbl 1328.62488
[11] Drovandi, CC; McGree, JM; Pettitt, AN, Sequential Monte Carlo for Bayesian sequentially designed experiments for discrete data, Comput Stat Data Anal, 57, 1, 320-335 (2013) · Zbl 1365.62318
[12] Drovandi, CC; McGree, JM; Pettitt, AN, A sequential Monte Carlo algorithm to incorporate model uncertainty in Bayesian sequential design, J Comput Graph Stat, 23, 1, 3-24 (2014)
[13] Durante F, Sempi C (2010) Copula theory: an introduction. In: Copula theory and its applications. Springer, New York, pp 3-31
[14] Genest, C.; MacKay, J., The joy of Copulas: bivariate distributions with uniform marginals, Am Stat, 40, 4, 280-283 (1986)
[15] Gueorguieva, I.; Aarons, L.; Ogungbenro, K.; Jorga, KM; Rodgers, T.; Rowland, M., Optimal design for multivariate response pharmacokinetic models, J Pharmacokinet Pharmacodyn, 33, 2, 97 (2006)
[16] Hill, WJ; Hunter, WG; Wichern, DW, A joint design criterion for the dual problem of model discrimination and parameter estimation, Technometrics, 10, 1, 145-160 (1968)
[17] Kullback, S.; Leibler, RA, On information and sufficiency, Ann Math Stat, 22, 1, 79-86 (1951) · Zbl 0042.38403
[18] Liu JS (2008) Monte Carlo strategies in scientific computing. Springer, New York (Incorporated) · Zbl 1132.65003
[19] Lopez-Fidalgo, J.; Tommasi, C.; Trandafir, PC, An optimal experimental design criterion for discriminating between non-normal models, J R Stat Soc Ser B (Stat Methodol), 69, 2, 231-242 (2007) · Zbl 1123.62056
[20] McGree, JM, Developments of the total entropy utility function for the dual purpose of model discrimination and parameter estimation in Bayesian design, Comput Stat Data Anal, 113, 207-225 (2017) · Zbl 1464.62130
[21] McGree, JM; Drovandi, CC; Thompson, M.; Eccleston, J.; Duffull, S.; Mengersen, K.; Pettitt, AN; Goggin, T., Adaptive Bayesian compound designs for dose finding studies, J Stat Plan Inference, 142, 6, 1480-1492 (2012) · Zbl 1242.62120
[22] McGree, JM; Drovandi, CC; White, G.; Pettitt, AN, A pseudo-marginal sequential Monte Carlo algorithm for random effects models in Bayesian sequential design, Stat Comput, 26, 5, 1121-1136 (2016) · Zbl 1505.62284
[23] McGree, JM; Eccleston, JA, Probability-based optimal design, Aust N Z J Stat, 50, 1, 13-28 (2008)
[24] Montgomery, DC, Design and analysis of experiments (2006), New York: Wiley, New York
[25] Nelsen, RB, An introduction to Copulas (2006), New York: Springer Series in Statistics, Springer-Verlag, New York · Zbl 1152.62030
[26] Overstall, AM; McGree, JM; Drovandi, CC, An approach for finding fully Bayesian optimal designs using normal-based approximations to loss functions, Stat Comput, 28, 2, 343-358 (2018) · Zbl 1384.62277
[27] Overstall, AM; Woods, DC, Bayesian design of experiments using approximate coordinate exchange, Technometrics, 59, 4, 458-470 (2017)
[28] Overstall AM, Woods DC, Adamou M (2018) acebayes: an R package for Bayesian optimal design of experiments via approximate coordinate exchange. arXiv:1705.08096
[29] Panagiotelis, A.; Czado, C.; Joe, H.; Stöber, J., Model selection for discrete regular vine Copulas, Comput Stat Data Anal, 106, 138-152 (2017) · Zbl 1466.62171
[30] Perrone, E.; Müller, WG, Optimal designs for Copula models, Statistics, 50, 4, 917-929 (2016) · Zbl 1353.62085
[31] Perrone E, Rappold A, Müller WG \((2016) {D}_s\)-optimality in Copula models. Stat Methods Appl 1-16 · Zbl 1384.62278
[32] Schoelzel, C.; Friederichs, P., Multivariate non-normally distributed random variables in climate research-introduction to the Copula approach, Nonlinear Process Geophys, 15, 5, 761-772 (2008)
[33] Sklar, A., Fonctions de répartition à n dimensions et leurs marges, Publications de l’Institut de Statistique de L’Université de Paris, 8, 229-231 (1959) · Zbl 0100.14202
[34] Song, PX-K; Li, M.; Yuan, Y., Joint regression analysis of correlated data using gaussian Copulas, Biometrics, 65, 1, 60-68 (2009) · Zbl 1159.62049
[35] Tao A (2010) A dose-finding method in joint modeling of efficacy and safety endpoints in phase II studies: an extension of the MCP-Mod method, PhD thesis, University of Medicine and Dentistry of New Jersey
[36] Tao, Y.; Liu, J.; Li, Z.; Lin, J.; Lu, T.; Yan, F., Dose-finding based on bivariate efficacy-toxicity outcome using archimedean Copula, PLoS ONE, 8, 11, 1-6 (2013)
[37] Tierney, L., Markov chains for exploring posterior distributions, Ann Stat, 22, 4, 1701-1728 (1994) · Zbl 0829.62080
[38] Waterhouse, T.; Woods, D.; Eccleston, J.; Lewis, S., Design selection criteria for discrimination/estimation for nested models and a binomial response, J Stat Plan Inference, 138, 1, 132-144 (2008) · Zbl 1144.62060
[39] Weir, CJ; Spiegelhalter, DJ; Grieve, AP, Flexible design and efficient implementation of adaptive dose-finding studies, J Biopharm Stat, 17, 6, 1033-1050 (2007)
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.