PyMC3 で複数のチェーンをサンプリングしようとしています。PyMC2 では、次のようにします。
for i in range(N):
model.sample(iter=iter, burn=burn, thin = thin)
PyMC3 で同じことを行うにはどうすればよいですか? 「sample」メソッドに「njobs」引数があるのを見ましたが、値を設定するとエラーがスローされます。「pymc.gelman_rubin」出力を取得するために、サンプリングされたチェーンを使用したいと考えています。