Stata可视化：biplot一图看尽方差、相关性和主成分
👇 连享会 · 推文导航 | www.lianxh.cn
🍎 Stata：Stata基础 | Stata绘图 | Stata程序 | Stata新命令 📘 论文：数据处理 | 结果输出 | 论文写作 | 数据分享 💹 计量：回归分析 | 交乘项-调节 | IV-GMM | 时间序列 | 面板数据 | 空间计量 | Probit-Logit | 分位数回归 ⛳ 专题：SFA-DEA | 生存分析 | 爬虫 | 机器学习 | 文本分析 🔃 因果：DID | RDD | 因果推断 | 合成控制法 | PSM-Matching 🔨 工具：工具软件 | Markdown | Python-R-Stata 🎧 课程：公开课-直播 | 计量专题 | 关于连享会
连享会寒假班
作者： 左祥太 (武汉纺织大学)
邮箱： Shutter_Z@outlook.com
目录
1. biplot 命令介绍
2. biplot 基本语法
3. Stata 实例
3.1 图形详解
3.2 图形美化
4. 结语
5. 参考资料
温馨提示： 文中链接在微信中无法生效。请点击底部「阅读原文」。或直接长按/扫描如下二维码，直达原文：
1. biplot 命令介绍
通过使用 biplot 可以将一份数据集的分布特征生动地展现在一副二维图上。biplot 将会同时显示观测值和变量的相对位置。标记符号（点）会显示指定的标签，而变量则用箭头表示。所有特征的观测值将被投射到两个维度上，从而使观测值之间的距离近似保留（压缩至两个维度的主成分）。箭头之间角度的余弦近似于变量之间的相关性，箭头的长短近似于表示该变量的方差。
2. biplot 基本语法
biplot 语法中的 [, options] 划分了主要选项、行选项与列选项，本文将着重介绍几个最基本的实用项，其余选项以及介绍可以通过 help biplot 查看。
biplot varlist [if] [in] [, options]
rowlabel: 用于指定图中点的标签 alpha: 默认值为 0.5，此时绘制对称比例双平面图；取 1 时绘制主成分双平面图 (也称 row-preserving biplot) stretch: 伸展列 (变量) 箭头，可简单记忆为调整图中箭头的长短 table: 用于设定是否需要返回图中点的坐标 rowover: 是否需要识别来自不同组别的变量 colopts(name()): 更改图例中箭头的名称 rowopts(name()): 更改图例中点的名称 title: 设置图片标题 autoaspect: 根据数据内容自动调整图片长宽比，默认为1
3. Stata 实例
3.1 图形详解
在使用 biplot 绘图前需要下载安装该命令至本地：
ssc install biplot
接下来调用该命令所提供的数据集进行绘图操作，为了方便后期图片中颜色的显示，我们采用连享会上期推文中的 schemes 对图片风格进行设置，具体参考：Stata 绘图：用 Stata 绘制一打精美图片-schemes| 连享会主页 (lianxh.cn) 。
首先我们调用命令集所提供的数据文件，并观察数据的基本结构与变量名称。
. webuse renpainters //调用数据
clear
input str15 painter byte(composition drawing colour expression)
"Del Sarto" 12 16 9 8
"Del Piombo" 8 13 16 7
"Da Udine" 10 8 16 3
"Giulio Romano" 15 16 4 14
"Da Vinci" 15 16 4 14
"Michelangelo" 8 17 4 8
"Fr. Penni" 0 15 8 0
"Perino del Vaga" 15 16 7 6
"Perugino" 4 12 10 4
"Raphael" 17 18 12 18
end
然后我们设置一副白色的背景风格，方便后续颜色的展示，最后通过 biplot 绘制一幅最基本的图片。
. set scheme white_tableau //设定绘图风格为white_tableau
//本推文的输出风格
. set scheme s2color // Stata 默认
. biplot composition-expression //简单绘制图片
/* 以下为 Stata 窗口显示内容 */
Biplot of 10 observations and 4 variables
Explained variance by component 1 0.6700
Explained variance by component 2 0.2375
Total explained variance 0.9075
Stata 的窗口展示了每个组成维度的解释性方差，不难发现降维后的两个维度对原样本的方差解释力度达到了 90.75% ，即在一个可以接受的范围内还原了原始数据 90.75 %的数据特征。
图形解读如下：
点： 上图中两点之间的距离近似于真实的多变量空间中两个观测点之间的欧氏距离。以第 4 和第 5 条样本为例，在原始数据中两个样本的取值一样，根据四维空间中欧氏距离的计算公式有:
可得两点之间在四维空间中的距离为 0 ，所以压缩至二维空间，在 bioplot 图中表现为重合状态，观察上图可知样本 4 和样本 5 重合。
箭头： 上图中线条的长度近似于变量的方差，即线条越长，方差越大。可以看到 colour 在所有变量中的方差最大，而 drawing 的方差最小。 夹角： 线条之间角度的余弦值，可以近似于它们所对应变量之间的关联性。角度越接近 90 度 () 或 270 度 () ，相关度就越小。同理， 0 度或 180 度的角度分别反映了 1 或 -1 的关联性。上图展现了 composition 和 expression 之间的强正相关关系， colour 和 drawing 之间的强负相关关系，而 composition 与 colour 以及 drawing 之间则表现为弱相关关系。
3.2 图形美化
为了方便对点标签的识别，引入 rowlabel 选项将变量 painter 中的文本作为上图中点的标签，同时通过 rowopts(name()) 和 colopts(name()) 分别改变图例中点和箭头的名称，并通过 title 选项为图加上标题。
biplot composition-expression, rowlabel(painter) ///
rowopts(name(Painters)) colopts(name(Attributes)) ///
title(Renaissance painters)
接下来通过 alpha 选项和 stretch 选项使得图片看起来更加舒适。
biplot composition-expression, rowlabel(painter) ///
rowopts(name(Painters)) colopts(name(Attributes)) ///
title(Renaissance painters) ///
alpha(1) stretch(10)
同时我们可以利用 table 选项来返回图中各点的坐标值。
biplot composition-expression, rowlabel(painter) ///
rowopts(name(Painters)) colopts(name(Attributes)) ///
title(Renaissance painters) ///
alpha(1) stretch(10) table
/* 以下为 Stata 窗口显示内容 */
Biplot coordinates
Painters | dim1 dim2
-------------+--------------------
Del Sarto | 1.2120 0.0739
Del Piombo | -4.5003 5.7309
Da Udine | -7.2024 7.5745
Giulio Rom~o | 8.4631 -2.5503
Da Vinci | 8.4631 -2.5503
Michelangelo | 0.1284 -5.9578
Fr Penni | -11.9449 -5.4510
Perino del~a | 2.2564 -0.9193
Perugino | -7.8886 -0.8757
Raphael | 11.0131 4.9251
----------------------------------
Attributes | dim1 dim2
-------------+--------------------
composition | 6.4025 3.3319
drawing | 2.4952 -3.3422
colour | -2.4557 8.7294
expression | 6.8375 1.2348
----------------------------------
最后，不知你是否觉得这个图片的风格很清爽，完全不同于 Stata 默认的粗笨风格。若想实现这个风格，只需在绘图之前执行如下命令，将绘图模板设定为 white_tableau 风格即可。
有关绘图模板的设定方法，参见往期推文
Stata绘图：一个干净整洁的-Stata-图形模板qlean Stata黑白图形模板：中文期刊风格的纯黑白图形
. ssc install schemepack, replace //安装white_tableau 模板
. set scheme white_tableau //设定绘图风格为white_tableau
若需重设为 Stata 默认绘图风格，可以按需执行如下两条命令之一：
. set scheme s2color // Stata 默认, 彩色
. set scheme s2mono // Stata 默认, 黑白
4. 结语
biplot 命令的出现很好地帮助我们在一张图内同时展现了变量之间的近似相关性、主成分等内容。与此同时，biplot 与 schemes 所提供的绘图背景相结合有助于我们画出更加惊艳的图片。
5. 参考资料
Ulrich Kohler, Magdalena Luniak, 2005, Data Inspection using Biplots, Stata Journal, 5(2): 208–223. -PDF- Stata manual, [MV] biplot GGE Biplot 双标图如何看？
连享会：因果推断专题
New！ Stata 搜索神器：
lianxh和songblGIF 动图介绍
搜： 推文、数据分享、期刊论文、重现代码 ……
👉 安装：
. ssc install lianxh
. ssc install songbl
👉 使用：
. lianxh DID 倍分法
. songbl all
🍏 关于我们
连享会 ( www.lianxh.cn，推文列表) 由中山大学连玉君老师团队创办，定期分享实证分析经验。 直通车： 👉【百度一下：连享会】即可直达连享会主页。亦可进一步添加 「知乎」,「b 站」,「面板数据」,「公开课」 等关键词细化搜索。