如果你已经用过PROCESS插件做中介、调节分析你大概率体会过那种选好模型、填好变量、点一下运行结果表格哗啦一下全出来的爽快感。但这里有个前提PROCESS默认你的样本是互相独立的观测。一旦数据结构变成学生嵌在班级里员工嵌在团队里病人嵌在医院里你再用PROCESS硬跑结果就很容易被审稿人质疑甚至直接打回。这时候你需要的是MLmed一个专门在SPSS里跑多水平中介分析的插件。这篇文章我就把MLmed在检验多水平中介的调节模型时的完整过程拆给你看包括模型设定、数据准备、界面操作、结果解读以及最让人头疼的0xc0000005报错怎么一步步排查。所有内容都基于我自己的实际使用经验适合正在处理嵌套数据的量化研究者也适合刚接触多水平模型、不知道从哪下手的同学。1. 为什么嵌套数据里的中介模型不能直接拿PROCESS硬跑1.1 你的数据什么时候已经嵌套了很多研究者一开始并没有意识到自己的数据有多水平结构。举个最直白的例子你想研究教师支持程度X是否通过学生学业自我效能感M影响学业成绩Y数据来自40个班级、每个班30个学生。这个数据天然就是两层的第一层是学生个体第二层是班级。同一个班的学生共享同一个班主任、同一种班级氛围他们之间的测验分数并不是相互独立的。类似的场景还非常多重复测量数据里同一被试在不同时间点的观测嵌套在个体内纵向追踪数据里每次测量嵌套在受访者内组织研究里员工嵌套在部门或团队里教育研究里学生嵌套在班级、学校甚至地区里临床研究里多次就诊记录嵌套在患者里。判定标准其实就一句话你的数据里有没有一个分组变量让组内的观测比组间的观测更相似如果有这个变量就是你的集群变量cluster variable也就是第二层单位。1.2 忽略组结构中介效应会发生什么偏差多水平数据最核心的问题是违背了回归分析里观测独立的假设。你可以这样理解跑中介模型时普通回归包括PROCESS用的OLS路径估计假设每个人的残差都是独立的。但嵌套数据里同一个班级的学生受到相同的班级环境影响他们的残差天然是相关的。IBM SPSS不报错PROCESS也不报错但你拿到的那一排标准误会偏小p值会偏小本来不显著的路径可能跑着跑着就显著了。这就是第一类错误膨胀。偏小多少这取决于组内相关性ICCIntraclass Correlation Coefficient。ICC越高同一个组里的人越像标准误被低估得越严重。有的模拟研究显示当ICC在0.10到0.30之间时忽略层级结构对中介路径显著性判断的影响就已经不能忽略了。更麻烦的是另一个问题嵌套数据的组内和组间效应可能是不同的。一个变量对中介变量的影响在个体层面within和组间层面between的方向、大小甚至性质都可能不一样。PROCESS干的是把两层混在一起估计一个平均关系这个关系在两个层面都不一定准确理论上还不容易解释。1.3 PROCESS和MLmed的分工一个处理单层一个处理多层MLmed的全称是Multilevel Mediation从名字就能看出它是专门为多层数据里的中介分析设计的。它底层采用的是多水平结构方程MSEM的思路先把每个变量分解成组均值between部分和个体偏差within部分然后在两层分别放中介路径。这样做有两个直接好处一是残差的独立性被正确处理了标准误可信二是你能看到中介效应到底是发生在个体层面比如我在这节课感觉自己更有信心所以考得好还是组间层面比如教学支持好的班级全班平均信心更高带动全班平均成绩上升。那是不是说PROCESS就没用了当然不是。如果你的数据确实是一层每个观测之间相互独立PROCESS依旧是效率高、结果清晰的好工具。两者的关系不是替代而是工具匹配问题。这篇文章后面每一节我都会拿PROCESS的常见习惯做参照方便你理解MLmed的特殊之处。2. 安装MLmed之前必须确认的版本与运行环境2.1 下载渠道与文件构成先说一个容易踩的坑你在搜索引擎里搜process插件会出来一堆名字里带process但完全不相干的东西比如某些游戏插件、系统进程工具。这里说的process插件指的是SPSS里的宏macro文件。PROCESS是Hayes开发的单水平中介调节宏而MLmed是处理多水平数据的另外一个宏文件两者不是同一个东西不要混淆。MLmed通常以几个文件的形式发布核心通常是一个.sps格式的语法文件可能还附带说明文档和示例数据。.sps文件本质上是一段SPSS语法代码它不像普通软件那样装完就有图标而是要放进SPSS里运行一次把宏定义加载到当前会话里。下载的时候务必看清楚你拿到的版本适配哪些SPSS版本、是否区分32位和64位。把宏文件解压到一个路径里建议路径里不要有中文、不要有空格比如C:\MLmed\这种。这个细节看起来不起眼但很多加载失败都是路径问题引起的。2.2 SPSS版本、系统位数与兼容策略在我用过的几个SPSS版本里MLmed的兼容性总体不错SPSS 22以上基本都能跑。但有一个经验必须分享能用64位SPSS就不要用32位。原因很简单多水平模型要估计的参数比单水平模型多数据量又不小32位进程在内存使用上很容易到瓶颈。真跑到一半弹出内存不足或者干脆崩溃你前面的工作量就全白搭了。另外MLmed在Windows环境下相对稳定Mac版SPSS上跑宏出诡异问题的概率更高。如果你的项目时间很紧我建议直接用Windows环境。还有一点很多人的SPSS是网上各种渠道装的版本被精简过有些组件缺失。MLmed在运行时要调用矩阵运算相关模块如果你的SPSS是精简版加载宏时可能报错说找不到某个过程或函数。这时候换一个完整版或者换SPSS官方试用版往往就好了。2.3 首次运行前的宏安全设置与路径权限SPSS对宏文件是有一套安全机制的。第一次运行宏文件时如果SPSS拦截了你需要去SPSS的选项设置里找到宏安全相关选项把它调整到允许执行本地宏文件的级别。每个版本SPSS的菜单名称不太一样但大致都能在Options或Preferences里找到。还有一个很容易被忽略的点如果你把宏文件放在C盘根目录或者Program Files下面而SPSS是以普通用户身份启动的它可能没有权限去读取或写入临时文件。解决办法很简单右键SPSS图标选择以管理员身份运行。如果你手上同时有PROCESS和MLmed记得两个宏文件不要放在同名同目录下覆盖因为它们的核心宏名称不同理论上可以共存但如果你下载的是奇怪的集成版保不齐里面的宏名冲突。我习惯的做法是把PROCESS放一个文件夹MLmed放另一个文件夹谁用就加载谁互不干扰。3. 建模型前先搭数据判断路径结构、中心化与ICC计算3.1 先判断你的中介模型属于哪条路径结构多水平中介的第一步不是打开插件而是拿起笔在白纸上写清楚X、M、Y、Cluster这四类变量到底各在哪个水平上。所谓水平就是数据的分析单位。如果变量值对同一个组里的所有人都一样那它就是第二层变量如果每个个体的取值可以不同那它就是第一层变量。常见的中介路径结构有这么几种模型结构X所在层M所在层Y所在层典型研究场景2-1-1第二层第一层第一层班级氛围通过个体效能感影响个体成绩2-2-1第二层第二层第一层班级氛围影响班级平均效能再影响个体成绩1-1-1第一层第一层第一层个体感知通过情绪状态影响行为控制组结构1-2-1第一层第二层第一层个体互动形成团队氛围再回落个体结果2-1-2第二层第一层第二层组织政策影响个体体验再影响组织绩效我这里8个H2不我计划6个H2足够了。3.2 中心化策略为什么默认建议组均值中心化多水平模型里中心化不是可有可无的数据清洗步骤它直接决定你估计出来的是组内效应还是组间效应。拿第1章的例子来说如果X是班级层面的变量每个学生的X值都是一样的那X只有组间变异怎么中心化都不影响组间路径。但如果M和Y是第一层变量它们自身的变异包含组内、组间两部分不处理就会混在一起。组均值中心化的操作很简单用个体原始值减去该个体所在组的均值。本质上是把我在这节课的自我效能感转化为我相对于本班平均水平的自我效能感。这么做可以去除组间差异带来的混淆让M在个体层面的变异更干净。但组均值中心化有个副作用它把组间关于M的信息全删掉了。所以MLmed这样的多水平中介工具在做组均值中心化后还需要把组均值作为一个显式变量回加进模型以便同时估计组间路径。这个技术细节是程序内部处理的但你在理解结果时需要知道输出里的within部分和between部分分别来自这两套变异来源。总均值中心化grand-mean centering则是用个体原始值减去全体平均值它不影响组间和组内的分解只是让截距和模型收敛更容易。如果MLmed提供了中心化选项我个人的建议是第一层变量做组均值中心化第二层变量做总均值中心化或者不中心化也问题不大具体以你下载版本里给的经验为指导不要瞎改。3.3 用SPSS自带Mixed Models手算ICC我一直建议在跑MLmed之前先算一下ICC这能帮你看清数据里到底有多少组间变异。手算ICC的步骤并不复杂。SPSS里打开Analyze - Mixed Models - Linear...把Y变量放进因变量框Fixed按钮里只勾选InterceptRandom按钮里勾选InterceptSubjects里放你的集群变量比如班级编号。其他的都不要动直接运行。输出结果里会有两个关键的协方差参数一个是Residual个体层面残差方差记作σ²另一个是Intercept [subject cluster]组间截距方差记作τ00。ICC τ00 / (τ00 σ²)。举例来说我跑过的一批学生数据Residual估计为0.72截距方差为0.28ICC算下来就是0.28/(0.280.72)0.28。意思是学生成绩的差异中有28%来自班级之间的差异。这个比例不低绝对不能忽略层级结构。如果ICC很接近0比如小于0.05那用单水平模型问题不大如果超过0.10甚至0.20老老实实走多水平路线。ICC这个数字以后写到论文的方法部分时也很有用很多审稿人第一个问题就问它。4. MLmed界面逐项操作从加载宏到跑出完整结果4.1 菜单入口与变量框的对应关系MLmed加载完成之后通常会出现一个图形界面对话框也可能是通过Analyze菜单下的某个入口打开不同版本入口略有差异。你拿到宏文件时附件里一般都会有一份用户指南第一次用之前花十分钟把指南里的界面截图看一遍绝对比瞎点效率高。对话框里的核心字段通常包括X自变量、M中介变量、Y因变量、W如果存在调节变量、Cluster集群变量。注意区分如果你做的是有调节的中介模型调节变量可能作用于中介路径aX→M也可能作用于路径bM→Y不同版本里这个选择的位置不一样有的直接在Model选项里选有的需要单独指定。选之前想清楚你的理论假设到底是什么是班级支持对效能感的影响被班级凝聚力调节还是效能感对成绩的影响被个体以往成绩调节这一步选错了后面读结果会非常拧巴。4.2 模型设置里每个选项的实际含义MLmed对话框里的选项大致可以分为三类。第一类是模型结构也就是第3节说的路径类型你选2-1-1还是2-2-1程序才知道哪些变量该放在哪个水平的路径上。第二类是估计方法。比较常见的是全信息极大似然ML和限制性极大似然REML。在固定效应是关注点的时候ML和REML对点估计差异不大但如果你的随机斜率比较多、组数又不太够调整估计方法往往是解决不收敛问题的常用手段。第三类是随机效应设置。默认一般会给随机截距也就是允许不同组在Y的基线水平上有所差异。要不要随机斜率则需要谨慎让某个路径的系数在不同组之间自由变化意味着你要估计更多方差参数对第二层样本量要求更高。我的经验是除非你有明确理论预期比如某个路径系数在组间差异确实值得研究否则先跑随机截距模型等结果稳定了再慢慢加随机斜率。另外输出选项里关于间接效应置信区间的设定也值得注意。多水平中介的间接效应分布通常不是正态的MLmed一般会提供Monte Carlo置信区间方法模拟次数可以自己填常见做法是设置成20000次。次数越多区间越稳定但运行时间越长。跑大样本时我第一次直接设了50000次结果等了很久才出结果。后来我摸到规律20000次已经足够稳定除非数据分析阶段反复被要求调整否则没有必要贪高。4.3 语法方式运行作为一种稳定替代图形界面虽然方便但如果你追求可重复性我更推荐把核心命令写成语法文件保存下来。MLmed本身以宏形式存在它的语法调用逻辑与PROCESS类似大致思路是指定数据、依次传入各变量、设定模型类型、再给输出选项。不同版本的参数名可能不一样以下是示意格式具体字段名以你下载的说明书为准* 加载宏文件. INSERT FILEC:\MLmed\MLmed.sps. * 调用宏具体命令名和参数因版本而异. MLmed XTeacherSupport /MSelfEfficacy /YAchievement /WCohesion /CLUSTERClassID /MODEL2-1-1 /CENTERGROUP /MC20000 /RANDOMINTERCEPT.把语法存成一个.sps文件的好处是以后改了数据或换了电脑只要运行同一份语法就能得到相同的结果。你在写论文、做复现报告时把这份语法作为附件或补充材料放出去可比贴几张截图有说服力多了。5. 输出表怎么读才不白跑直接效应、间接效应与条件中介5.1 固定效应路径系数a、b、c以及交互项MLmed结果输出的第一大部分是固定效应表直接对应你模型里的核心路径。a路径是X对M的效应b路径是M对Y的效应c是X对Y的直接效应。如果模型里放了调节变量W表格里还会出现交互项比如a路径的交互项反映W是否调节了X→M。读这一部分的重点是看系数方向和显著性但不要只盯着p值。我建议你把未标准化系数、标准误、t值或z值、p值一起抄下来写论文时这些都是必需的。由于MLmed把变异分成within和between两层固定效应表里有些路径可能会有两个版本。比如在2-1-1模型里X是第二层变量它到M的路径通常只会出现在between部分但M到Y的路径可能同时有within和between两个估计——个体层面M的变化和组间层面M的平均水平变化对Y的影响可能是不同的。很多第一次跑的人看到结果里同一个M→Y路径出现两行就懵了其实这正是多水平模型的优势它允许你回答个体层面的中介和组间层面的中介分别是否显著。5.2 方差组分和ICC组间效应占比到底多少固定效应之外输出里还会有随机效应表也就是方差组分。里面通常包含第二层的截距方差和第一层的残差方差。用这两个值你就能回算模型允许随机效应变化后的ICC检查加入X、M和W之后组间方差的解释程度。如果加入第二层X之后截距方差明显下降说明X确实解释了组间差异的一部分。这一部分还能帮你发现一个常见问题如果第二层方差组分不显著说明你的数据组间差异本身就很小多水平中介可能并没有比单水平中介提供更多有效信息。这不是模型跑错了而是数据本身不支持复杂的组间假设。方差组分还能用来判断一个更隐蔽的问题你的模型是否过参数化。当你发现某个随机斜率的方差估计是0或者非常接近0并伴随收敛警告时通常意味着这个随机斜率在你的数据里根本不存在应该把它从模型里去掉。5.3 间接效应的分解与Monte Carlo置信区间这一部分是整个输出里最有价值的区域。MLmed会把间接效应indirect effect按路径来源拆分出来你可能看到within部分的a×b、between部分的a×b以及总的间接效应。对于调节的中介模型还会看到在不同调节变量水平下的条件间接效应conditional indirect effect。判断间接效应是否显著的做法现在的主流已经不看简单的p值了而是看置信区间是否包含0。MLmed提供的Monte Carlo置信区间会给出一个区间估计比如[0.02, 0.15]不包含0就说明对应的间接效应在0.05水平上显著。需要牢记的规律是即使一条路径的p值略大于0.05只要间接效应置信区间不包含0中介效应依然可能显著反过来也一样路径各自显著并不代表间接效应一定显著。所以一切以间接效应区间的结论为准。5.4 调节效应的实操解释当模型里有调节变量W时条件间接效应表通常会列出W取低值均值减一个标准差、中值均值、高值均值加一个标准差时间接效应的大小和置信区间。假如你的W调节的是a路径那这个表就是在告诉你在不同水平的W之下X通过M影响Y的中介作用是否成立。如果低值下置信区间包含0高值下不包含0那就能讲出一个有层次的故事中介效应只在某种条件下存在。这里我最想提醒的是调节效应显著不显著不能只看交互项p值。交互项只是斜率差异的检验条件间接效应才是你研究里的实质证据。有的情况下交互项边缘不显著但条件间接效应呈现清晰的趋势这时候最好多结合理论解释而不是机械地说因为没有交互效应所以模型不支持调节。6. 0xc0000005报错排查链路以及其它常见翻车点6.1 3221225477到底代表什么很多人在网上搜process插件时会看到这样一条报错信息process exited with code 3221225477 / 0xc0000005 (memory access violation)。如果你在跑MLmed时撞到这个错先别慌看明白它的含义0xc0000005是Windows系统层面的STATUS_ACCESS_VIOLATION也就是访问违规——某个进程试图读取或写入它没有权限访问的内存地址。它不是SPSS语法错误也不是你的公式写错了说明SPSS在调用插件或底层组件时撞上了操作系统的内存保护机制。这种报错在涉及大量矩阵运算或第三方宏加载时偶尔会出现尤其在高负荷运算或插件与SPSS版本不兼容时更容易触发。它不一定是机器坏了但确实需要你按顺序排查环境问题。6.2 从最外层到最内层的六步排查我自己遇到过一次当时模型跑了两遍就崩后来一步步排查才发现是杀毒软件把宏文件运行时的临时进程给盯上了。下面这六步是我踩坑之后总结的排查顺序从成本最低的开始。第一步重启SPSS并关闭所有不相关的程序。这一步看着简单但多水平模型的运算会吃大量内存内存被占满时系统就可能触发访问违规。第二步用管理员身份重新启动SPSS。右键图标选以管理员身份运行再重新加载宏。这一步能解决权限不足导致的写临时文件失败。第三步核对SPSS位数与插件版本。64位SPSS就配64位的宏文件如果你装的SPSS实际上是X86版本而系统是64位Windows趁早换掉。第四步检查宏文件路径。文件名和路径里不要有中文、空格、特殊字符。把宏文件放到干净路径并确认当前用户对该文件夹有完全读写权限。第五步暂时关闭杀毒软件或安全卫士的实时防护重新跑一次。如果问题消失那就是安全软件在拦截插件进程。跑完以后再打开防护就行。公司电脑如果有统一安全策略这一步要跟管理员确认不要自己乱关。第六步清理SPSS临时目录或者换一个SPSS版本重装。这一步成本最高但如果前五步都没解决那大概率是SPSS安装本身出现了组件缺失或损坏。如果以上都试完还报同样错误换个电脑跑一次是最快的排除法。如果换电脑就正常那么问题就锁定在你原机器的系统环境上了。6.3 宏没有被识别、收敛失败、结果全是NaN等其它常见报错除了0xc0000005还有几个翻车点经常出现。第一种是宏没有被找到或提示未知命令。原因几乎都是没有先加载宏文件。你要先通过File - Open - Syntax打开MLmed的.sps文件全选运行把宏定义加载进当前会话再执行你的分析语法。如果直接输入主命令而不先加载宏SPSS当然不认识它。第二种是收敛失败或输出警告信息。多水平模型的自带难点是迭代收敛。常见解法有把随机斜率改成随机截距、把估计方法从ML换成REML、增加最大迭代次数、对变量重新做中心化。如果这些都不行检查一下你的第二层样本量。我个人的经验门槛是少于30个组就要非常谨慎地使用复杂随机效应少于20个组可能连随机截距模型都很难稳定。第三种是结果全是NaN或出现极端系数。这通常意味着某个变量有大量缺失值、某个变量几乎没有变异比如所有人分数都一样或者类别变量没有正确进行虚拟编码。还有一个特别容易被忽视的情况你的调节变量和自变量高度相关导致交互项和主效应之间出现多重共线性系数和标准误就会变得极不稳定。这时可以考虑把调节变量做总均值中心化再生成交互项。最后再说一个很多人忽略的细节跑MLmed时用的数据扫描、处理、清洗的每一个步骤都要在语法或操作日志里留痕。多水平模型的变量结构复杂中间任何一个环节出了偏差最直观的表现就是结果明显不合常理。我一直把先跑一遍空模型、再跑一遍只含直接效应的模型、最后再加入中介和调节当作固定流程每加一环就对比一次结果。这样既能防止模型设置错误也能在论文答辩或审稿问为什么你的结果长这样的时候拿出一个清晰的逻辑链。我自己的体会是MLmed这类插件并不是用来取代PROCESS的而是把中介调节分析从单层世界带进多层世界的桥。单水平数据PROCESS依然是利器一旦数据有嵌套结构MLmed能让你的结果经得住审稿人追问。跑多水平中介最耗时间的地方其实不是点鼠标而是前期把变量层级、中心化方式、随机效应结构想明白。这几件事想清楚了后面不管报错也好、迭代不收敛也好都知道往哪个方向去调整。希望这篇文章能帮你少走一点我当年走过的弯路。
