马尔可夫链蒙特卡罗（MCMC）方法是一类基于构造马尔可夫链的概率分布抽样算法，该马尔可夫链具有所需的分布作为其均衡分布。然后将多个步骤之后的链状态用作所需分布的样品
我在 MCMCglmm 中运行混合模型时遇到一些麻烦。当我尝试将其作为线性回归运行时，它工作得很好： Model.Fit = MCMCglmm(Y~X1, 随机 = ~Cluster, 数据 = 数据集) 但我不能...
无法让 DiscreteMarkovChain 在 PyMC 中工作
我正在尝试使用 PyMC 拟合两国政权切换模型。 完整模型如下所示。 模型的具体细节并不重要（如果你还是好奇的话，我正在尝试拟合地理......
如何在Mathematica中写出Profile分布（Profile Likelihood）？
我是 Mathematica 新手。我想使用从 MCMC 链获得的结果并使用此方法 2203.16285 来获取配置文件分布。我什至找到了作者演示文稿中的步骤...
有没有办法在一组值上并行化 scipy.integrate.quad？
我正在使用计算速度较慢的似然函数计算某些参数 θ 的后验概率，我想知道是否有办法加快速度。 我的代码中最慢的部分是
有没有办法在整组值上并行化 scipy.integrate 函数？
我正在计算某些参数 $ heta$ 的后验概率，其似然函数计算速度很慢，我想知道是否有办法加快速度。 我的代码中缓慢的部分...
我正在尝试编写一个预订模型，但遇到了一些问题： 当我按“加载初始化”时，我收到消息 该链包含统一变量 虽然我已经初始化了所有的
如何在使用自定义对数似然的 PyMC V5 模型上应用 Blackjax 采样器
我正在使用 PyMC v5 在模型中执行哈密顿蒙特卡罗。我可以运行下面的代码，但即使有多个核心，它也非常慢。我有一个 applyMCMC 函数用于此目的： # 定义一个
blavaan 错误：从 lavaan 到 MCMC 语法的翻译出现问题
我正在使用 blavaan 运行测量不变性测试，但出现下一个错误： [1] “lamsign[l1, 1] 中的错误：下标超出范围 ” attr(,"类") [1] &
我正在研究 SAS 中 MCMC 的一种实现 https://support.sas.com/documentation/cdl/en/statug/67523/HTML/default/viewer.htm#statug_mcmc_examples13.htm 它实现了 proc mcmc 数据=种子
我有一个有很多参数/变量的函数。为了简单起见，我们可以假设它是 f(a,b,c)=a + b * x + c * x**2，其中参数“a”、“b”和“c”可以...
无法使用 Pymc 3.11.4 加载之前保存的跟踪以供将来推理
我正在尝试实现一个简单的分层贝叶斯推理。我使用 Pymc 推断简单线性模型 y = mx 的参数“m”。我想保存跟踪以便稍后加载。
R - 仅在一台设备上出现错误：chol.default() 中的错误：顺序的前导次要不是正定的，BVAR
我在 R 中估计 BVAR 模型并得到以下错误： chol.default(V_post) 中的错误：1 阶的前导次要不是正定的 这实际上不是问题，因为我工作......
我想在pymc2中为iris数据集构建BNN模型。 我定义了模型并尝试进行训练，但训练和测试数据的准确度仅为 0.33。 这是我当前的代码 迭代次数 = 2000 伊里...
arviz.from_pyjags() ValueError：需要解包的值太多（预计为 3 个）
当我使用 az.from_pyjags 时，它引发了错误。 -------------------------------------------------- ------------------------ ValueError Traceback（最近一次调用最后一次）...
我想知道是否有办法在 Google Colab 应用程序上并行化 Jupyter 笔记本。我使用 cobaya 包进行宇宙学分析，并且执行了许多蒙特卡洛马尔可夫 Ch...
MCMCglmm 具有生态因素和系统发育树：“某些级别的动物在 ginverse 中没有行条目”
我正在尝试针对数量性状和生态因素运行 MCMCglmm，同时考虑系统发育。我正在研究猫头鹰物种的头骨不对称性，所以我使用系统发育树
R 中的 JAGS：生成的 MCMC 链不会在初始值附近开始
我正在运行这个简单的 Rjags 教程中的代码。请注意，参数的初始值由 model.inits 指定<- list(tau=1, beta0=1, beta1=0, beta2=0). I wanted to test tha...
基于“R2jags”和“bridgesampling”验证贝叶斯因子计算
我正在尝试通过桥接采样计算单向方差分析的 Jeffreys-Zellner-Siow (JZS) 贝叶斯因子。然而，结果远不如 anovaBF 函数计算的值可靠...
我正在尝试构建一个吉布斯抽样算法来估计线性回归的系数和方差。由于假设方差是常数，我尝试使用 MVNormal-conjugate
使用 MCMCpack 作为潜在变量并从 mcmc 对象中提取潜在分数
我正在使用 R 包 MCMCpack 中的函数 MCMCordfactanal() 和 MCMCirtKd() 生成后验分布并创建潜在因子变量。 我可以使用 summary() 和