Maybe you have already seen this, but the new version of modeest (2.3.2) can lead to warnings like Warning message: In runjags.summaries(fullmcmclist = mcmc, thinnedmcmclist = thinnedmcmc, : An unexpected error occured while calculating the mode I think that a small piece of line 280 of runjags.summaries.R needs to be deleted/commented out: return(modeest::mlv(x, ...)) #$M[1]) Best, Ed
Thanks, that's exactly what I was looking for, and exactly what I would have never figured out by myself! Ed
Is there a way to modify the axis labels on trace plots? I ran the first regression model in the examples of the run.jags() help file. Then I tried this command: plot(results, "trace", trace.options=list(ylab=c("1","2","3"))) which gives me the label "1 2 3" on all three panels of the plot. I also tried supplying a list to the ylab argument instead of a character vector, and I tried using the panel argument instead of ylab with no luck (but reading the lattice documentation brought back fond memories)....
All, There exists a handful of methods to calculate effective number of parameters...
Matt, I have enjoyed the runjags package and use it heavily in my "blavaan" package....