[R-sig-ME] lme4 question

Thierry Onkelinx thierry.onkelinx at inbo.be
Wed Jul 19 16:19:23 CEST 2017


Dear Ahreum,

Keep the mailinglist in cc.

Your coding is wrong. Look at the model output. You'll see 5 groups for A
and B instead of 4. You need all stimuli in a single variable (e.g. A1, A2,
A2, A4, B1, ...). The model becomes DV ~ type*IV+(1|id)+(0+type|id)+(
1|stimulus)

Best regards,

ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Institute for Nature and
Forest
team Biometrie & Kwaliteitszorg / team Biometrics & Quality Assurance
Kliniekstraat 25
1070 Anderlecht
Belgium

To call in the statistician after the experiment is done may be no more
than asking him to perform a post-mortem examination: he may be able to say
what the experiment died of. ~ Sir Ronald Aylmer Fisher
The plural of anecdote is not data. ~ Roger Brinner
The combination of some data and an aching desire for an answer does not
ensure that a reasonable answer can be extracted from a given body of data.
~ John Tukey

2017-07-19 15:58 GMT+02:00 Ahreum Maeng <amaengwork op gmail.com>:

> No it is not a complete data. I show just part of it as an example.
> The grouping variable is "type": there are two types A (coded as 1) and B
> (coded as 0).
>
> id: individuals
> type: 2 types of stimuli nested within id. repeated within individual
> A_id: there are 4 different stimuli within each type (A and B). This
> variable indicates 4 different stimuli within type A. coded as 1, 2, 3, 4.
> B_id: there are 4 different stimuli within B.coded as 1, 2, 3, 4.
> IV: this variables is characteristics of each stimuli.
> DV: responses of individuals for stimuli.
>
> Please let me know if you need further information about the data set.
>
> Thank you so much for your kind help in advance,
> Ahreum
>
>
> On Wed, Jul 19, 2017 at 7:52 AM, Thierry Onkelinx <
> thierry.onkelinx op inbo.be> wrote:
>
>> Dear Ahream,
>>
>> You need to tell us more about the data set. Is this the complete data?
>> What variable indicates the grouping?
>>
>> Best regards,
>>
>> ir. Thierry Onkelinx
>> Instituut voor natuur- en bosonderzoek / Research Institute for Nature
>> and Forest
>> team Biometrie & Kwaliteitszorg / team Biometrics & Quality Assurance
>> Kliniekstraat 25
>> 1070 Anderlecht
>> Belgium
>>
>> To call in the statistician after the experiment is done may be no more
>> than asking him to perform a post-mortem examination: he may be able to say
>> what the experiment died of. ~ Sir Ronald Aylmer Fisher
>> The plural of anecdote is not data. ~ Roger Brinner
>> The combination of some data and an aching desire for an answer does not
>> ensure that a reasonable answer can be extracted from a given body of data.
>> ~ John Tukey
>>
>> 2017-07-19 14:21 GMT+02:00 Ahreum Maeng <amaengwork op gmail.com>:
>>
>>> Hello,
>>>
>>> I am trying to run the following model:
>>>
>>> DV ~ type*IV+((1|id)+(0+type|id)+(1|A_id)+(1|B_id))
>>>
>>> As you see on the following sample data structure, "type" is repeated
>>> measure where 0=A, 1=B. There are 4 ids within each type A and B.
>>> Thus, I coded "B_id" as "0" when the type is 1 (A) and coded "A_id" as
>>> "0"
>>> when type is 0 (B).
>>> Would it be a right way to deal with this repeated measure issue?
>>>
>>> ​id type    A_id B_id IV    DV
>>> 1  1   0   1 1     3.14
>>> 2  1   0   2 2     4.67
>>> 3  1   0   3 3     4.23
>>> 4  1   0   4 1     7.00 ​
>>> ​1  0   1   0 2     3.00
>>> 2  0   2   0 3     4.77
>>> 3  0   3   0 1     4.25
>>> 4  0   4   0 2     7.12 ​​
>>>
>>> T​hank you so much for your help in advance!
>>>
>>> _______________________________________________
>>> R-sig-mixed-models op r-project.org mailing list
>>> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>>>
>>
>>
>
>
> --
> Ahreum Maeng
> Assistant Professor in Marketing
> KU School of Business
> University of Kansas
> #2183 Capitol Federal Hall
> 1654 Naismith Drive.
> Lawrence, KS 66045
> https://sites.google.com/site/ahreummaeng1/research
> <https://sites.google.com/site/ahreummaeng1/home>
>

	[[alternative HTML version deleted]]



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