腾讯云
开发者社区
文档
建议反馈
控制台
登录/注册
首页
学习
活动
专区
圈层
工具
文章/答案/技术大牛
搜索
搜索
关闭
发布
精选内容/技术社群/优惠产品,
尽在小程序
立即前往
文章
问答
(9999+)
视频
沙龙
1
回答
R
中
的
MCMC
Gibbs
采样器
、
、
、
出于个人目的,我想使用
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
2
回答
Matlab或Python
中
的
MCMC
图像处理示例
、
、
我正在寻找一个示例代码,使用马尔可夫链蒙特卡罗方法进行图像处理,最好是用于分割,在Matlab或Python
中
。除了MRP ICM.m示例之外,我找不到任何示例。有什么建议吗?
浏览 2
提问于2011-11-03
得票数 2
回答已采纳
3
回答
使Rcout输出立即出现在
R
控制台上
、
、
我已经使用Rcpp (特别是RcppArmadillo)编写了一个
MCMC
采样器
。我希望我
的
采样器
将周期性
的
状态消息打印到
R
控制台,但当
采样器
终止时,这些消息会立即出现。<< double(iterations) / double(maxit) * 100 << std::endl; 可以扩充这段代码,使其输出立即出现在
R
控制台上吗
R</
浏览 4
提问于2012-11-11
得票数 5
1
回答
从
MCMC
采样器
输出
中
查找最可能
的
值
、
、
、
我对PyMC比较陌生,关于
MCMC
采样器
的
输出我有一个简短
的
问题。我希望找到
MCMC
采样器
找到
的
变量
的
最可能值(后验最大值)。有没有一种快速
的
方法可以做到这一点?据推测,
MCMC
采样器
发现
的
后验最大值
的
变量值可能与PyMC
的
MAP方法发现
的
变量值有很大不同。 感谢开发者提供了PyMC。它对我
的
工作非常有用。这个
浏览 0
提问于2013-09-10
得票数 3
1
回答
如何在BRugs中找到后验模式
、
我正在尝试使用
R
包"BRugs“来实现
Gibbs
采样器
,但是生成汇总后验统计信息
的
函数,如samplesStats(),只返回平均值和中位数。是否有可能提取后验模式?
浏览 2
提问于2014-02-19
得票数 3
3
回答
主题建模(LDA)
中
的
可复制性/重复性
、
如果我没有错,主题建模(LDA)是不可复制
的
,也就是说,它在不同
的
运行
中
给出了不同
的
结果。这是从哪里来
的
(这种随机性是从哪里来
的
?为什么有必要?)怎样才能解决这个问题,或者在运行
中
获得更多
的
稳定呢? 谢谢你
的
帮助。
浏览 0
提问于2015-09-24
得票数 2
回答已采纳
1
回答
如何在TensorFlow中进行混合采样
、
、
、
我正在用TensorFlow概率
的
MCMC
包实现一个复杂
的
采样算法。假设我用
Gibbs
抽样来表示一个分量,哈密顿蒙特卡罗对另一个分量使用,最后使用随机漫步大都会。我
的
自定义抽样迭代大致如下: # Sampling `a`: # part o
浏览 4
提问于2020-05-05
得票数 0
1
回答
在使用对数似然法时,如何评估metropolis-hastings建议值
的
可接受性?
、
我目前正在用
R
编写一个
MCMC
过程,用于估计Rasch模型参数。为此,我在
Gibbs
采样器
中使用metropolis hastings算法。 d.l.c[i] <- d.l.p[i] 我
的
问题是我知道在非对数似然
的
情况下,您可以使用if(r
浏览 1
提问于2015-03-12
得票数 0
1
回答
mcmc
中
从C到
R
调用
R
、
、
所以我必须运行大型
的
mcmc
,这在
R
中
是不可行
的
,所以我用C写了循环,然后用.C从
R
调用它。但是,在这个循环中,我偶尔会有采样步骤或建模步骤,在
R
中有标准
的
包,用C语言从头编码这些东西对我来说太过分了。因此,
R
嵌入到C
中
,这本
浏览 0
提问于2013-08-28
得票数 1
回答已采纳
1
回答
用螺母(
MCMC
)
采样器
时离散潜变量
的
变维
、
感谢您抽出时间阅读我
的
问题如下。
mcmc
.run(model, data) 当我运行
M
浏览 3
提问于2020-11-04
得票数 1
回答已采纳
2
回答
PyMC:一步一步取样?
、
我想知道为什么
采样器
在一步一步地取样时速度非常慢。例如,如果我运行:
mcmc
.sample(1000)有办法加快速度吗?在这里,我更详细地介绍具体
的
问题: 我在竞争中有两个模型,它们是一个更大
的
模型
的
一部分,它有一个分类变量,作为两者之间
的
“切
浏览 2
提问于2014-03-22
得票数 0
回答已采纳
4
回答
SciPy
中
截断
的
多元正态?
、
我正在尝试自动化一个过程,这个过程在某个时候需要从截短
的
多元正常值
中
抽取样本。也就是说,它是正态分布
的
多元正态分布(即高斯分布),但变量被限制在长方体上。我给定
的
输入是完全多元正态
的
均值和协方差,但我需要在我
的
盒子里
的
样本。所以我想要做
浏览 12
提问于2013-11-21
得票数 8
回答已采纳
2
回答
如何在Haskell
中
从复杂或复合分布
中
采样?
、
、
、
、
我正在尝试为Haskell
中
的
假想行星生成随机质量。我想通过抽样双峰分布(理想情况下是两个正态分布
的
叠加:一个对应于小行星,一个对应于气体巨星)来产生这些质量。我看过,它提供了quantile函数,该函数可以在许多发行版上将统一分布
的
Double转换为Double。但是似乎没有任何对组合发行版
的
支持。这种特殊
的
情况可以通过预先选择一个发行版或另一个发行版进行采样来解决,但我希望使用单一发行版,特别是因为我可能需要在以后调整整个发行版。最终,我可能会用来自天空测量
的
真实数
浏览 0
提问于2012-05-31
得票数 6
回答已采纳
1
回答
如何在PyMC
中
定义通用确定性函数
、
、
、
、
在我
的
模型
中
,我需要使用一个复杂
的
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
回答
为什么下面的for循环使用
R
中
我机器上
的
所有内核?
、
、
、
、
我有以下未显式启用并行化
的
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
1
回答
如何对具有TensorFlow概率
的
自定义函数使用
MCMC
抽样
、
我刚刚开始研究TensorFlow,我不确定如何从自定义概率分布
中
采样,该概率分布不容易表达为股票分布
的
组合。 如何为
MCMC
采样器
的
target_log_prob输入使用自定义函数?
浏览 14
提问于2020-09-19
得票数 0
回答已采纳
2
回答
模板化
的
Barton和Nackman技巧问题
、
、
、
(我正在编写性能非常重要
的
MCMC
代码。)我不是C++专家,但基本
的
技巧在其他地方对我有效。然而,我现在遇到了第二个派生类需要模板化
的
情况。这似乎会带来问题。我
的
代码概要是:template<class DerivedStepType>public: } void
浏览 1
提问于2009-06-05
得票数 4
点击加载更多
扫码
添加站长 进交流群
领取专属
10元无门槛券
手把手带您无忧上云
热门
标签
更多标签
云服务器
ICP备案
对象存储
云直播
腾讯会议
活动推荐
运营活动
广告
关闭
领券