I am post doctoral fellow in neurosciences. I have a quite basic question that I failed to solve by myself.
In a pretest experiment for an fMRI study, we have asked to 3 age groups (~30 participants in each group) to listen 6 different musics. Each of the 6 musics is supposed to elicit only one particular emotion. After listening each music, we collected a binary outcome.
To sum up, each participant has 6 binary responses, one for each modality of the music factor.
Intuitively, because the DV is binary and participants answered several "trials", I ran a GEE model. Since I found a significant Music*Age interaction, I ran post-hoc tests sorted by the music factor (via the emmeans package in R) to compare all age groups at each music modality.
After running these analyses, I felt that something was wrong but I am still not clearly able to explain why. Is it correct or are some analyses more appropriate with this design?
Let me know if something should be clarified and thanks for your help!