[R-sig-ME] Are there sampling methods for glmer objects?

Joshua Wiley jwiley.psych at gmail.com
Mon Dec 17 18:19:57 CET 2012

 Hi Ben,

Thanks a lot!  I appreciate your comments.



On Mon, Dec 17, 2012 at 8:09 AM, Ben Bolker <bbolker at gmail.com> wrote:
> Joshua Wiley <jwiley.psych at ...> writes:
>> Hi All,
>> The subject pretty much says it all.  Using: lme4_0.999999-0
>> >From the archives, I have seen Ben Bolker mention a function "glmersim".
>> ?mcmcsamp seems to indicate it can take model objects from glmer, but
>> when I try it with a simple random intercept logit, I get:
>> Error in validObject(.Object) :
>>   invalid class "merMCMC" object: Slot sigma must have length 1.
>> I also seem to recall something about a "bootMer" function to do
>> (faster) bootstrapping.
>   See ?"simulate-mer", which replaces glmersim.
>   I've just added a note on the mcmcsamp status to the FAQ:
> http://glmm.wikidot.com/faq#mcmcsamp_status
> (That said, the error message you got does seem a bit odd -- but
> I'm not motivated to pursue it since I don't think it will work
> properly anyway).
>   Note that arm::sim is also available (although I'm not sure
> what algorithm it uses, and it only samples beta, not theta [variance
> component] parameters).
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models

Joshua Wiley
Ph.D. Student, Health Psychology
Programmer Analyst II, Statistical Consulting Group
University of California, Los Angeles

More information about the R-sig-mixed-models mailing list