I have been running several linear mixed effects models for some data of my current project, and now I'm moving on to different data I have. I say that because I'm in a mindset to use LME, and didn't think about good ole ANOVA, though I don't think it's appropriate here.
Here's my design:
For this specific group, there's 30 subjects. Each subject reads 80 sentences in Spanish, and 80 sentences in English.
My first instinct was to run lmer, thinking that Sentence Language (Spa or Eng) is nested within Subject.
mod1 = lmer(totreadtime ~ langcode + (1|subject), REML=F, data=data) > anova(mod1) Type III Analysis of Variance Table with Satterthwaite's method Sum Sq Mean Sq NumDF DenDF F value Pr(>F) langcode 19268687 19268687 1 4740 3.8996 0.04835 * --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Great, there's a difference between Spanish Sentences and English sentences!
Am I correct in using lmer? Just want to verify.