poissonconsulting / mcmcr

An R package to manipulate MCMC samples
https://poissonconsulting.github.io/mcmcr/
Other
17 stars 2 forks source link

incorrect coercion from mcmc when missing terms! #27

Closed joethorley closed 5 years ago

joethorley commented 5 years ago

mcmclist.RData.zip

library(mcmcr)
mcmclist <- readRDS("mcmclist.RData")
class(mcmclist)
mcmclist[[1]][1,]
names(mcmclist[[1]][1,])
#if you look at rhat[,3] for example, it only exists
#for rhat[1,3], rhat[2,3], rhat[16,3] and rhat[21,3]
mcmcr.object <- as.mcmcr(mcmclist)
mcmcr.object$rhat
mcmcr.object$rhat[1,1,1:21,1:3]
mcmcr.object$rhat[1,2,1:21,1:3]
#now when we look at rhat, it has numerical values everywhere...
#whereas it should not exist for rhat[3,3], rhat[4,3], etc.
joethorley commented 5 years ago

fixed!