haskell - 多次重复 monadic 计算并打印出结果？ (MonadRandom)

Question

现在我正在使用 MonadRandom 库。我有一个计算:

metroChain :: (RandomGen g) => Rand g Double

我想执行多次，然后依次打印出结果。或者更确切地说，我想创建某种多重计算列表。

做一次，我会用

main = do
  result <- evalRandIO metroChain
  print result

或

main = evalRandIO metroChain >>= (\result -> print result)

但是，我在打印出任意 (n) 个 metroChain 结果时遇到了很多麻烦。

每个结果都应该使用最后一个结果末尾给出的 RandomGen...这就是 MonadRandom 应该如何工作，对吧？

我研究过 replicateM、fmap 和一些变形金刚(尽管我承认我似乎无法充分理解它们以掌握它们在我的应用程序中的应用问题)。

任何人都可以帮助我实现我正在寻找的功能吗？我觉得我错过了一些非常简单的东西。但我对 Haskell 还很陌生。

Answer 1

我要进行一次飞跃，并假设 metroStep 是 MCMC Metropolis-Hastings 迭代。

您遇到的问题是您希望 MH 步骤是马尔可夫的，但仅共享 RandomGen 状态(这正是 replicateM n metroStep 所做的)是不够的。这只会使每个步骤都能够基于独立的随机变量。相比之下，如果未共享 RandomGen 状态，则不变性将保证每个 metroStep 都是相同的。

因此，您真正需要的是同时具有 RandomGen 状态的东西，以便提供用于生成自变量样本的伪随机数链 和 固定状态，以便在每一步你都可以有 P(x_i | theta, x_(i-1))。我们构建了一个转换器堆栈来执行此操作——我将使用 mtl 库和 random-fu，因为几天前我刚刚使用这些库编写了一个 MCMC。

metroStep :: (MonadRandom m, MonadState StateSpace m) => m StateSpace

其中 StateSpace 是状态空间中的一个点，包括观察到的和未观察到的变量——它是似然函数右侧的每个参数。现在，replicateM n metroStep::(MonadRandom m, MonadState StateSpace m) => m [StateSpace] 是马尔可夫序列 StateSpace 点的列表。

然后我们像这样“运行”这个 monad 堆栈的具体版本

do steps <- (`runRVar` StdRandom) . (`evalStateT` ss0) $ (replicateM n metroStep)
   mapM_ print steps