Estimating boundary curves has many applications such as economics, climate science, and medicine. Bayesian trend filtering has been developed as one of locally adaptive smoothing methods to estimate the non-stationary trend of data. This paper develops a Bayesian trend filtering for estimating boundary trend. To this end, the truncated multivariate normal working likelihood and global-local shrinkage priors based on scale mixtures of normal distribution are introduced. In particular, well-known horseshoe prior for difference leads to locally adaptive shrinkage estimation for boundary trend. However, the full conditional distributions of the Gibbs sampler involve high-dimensional truncated multivariate normal distribution. To overcome the difficulty of sampling, an approximation of truncated multivariate normal distribution is employed. Using the approximation, the proposed models lead to an efficient Gibbs sampling algorithm via P\'olya-Gamma data augmentation. The proposed method is also extended by considering nearly isotonic constraint. The performance of the proposed method is illustrated through some numerical experiments and real data examples.
翻译:边界曲线估计在经济学、气候科学和医学等领域具有广泛应用。贝叶斯趋势滤波作为一种局部自适应平滑方法,被用于估计数据的非平稳趋势。本文提出了一种用于估计边界趋势的贝叶斯趋势滤波方法。为此,引入了基于截断多元正态分布的工作似然函数,以及基于正态分布尺度混合的全局-局部收缩先验。特别是,针对差分项采用著名的马蹄先验,可实现边界趋势的局部自适应收缩估计。然而,吉布斯采样器的全条件分布涉及高维截断多元正态分布。为克服采样困难,采用截断多元正态分布的近似方法。借助该近似,所提模型通过Pólya-Gamma数据增广实现了高效的吉布斯采样算法。此外,该方法进一步扩展至近似等序约束情形。通过数值实验和真实数据实例验证了所提方法的有效性。