Just out of curiosity, why might this be occuring: > class(x6) [1] "mcmc" > crosscorr.plot(x6) NULL # Replicable code example(lmer) x6 <- mcmcsamp(fm1, n=1000) crosscorr.plot(x6)