I have a simple dataset from a within-subject design. Each participant provided a verbal description of 3 stimuli. The descriptions were coded so that they consist from objects each belonging to 1 out of 3 classes. I am interested in the proportion of these objects between the classes. I want to know if some classes were occurring more often than others and if some stimuli caused more often occurrence than others. So we have a following dataset of counts:
Participant stimulus class_a class_b class_c 1 stim1 8 0 3 1 stim2 7 9 5 1 stim3 1 2 4 2 stim1 3 4 6 ...
The simplest way to analyse this data would be a Repeated-Measures ANOVA. However, I would have to run them 3 times - once for each 'class' as a dependent variable. And it doesn't give me the relations between the classes.
Another solution might be MANOVA. But I am a bit confused whether it can be used when there is no between-subject independent variable (like an experimental condition). And how does it deal with the within-subject factor? I searched for R tutorials on that, but they either miss the repeated-measures component or deal with more than 1 experimental condition.
One more idea that I had was to build a mixed-effect model. Since mixed-effect models require a single output variable though, I would need to transform the data to the following format:
Participant stimulus class response 1 stim1 a 8 1 stim1 b 0 1 stim1 c 3 1 stim2 a 7 ...
My concern is, that any model I construct from such a dataset will treat both 'stimulus' and 'class' as Independent Variables / Predictors. In other words, they will lie on the right hand side of the formula: response ~ stimulus*class + (1|Participant).
At the same time, mixed-models seem to offer the highest power, require less assumptions, and provide the most flexibility. Can they be used for this multivariate analysis? Or can multiple models be constructed (one for each class) and somehow related to each other to show the differences across classes?