中介效应模型
一、什么是中介效应模型？
https://www.zhihu.com/question/599491767?utm_id=0
考虑自变量X对因变量Y的影响, 如果X通过影响变量M而对Y产生影响, 则称M为中介变量。路径图如下图所示：
①方程(1)的系数c为自变量X对因变量Y的总效应;
②方程(2)的系数a为自变量X对中介变量M的效应;
③方程(3)的系数b是在控制了自变量X的影响后, 中介变量M对因变量Y的效应;
④系数c’是在控制了中介变量M的影响后, 自变量X对因变量Y的直接效应;
⑤e1~e3是回归残差。
对于这样的简单中介模型 , 中介效应等于间接效应即等于系数乘积ab, 它与总效应和直接效应有方程（4）这样的关系。
检验中介效应最常用的方法是逐步检验回归系数，即通常说的逐步法:
(i)检验方程(1)的系数c (即检验H0 : c = 0,检验的是X对Y的总效应);
(ii)依次检验方程(2)的系数a (即检验H0 : a = 0) 和方程(3)的系数b (即检验H0 : b = 0)，第二步实际上是检验系数乘积的显著性；
如果(i)系数c显著, (ii)系数a和b都显著, 则中介效应显著。
(iii)如果方程(3)的系数c’不显著，则是完全中介效应（第三步检验用来区分完全中介还是部分中介）。
具体检验步骤如下：
1.逐步法
①quietly reg z_P z_T z_age z_job z_M_salary z_sex z_mar z_degree z_Party z_Q1 z_Q2 z_Q3
outreg2 using abc,word
z_P为被解释变量， z_T 为解释变量，后面的一串为控制变量
②quietly reg z_R1 z_T z_age z_job z_M_salary z_sex z_mar z_degree z_Party z_Q1 z_Q2 z_Q3
outreg2 using abc,word
z_R1为中介变量， z_T为解释变量，后面一串为控制变量
③quietly reg z_P z_R1 z_T z_age z_job z_M_salary z_sex z_mar z_degree z_Party z_Q1 z_Q2 z_Q3
outreg2 using abc,word
z_P为被解释变量， z_R1为中介变量， z_T为解释变量，后面一串为控制变量
stata检验回归一般用reg，以上代码可以另外输出word形式，方便论文使用。
2.sobel检验法
Sobel法是检验系数乘积法，其检验力高于依次检验，但这个检验统计量的推导需要假设a的估计*b的估计服从正态分布,具有局限性。其代码如下：net install sgmediation2, from (“https://tdmize.github.io/data/sgmediation2”)
sgmediation z_P, mv(z_R1) iv(z_T) cv(z_age z_job z_M_salary z_sex z_mar z_degree z_Party)
z_P为被解释变量，mv( )括号里为中介变量，iv( )括号里为解释变量，cv（ ）括号里为控制变量。
如何看sobel检验结果：看Goodman-1(Aroian)那一行的p值，如果p值<0.05，则中介效应存在。如果中介效应存在，那么具体的直接效应值（direct effect）、间接效应值(indirect effect)和总效应值(total effect)可以看下面那个表的Coef.值。
3.bootstrap检验法（直接检验系数乘积的显著性）
Bootstrap法是一种从样本中重复取样的方法, 前提条件是样本能够代表总体(当然这也是通常取样进行统计推论的要求)。检验力高于Sobel检验。代码如下：
bootstrap r(ind_eff) r(dir_eff), reps(500) : sgmediation z_P, mv(z_R1) iv(z_T) cv(z_age z_job z_M_salary z_sex z_mar z_degree z_Party)
estat bootstrap, percentile bc
同上，z_P为被解释变量，mv( )括号里为中介变量，iv( )括号里为解释变量，cv（ ）括号里为控制变量。reps(500) 表示重复抽样500次。
bootstrap检验法如何看结果：不是看P值，是看[BootLLCl , BootULCl]这个置信区间，如果这个区间包含0，则不显著，如果不包含0，则中介效应显著。
综上，检验中介效应可以先进行依次检验, 如果不显著就使用效力更高的Bootstrap检验法，总体检验流程如下：
①检验方程(1)的系数c, 如果显著, 按中介效应立论, 否则按遮掩效应立论。但无论是否显著, 都进行后续检验。
②依次检验方程(2)的系数a和方程(3) 的系数b, 如果两个都显著, 则间接效应显著, 转到第四步; 如果至少有一个不显著, 进行第三步。
③用Bootstrap法直接检验H0 : ab = 0。如果显著, 则间接效应显著, 进行第四步; 否则间接效应不显著, 停止分析。
④检验方程(3)的系数c’,如果不显著, 即直接效应不显著, 说明只有中介效应。如果显著, 即直接效应显著, 进行第五步。
⑤比较ab和c’的符号, 如果同号, 属于部分中介效应, 报告中介效应占总效应的比例ab/c。如果异号, 属于遮掩效应, 报告间接效应与直接效应的比例的绝对值|ab/c’|。
注：①用stata做中介效应检验还需要一个外部命令sgmediation.ado文件，如有需要，可分享。
②如果在做回归时，想要显著性水平上面有星号，需要安装一个pwcorr_a外部命令，即把解压好的pwcorr_a文件中的pwcorr_a.ado和pwcorr_a.help两个文件放到stata/ado/base/p目录下，如有需要，可分享。
内生性