runjags

Inserting jagam code into runjags (JAGS) model

时光毁灭记忆、已成空白 提交于 2021-02-10 18:55:19
问题 I’ve been trying to incorporate smoothing into a runjags model that I’ve created to model seabird burrow numbers and distribution across an island. I’ve managed to generate some smoothing code by extracting the count data and x and y coordinates from the model output and using the JAGAM tutorial on this page http://www.petrkeil.com/?p=2385 I think I might be able to improve model performance by incorporating the smoothing into the jags model but I’m at a loss as to how to do this. Can you

Inserting jagam code into runjags (JAGS) model

与世无争的帅哥 提交于 2021-02-10 18:55:16
问题 I’ve been trying to incorporate smoothing into a runjags model that I’ve created to model seabird burrow numbers and distribution across an island. I’ve managed to generate some smoothing code by extracting the count data and x and y coordinates from the model output and using the JAGAM tutorial on this page http://www.petrkeil.com/?p=2385 I think I might be able to improve model performance by incorporating the smoothing into the jags model but I’m at a loss as to how to do this. Can you

Inserting jagam code into runjags (JAGS) model

独自空忆成欢 提交于 2021-02-10 18:54:37
问题 I’ve been trying to incorporate smoothing into a runjags model that I’ve created to model seabird burrow numbers and distribution across an island. I’ve managed to generate some smoothing code by extracting the count data and x and y coordinates from the model output and using the JAGAM tutorial on this page http://www.petrkeil.com/?p=2385 I think I might be able to improve model performance by incorporating the smoothing into the jags model but I’m at a loss as to how to do this. Can you

Fitting a multivariate dirlichet model in JAGS for R

家住魔仙堡 提交于 2019-12-14 01:28:35
问题 I am trying to fit a a multivariate model to species composition data using JAGS, implemented in R. I have data on 3 species relative abundances (bounded between [0,1]), two of which are correlated. Here is code to generate similar data. #generate some correlated fractional composition data. y1 <- runif(100,10,200) y2 <- y1*1.5 + rnorm(100, sd = 5) y3 <- runif(100,10,200) total <- y1+y2+y3 y1 <- y1/(total) y2 <- y2/(total) y3 <- y3/(total) y <- data.frame(y1,y2,y3) y is a data.frame of my

runjags object is too big

夙愿已清 提交于 2019-12-07 17:55:40
问题 I am trying to run JAGS using a new package runjags, because R2jags has a bug (the full model code is in the question https://stats.stackexchange.com/q/62006/5509): require("runjags") out <- run.jags("Poisson.OD.t.test.txt", params, win.data, nc, inits, nb*4/5, ni, nb*1/5) plot(out2, layout = c(4, 2)) It works as a charm, but the drawbacks of this package is that the runjags object returned by the run.jags function is already bundled with prepared charts and outputs and is too big. Just for