[R-sig-ME] Specifying the correct LMM for 'unsual' data
Maarten.Jung at mailbox.tu-dresden.de
Thu Jan 25 15:35:06 CET 2018
a colleague of mine asked me to help her planing a linear mixed models
analysis and, as handling her data and the corresponding research questions
with lmer seems kind of tricky to me, I hope one of you can help me along.
The experiment is as follows:
Participants (46 younger and 45 older children) looked at a series of
pictures (one picture per trial) and had to solve two tasks consecutively:
- Task block 1: Prospective memory (PM) task: while doing other tasks,
participants had to remember to press a specified button when they saw a
- Task block 2. Visual search: participants had only this one task –
pressing a button as soon as possible when seeing a certain object
Each child saw the same pictures in the same task block – pictures 1-6 in
task block 1 and pictures 7-18 in task block 2. Each picture was presented
only once, so there were different pictures in the task blocks.
Trials with target object in task 1 are allocated regarding the
participant’s reactions in PM hits (participants did press the button) and
PM misses (participants did not press the button). (Therefore, a certain
picture can be a PM hit trial for one child and a PM miss trial for the
other.) As there were six trials (= pictures), which contained the target
object, each participant can have a minimum of zero and a maximum of six PM
hits with the according number of PM misses.
Here is the number of PM hits per age group:
- 2 children: 0 hits
- 9 children: 1 hit
- 8 children : 2 hits
- 12 children: 3 hits
- 4 children: 4 hits
- 4 children: 5 hits
- 7 children: 6 hits
- 2 children: 0 hits
- 3 children: 1 hit
- 4 children: 2 hits
- 6 children: 3 hits
- 7 children: 4 hits
- 11 children: 5 hits
- 12 children: 6 hits
(In the visual search task almost all children have pressed the button
correctly in all 12 visual search target trials).
She is interested in how long participants looked at the PM and visual
search target, respectively, depending on if it was a PM hit, a PM miss or
a visual search hit and how this is influenced by the age group. Therefore,
she has got only one data point per trial. And if a participant has no PM
misses there is no data point at all in this condition for this participant.
The variables are defined as follows:
- age_group: categorical predictor with 2 levels (younger and older
- condition: categorical predictor with 3 levels (PM hit, PM miss, visual
My suggestion for the maximal linear mixed model would be:
lmer(dwell_time ~ age_group*condition + (1 + condition|participant) +
I intentionally didn`t use (1 + condition|picture) here because there are
different pictures in the task blocks (see above) - hope this makes sense.
I have two questions:
1. Am I correct with the maximal linear mixed model specifications?
2. I think that the data points in the PM-miss-condition (or
PM-hit-condition) are not missing at random because they are missing if
(and only if) there are 6 data point for the same participant in the
PM-hit-condition (and vice versa). Do you think one has to worry about this
and are there any suggestions how to deal with it?
[[alternative HTML version deleted]]
More information about the R-sig-mixed-models