查看原文
其他

Stata绘图:精美的散点图和线性拟合图-scatterfit-slopefit

连享会 连享会 2023-10-24

👇 连享会 · 推文导航 | www.lianxh.cn

连享会课程 · 2023 暑期班

温馨提示: 文中链接在微信中无法生效。请点击底部「阅读原文」。或直接长按/扫描如下二维码,直达原文:

作者:刘东 (中国农业大学)
E-mail:ld01@cau.edu.cn

Source:https://github.com/leojahrens/scatterfit


目录

  • 1.简介

  • 2.安装命令

  • 3. scatterfit 实例

  • 4.slopefit 实例



1.简介

scatterfit 可生成多种散点图并在其上叠加拟合线。它包含两个命令:scatterfit 用于可视化两个变量 x 和 y 之间的关系,slopefit 用于可视化在另一个连续变量 z 条件下 x 和 y 之间的关系。

2.安装命令

. ssc install scatterfit, replace
. help scatterfit
. help slopefit

要求: Stata version 13.1 及以上版本,并且需要安装 gtools, ftools, palettes, colrspace, blindschemes, labutil 命令。

. ssc install gtools
. ssc install ftools
. ssc install palettes
. ssc install colrspace
. ssc install blindschemes
. ssc install labutil

3. scatterfit 实例

调入数据

. sysuse "auto.dta", clear

画出表示 weight 和 length 之间的散点图

. scatterfit weight length

在散点图的拟合线上添加置信区间

. scatterfit weight length, fit(lfitci)

局部多项式拟合

. scatterfit weight length, fit(polyci) bw(10)

bwidth(num)表示用于控制局部多项式或 LOWESS 平滑拟合线平滑程度的带宽参数。

按照不同的分类 ( 是否为 foreign)展示 weight 和 length 之间的散点图和拟合线

. scatterfit weight length, by(foreign)

将 X 轴变量基于分位点切分成大小相等的区间,并在这些区间内绘制 Y 变量和 X 变量的均值

. scatterfit weight length, binned nq(20)

nq(num) 表示选择等大小的分组 / 分位数的数量。参数 num 用于指定要将变量分成多少个等大小的组 / 分位数。具体来说,该命令将变量的值分为 num 个等大小的分组 / 分位数,以进行后续的分析或可视化。

scatterfit 可以对数据进行预处理,以便在考虑控制变量后,仅显示 x 和 y 变量之间的残差协变关系。这是通过首先在控制变量上回归 x 和 y 变量,然后使用残差进行绘图来实现的

. scatterfit weight length, controls(trunk) ///
fcontrols(foreign) ///
regparameters(coef sig pval)

4.slopefit 实例

Slopefit 可以在一个连续变量 z 的条件下,可视化 x 和 y 之间的关系。它也可以创建带有拟合线的散点图。散点图和拟合线都显示了在不同的 z 值下 x 的斜率(即效应)。拟合线是从一个简单的交互模型中得出的,该模型假设 x 的效应随 z 线性变化。相比之下,散点图显示了在用户定义的 z 分组中(例如,z 的第一十分位数等)x 的效应。这些单独的斜率是从一个交互回归模型中得出的,该模型使用 z 的分解版本来获取不同子集的 z 的分别斜率系数。

有几种方法可以创建 z 的分段。与 scatterfit 一样,标准设置是根据分位数截断点将 z 分类为等大小的分段。

. slopefit weight length turn, ///
indslopes(quantiles) nquantiles(20)

此时,拟合线显示了连续交互模型的结果,最底部的散点标记显示的是在 z 分布的最低 5% 范围内 x 的效应。

scatterfit 类似, 也可以通过将 z 排序为等间距的分段来创建这些分段。 the bins can also be created by sorting z into equally spaced bins.

. slopefit weight length turn, ///
indslopes(unibin) nunibin(20)

或者将 z 本身视为离散变量, 并使用 z 的每个不同的值来创建分段。这意味着将 z 的每个不同取值单独作为一个分段,并在每个分段中计算 x 的效应。

. slopefit weight length turn, ///
indslopes(discrete)

也可以附加 indslopesci 选项,以绘制个别斜率的置信区间,对于每个不同的左右自我定位值,都可以进行绘制。indslopesci 的作用是计算和绘制个别斜率的置信区间。该选项通常用于多元回归分析中,用于检验每个自变量的效应是否显著。

. slopefit weight length turn, ///
indslopes(discrete) ci indslopesci

课程推荐:2023 暑期班
主讲老师:连玉君,王群勇
🍓 课程主页https://www.lianxh.cn/news/fdc69c3695aec.html

New! Stata 搜索神器:lianxhsongbl  GIF 动图介绍
搜: 推文、数据分享、期刊论文、重现代码 ……
👉 安装:
. ssc install lianxh
. ssc install songbl
👉  使用:
. lianxh DID 倍分法
. songbl all

🍏 关于我们

  • 连享会 ( www.lianxh.cn,推文列表) 由中山大学连玉君老师团队创办,定期分享实证分析经验。
  • 直通车: 👉【百度一下: 连享会】即可直达连享会主页。亦可进一步添加 「知乎」,「b 站」,「面板数据」,「公开课」 等关键词细化搜索。


您可能也对以下帖子感兴趣

文章有问题?点此查看未经处理的缓存