复现丨带置信区间的调节效应图
不知道大家在阅读文献时有没有感觉，现在的图越来越卷了。
最新录用中，很多调节效应图都有置信区间。
秉承着打不过就加入的原则，萜妹前阵子去研究了下如何画出这种图。
正好看到了篇图画的很漂亮的研究，那我们就一起来复现吧。
Madan, S., Savani, K., & Johar, G. V. (2022). How you look is who you are: The appearance reveals character lay theory increases support for facial profiling. Journal of Personality and Social Psychology. https://doi.org/10.1037/pspa0000307
整体介绍
这是一篇偏心理学的文章，主要结论是：人们越相信外貌可以揭示性格，他们对基于外貌的判断就越有信心，因此他们越支持在执法、教育和商业中使用面部分析技术。
作者做了十个研究来佐证它。但考虑到篇幅问题（主要是写长了你们也不看！），我们这次只复现它的Study 5a，即涉及调节效应和带置信区间的调节效应图的部分。
公开数据
数据及代码：https://osf.io/aj86v
补充材料：https://doi.org/10.1037/pspa0000307.supp
作者的公开材料里没有作图部分的代码，萜妹自学失败后，发邮件问了作者，她超大方的分享了（大家真的都好nice）。我也会在接下来的示例中展示。
研究概述
Study 5a就是一个简简单单的调节效应检验。
作者通过李克特量表获得了被试在自变量、调节变量的得分。
然后向被试描述了四种使用面部分析技术的政策及其利弊，并要求他们评价在多大程度上支持使用面部分析技术。
数据分析
假设检验
作者使用的是Stata软件，公开文件里给出了stata语句与结果。
我们只需要打开Stata，将代码复制进Stata的命令窗口即可。
regress supportavg laycen entcen entcenXlaycen if gibberish==0
supportavg是因变量，支持程度的平均值。 laycen是中心化后的自变量，是先算了各题项平均值再进行的中心化。 entcen是中心化后的调节变量，也是先算了均值再中心化。
按【Enter】，软件就能运行出结果。
文中的汇报也比较简单，不过我没有找到β的来源。
简单斜率分析
作者提前计算出了两个变量，分别为调节变量的均值加1个标准差（eplus1sd）和减1个标准差（eminus1sd）。然后调整了下上面的语句，又进行了两次回归。
regress supportavg laycen eplus1sd eplus1sdXlaycen if gibberish==0
regress supportavg laycen eminus1sd eminus1sdXlaycen if gibberish==0
不过要注意，此时是简单斜率分析，所以我们要看主效应，而不再是调节效应。
PS：软件中结果与文中汇报是反的。
带置信区间的调节效应图
接下来是我们本次的重点啦！
重新进行回归
首先要重新进行调节效应的回归，需要告诉软件，我们这里的调节变量和自变量是连续变量。如果是类别变量则改为i.variable。
regress supportavg c.entavg##c.layavg if gibberish==0
这个回归的结果会在后续作图时调用。需要注意：这里的结果和我们前面的回归结果是有一些出入的。
调节变量的描述性统计
然后对调节变量进行描述性统计。
summ entavg
计算可得调节变量的均值和标准差。
作图
首先要利用下面的语句，得到作图的数据。entavg括号内的数即为均值加减一个标准差，
margins, at(layavg =(1 2 3 4 5 6) entavg =(1.94 3.28 4.62)) vsquish
这里的结果不重要，所以就不截图上来啦。
然后我们再进行作图。
marginsplot, recastci(rarea) x(layavg ) recast(line) xlabel(1.94 3.28 4.62)
这里是最基础的图形结果。我们可以在Stata的图形编辑里改变样式，或者也可以利用语句来优化。不过我觉得图形编辑好用一些。
作者原文中的图形长这样：
今天的推送就到这里啦。这次的复现其实非常简单，但就是这么简单的内容，我自己摸索了一个月都没摸索出来，所以我觉得还是非常有必要和小可爱们分享的。
如果小可爱们知道其他能在调节效应图上添加置信区间的方法也欢迎与我交流呀。希望能和小可爱们，一起学习，共同进步~
另外，如果小可爱们对Stata作图有兴趣也可以留言，虽然我还在摸索之中，但是如果大家有需求，我会学得快一点吧，hhh~
下期预告：《范文丨顶刊引言精读（四）》
往期推送：