Re: [R] Crosscorr.plot

From: Prof Brian Ripley <>
Date: Tue, 10 Jun 2008 12:15:37 +0100 (BST)

I presume that this is the function from package coda (not mentioned and we have no sessionInfo() and the 'replicable code' does not work).

On Mon, 9 Jun 2008, Doran, Harold wrote:

> Just out of curiosity, why might this be occuring:
>> class(x6)
> [1] "mcmc"
>> crosscorr.plot(x6)
And what disturbed you about this? coda::crosscorr.plot ends


when I think it would better end invisible(), but that is something to raise with the maintainer (Cc:ed here), not the list.

return() does return NULL, visibly (as its help page says).

> # Replicable code

It needs 'library(lme4)' to be replicable!

> example(lmer)
> x6 <- mcmcsamp(fm1, n=1000)
> crosscorr.plot(x6)

which does give a credible plot.

Brian D. Ripley,        
Professor of Applied Statistics,
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272866 (PA)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595

______________________________________________ mailing list
PLEASE do read the posting guide
and provide commented, minimal, self-contained, reproducible code.
Received on Tue 10 Jun 2008 - 12:02:55 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Tue 10 Jun 2008 - 13:30:43 GMT.

Mailing list information is available at Please read the posting guide before posting to the list.

list of date sections of archive