使用Stan估计集群分配

文章/答案/技术大牛

发布

1回答

mcmc、stan、rstan

我正在尝试使用Stan拟合一个具有未知均值和协方差的有限高斯混合模型。我知道，由于HMC不能应用于离散分布的样本，因此边缘化技术通常用于使用Stan来推断混合参数。但是，对于我的应用程序，我需要数据集群分配。在Stan中推断它们的最好方法是什么？我们将非常感谢您的建议。

浏览 6提问于2018-01-15得票数 2

回答已采纳

1回答

从.to_dataframe中提取() StanFit4Model ()在Pystan2.18上不起作用

python-3.x、pystan

我用pystan估计了一个模型samples =stan_model.sampling(data = sdata, iter = 10, chains = 1, seed = 42) model_results = {'mdata' : model_data, 'sdata&#

浏览 0提问于2018-09-26得票数 2

回答已采纳

1回答

评估未见数据的对数似然率

r、bayesian、stan、rstan、log-likelihood

[N] log_lik; log_lik[i] = bernoulli_lpmf((y[i] - 1) | p[i, 2]);拟合后，我可以使用loo包提取对数似然：但我想评估未见数据的对数可能性。作为参考，我正在尝试使用估计的LFO-CV来评估和比较我的时间序列模型。(例如) ()

浏览 6提问于2019-03-28得票数 1

回答已采纳

2回答

在R中保存来自贝叶斯估计函数的表(例如latex)

r、bayesian、stan

我想将贝叶斯估计函数的结果(例如brm)保存到文件中(例如以Latex格式)，以便共享结果和准备出版物。我已经找到了从使用JAGS/BUGS 估计的模型制作表格的工具，但是brm使用stan来估计模型。

浏览 24提问于2018-02-28得票数 2

回答已采纳

1回答

当我使用例如doSNOW和foreach并行运行多个Stan进程时，我希望我的代码在两个级别上都并行运行*。取而代之的是，Stan进程被委托给我的工作人员，并且似乎在那里按顺序运行它们的链，就好像一旦它们被分配到一个内核，它们就看不到机器的其他内核，并且认为它们在单核机器上。有没有一种方法可以创建4核节点的集群，我可以将其传递给R中的某个并行化包，这样我就可以从我的机器中获得最大的效率？ *假设我有一台36核的机器，9个Stan场景分别运行4个链。现在，我一次使用9个内

浏览 5提问于2018-03-23得票数 1

1回答

解析R Mclust软件包的结果

r、cluster-analysis

我使用R包mclust来估计数据中的集群数量，并得到以下结果： 2 7 8 9 1: In map(z) : no assignment to 1,3,4,5,6我有9个集群作为最好的，但它没有分配给5个集群。这是否意味着我想使用9或

浏览 3提问于2014-03-10得票数 2

回答已采纳

1回答

独立随机变量的STAN定义先验

bayesian、calibration、stan、rstan

我使用人工神经网络(ANN)创建了一个元模型，因为仿真模型的运行时间非常长。我试图使用贝叶斯校准来估计未观测的参数，其中先验是基于当前的知识，观测数据的可能性是从元模型中估计出来的。查询：我有两个随机变量X和Y，我试图使用STAN获得后验分布。X的先验分布是一致的，U(0,2)。Y的先验也是一致的，但它总是超过X，即，Y~ U(X,2)。既然Y与X相连，我如何在STAN中定义Y的先验分布，以便约束Y>X保持？我是新来的斯坦，所以我希望有任何建议或指导如何进行。非常

浏览 6提问于2021-11-08得票数 0

2回答

为什么斯坦抽样不符合理论值？

stan、rstan

我正在学习stan，并且尝试了一个非常简单的模型(bernoulli)，如下所示，我期望后验抽样的平均值为0.3，因为先验只是一个均匀分布，但stan实际上给出了一个平均值0.33。y ~ bernoulli(theta); N = 10, ) print(fit) model =

浏览 0提问于2019-09-14得票数 1

回答已采纳

1回答

python中项目反应理论的好模块

python、r

下面是我想做的事情的清单：使用GRM等价物使用多维GRM的方程谢谢!

浏览 0提问于2020-02-21得票数 2

2回答

Python中循环的累加和和标准差

python、finance、standard-deviation

)^2 + 0.3*(0.3 - 0.14)^2 = 0.0114Var^(1/2) = 0.10677078 rounded to 0.10677import math Pro = 0 Var = 0 Pro += prob[i]*ret[i] Var +=

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

1回答

测试NATS-以最小的努力在Kubernetes中进行流

go、kubernetes、nats.io、kubernetes-custom-resources、nats-streaming-server

在一个能够访问集群的终端中(在我的例子中，它是一个-cluster )。我使用作为NATS流客户端.下面是我试图连接到NATS流服务器的代码： sc, err := stan.Connect("stan", "test-client") if err !

浏览 3提问于2021-05-07得票数 0

回答已采纳

2回答

我如何在R中将一个字符串拆分成一个I()？

r、string、list

我使用的是包，它有一个函数map()和map2stan()来执行贝叶斯估计(后者作为rstan的接口)。map()和map2stan()采用公式的alist()来定义似然和先验，例如： height ~ dnorm(mu, sigma),s <- "height ~ dnorm(mu, sigma); mu ~ dnorm(178, 20); sigma ~ dunif(0, 50)" 我认为我应该能够对;进行拆分，并将

浏览 11提问于2016-08-30得票数 1

回答已采纳

1回答

如何在R中的stan_glm()函数中更改间隔？

r、rstanarm

当我使用stan_glm()时，我不知道如何改变我的间隔。我的代码： mod_rstan <- stan_glm( data = dat, prior_intercept

浏览 54提问于2020-12-15得票数 0

1回答

最大SQL内存包括SQL内存要求吗？

sql-server、availability-groups、clustering、sql-server-2019

其中94 as作为Max sql server内存分配。有15个数据库的windows集群和AG设置。我认为集群和AG的内存与94 AG的最大sql内存设置是分开的。是否有任何指南来估计内存集群和AG需要多少内存？

浏览 0提问于2022-05-19得票数 -3

1回答

Stan中的多元偏态正态

distribution、stan、rstan

我正在尝试估计一个贝叶斯分层模型，并希望在我的模型中包含一个多变量的偏态正态分布。这不是一个已经在Stan中定义的发行版，但文档似乎表明可以使用Cholesky因子来实现它。例如，Stan 2.15.0文档在pp.333-334上写道：根据标准正态变量对多变量正态分布的重新参数化可以扩展到其他可以概念化为多变量正态污染的多变量分布，如多变量学生t分布和偏多变量正态分布。我考虑过自己在Stan中实现倾斜多变量正态分布，但它看起来没有一个很好的封闭形式的分布，这将使自己直接实现……

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

1回答

如何连接nats流媒体集群

kubernetes、nats.io、kustomize、nats-streaming-server

我是kubernetes的新手，正在尝试设置nats流媒体集群。我正在使用下面的清单文件。但是我对如何在我的应用程序中访问nats流媒体服务器感到困惑。我正在使用azure kubernetes服务。---kind: ConfigMap name: stan-config stan.conf: | clus

浏览 63提问于2020-05-08得票数 1

1回答

Rstudio MCMC上的Rstan运行时间过长(可用CPU和RAM使用受限)

r、rstudio、mcmc、stan、rstan

我实际上是在使用脚本和来自纽约大学的一个类似的数据集，他报告了一个大约18小时的类似DS的估计时间。然而，当我尝试运行我的模型时，它在18小时内不会超过10%。我正在运行一个500ITER，100个预热2链模型，使用超过5个参数的Bernoulli_logit函数，试图通过No U Turn MC过程估计其中的2个。(在每一步，它从随机正态a中提取每个参数，然后估计y，并将其与实际数据进行比较，以查看新参数是否更适合数据) y[n] ~ bernoulli_logit( alpha[kk[n]] + b

浏览 2提问于2015-04-26得票数 1

1回答

如何继承斯坦的S4类

stan、rstan

我使用rstan::stan()开发了一些软件包。我创建了一个函数，其返回值是由S4生成的一个rstan::stan()类对象。为了方便地访问估计值或为数据添加信息，我希望创建一个新的S4类对象，它继承了rstan::stan()的S4类，这样就有了新的插槽。fit <- rstan::stan( model_name=scr, data=data) # This is a fictitious code.character",

浏览 2提问于2018-12-05得票数 0

回答已采纳

2回答

标准中的嵌套模型？

r、stan

假设我想在两个级别上模拟随机效果，即我有两个嵌套级别:父组中的个人和祖父母组中的父组。我知道如何从等示例中编写单个随机效果的基本模型(如下所示)，但我不知道如何编写与在伊默尔。data{ int<lower=0> K; vector[N] price; int<lower=1,upper=J> re[N];parameters{

浏览 0提问于2015-04-01得票数 6

点击加载更多