jags

    0熱度

    2回答

    我嘗試加載rjags時收到一些錯誤。我使用標準的install.packages(「rjags」),這似乎很好。但是當我加載包時,它不起作用。這是我得到: > library(rjags) Loading required package: coda Error : .onLoad failed in loadNamespace() for 'rjags', details: cal

    2熱度

    1回答

    我是新來JAGS(只是另一個吉布斯採樣),我想知道是否有可能延長它從隨機排列的空間來樣? 我問的原因是,我遇到了關於擴展JAGS的本教程,其中一個所需函數「logDensity」具有所需的輸入作爲double,但要使用排列我需要輸入double。 http://www.cidlab.com/prints/wabersich2013extending.pdf 我也沒有隻JAGS工作,如果任何人有另一

    0熱度

    1回答

    我正在處理實驗設計問題,並嘗試通過R和r2jags來匹配JAGS模型。 要測量殘留效應,我必須訪問其中一個變量列表中的i-1元素。當i=1時,此變量必須返回值列表中的最後一項。我試圖使用ifelse(),但沒有奏效。 我的嘗試: for (i in 1:Ntotal){ j <- ifelse(i==1,Ntotal,j) y[i] ~ dnorm(y.hat[i], tau)

    0熱度

    1回答

    我已經用rjags實現了LDA模型。我順利地拿到了最後的樣本: jags <- jags.model('../lda_jags.bug', data = data, n.chains = 1, n.adapt = 100) update(jags, 2000) samples <- jags.samples(jags, c('th

    0熱度

    1回答

    我組建了一個模型中JAGS爲UScrime數據的層次迴歸(從庫(MASS)包運行代碼。犯罪率作爲響應,用15個預測因子。 我在我的代碼中的一些錯誤,我無法彌補 JAGS型號: model{ for (i in 1:m){ for (j in 1:47){ y[j]~dnorm(beta[i,1]+beta[i,2]*x[j]+beta[i,3]*x[j]+beta

    0熱度

    1回答

    我嘗試估計我的數據的衰減。讓我解釋。想象一下,你觀察事件X,並且每個事件都有一個影響y,隨着時間exp(-t/Tau)衰減。你觀察時間t和事件x以及預測它的影響y。讓我給你看看我的JAGS代碼。 model{ for(j in 1:N){ for(i in 1:p){ td[j,i] <- exp(- t[j,i]/Tau[i]) } mu[j] <-

    3熱度

    1回答

    我正在使用rJAGS構建一個多層次的貝葉斯模型,我想指定一個Cauchy先前的幾個參數。有沒有辦法在JAGS中執行此操作,還是需要切換到STAN?我的JAGS模型如下。我想用Cauchy替換dnorm分佈,但JAGS找不到標準R Cauchy分佈,例如, dcauchy,pcauchy model_string <- "model{ for (i in 1:n){ y[i] ~ dbi

    1熱度

    1回答

    我正在編寫pdf,使用sweave和latex在r。對於MCMC方法,我正在運行rjags。運行這些模型大約需要一個小時才能收斂。每次運行我的sweave代碼編譯PDF時,它也會再次運行所有jags型號。這使編輯和發現,如果我需要花費一個小時來編譯PDF,那麼是否會造成小的語法錯誤。我如何保留我的jags代碼生成的所有變量,但不必每次都要評估sweave? 數據可以在這裏找到:https://uw

    2熱度

    1回答

    我試圖使用索引作爲響應(D47),溫度作爲預測因子(Temp)並考慮離散變量(材料)的隨機效應來建模貝葉斯迴歸。我發現了關於非等級迴歸的非常好的信息,有些帖子甚至包括這些模型的預測策略。儘管如此,在我的模型中預測D47值時,我發現了一個顯着的問題,主要是因爲隨機截取。 在預測JAGS迴歸期間有沒有辦法處理隨機截距? 謝謝您的回答, 克里斯蒂安 model1<-"model { # Priors

    3熱度

    1回答

    我試圖通過標記選票複製模擬一塊jags代碼,但是jags發送給我一個錯誤消息。 如果我正確理解它,它應該與索引的房子有問題影響每個黨的某個地方,但我無法找到它,因爲該節點似乎已經索引。有沒有人有一個想法是什麼錯誤? model <- jags.model(textConnection(model), data = data, n.chains=4,