[R-sig-ME] Simulating Data from a Linear Mixed Model

Perry de Valpine pdevalpine at berkeley.edu
Tue Nov 10 01:38:28 CET 2015


If you are willing to write the model in BUGS, you can use NIMBLE (R-nimble.org <http://r-nimble.org/>) to simulate from it.  One of the examples from the web site illustrates how to set up a GLMM in BUGS, and you can do it more compactly using linear algebra if you’d like.
Perry


> Date: Fri, 06 Nov 2015 19:21:25 +0000
> From: Douglas Bates <bates at stat.wisc.edu>
> To: "Hintz, F. (Florian)" <F.Hintz at let.ru.nl>,
> 	"r-sig-mixed-models at r-project.org" <r-sig-mixed-models at r-project.org>
> Subject: Re: [R-sig-ME] Simulating Data from a Linear Mixed Model
> Message-ID:
> 	<CAO7JsnRyTD1w=F3zk+ngyCnUuXDzW4xAAV7tPPJT_b+nmvpzsQ at mail.gmail.com>
> Content-Type: text/plain; charset="UTF-8"
> 
> On Thu, Nov 5, 2015 at 8:56 AM Hintz, F. (Florian) <F.Hintz at let.ru.nl>
> wrote:
> 
>> Hi,
>> 
>> I have a question that is very much related to an already existing post (
>> https://stat.ethz.ch/pipermail/r-sig-mixed-models/2007q3/000293.html),
>> however, I don?t seem to be able to get it to run for my purposes.
>> 
>> I would like to simulate additional data based on a linear mixed effect
>> model that has the following structure:
>> 
>> model.full = lmer(rt_log ~ condition + (1+condition|item) +
>> (1+condition|subj), data = data)
>> 
>> The dependent variable is continuous. The fixed factor ?condition? has two
>> levels. Both random factors have random intercepts and random slopes by
>> ?condition?.
>> 
>> Best,
>> Florian.
>> 
>> --
>> Florian Hintz
>> Centre for Language Studies
>> Radboud University
>> Nijmegen (The Netherlands)
>> 
> 
> It happens that http://arxiv.org/abs/1511.01864, which was uploaded
> yesterday, happens to deal with exactly that same model.
> 
> 	[[alternative HTML version deleted]]
> 
> 


	[[alternative HTML version deleted]]



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