Stata U型关系检验、调节效应与绘图案例

'
///////////////////////////////////////////////////////////////////////////////////////
                *U型/倒U型关系的验证及调节效应分析绘图Stata实现代码*
///////////////////////////////////////////////////////////////////////////////////////


//导入原始数据
use 案例数据,replace
//为了减少极端值对模型估计的影响，我们进行缩尾处理
winsor2 Size Lev Growth, replace cuts(1 99)  

///案例说明：我们假设我们要考察资产负债率Lev和企业规模Size之间是否存在非线性关系，随着资产负债率的增加企业的规模是上升还是下降呢？是会一直上升还是会在一个节点转变趋势呢？我们带着这个问题进行接下来的讨论。

///////////////////////////////////////////////////////////////////////////////////////
//第一步：我们需要先验地认为资产负债率Lev和企业规模Size这两个变量之间存在非线性关系，我们建立非线性模型，并通过回归等操作来证明该模型的合理性即可证明资产负债率Lev和企业规模Size这两个变量之间的非线性关系确实存在。/////////////////////////////////////////////////
///////////////////////////////////////////////////////////////////////////////////////

//我们的非线性模型是这样的： Size=a+bLev+cLev^2+dControl+residual   

//然后我们可以使用twoway生成Lev与Size之间的散点图和拟合曲线并初步观察下Lev与Size之间是否存在非线性关系，scatter绘制散点图，qfit绘制拟合曲线（蓝色线）。

twoway (scatter Size Lev) (qfit Size Lev, lcolor(red))
/////////////////////////////////////////////////////////////////////////////////////////第二步：刚才第一步中绘制的拟合曲线与散点图如果能初步看出Size与Lev之间可能存在U型关系，那么我们就可以进行下一步回归分析，并计算U型曲线的转折点的数值。//////////////////////
///////////////////////////////////////////////////////////////////////////////////////

//1.回归拟合

gen Lev2= Lev^2 //生成自变量的平方项
reg Size Lev Lev2 
  //Lev2为Lev的平方项，回归中的控制变量可以加也可以不加，初步判断Lev2的系数显著性即可，显著即存在为U型的可能，可继续进行计算转折点。