Menu
  • HOME
  • TAGS

Is there a changelog for JAGS?

jags

There is no changelog on web or in the release, one must look directly into the source code: http://sourceforge.net/p/mcmc-jags/code-0/ci/default/tree/NEWS...

Similar function to R “rep” in jags to create array?

r,jags

I have added a rep function to the development version (future JAGS 4.0.0) as Matt and John have alluded to, this requires the second argument to be fixed so that the length of the resulting vector can be determined at compile time.

Step by step right-censored survival analysis in JAGS

survival-analysis,mcmc,jags

I don't do a lot of survival analysis (and you don't state which distribution you would like to use for this part - there are different options), but this code should get you started for the interval censoring part: library("runjags") # Your data t1 <- c(1.73, NA, NA, NA, NA,0.77,...

Too low psrf values in runjags summary?

r,jags,runjags

It appears (from information sent to me off-line) that the psrf is being calculated correctly, but is being reported out of sequence because some information is not available for some semi-stochastic monitored variables. The fact that this isn't being picked up by the software is a bug that I will...

Zeros-one trick for Bayesian Generalized Poisson using JAGS

r,glm,jags

Your lines zeros[i] <- 0 zeros[i] ~ dpois(zeros.mean[i]) cause a problem. In JAGS you can't redefine the given value of a variable. I think you should drop the line: zeros[i] <- 0 from your code...

Logistic regression when response is a proportion (using JAGS)

r,logistic-regression,jags

You don't need to compute p in your data set at all. Just let it be a logical node in your model. I prefer the R2jags interface, which allows you to specify a BUGS model in the form of an R function ... jagsdata <- data.frame(y=rbinom(10, 500, 0.2), n=sample(500:600, 10),...

Inserting estimates from R to jags

r,jags

If you write your JAGS file as a brew package template you can do something like this: est.sd.mu <- <%= Xmu %> est.sd.beta <- <%= Xbeta %> est.sd.phi <- <%= Xphi %> [etc] Then then you process this with the brew function to generate the JAGS file, with the R...

Observed node inconsistent with unobserved parents at initialization in JAGS

r,initialization,jags

Reason of th error: In my data, the response time is in millisecond and when the RT is divided by 1000 (seconds) it works. However, I haven't seen any documentation that limits the use of ms in dwiener module. But I can guess that somthing like NaN in R is...

JAGS: unit-specific time trends

r,bayesian,jags

Your lm example can also be written: m1 <- lm(gdp ~ -1 + rain + factor(ccode) + factor(ccode):year) The equivalent JAGS model would be: M <- function() { for(i in 1:N) { gdp[i] ~ dnorm(yhat[i], sd^-2) yhat[i] <- b0[ccode[i]] + b1*rain[i] + b2[ccode[i]]*year[i] } b1 ~ dnorm(0, 0.001) for (j...

The reason for using update() on a model created with jags

r,jags

The function update performs burn-in iterations for your Markov chain. If you fail to run those burn-in iterations, then your chain will not have converged to its stationary distribution and any samples you generate will not actually be from the correct posterior distribution. In short, you absolutely need to update...