1回答
出于个人目的，我想使用MCMC Gibbs采样，我找到了以下MATLAB代码： https://theclevermachine.wordpress.com/2012/11/05/mcmc-the-gibbs-sampler/ 然而，与MATLAB相比，我更喜欢R。我认为我很好地转换了代码的最大部分： library("phonTools")mu<-c(0
浏览 33提问于2020-12-26得票数 0
我正在寻找一个示例代码，使用马尔可夫链蒙特卡罗方法进行图像处理，最好是用于分割，在Matlab或Python中。除了MRP ICM.m示例之外，我找不到任何示例。有什么建议吗？
浏览 2提问于2011-11-03得票数 2
回答已采纳
3回答
我已经使用Rcpp (特别是RcppArmadillo)编写了一个MCMC采样器。我希望我的采样器将周期性的状态消息打印到R控制台，但当采样器终止时，这些消息会立即出现。<< double(iterations) / double(maxit) * 100 << std::endl;
可以扩充这段代码，使其输出立即出现在R控制台上吗R</
浏览 4提问于2012-11-11得票数 5
1回答
我对PyMC比较陌生，关于MCMC采样器的输出我有一个简短的问题。我希望找到MCMC采样器找到的变量的最可能值(后验最大值)。有没有一种快速的方法可以做到这一点？据推测，MCMC采样器发现的后验最大值的变量值可能与PyMC的MAP方法发现的变量值有很大不同。
感谢开发者提供了PyMC。它对我的工作非常有用。这个
浏览 0提问于2013-09-10得票数 3
1回答
我正在尝试使用R包"BRugs“来实现Gibbs采样器，但是生成汇总后验统计信息的函数，如samplesStats()，只返回平均值和中位数。是否有可能提取后验模式？
浏览 2提问于2014-02-19得票数 3
如果我没有错，主题建模(LDA)是不可复制的，也就是说，它在不同的运行中给出了不同的结果。这是从哪里来的(这种随机性是从哪里来的?为什么有必要？)怎样才能解决这个问题，或者在运行中获得更多的稳定呢？
谢谢你的帮助。
浏览 0提问于2015-09-24得票数 2
回答已采纳
1回答
我正在用TensorFlow概率的MCMC包实现一个复杂的采样算法。假设我用Gibbs抽样来表示一个分量，哈密顿蒙特卡罗对另一个分量使用，最后使用随机漫步大都会。我的自定义抽样迭代大致如下： # Sampling `a`: # part o
浏览 4提问于2020-05-05得票数 0
我目前正在用R编写一个MCMC过程，用于估计Rasch模型参数。为此，我在Gibbs采样器中使用metropolis hastings算法。 d.l.c[i] <- d.l.p[i]
我的问题是我知道在非对数似然的情况下，您可以使用if(r
浏览 1提问于2015-03-12得票数 0
1回答
所以我必须运行大型的mcmc，这在R中是不可行的，所以我用C写了循环，然后用.C从R调用它。但是，在这个循环中，我偶尔会有采样步骤或建模步骤，在R中有标准的包，用C语言从头编码这些东西对我来说太过分了。因此，R嵌入到C中，这本
浏览 0提问于2013-08-28得票数 1
回答已采纳
2回答
我想知道为什么采样器在一步一步地取样时速度非常慢。例如，如果我运行：mcmc.sample(1000)有办法加快速度吗？在这里，我更详细地介绍具体的问题：
我在竞争中有两个模型，它们是一个更大的模型的一部分，它有一个分类变量，作为两者之间的“切
浏览 2提问于2014-03-22得票数 0
回答已采纳
4回答
我正在尝试自动化一个过程，这个过程在某个时候需要从截短的多元正常值中抽取样本。也就是说，它是正态分布的多元正态分布(即高斯分布)，但变量被限制在长方体上。我给定的输入是完全多元正态的均值和协方差，但我需要在我的盒子里的样本。所以我想要做
浏览 12提问于2013-11-21得票数 8
回答已采纳
2回答
我正在尝试为Haskell中的假想行星生成随机质量。我想通过抽样双峰分布(理想情况下是两个正态分布的叠加:一个对应于小行星，一个对应于气体巨星)来产生这些质量。我看过，它提供了quantile函数，该函数可以在许多发行版上将统一分布的Double转换为Double。但是似乎没有任何对组合发行版的支持。这种特殊的情况可以通过预先选择一个发行版或另一个发行版进行采样来解决，但我希望使用单一发行版，特别是因为我可能需要在以后调整整个发行版。最终，我可能会用来自天空测量的真实数
浏览 0提问于2012-05-31得票数 6
回答已采纳
1回答
在我的模型中，我需要使用一个复杂的python函数从一组父变量中获得确定性变量的值。下面是一个pyMC3代码，它显示了我试图在一个简化的情况下所做的事情。= pm.find_MAP() # Find starting value by optimization
step = pm.NUTS(state=start) # Instantiate MCMC从一个预先计算的网格中找到y_hat是很重要的<
浏览 3提问于2014-10-16得票数 1
回答已采纳
1回答
R拆分包含括号的公式
、、、
我目前正在编写一个软件包，用MLE或Gibbs采样器来估计随机前沿模型。终端用户应该能够输入如下公式：我的问题是:我不知道如何正确地将公式分成两部分。分裂后的结果应是：2："z1 + . + zn“最好的
浏览 1提问于2018-07-16得票数 0
回答已采纳
1回答
电阻抗断层成像与MCMC
、、、
我想用大都会黑斯廷斯采样器来估计身体的内部空洞。我已经使用了正演模型，这是一种基本解的方法来估计电势和电流的通量值(数据)。我很难使用贝叶斯方法(mcmc)来估计空洞或半径。有人能帮我做点什么吗？之所以这样做，是因为我是MCMC的新手
浏览 2提问于2020-04-08得票数 1
1回答
我有以下未显式启用并行化的R代码： matrix <- matrix(rnorm(1000^2), ncol = 1000)
for (i in 1:10000我的理解是，R中的for循环总是串行的。我确实注意到，对于单个大型矩阵乘法，只利用了一个核心，因此我不相信在矩阵乘法中会发生并行化。这里更大的问题是，我编写了一个MCMC采样器，它需要作为马尔科夫链串行运行，
浏览 8提问于2019-09-27得票数 0
回答已采纳
1回答
小鼠R模拟与回归()
、、、、
我正在使用R中的鼠标包进行多重计算，并试图理解其背后的算法。
从它的文档来看，老鼠算法被认为是被使用的。据我理解，它使用吉布斯采样器来执行MCMC，其中模拟参数β，它定义了给定Y-(所有其他变量都没有Y)的Y(带有缺失值的变量)的条件分布。利用模拟β定义了相应的条件分布。然后从条件分布中提取值，并用它替换缺失的值。它在所有缺少值的变量中<
浏览 1提问于2015-01-28得票数 1
我刚刚开始研究TensorFlow，我不确定如何从自定义概率分布中采样，该概率分布不容易表达为股票分布的组合。 如何为MCMC采样器的target_log_prob输入使用自定义函数？
浏览 14提问于2020-09-19得票数 0
回答已采纳
2回答
(我正在编写性能非常重要的MCMC代码。)我不是C++专家，但基本的技巧在其他地方对我有效。然而，我现在遇到了第二个派生类需要模板化的情况。这似乎会带来问题。我的代码概要是：template<class DerivedStepType>public: }
void
浏览 1提问于2009-06-05得票数 4
云服务器
ICP备案
对象存储
云直播
腾讯会议
腾讯云开发者
