{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2025,2,26]],"date-time":"2025-02-26T16:40:06Z","timestamp":1740588006660,"version":"3.38.0"},"reference-count":32,"publisher":"Springer Science and Business Media LLC","issue":"2","license":[{"start":{"date-parts":[[2025,2,5]],"date-time":"2025-02-05T00:00:00Z","timestamp":1738713600000},"content-version":"tdm","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0"},{"start":{"date-parts":[[2025,2,5]],"date-time":"2025-02-05T00:00:00Z","timestamp":1738713600000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0"}],"funder":[{"DOI":"10.13039\/501100001659","name":"Deutsche Forschungsgemeinschaft","doi-asserted-by":"publisher","award":["KL 3037\/1-1","KL 3037\/1-1"],"award-info":[{"award-number":["KL 3037\/1-1","KL 3037\/1-1"]}],"id":[{"id":"10.13039\/501100001659","id-type":"DOI","asserted-by":"publisher"}]},{"name":"Karlsruher Institut f\u00fcr Technologie (KIT)"}],"content-domain":{"domain":["link.springer.com"],"crossmark-restriction":false},"short-container-title":["Stat Comput"],"published-print":{"date-parts":[[2025,4]]},"abstract":"<jats:title>Abstract<\/jats:title>\n          <jats:p>We introduce a highly efficient fully Bayesian approach for anisotropic multidimensional smoothing. The main challenge in this context is the Markov chain Monte Carlo (MCMC) update of the smoothing parameters as their full conditional posterior comprises a pseudo-determinant that appears to be intractable at first sight. As a consequence, existing implementations are computationally feasible only for the estimation of two-dimensional tensor product smooths, which is, however, too restrictive for many applications. In this paper, we break this barrier and derive closed-form expressions for the log-pseudo-determinant and its first and second order partial derivatives. These expressions are valid for arbitrary dimension and very fast to evaluate, which allows us to set up an efficient MCMC sampler with derivative-based Metropolis\u2013Hastings (MH) updates for the smoothing parameters. We derive simple formulas for low-dimensional slices and averages to facilitate visualization and investigate hyperprior sensitivity. We show that our new approach outperforms previous suggestions in the literature in terms of accuracy, scalability and computational cost and demonstrate its applicability through an illustrating temperature data example from spatio-temporal statistics.\n<\/jats:p>","DOI":"10.1007\/s11222-025-10569-y","type":"journal-article","created":{"date-parts":[[2025,2,5]],"date-time":"2025-02-05T16:43:02Z","timestamp":1738773782000},"update-policy":"https:\/\/doi.org\/10.1007\/springer_crossmark_policy","source":"Crossref","is-referenced-by-count":0,"title":["Anisotropic multidimensional smoothing using Bayesian tensor product P-splines"],"prefix":"10.1007","volume":"35","author":[{"given":"Paul","family":"Bach","sequence":"first","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Nadja","family":"Klein","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"297","published-online":{"date-parts":[[2025,2,5]]},"reference":[{"key":"10569_CR1","unstructured":"Belitz, C., Brezger, A., Klein, N., Kneib, T., Lang, S., Umlauf, N.: BayesX\u2014Software for Bayesian inference in structured additive regression models. http:\/\/www.bayesx.org. Version 3.0.2 (2015)"},{"key":"10569_CR2","doi-asserted-by":"crossref","unstructured":"B\u00fcrkner, P.-C., Gabry, J., Kay, M., Vehtari, A.: posterior: Tools for Working with Posterior Distributions. R package version 1.3.1 (2022)","DOI":"10.32614\/CRAN.package.posterior"},{"issue":"1","key":"10569_CR3","doi-asserted-by":"publisher","first-page":"1","DOI":"10.18637\/jss.v076.i01","volume":"76","author":"B Carpenter","year":"2017","unstructured":"Carpenter, B., Gelman, A., Hoffman, M.D., Lee, D., Goodrich, B., Betancourt, M., Brubaker, M., Guo, J., Li, P., Riddell, A.: Stan: A probabilistic programming language. J. Stat. Softw. 76(1), 1\u201332 (2017)","journal-title":"J. Stat. Softw."},{"issue":"2","key":"10569_CR4","doi-asserted-by":"publisher","first-page":"159","DOI":"10.1016\/S0169-7439(03)00029-7","volume":"66","author":"PHC Eilers","year":"2003","unstructured":"Eilers, P.H.C., Marx, B.D.: Multivariate calibration with temperature interaction using two-dimensional penalized signal regression. Chemom. Intell. Lab. Syst. 66(2), 159\u2013174 (2003)","journal-title":"Chemom. Intell. Lab. Syst."},{"issue":"1","key":"10569_CR5","first-page":"1","volume":"19","author":"JH Friedman","year":"1991","unstructured":"Friedman, J.H.: Multivariate adaptive regression splines. Ann. Stat. 19(1), 1\u201367 (1991)","journal-title":"Ann. Stat."},{"issue":"4","key":"10569_CR6","doi-asserted-by":"publisher","first-page":"775","DOI":"10.1081\/STA-120018828","volume":"32","author":"J Geweke","year":"2003","unstructured":"Geweke, J., Tanizaki, H.: Note on the sampling distribution for the Metropolis-Hastings algorithm. Commun. Stat. Theory Methods 32(4), 775\u2013789 (2003)","journal-title":"Commun. Stat. Theory Methods"},{"issue":"2","key":"10569_CR7","doi-asserted-by":"publisher","first-page":"359","DOI":"10.1214\/17-BA1051","volume":"13","author":"J Ghosh","year":"2018","unstructured":"Ghosh, J., Li, Y., Mitra, R.: On the use of Cauchy prior distributions for Bayesian logistic regression. Bayesian Anal. 13(2), 359\u2013383 (2018)","journal-title":"Bayesian Anal."},{"key":"10569_CR8","unstructured":"Goodrich, B., Gabry, J., Ali, I., Brilleman, S.: rstanarm: Bayesian applied regression modeling via Stan. R package version 2.21.3. (2022)"},{"key":"10569_CR9","doi-asserted-by":"publisher","DOI":"10.1007\/978-1-4614-5369-7","volume-title":"Smoothing Spline ANOVA Models","author":"C Gu","year":"2013","unstructured":"Gu, C.: Smoothing Spline ANOVA Models, 2nd edn. Springer, New York, NY (2013)","edition":"2"},{"key":"10569_CR10","doi-asserted-by":"publisher","DOI":"10.1007\/978-1-4939-8853-2","volume-title":"Semiparametric Regression with R","author":"J Harezlak","year":"2018","unstructured":"Harezlak, J., Ruppert, D., Wand, M.P.: Semiparametric Regression with R. Springer, New York, NY (2018)"},{"issue":"1","key":"10569_CR11","doi-asserted-by":"publisher","first-page":"124","DOI":"10.1214\/19-STS755","volume":"36","author":"M Katzfuss","year":"2021","unstructured":"Katzfuss, M., Guinness, J.: A general framework for Vecchia approximations of Gaussian processes. Stat. Sci. 36(1), 124\u2013141 (2021)","journal-title":"Stat. Sci."},{"issue":"4","key":"10569_CR12","doi-asserted-by":"publisher","first-page":"1071","DOI":"10.1214\/15-BA983","volume":"11","author":"N Klein","year":"2016","unstructured":"Klein, N., Kneib, T.: Scale-dependent priors for variance parameters in structured additive distributional regression. Bayesian Anal. 11(4), 1071\u20131106 (2016)","journal-title":"Bayesian Anal."},{"issue":"509","key":"10569_CR13","doi-asserted-by":"publisher","first-page":"405","DOI":"10.1080\/01621459.2014.912955","volume":"110","author":"N Klein","year":"2015","unstructured":"Klein, N., Kneib, T., Lang, S.: Bayesian generalized additive models for location, scale, and shape for zero-inflated and overdispersed count data. J. Am. Stat. Assoc. 110(509), 405\u2013419 (2015)","journal-title":"J. Am. Stat. Assoc."},{"issue":"1","key":"10569_CR14","doi-asserted-by":"publisher","first-page":"1","DOI":"10.1007\/s11749-019-00631-z","volume":"28","author":"T Kneib","year":"2019","unstructured":"Kneib, T., Klein, N., Lang, S., Umlauf, N.: Modular regression\u2014a Lego system for building structured additive distributional regression models with tensor product interactions. TEST 28(1), 1\u201339 (2019)","journal-title":"TEST"},{"issue":"30","key":"10569_CR15","doi-asserted-by":"publisher","first-page":"4771","DOI":"10.1002\/sim.7967","volume":"37","author":"M K\u00f6hler","year":"2018","unstructured":"K\u00f6hler, M., Umlauf, N., Greven, S.: Nonlinear association structures in flexible Bayesian additive joint models. Stat. Med. 37(30), 4771\u20134788 (2018)","journal-title":"Stat. Med."},{"issue":"1","key":"10569_CR16","doi-asserted-by":"publisher","first-page":"183","DOI":"10.1198\/1061860043010","volume":"13","author":"S Lang","year":"2004","unstructured":"Lang, S., Brezger, A.: Bayesian P-splines. J. Comput. Graph. Stat. 13(1), 183\u2013212 (2004)","journal-title":"J. Comput. Graph. Stat."},{"issue":"2","key":"10569_CR17","doi-asserted-by":"publisher","first-page":"223","DOI":"10.1007\/s11222-012-9366-0","volume":"24","author":"S Lang","year":"2014","unstructured":"Lang, S., Umlauf, N., Wechselberger, P., Harttgen, K., Kneib, T.: Multilevel structured additive regression. Stat. Comput. 24(2), 223\u2013238 (2014)","journal-title":"Stat. Comput."},{"issue":"1","key":"10569_CR18","doi-asserted-by":"publisher","first-page":"49","DOI":"10.1177\/1471082X1001100104","volume":"11","author":"D-J Lee","year":"2011","unstructured":"Lee, D.-J., Durb\u00e1n, M.: P-spline ANOVA-type interaction models for spatio-temporal smoothing. Stat. Model. 11(1), 49\u201369 (2011)","journal-title":"Stat. Model."},{"issue":"3","key":"10569_CR19","doi-asserted-by":"publisher","first-page":"705","DOI":"10.1214\/aos\/1056562461","volume":"31","author":"RM Neal","year":"2003","unstructured":"Neal, R.M.: Slice sampling. Ann. Stat. 31(3), 705\u2013767 (2003)","journal-title":"Ann. Stat."},{"key":"10569_CR20","volume-title":"Numerical Optimization","author":"J Nocedal","year":"2006","unstructured":"Nocedal, J., Wright, S.: Numerical Optimization, 2nd edn. Springer, New York, NY (2006)","edition":"2"},{"key":"10569_CR21","unstructured":"Plummer, M.: JAGS: a program for analysis of Bayesian graphical models using Gibbs sampling. In: Hornik, K., Leisch, F., Zeileis, A. (eds.) Proceedings of the 3rd International Workshop on Distributed Statistical Computing. Technische Universit\u00e4t Wien, Vienna (2003)"},{"key":"10569_CR22","volume-title":"Gaussian Processes for Machine Learning","author":"CE Rasmussen","year":"2006","unstructured":"Rasmussen, C.E., Williams, C.K.I.: Gaussian Processes for Machine Learning. MIT Press, Cambridge, MA (2006)"},{"key":"10569_CR23","doi-asserted-by":"crossref","unstructured":"Rodr\u00edguez-\u00c1lvarez, M.X., Lee, D.-J., Kneib, T., Durb\u00e1n, M., Eilers, P.: Fast smoothing parameter separation in multidimensional generalized P-splines: the SAP algorithm. Stat. Comput. 25(5), 941\u2013957 (2015)","DOI":"10.1007\/s11222-014-9464-2"},{"key":"10569_CR24","doi-asserted-by":"publisher","DOI":"10.1201\/9780203492024","volume-title":"Gaussian Markov Random Fields: Theory and Applications","author":"H Rue","year":"2005","unstructured":"Rue, H., Held, L.: Gaussian Markov Random Fields: Theory and Applications. Chapman & Hall\/CRC Press, Boca Raton, FL (2005)"},{"issue":"3","key":"10569_CR25","doi-asserted-by":"publisher","first-page":"612","DOI":"10.1080\/10618600.2017.1407325","volume":"27","author":"N Umlauf","year":"2018","unstructured":"Umlauf, N., Klein, N., Zeileis, A.: bamlss: Bayesian additive models for location, scale, and shape (and beyond). J. Comput. Graph. Stat. 27(3), 612\u2013627 (2018)","journal-title":"J. Comput. Graph. Stat."},{"issue":"2","key":"10569_CR26","doi-asserted-by":"publisher","first-page":"667","DOI":"10.1214\/20-BA1221","volume":"16","author":"A Vehtari","year":"2021","unstructured":"Vehtari, A., Gelman, A., Simpson, D., Carpenter, B., B\u00fcrkner, P.-C.: Rank-normalization, folding, and localization: an improved R for assessing convergence of MCMC (with discussion). Bayesian Anal. 16(2), 667\u2013718 (2021)","journal-title":"Bayesian Anal."},{"issue":"4","key":"10569_CR27","doi-asserted-by":"publisher","first-page":"1025","DOI":"10.1111\/j.1541-0420.2006.00574.x","volume":"62","author":"SN Wood","year":"2006","unstructured":"Wood, S.N.: Low-rank scale-invariant tensor product smooths for generalized additive mixed models. Biometrics 62(4), 1025\u20131036 (2006)","journal-title":"Biometrics"},{"issue":"1","key":"10569_CR28","doi-asserted-by":"publisher","first-page":"3","DOI":"10.1111\/j.1467-9868.2010.00749.x","volume":"73","author":"SN Wood","year":"2011","unstructured":"Wood, S.N.: Fast stable restricted maximum likelihood and marginal likelihood estimation of semiparametric generalized linear models. J. R. Stat. Soc. Ser. B (Stat. Methodol.) 73(1), 3\u201336 (2011)","journal-title":"J. R. Stat. Soc. Ser. B (Stat. Methodol.)"},{"key":"10569_CR29","unstructured":"Wood, S.N.: mgcv: Mixed GAM Computation Vehicle with GCV\/AIC\/REML smoothness estimation (2012)"},{"issue":"7","key":"10569_CR30","doi-asserted-by":"publisher","first-page":"1","DOI":"10.18637\/jss.v075.i07","volume":"75","author":"SN Wood","year":"2016","unstructured":"Wood, S.N.: Just another Gibbs additive modeler: interfacing JAGS and mgcv. J. Stat. Softw. 75(7), 1\u201315 (2016)","journal-title":"J. Stat. Softw."},{"issue":"4","key":"10569_CR31","doi-asserted-by":"publisher","first-page":"1071","DOI":"10.1111\/biom.12666","volume":"73","author":"SN Wood","year":"2017","unstructured":"Wood, S.N., Fasiolo, M.: A generalized Fellner-Schall method for smoothing parameter optimization with application to Tweedie location, scale and shape models. Biometrics 73(4), 1071\u20131081 (2017)","journal-title":"Biometrics"},{"issue":"3","key":"10569_CR32","doi-asserted-by":"publisher","first-page":"341","DOI":"10.1007\/s11222-012-9314-z","volume":"23","author":"SN Wood","year":"2013","unstructured":"Wood, S.N., Scheipl, F., Faraway, J.J.: Straightforward intermediate rank tensor product smoothing in mixed models. Stat. Comput. 23(3), 341\u2013360 (2013)","journal-title":"Stat. Comput."}],"container-title":["Statistics and Computing"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/link.springer.com\/content\/pdf\/10.1007\/s11222-025-10569-y.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"text-mining"},{"URL":"https:\/\/link.springer.com\/article\/10.1007\/s11222-025-10569-y\/fulltext.html","content-type":"text\/html","content-version":"vor","intended-application":"text-mining"},{"URL":"https:\/\/link.springer.com\/content\/pdf\/10.1007\/s11222-025-10569-y.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2025,2,26]],"date-time":"2025-02-26T16:15:20Z","timestamp":1740586520000},"score":1,"resource":{"primary":{"URL":"https:\/\/link.springer.com\/10.1007\/s11222-025-10569-y"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2025,2,5]]},"references-count":32,"journal-issue":{"issue":"2","published-print":{"date-parts":[[2025,4]]}},"alternative-id":["10569"],"URL":"https:\/\/doi.org\/10.1007\/s11222-025-10569-y","relation":{},"ISSN":["0960-3174","1573-1375"],"issn-type":[{"type":"print","value":"0960-3174"},{"type":"electronic","value":"1573-1375"}],"subject":[],"published":{"date-parts":[[2025,2,5]]},"assertion":[{"value":"7 April 2024","order":1,"name":"received","label":"Received","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"13 January 2025","order":2,"name":"accepted","label":"Accepted","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"5 February 2025","order":3,"name":"first_online","label":"First Online","group":{"name":"ArticleHistory","label":"Article History"}}],"article-number":"43"}}