如何使用具有自定义对数概率的MCMC及其矩阵的求解

文章/答案/技术大牛

发布

1回答

pymc3、pymc、mcmc

代码是用PyMC3编写的，但这是一个一般性问题。我想找出哪个矩阵(变量的组合)给我的概率最高。取每个元素的迹线的平均值是没有意义的，因为它们相互依赖。下面是一个简单的例子；为简单起见，代码使用向量而不是矩阵。目标是找到一个长度为2的向量，其中每个值都在0和1之间，因此总和为1。

浏览 16提问于2019-03-27得票数 1

1回答

如何对具有TensorFlow概率的自定义函数使用MCMC抽样

tensorflow、probability

我刚刚开始研究TensorFlow，我不确定如何从自定义概率分布中采样，该概率分布不容易表达为股票分布的组合。如何为MCMC采样器的target_log_prob输入使用自定义函数？

浏览 14提问于2020-09-19得票数 0

回答已采纳

1回答

使用随机矩阵变量定义自定义PyMC分布

python、pymc

我正在尝试定义以下分布：其中R是2x2速率矩阵，我想用数据来求解(它的每一行的和必须是0)。这是对数似然函数，以及一些数据：# The data is in the format: (row index, column index, t) # == sum( log( exp(r*t) ) ) return np.su

浏览 5提问于2018-08-16得票数 1

1回答

TensorFlow概率条件分布的形成

tensorflow、machine-learning、autoencoder、tensorflow-probability

我正在使用Tensorflow概率建立一个VAE，其中包括图像像素以及其他一些变量。VAE的输出：我试图了解如何在此基础上形成其他条件分布，可以与推理方法(MCMC或VI)一起使用。假设上面的输出是P(A，B，C= Z)，我将如何利用这个分布形成一个后

浏览 0提问于2018-11-21得票数 0

回答已采纳

1回答

在pymc中执行两次MAP会产生不同的值

python、distribution、probability、pymc、mcmc

在使用mcmc进行采样后，我得到了以下lambda的后验分布显然，在λ=0.20时后验概率最大，95%的区间是据我所知，MAP给出了一个点估计(具有最大后验概率的lambda的值)，但是当我运行两次MAP时，我得到了不同的值，这不应该是这种情况。我在d相同程序的两次执行中使用map之后，在nd之前打印lambda的值。在使用MAP 0.

浏览 0提问于2014-06-12得票数 0

1回答

如何在转移概率矩阵中创建转移概率矩阵？

matlab、matrix、multidimensional-array、probability、markov

我必须创建一个3D矩阵使用马尔科夫链为基础的方法。为了理解我的问题，我请求您先查看。这种方法有一个大的三维速度和加速度矩阵(即速度和加速度的不同状态)表，表中的.Each单元具有2*2的概率矩阵。在转移概率矩阵的每个单元中，下一次速度和加速度的概率矩阵都是概率</em

浏览 15提问于2018-01-11得票数 0

3回答

在openCV中执行SparseMat (稀疏矩阵)运算

opencv、sparse-matrix

我需要在OpenCV中对稀疏矩阵SparseMat进行矩阵运算(主要是乘法和求逆)。有没有我可以使用的外部代码？(或者我错过了什么？)

浏览 0提问于2011-07-29得票数 9

回答已采纳

2回答

在PyMC中学习离散HMM参数

python、time-series、hidden-markov-models、pymc、mcmc

我正在尝试使用PyMC学习一个简单的离散HMM的参数。我在HMM上的页面上模拟了晴雨如雨的模型。该模型如下所示：我使用的是以下先例。theta_emission_rainy) y_train[i] = pm.Categorical('y_train_%d' %i, theta_emission_sunny) 但是，我不能理解如何使用我做了如下的开始。

浏览 0提问于2014-05-04得票数 1

1回答

是否有一个fortran求解器能够解出下列特征值？

fortran、eigenvalue

A是N乘N矩阵。I是(N2)与(N2)的单位矩阵.B是另一个N乘以N的矩阵，定义为 0 0 0;。x是一个具有N个元素的数组。如何使用特征值求解器求解如下形式的特征值：x=cB x，其中c是特征值？

浏览 1提问于2017-09-22得票数 1

回答已采纳

1回答

VarianceCovariance矩阵的行列式

matlab、correlation、matrix-inverse、determinants、fminsearch

在求解自回归模型的对数似然表达式时，我跨越了在幻灯片9 下给出的方差协方差矩阵。现在，为了使用为了使似然函数表达式最大化，我需要表示方差协方差矩阵出现的似然函数。有人能给我举个例子说明我如何实现(determinant of Gamma)^-1/2吗？除了自回归模型之外，任何其他的例子也可以。

浏览 4提问于2014-11-21得票数 0

回答已采纳

1回答

方程组稀疏矩阵索引的MATLAB内存管理

matlab、memory-management、sparse-matrix、matrix-multiplication

我想知道MATLAB是如何对下面的方程系统进行内存管理的，包括稀疏矩阵的索引。x = A(indices,indices) \ b(indices);x = Aindexed \ bindexed; MATLAB有没有特殊的

浏览 0提问于2018-10-14得票数 1

1回答

如何在MATLAB中创建自定义的“概率分布对象”

matlab、oop、machine-learning、statistics、distribution

我想在MATLAB (使用统计工具箱)中包括一些简单的方法来从现有的概率分布中创建新的概率分布。例如:有限的混合物或复合分布。我的目标是实现与现有概率分布及其周围环境的无缝集成。是否有一些文档或例子应该如何做到这一点？文档只包含预定义发行版的描述。对于“概率分布对象”的类结构，以及如何编写自定义的</

浏览 7提问于2016-03-30得票数 6

回答已采纳

1回答

如何从一个PETSc (用C++ )代码链接到它？

c++、armadillo、petsc

我有一段C++代码(main.cpp)，它使用armadillo来编写一个非常大的稀疏矩阵及其r.h.s，以便像x= A\b那样进行求解。问题：在armadillo中，不存在并行迭代线性稀疏求解器.我想考虑一下PETSc，它有很多有能力的解决者。但是，我不知道如何将它链接到我的main.cpp代码。提示: 1-我的代码是串行的，然而，armadillo使用线程进程。2- PETSc

浏览 1提问于2021-03-01得票数 0

1回答

用特征求解稀疏矩阵

matrix、constraints、sparse-matrix、eigen、particle-system

我正在写一个基于论文的基于约束的粒子系统。约束的雅可比C与系统中粒子的数目和约束的数目成比例。因为每个约束通常只有几个依赖于它的粒子，所以矩阵将非常稀疏。我认为使用本征中的稀疏矩阵来求解这个系统可能是有益的。根据，似乎有几种不同的方法来解决这些稀疏矩阵方程。我的问题是：这些矩阵的大小是否

浏览 3提问于2016-02-24得票数 3

1回答

PyMC -方差-协方差矩阵估计

bayesian、covariance、pymc

我阅读了下面的论文()，其中他们将方差-协方差矩阵Σ建模为：S是标准差的k×1向量，diag(S)是具有对角元素S的对角矩阵，R是k×k相关矩阵我如何使用PyMC实现这一点？下面是我写的一些初始代码：import pandas as pdprior_mu=np.one

浏览 0提问于2014-02-11得票数 0

1回答

对logistic回归结果的认识

excel、regression、logistic-regression

假设我们有以下具有二进制响应输出的数据(优惠券)年消费是1000单位，我的目标是估计如果客户消费超过2000并且有西蒙斯卡，也会有优惠券，首先我已经根据响应数据对数据进行排序，我得到了下面的图片在下一阶段，我计算了每个数据的logit，对于那些最初我选择了以下系数的数据B1 0.1我根据下面的公式计算了L在下一阶段，我计算了e^L (在excel=EXP(D2)=E2/(1+

浏览 2提问于2017-09-08得票数 1

回答已采纳

1回答

我有一个生成系统矩阵的Python脚本。这是连续发生的，在一个处理器上，在一个进程中，没有并行化。我也有一个解算码。代码使用MPI在许多处理器上运行。目前，Python脚本创建矩阵，将其写入文件，通过subprocess.call(["mpirun ....."])调用求解器，求解器从文件中读取矩阵，进行求解，将其写回文件，最后Python脚本从文件中读取结果现在我正在寻找更有效的方法，避免文件<

浏览 5提问于2017-05-16得票数 0

1回答

使用对数比例分布和0绘制

python、matplotlib、data-visualization

我正在尝试绘制一个概率分布(比如k个事件的概率)。它应该绘制为水平轴上的对数刻度，因为k值较大时的行为看起来类似于k^{-alpha}。所以这是一条直线，表示对数图上的大k。但是0发生了。我想以一种易于理解的方式来绘制这个图。举个例子，考虑定义一个概率，这样p_0 = 0.5和对于k= 1, 2, 3, ...，我们设置p_k = Ck^{-2}，如果我正确地计算了C=3/pi^2。这应该等于1，并为k

浏览 21提问于2020-10-03得票数 0

回答已采纳

1回答

在NumPy中按非标准化对数概率从集合中采样

python、numpy、random

我有一个一维np.ndarray，里面充满了定义分类分布的非标准化对数概率。我想从这个分布中抽取一个整数索引。由于许多概率都很小，对对数概率进行归一化和指数化会带来很大的数值误差，因此我不能使用。实际上，我正在寻找一个与TensorFlow的等效的NumPy，它处理非标准化的对数概率。如果NumPy中没有直接实现此目的的函数，那么实现这种采样<

浏览 0提问于2021-01-24得票数 2

1回答

Logistic回归框架

python-2.7、scikit-learn、classification、logistic-regression

培训集-约占整个数据集测试集的80%，其余20%我们是否只使用决策边界图，并根据新的数据点是在图中还是不在图中作出决定？我可以得到系数，但是计算目标的公式是什么？它是乙状结肠伞下的线状聚光体吗？是否有一个函数来计算我们的决定是正确的概率(是或否)？我可以使用</em

浏览 6提问于2015-10-07得票数 0

点击加载更多