[R] Using the mcmcsamp function

Hank Stevens HStevens at MUOhio.edu
Sat Sep 27 15:18:37 CEST 2008


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

> 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
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.



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)



More information about the R-help mailing list