mcmc 相关问题

马尔可夫链蒙特卡罗(MCMC)方法是一类基于构造马尔可夫链的概率分布抽样算法,该马尔可夫链具有所需的分布作为其均衡分布。然后将多个步骤之后的链状态用作所需分布的样品

如何在Mathematica中写出Profile分布(Profile Likelihood)?

我是 Mathematica 新手。我想使用从 MCMC 链获得的结果并使用此方法 2203.16285 来获取配置文件分布。我什至找到了作者演示文稿中的步骤...

回答 1 投票 0

有没有办法在一组值上并行化 scipy.integrate.quad?

我正在使用计算速度较慢的似然函数计算某些参数 θ 的后验概率,我想知道是否有办法加快速度。 我的代码中最慢的部分是

回答 1 投票 0

有没有办法在整组值上并行化 scipy.integrate 函数?

我正在计算某些参数 $ heta$ 的后验概率,其似然函数计算速度很慢,我想知道是否有办法加快速度。 我的代码中缓慢的部分...

回答 1 投票 0

gamma 全条件 tau 的形状参数太小,WinBUGS

我正在尝试编写一个预订模型,但遇到了一些问题: 当我按“加载初始化”时,我收到消息 该链包含统一变量 虽然我已经初始化了所有的

回答 1 投票 0

如何在使用自定义对数似然的 PyMC V5 模型上应用 Blackjax 采样器

我正在使用 PyMC v5 在模型中执行哈密顿蒙特卡罗。我可以运行下面的代码,但即使有多个核心,它也非常慢。我有一个 applyMCMC 函数用于此目的: # 定义一个

回答 1 投票 0

blavaan 错误:从 lavaan 到 MCMC 语法的翻译出现问题

我正在使用 blavaan 运行测量不变性测试,但出现下一个错误: [1] “lamsign[l1, 1] 中的错误:下标超出范围 ” attr(,"类") [1] &

回答 1 投票 0

了解mcmc实现过程的一些部分

我正在研究 SAS 中 MCMC 的一种实现 https://support.sas.com/documentation/cdl/en/statug/67523/HTML/default/viewer.htm#statug_mcmc_examples13.htm 它实现了 proc mcmc 数据=种子

回答 1 投票 0

如何仅使用函数变量的子集进行 MCMC 拟合?

我有一个有很多参数/变量的函数。为了简单起见,我们可以假设它是 f(a,b,c)=a + b * x + c * x**2,其中参数“a”、“b”和“c”可以...

回答 1 投票 0

无法使用 Pymc 3.11.4 加载之前保存的跟踪以供将来推理

我正在尝试实现一个简单的分层贝叶斯推理。我使用 Pymc 推断简单线性模型 y = mx 的参数“m”。我想保存跟踪以便稍后加载。

回答 1 投票 0

R - 仅在一台设备上出现错误:chol.default() 中的错误:顺序的前导次要不是正定的,BVAR

我在 R 中估计 BVAR 模型并得到以下错误: chol.default(V_post) 中的错误:1 阶的前导次要不是正定的 这实际上不是问题,因为我工作......

回答 1 投票 0

如何在pymc2中制作一个简单的贝叶斯神经网络进行多类分类

我想在pymc2中为iris数据集构建BNN模型。 我定义了模型并尝试进行训练,但训练和测试数据的准确度仅为 0.33。 这是我当前的代码 迭代次数 = 2000 伊里...

回答 1 投票 0

arviz.from_pyjags() ValueError:需要解包的值太多(预计为 3 个)

当我使用 az.from_pyjags 时,它引发了错误。 -------------------------------------------------- ------------------------ ValueError Traceback(最近一次调用最后一次)...

回答 1 投票 0

如何并行化 Google Colab 笔记本以使其更快?

我想知道是否有办法在 Google Colab 应用程序上并行化 Jupyter 笔记本。我使用 cobaya 包进行宇宙学分析,并且执行了许多蒙特卡洛马尔可夫 Ch...

回答 1 投票 0

MCMCglmm 具有生态因素和系统发育树:“某些级别的动物在 ginverse 中没有行条目”

我正在尝试针对数量性状和生态因素运行 MCMCglmm,同时考虑系统发育。我正在研究猫头鹰物种的头骨不对称性,所以我使用系统发育树

回答 1 投票 0

R 中的 JAGS:生成的 MCMC 链不会在初始值附近开始

我正在运行这个简单的 Rjags 教程中的代码。请注意,参数的初始值由 model.inits 指定<- list(tau=1, beta0=1, beta1=0, beta2=0). I wanted to test tha...

回答 0 投票 0

基于“R2jags”和“bridgesampling”验证贝叶斯因子计算

我正在尝试通过桥接采样计算单向方差分析的 Jeffreys-Zellner-Siow (JZS) 贝叶斯因子。然而,结果远不如 anovaBF 函数计算的值可靠...

回答 0 投票 0

Gibbs 采样器:数组的第 n 个主要次要不是正定的

我正在尝试构建一个吉布斯抽样算法来估计线性回归的系数和方差。由于假设方差是常数,我尝试使用 MVNormal-conjugate

回答 0 投票 0

使用 MCMCpack 作为潜在变量并从 mcmc 对象中提取潜在分数

我正在使用 R 包 MCMCpack 中的函数 MCMCordfactanal() 和 MCMCirtKd() 生成后验分布并创建潜在因子变量。 我可以使用 summary() 和

回答 0 投票 0

从glmmTMB输出中提取后置模式和可信区间。

我通常使用lme4包,但glmmTMB包越来越成为更好的工作与高度复杂的数据(认为过度分散或零膨胀)。有没有一种方法...

回答 1 投票 1

拟合多项式模型时的运行误差。"部分观测,部分缺失"

我是Jags的新手,我试图对我的数据拟合一个多项式模型。当我运行代码时,我得到以下错误。"positive.counts[1,1:9] 部分观察到,部分缺失"。我上网查了一下, ...

回答 1 投票 0

© www.soinside.com 2019 - 2024. All rights reserved.