Input: 𝑌 , 𝑛 m c , 𝜔 i n i t
Output:   𝜔
𝜔 = 𝜔 𝜔 i n i t
for   𝑛 = 1 to 𝑛 m c   do
    propose 𝜔 based on 𝜔 (see [9])
    sample 𝑈 from Unif[ 0 , 1 ]
    𝜔 𝜔 if 𝑈 < 𝐴 ( 𝜔 , 𝜔 )
      𝜔 𝜔 if 𝑝 ( 𝜔 𝑌 ) / 𝑝 ( 𝜔 𝑌 ) > 1
end for
Algorithm 1: MCMCDA.