吉布斯采样matlab代码(回收)No-U-Turn-Sampler:Matlab实现
该存储库包含Hoffman和Gelman(2014)的No-U-Turn-Sampler(NUTS)的Matlab实现以及Nishimura和Dunson(2016)的扩展Recycled
NUTS。
脚本“
getting_started_with_NUTS_and_dual_averaging_algorithm.m”说明了主要功能“
NUTS”和“
dualAveraging”的用法。
其他示例可以在“示例”文件夹下找到。
回收的NUTS实现“
ReNUTS”位于“回收”文件夹下,该功能通过回收NUTS轨迹的中间状态,提供了改进的统计效率,并且几乎没有额外的计算时间。
这里的代码适合于研究目的,因为它提供了对NUTS内部工作的访问,并且是可自定义的。
作为一个示例,此处的实现允许人们将NUTS用作Gibbs步骤。
对于希望更好地了解NUTS和HMC如何工作(以及何时可能表现不佳)的人员,该代码也应该有用。
但是,对于应用贝叶斯建模,使用Stan将是利用NUTS和HMC通用性的最简单方法。
此外
2025-06-25 14:29:39
35KB
系统开源
1