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@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)
______________________________________________
R-help@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.