|
在本附录的以下章节中,我们描述了如何对第5节中描述的四个模型中的每个模型的参数进行采样。请注意,在以下部分中,我们使用Θ来描述模型参数集和符号Θ-α来描述不包括参数α的一组模型参数。D、 1。Metropolis StepLetαide的详细信息请注意参数块α的第i个样本。在Metropolis Hastings随机游走的每一步中,样本(i+1)的建议值都是从具有平均α和方差的正态分布中提取的. 在参数块α是向量的情况下,我们使用平均α和对角协方差矩阵相等的多元正态建议分布一、 Graves(2011)发现,调整MCMC算法的步长有助于在采样期间保持较高的接受概率。因此,我们调整步长 根据PyMC3框架默认实施(Salvatier et al.2016)提出的时间表,基于先前样品的验收率,如下所示:1。将步长初始化为1.2。每100个样品,执行以下操作:o如果验收率小于。05,乘法 通过5o如果验收率低于。2,乘法 通过9o如果验收率大于。5,乘法 如果验收率大于,则乘以1.1。75,乘法 如果验收率大于,则增加2。95,乘法 10之前o否则不要更改D、 2。同构模型参数第5.1节中描述的同构模型中的模型参数集为Θ={βyf,对于我们的数据中的所有y,σ,δ,ν,τ}。对于该模型,我们将超参数设置为uβ=0,σβ=100,aσ=3,bσ=1,λ=1,aν=0,bν=50,aτ=-10,且bτ=10。一般来说,我们选择的超优先级是fl-at和noninformation。对于ν和τ的先验值,我们将数据分析中的信息纳入第3节。
|