We introduce a highly efficient fully Bayesian approach for anisotropic multidimensional smoothing. The main challenge in this context is the Markov chain Monte Carlo 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, most 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 efficient to evaluate, which allows us to set up an efficient MCMC sampler with adaptive Metropolis-Hastings updates for the smoothing parameters. We investigate different priors for the smoothing parameters and discuss the efficient derivation of lower-dimensional effects such as one-dimensional main effects and two-dimensional interactions. We show that the suggested approach outperforms previous suggestions in the literature in terms of accuracy, scalability and computational cost and demonstrate its applicability by consideration of an illustrating temperature data example from spatio-temporal statistics.
翻译:我们引入了一种高效的全巴伊西亚方法,用于厌食性多层滑动。在这方面,主要的挑战在于Markov链链Monte Carlo更新了光滑参数,因为其完全有条件的后背体由初步看来似乎难以解决的假决定物组成。因此,大多数现有实施方法在计算上都是可行的,仅用于估算二维抗拉产品平滑,但对于许多应用而言,这种光滑过于严格。在本文中,我们打破了这一屏障,并产生了对正对流-假行-确定及其第一和第二顺序部分衍生物的闭式表达式。这些表达式对任意性层面是有效的,而且非常高效地进行评估,使我们能够建立一个高效的MCMC取样器,对通融参数进行适应性Metropolis-Hasting更新。我们调查了平滑的参数的不同前科,并讨论了一维主要效应和二维互动等低维效应的有效衍生法。我们表明,所建议的方法在精确性、可缩性和计算成本方面,并展示了从数据温度统计中展示出其适用性。