在JAGS中以“计数过程”形式表示参数生存模型
我正在尝试在JAGS中构建一个生存模型，允许时间变量的协变量。我希望它是参数模型 - 例如，假设生存遵循Weibull分布（但我想允许危险率变化，因此指数分布过于简单）。因此，这本质上是可以在flexsurv包中完成的内容的贝叶斯版本，该包允许参数模型中的时间变量协变量。 因此，我希望能够以...
Stan. 使用 target += 语法
我开始学习Stan。有人能解释一下如何使用...这样的语法吗？ target += 不是只有： y ~ normal(mu, sigma) 例如，在Stan手册中，您可以找到以下示例。 model { real ps[K]; // temp for log component d...
如何在rstan中表示分类预测变量？
如何正确格式化分类预测变量以在 STAN 中使用？我似乎无法将分类预测变量作为普通因子变量输入，那么将普通分类变量转换为 Stan 能够接受的最快方式是什么？ 例如，假设我有一个连续预测变量和一个分类预测变量。your_dataset = data.frame(income = c(62085...
如何控制来自fbprophet的输出？
我希望在拟合预测模型时能够抑制来自fbprophet的某些输出。这些输出（“Initial log joint probability...”，“Optimization terminated normally:”，“Convergence detected:...”等等）显然来自Stan的c...
无法在Ubuntu 18.04 LTS上安装rstanarm。
我无法在运行Ubuntu 18.04 LTS的R 3.5.3上安装rstanarm。 我使用了以下命令：install.packages("rstanarm") anddevtools::install_github("stan-dev/rstanarm", build_vignettes =...
超过了R中DLL的最大数量
我正在使用RStan从大量高斯过程（GPs）中进行采样，即使用函数stan()。对于我拟合的每个GP，另一个DLL会被加载，可以通过运行R命令来查看。 getLoadedDLLs() 我遇到的问题是，由于需要适应如此多独特的GPs，我超出了可以加载的DLL的最大数量，此时我会收到以下错误...
抑制RStan在Rmarkdown HTML输出中的运行消息
以下是我的 R 代码。 ```{r message=FALSE, warning=FALSE, cache=0,eval=TRUE, error=FALSE} stan_m1 <- rethinking::map2stan( alist( y ~ dbinom(n, p),...
我如何在STAN中获得最大似然估计的标准误差？
我正在使用Stan中的最大似然优化，但不幸的是optimizing()函数没有报告标准误差： > MLb4c <- optimizing(get_stanmodel(fitb4c), data = win.data, init = inits) STAN OPTIMIZATI...
RStan能在超级计算机上运行吗？
Stan是由Gelman等人开发的新的贝叶斯分析软件。 RStan是一种从R中调用Stan的方式（猜测）。 Stan / RStan是否可以在带有Linux操作系统的超级计算机上运行，如果可以，是否能够利用超级计算机的多处理器？ 我被告知WinBUGS无法在Linux机器上运行和/或无法利...
如何在rmarkdown中使用stan
我希望在R Notebook中使用rstan来获得模型的估计系数。 我有以下stan代码块: ```{stan output.var="rats"} data { int<lower=0> N; int<lower=0> T; real x[T]; ...