[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