[R] mcmcsamp() in lmer

2006-01-18 Thread Andrew Gelman
I am working with lmer() in the latest release of Matrix, doing various things including writing a function called mcsamp() that acts as a wrapper for mcmcsamp() and automatically runs multiple chains, diagnoses convergence, and stores the result as a bugs object so it can be plotted. I recogn

Re: [R] mcmcsamp() in lmer

2006-01-21 Thread Spencer Graves
Dear Prof. Gelman: Thanks for providing such a simple, replicable example. I got the same results you describe under Windows XP Pro with: > sessionInfo() R version 2.2.1, 2005-12-20, i386-pc-mingw32 attached base packages: [1] "methods" "stats" "graphics" "grDevices" "utils"

Re: [R] mcmcsamp() in lmer

2006-01-21 Thread Douglas Bates
On 1/21/06, Spencer Graves <[EMAIL PROTECTED]> wrote: > Dear Prof. Gelman: > > Thanks for providing such a simple, replicable example. I got the > same results you describe under Windows XP Pro with: > > sessionInfo() > R version 2.2.1, 2005-12-20, i386-pc-mingw32 > > attached base pack

Re: [R] mcmcsamp() in lmer

2006-01-21 Thread Spencer Graves
Hi, Doug: Thanks for all your hard work and creativity with mixed-effects models -- and for your quick reply to my question about how to use 'showMethods' to view the source. Best Wishes, spencer graves Douglas Bates wrote: > On 1/21/06, Spencer Graves <[EMAIL PR

Re: [R] mcmcsamp() in lmer

2006-01-21 Thread Martin Maechler
> "DougB" == Douglas Bates <[EMAIL PROTECTED]> > on Sat, 21 Jan 2006 16:13:37 -0600 writes: DougB> On 1/21/06, Spencer Graves <[EMAIL PROTECTED]> DougB> wrote: >> Dear Prof. Gelman: >> >> Thanks for providing such a simple, replicable example. >> I got the same