Im using lme to calculate a mixed factors ANOVA according to:
px_anova = anova(lme(dep~music*time*group, random = ~1|id, data = px_data))
dep is a threshold,
time is a repeated measures variable (2 levels)
group is a between subjects variable (2 levels)
id is a random factor (subject id)
music is a between subjects variable (2 levels) indicating if a person has a musical experience, or not
Musical experience is now decided by categorizing depending on the number of years practicing playing an instrument.
I would like to use the years of playing an instrument as a covariate instead of creating categories.