我希望JAGS模型中的變量被分配一個概率爲p [i]的值,另一個值爲概率1-p [i]。 p [i]是模型中的變量。 我試過以下,但我不確定它是否合適。恐怕單個值(0或1)會歸因於主體i,並且此值將始終保持不變:對於給定主體i,向量a [i]中的值將始終等於20或30我想得到的是一個包含20s和30s的矢量,其中20的概率等於p [i]。 for (i in 1:n) {
m[i] ~
我想用JAGS來推斷(隨機)純出生過程中的出生率。 在化學的語言,這種模式是相當於反應:X-> 2X與速率的α* X(也可以被看作是一個鏈式反應的模型) 這是R代碼我用於生成過程(在固定時間)以及用於對參數α進行推斷的鋸齒代碼。 library(rjags)
y <- 1; # Starting number of "individuals"
N <- 25 # number of time
我想使用OpenBUGS code來估計cox regression model的beta值。 我修改了示例代碼,因爲在該示例中它只有一個測試版,我需要各種不同數量的測試版,具體取決於我提供的模型。 這是我的openBUGS模型;它運行正常: bugsmodel <- function(){
# Set up data
for(i in 1:N) {
for(j in 1:T) {