stan

1熱度

1回答

可以說我有一個表格y=a_{i} + b_{i,1}*x_{1} + b_{2}*x_{2}模型，其中i=1,2,...,12和我想估計這個模型使用rstanarm。是否有可能爲每個攔截a_{i}設置不同的先驗（因此可以說，前4個具有normal(location = 0, scale = 1, autoscale = TRUE)，在接下來的4有normal(location = 1, scal

1熱度

1回答

斯坦 - 難度矢量化

我： mu_x - 實際值的10×10矩陣 p型 - 10000長的矢量 sender_name_type - 10000長的矢量我想找到的價值 x_real - 10000長矢量我現在有代碼，和作品，是： for(i in 1:N){ if(x_available[i]){ x_real[i]~normal(x[i],0.01);} else{ x_real[i

0熱度

1回答

斯坦代碼適用於Windows，但不適用於Linux

在Windows上與rstan一起使用時，我有以下Stan代碼工作得非常好。然而，使用Linux（CentOS 6的）集羣上運行時，它引發了一個很長的錯誤，其中包括〜500行的，我想，RCPP碼和最後一個塊如下： Error in compileCode(f, code, language = language, verbose = verbose) : Compilation ERROR, fu

0熱度

1回答

運行map2stan功能錯誤127

我安裝rstan，當我打電話map2stan我得到的錯誤如下： > library(rstan) > packageVersion("rstan", lib.loc = "C://TOOLS//R-3.4.0patched//library" ) [1] ‘2.16.2’ 這裏的調用map2stan library(rethinking) data(rugged) d <- rugge

2熱度

2回答

與Stan有關的分層模型的超級啓發

我正在尋找一個模型來估計與Stan有關的二項數據的多個概率。我對每種概率都使用了Beta測試版，但我一直在閱讀關於使用hyperpriors來收集信息並鼓勵估計收縮的文章。我已經看到了這個例子來定義pymc的hyperprior，但我不知道怎麼做類似的事情斯坦 @pymc.stochastic(dtype=np.float64) def beta_priors(value=[1.0, 1.0]

3熱度

1回答

如何使用rstanarm以APA風格報告貝葉斯線性（混合）模型？

我目前正在努力遵循APA-6建議，如何報告rstanarm::stan_lmer()的輸出。首先，我將在頻率方法中使用混合模型，然後嘗試使用貝葉斯框架來做同樣的工作。下面是可重複的代碼來獲取數據： library(tidyverse) library(neuropsychology) library(rstanarm) library(lmerTest) df <- neuropsy

1熱度

1回答

斯坦的多項式模型 - 如何擬合狄利克雷分佈參數？

我想了解如何在stan中使用Dirichlet分佈。我得每六個級別的因子變量的觀測總數的表： counts n factor_var -------- ------ ------------ 3710 4386 level 1 252 4386 level 2 332 4386 level 3 59 4386 level 4

0熱度

1回答

貝葉斯線性（混合）模型的總體預測能力（例如R2）

正在關注this關於期刊報告的問題，我想知道是否對使用stan_lmer擬合的貝葉斯模型有任何總體效應量指標？在頻率論框架，有例如僞R2（由this包計算的），它返回邊際和條件（方差由固定因子（一個或多個單獨的）解釋的比例）（由固定和隨機因素解釋的方差比例）R2。是否有任何等價物可以幫助我們量化和限定效應/預測功率大小？謝謝。

1熱度

1回答

R STAN - 如何制定多元協方差矩陣

我目前使用R斯坦來擬合多元正態分佈。當前模型是 b〜MVN（0，Sigma公司）其中 B =（X1，X2，X3） 0 =（0,0，0）西格瑪= 我能夠使用建立的共方差矩陣如下： parameters { row_vector b[3]; real<lower=0> b_sigma[3]; real<lower=-1, upper=1> b_rho[3]; }

0熱度

1回答

brm迴歸參數的含義

我正在使用brms包在預測變量x上構建一個帶有高斯過程的多級模型。這個模型看起來像這樣：make_stancode（y〜gp（x，cov =「exp_quad」，by = groups）+（1 | groups），data = dat）這樣一個gp對x預測變量和一個多級組變量。在我的情況下，我有5個組。我一直在尋找代碼（下面），我試圖找出一些參數的含義和尺寸。我看到M_1是組數我的問題是：什