• *----------------------------------------------* *==== stata绘制人口金字塔 =====* *---------------------------------------------* sysuse pop2000,clear replace maletotal = -maletotal/1000000 ...
*----------------------------------------------*
*====         stata绘制人口金字塔         =====*
*---------------------------------------------*
sysuse pop2000,clear
replace maletotal = -maletotal/1000000
replace femtotal = femtotal/1000000
generate zero = 0
#delimit ;
twoway (bar maletotal agegrp,horizontal)||(bar femtotal agegrp,horizontal)  //horizontal是使条形图水平放置
(scatter agegrp zero,msymbol(i) mlabel(agegrp) mlabcolor(black)),   //添加年龄组的标签
xlabel(-12 "12" -8 "8" -4 "4" 4 8 12) xtitle("Population in millions",height(0))  //给x轴添加标题，修改标签，height(0)是把标签和标题放在一行
yscale(off) ylabel(,nogrid) legend(order(1 "Male" 2 "Female"))  //去掉y轴，去掉网格线，完善图例
title("US Male and Female Population by Age,2000")
note("Source:U.S. Census Bureau,Census 2000,Tables 1,2 and 3")
;
#delimit cr

*----------------------------------------------*
*====             stata绘制饼图           =====*
*---------------------------------------------*
webuse census,clear
#delimit ;
graph pie pop,over(region) pie(2,explode) plabel(2 percent,gap(-10))  //突出显示顺时针第二个扇区，并标注第二个扇区的百分比，可以将percent改为sum
plabel(_all name,size(*1.5) color(white))  //显示每个扇区名称
legend(off)
plotregion(lstyle(none))
title("Polulation by region")
subtitle("1980")
note("Source:1980 Census data by state")
;
#delimit cr



展开全文
• stata绘图
• 应用stata学习计量经济学原理 Practice 3 stata绘图进阶 Practice 3 stata绘图进阶 log using homework3 ,replace text * Q1 sysuse uslifeexp.dta , clear * 1 twoway line le year , title(1900-2000 US life ...
应用stata学习计量经济学原理 Practice 3 stata绘图进阶
Practice 3 stata绘图进阶

log using homework3 ,replace text
* Q1
sysuse uslifeexp.dta , clear
* 1
twoway line le year , title(1900-2000 US life expectancy)
* 2
generate diff = le_male - le_female
label variable diff "male-female"

twoway (lfitci diff year,yaxis(1) level(95)) (line diff year,yaxis(1)) (line le_male year,yaxis(2)) (line le_female year,yaxis(2)), ///
xlabel(1900(10)2000) ///
title(1900-2000 US life expectancy)  ///
note("Source: U.S. life expectancy, 1900-1999")  ///
xtitle(year) ytitle(life expectancy,axis(2)) ytitle(diff between male and female,axis(1))

* Q2
use wage.dta,clear
*1
qnorm wage

sktest wage
swilk wage
sfrancia wage
sktestdc wage

* 不满足正态分布

*2
generate lnw = log(wage)

qnorm lnw

sktest lnw
swilk lnw
sfrancia lnw
sktestdc lnw

* 好像还是不满足

log close




展开全文
•   作者：王月琴 (中国农科院 ) Stata 连享会： 知乎 | 简书 | 码云 ...Stata连享会 精品专题 || 精彩推文 ...如果用 Excel 绘制，似乎并不复杂，但后续修图比较麻烦，且...不过，经过此番折腾，对于 Stata 绘图命...


🍎 连享会     主页  || 视频 || 推文
扫码查看最新推文和分享

NEW！连享会·推文专辑：
Stata资源 | 数据处理 | Stata绘图 | Stata程序
结果输出 | 回归分析 | 时间序列 | 面板数据 | 离散数据
交乘调节 | DID | RDD  |  因果推断 |  SFA-TFP-DEA
文本分析+爬虫 | 空间计量 | 学术论文 | 软件工具

连享会 - Stata 暑期班
线上直播 9 天：2020.7.28-8.7
主讲嘉宾：连玉君 (中山大学) | 江艇 (中国人民大学)
课程主页：https://gitee.com/arlionn/PX | 微信版

作者：王月琴 (中国农科院 )
Stata 连享会： 知乎 | 简书 | 码云

Stata连享会 精品专题  || 精彩推文

任务
想绘制一幅多维柱状图。如果用 Excel 绘制，似乎并不复杂，但后续修图比较麻烦，且绘图工作不具有可重复性。
为此，想用 Stata 绘制这幅图形。由于对 Stata 绘图不熟悉，前后折腾了好几个小时。不过，经过此番折腾，对于 Stata 绘图命令的语法特征，各个要素的选项设定有了比较清晰的了解。

Stata 绘图过程
原始数据存储于 Excel 表格中，可以用 import excel 命令导入。这里，为了便于展示，使用更原始的 input 命令导入：
clear
input year	str30 region level
2008	"China"	5.368679
2009	"China"	10.27191
2010	"China"	11.32058
2011	"China"	11.88313
2012	"China"	18.64346
2013	"China"	21.03853
2014	"China"	14.79839
2015	"China"	20.26556
2016	"China"	24.42027
2017	"China"	25.32076
2018	"China"	26.54098
2008	"Manitoba"	40.45824
2009	"Manitoba"	59.37749
2010	"Manitoba"	60.37274
2011	"Manitoba"	63.62734
2012	"Manitoba"	52.39755
2013	"Manitoba"	61.7517
2014	"Manitoba"	70.32643
2015	"Manitoba"	65.98277
2016	"Manitoba"	61.32922
2017	"Manitoba"	48.68953
2018	"Manitoba"	52.63895
2008	"United States of America"	52.33519
2009	"United States of America"	93.42258
2010	"United States of America"	51.09012
2011	"United States of America"	72.44342
2012	"United States of America"	61.02513
2013	"United States of America"	80.46352
2014	"United States of America"	77.81561
2015	"United States of America"	82.49318
2016	"United States of America"	73.81653
2017	"United States of America"	72.35381
2018	"United States of America"	74.25181
end
* save  "Fig_Data.dta", replace  // 保存数据

绘图语句如下。几点说明：

由于语句比较长，故用 #delimit ; 将 Stata 的换行方式由默认的回车键更改为 ; （#delimit ; 可以简写为 #d ;）；待绘图命令结束后，用 #delimit cr 恢复默认设定。注意： 绘图命令的最后一行只有一个结束符 ;。
命令中涉及各个选项的填写，均已做了简要说明，详情参见 help graph bar

* use  "Fig_Data.dta",  clear
#delimit ;
graph bar (asis) level
,  // 后面都是选项设定
over(region, gap(0) label(nolabel)) //首先按照地区分类, 无横坐标标签
over(year,label(labsize(small)))    //其次按照年份分类，横坐标标签字体小号
asyvars                             //按照地区在y轴分类
ylabel(0(20)100, tposition(inside) labsize(*0.8) angle(0))
// 纵轴刻度和标签设定
// 纵坐标从0-100，刻度线朝内, y轴标签字号为默认字号的 0.8 倍
// 标签为纵向放置 (默认是横向)
ytitle("%", tstyle(smbody)) // y 轴标题为小号字体
blabel(bar, size(vsmall) format(%3.1f))
//产生数字标签，字体小号，固定格式保留1位小数
bar(1, color(red))
bar(2, color(green))
bar(3, color(orange)) //中国条形图为红色，加拿大为绿色，美国为橙色
legend(row(1) ring(1) position(6)
label(1 "中国")
label(2 "加拿大（曼）")
label(3 "美国")
order(1 2 3)
size(small))
//图例: 呈一列放在图表外6点钟方向，标签按照地区排列，图例字体小号
graphregion(color(white)) //图片底色为白色
;
#d ;
. graph export Figure.png, replace //保存图片

输出 Stata 图形如下：

关于我们

Stata连享会 由中山大学连玉君老师团队创办，定期分享实证分析经验。直播间 有很多视频课程，可以随时观看。
你的颈椎还好吗？ 您将 ::连享会-主页:: 和 ::连享会-知乎专栏:: 收藏起来，以便随时在电脑上查看往期推文。
公众号推文分类： 计量专题 | 分类推文 | 资源工具。推文分成  内生性 | 空间计量 | 时序面板 | 结果输出 | 交乘调节 五类，主流方法介绍一目了然：DID, RDD, IV, GMM, FE, Probit 等。
公众号关键词搜索/回复 功能已经上线。大家可以在公众号左下角点击键盘图标，输入简要关键词，以便快速呈现历史推文，获取工具软件和数据下载。常见关键词：

课程, 直播, 视频, 客服, 模型设定, 研究设计, 暑期班
stata, plus，Profile, 手册, SJ, 外部命令, profile, mata, 绘图, 编程, 数据, 可视化
DID，RDD, PSM，IV，DID, DDD, 合成控制法，内生性, 事件研究
交乘, 平方项, 缺失值, 离群值, 缩尾, R2, 乱码, 结果
Probit, Logit, tobit, MLE, GMM, DEA, Bootstrap, bs, MC, TFP, 面板, 直击面板数据, 动态面板, VAR, 生存分析, 分位数
空间, 空间计量, 连老师, 直播, 爬虫, 文本, 正则, python
Markdown, Markdown幻灯片, marp, 工具, 软件, Sai2, gInk, Annotator, 手写批注, 盈余管理, 特斯拉, 甲壳虫, 论文重现, 易懂教程, 码云, 教程, 知乎


展开全文
• 这是Stata绘图命令twoway area的第三弹。我们已经在之前的文章中铺垫了很多关于twoway area的用法。它看起来是那么酷炫，但似乎毫无用处。这一篇我们将使用twoway area命令绘制地图，同时解密Stata中用于绘制地图...
这是Stata绘图命令twoway area的第三弹。我们已经在之前的文章中铺垫了很多关于twoway area的用法。它看起来是那么酷炫，但似乎毫无用处。这一篇我们将使用twoway area命令绘制地图，同时解密Stata中用于绘制地图命令——spmap——的原理。在开始之前，你可参考一下文章：张大神气：Stata绘图之twoway area（一）：绘制芒星​zhuanlan.zhihu.com张大神气：Stata绘图之twoway area（二）：绘制国旗​zhuanlan.zhihu.com1. 地图坐标数据我们在之前的推文中已经介绍过了Stata中地图的绘制过程，相信大家已有所涉猎。在Stata中绘制地图，需要两份数据：一是地图的坐标数据。它用来标定地图的各个地区的坐标或者经纬度。二是你所要要绘制的一些数据，比如各个省份的人口密度、各个城市的平均气温等。这两份数据通过某一个变量（如_ID）关联。从而使人们可以方便地绘制出地图。我已经将一份中国各省份地图的坐标数据上传至我的github上了，通过设置以下参数，可以从Stata中直接调用该数据。// 设置参数
global github "https://raw.githubusercontent.com/zhangdashenqi"
webuse set "\${github}/the_zen_of_stata/master/data"

//载入数据
webuse china_map.dta, clear
br这份坐标数据就是这个摸样：想一下，与我们之前绘制芒星和国旗时的数据是不是一个样？！常规操作，绘制散点：scatter _Y _X, msize(vtiny)2. 绘制一个省在这份坐标数据中，_ID = 1对应为黑龙江省，如果你看过之前的推文，相信可以对于下面的操作并不陌生：twoway (area _Y _X if _ID==1, nodropbase)这个就是我们使用twoway area绘制一个省份，与之前绘制芒星和国旗并无不同。一切都是水到渠成！结果如下：那么相对应地，绘制多个省份的代码如下:#delimit ;
twoway (area _Y _X if _ID==1, nodropbase)
(area _Y _X if _ID==2, nodropbase)
(area _Y _X if _ID==3, nodropbase)
(area _Y _X if _ID==4, nodropbase)
(area _Y _X if _ID==5, nodropbase)
(area _Y _X if _ID==6, nodropbase)
(area _Y _X if _ID==7, nodropbase)
(area _Y _X if _ID==8, nodropbase)
(area _Y _X if _ID==9, nodropbase)
;
#delimit cr结果如下：3. 绘制一份完整的地图在这份坐标数据中，_ID从1到66，你当然可以将上述代码重复66行，但是为了简洁美观起见，我们可以采用循环语句：local cmd ""
forvalues i = 1(1)66{
local cmd0 "(area _Y _X if _ID==i', nodropbase)"
local cmd cmd' cmd0'
}
twoway cmd', legend(off)第2~5行，我们使用循环语句将那些重复的命令全部放到cmd这个局部暂元中。 第6行，绘制出全部地图，并且不显示图例。结果如下：可以看出，这已经是一个完整的地图了。只不过每个省份的颜色比较随机，没有具体含义。事实上，Stata绘制地图的命令spmap就是基于twoway area命令进行封装的，感兴趣的同学可以学习一下它的源码。spmap的作者Maurizio Pisati在代码中花费了大量的篇幅来定义地图上填充的颜色，以使不同的颜色表征不同范围的数值。如果你有自定义地图的需要，完全可以使用twoway area绘制一幅你想要的地图！
展开全文
• ## 一文看尽 Stata绘图

万次阅读 多人点赞 2019-07-17 23:27:07
作者：万莉 (北京航空航天大学) ...基于由浅入深的原则，本推文结合 Stata 范例，梳理了 Stata绘图命令，便于读者掌握其中规律。用图说话，就从这里开始！ 1. 绘图简介 Stata 提供各种绘图命令，包括点...
• 应用stata学习计量经济学原理 Practice 2 Practice 2 * Q1 clear all * 1 sysuse uslifeexp.dta * 2 scatter le year scatter le year,connect(1) scatter le year,connect(1) msymbol(i) * 3 scatter le year, /...
• 作者：连玉君 (知乎 | 简书 | 码云) 连享会   最新专题   直播简介有些中文期刊 (如《统计研究》) 要求图片中的网格线为虚线，而不是 Stata 默认的浅灰色实线。我们可以使用 yline() 选项自行修改。Stata 范例代码...
• 1. Stata绘制双轴图之前我们介绍过很多Stata绘图的小技巧（没看过的可以翻一翻公众号的历史记录），今天我们介绍如何将不同x轴的数据绘制到一幅图中。先看数据：这是人民币兑美元的日度数据：这是中国对美出口的年度...
• 2019连享会暑期“实证研究方法与经典论文”专题班编者按：画图时，Stata 在默认情况下的图形输出通常为浅蓝色带横向网格底纹的图形，且不同输出图形通常以差异化的颜色加以标注，但很多时候，期刊排版要求图片为可...
• 作者：连玉君 (知乎 | 简书 | 码云)连享会   最新专题   直播连享会-知乎推文列表连享会 - Stata 暑期班：直播 (不必舟车劳顿了) 时间： 2020.7.28-8.7 主讲嘉宾：连玉君 (中山大学) | 江艇 (中国人民大学) 课程...
• 郑连虎，在数学学院取得理学学位的文科生，中国人民大学硕博连读生在读，山东大学管理学学士、理学学士， R语言中文社区专栏作者个人公众号：阿虎定量笔记往期回顾：R语言绘图：28个实用程序包全文注释* 初识Stata，...
• 不过，幸运的是，Stata 允许用户自定义绘图主题，本文尽可能地收集了这些第三方绘图主题，合计 82 个。最后我使用 putdocx + 循环创建了一份文档展示这些 Stata 主题。安装 tssctssc 命令是我编写的一个小命令，能够...