采用R软件进行数据统计分析,采用lme4包拟合混合效应模型。依据前人标准删除正确率低于90%的被试(
孙桂芹 等, 2021),低水平对照组中有一人数据被排除。由于三组被试的平均正确率较高,产生了天花板效应(高水平二语组: 98%; 高水平音乐组: 97%; 低水平对照组: 96%),故不做正确率分析。反应时数据剔除错误反应及超出均值±2个标准差的数据(
Li et al., 2017),约4%的数据被删除。为减少异常值影响并使数据更接近正态分布,对反应时数据进行对数转换后再作为因变量(
马拯 等, 2022)。根据Stroop任务语言类型的不同,分别建立了中文和英文任务的混合效应模型。首先,将组别(高水平音乐组、高水平二语组、低水平对照组)、刺激条件(色词一致、色词不一致、无关中性词)及其交互作用纳入固定效应。对于显著的交互作用,采用emmeans包进行简单效应分析,并采用Tukey方法校正事后检验的
p值。在考察混合模型的随机结构时,同时包括被试的随机斜率和随机截距,其中被试和项目作为随机因子。另外,为控制被试二语水平以及音乐经验对实验结果造成的影响,将标准化后的LexTALE和MET得分作为协变量放入模型中。通过比较不同随机效应结构的模型(如是否包含随机斜率),并基于赤池信息量准则(Akaike information criterion, AIC)选择了AIC值最小的模型作为最优模型。