Skip to content

Using the mcmcsamp function

2 messages · Parry Clarke, Martin Henry H. Stevens

#
Hello,

I'm building a couple of mixed models using the lmer function.
The actual modelling is going well, but doing some reading on the use of
crossed random effects and the comparison of models with and without
random effects it is clear that I need to generate some Markov Chain Monte
Carlo samples. However, I'm struggling because everyone time I go to
generate a sample I get the following error:

'Error in .local(object, n, verbose, ...) : Update not yet written'

Can anyone tell me what this means and how I can resolve it, I'm using R.
2.7.2 and updated 'lme4' yesterday

The code I use is as follows:

model<-lmer(FemApp~MaleWith+(1|MaleID)+(1|FemID),poisson,data=Compliance)

mcmp<-mcmcsamp(model3, n= 1000)

Thanks in advance

Parry
#
mcmcsamp does not yet operate on non-normal distributions.
Hank
On Sep 27, 2008, at 9:06 AM, Parry Clarke wrote:

            
Dr. Hank Stevens, Associate Professor
338 Pearson Hall
Botany Department
Miami University
Oxford, OH 45056

Office: (513) 529-4206
Lab: (513) 529-4262
FAX: (513) 529-4243
http://www.cas.muohio.edu/~stevenmh/
http://www.cas.muohio.edu/ecology
http://www.muohio.edu/botany/

"If the stars should appear one night in a thousand years, how would men
believe and adore." -Ralph Waldo Emerson, writer and philosopher  
(1803-1882)