在绘制核密度图的过程中，曲线平滑与否和带宽关系较大，不同的带宽会带来不同的估计效果和图形呈现。本文以Stata中的核密度图为例展示不同带宽的绘图效果。
0、演示数据
以国家统计局[1]官方网站上关于2017-2020年我国31个省份人均GDP（未包含港澳台数据）为例。数据可在国家统计局官网获得，或在推文Stata绘图 | Stata绘制多期（年）核密度图中获得。
cd "C:\Users\Acer\Desktop" \\工作路径
import delimited "density.csv", encoding(GB18030) \\数据导入
* 对数化
gen ln2020 = ln(y2020)
gen ln2019 = ln(y2019)
gen ln2018 = ln(y2018)
gen ln2017 = ln(y2017)
1、单期核密度
以默认带宽和带宽分别为0.05、0.5、1为例，展示不同的带宽效果
kdensity ln2020, note("默认bwidth") //默认带宽
kdensity ln2020,bw(0.05) note("bwidth=0.05") //带宽为0.05
kdensity ln2020,bw(0.5) note("bwidth=0.5") //带宽为0.5
kdensity ln2020,bw(1) note("bwidth=1") //带宽为1
2、多期核密度
以默认带宽和带宽分别为0.05、0.5、1为例，展示不同的带宽效果
默认带宽
kdensity ln2017,addplot((kdensity ln2018) (kdensity ln2019) (kdensity ln2020)) ///
legend( label (1 "2017 kernel") label (2 "2018 kernel") label (3 "2019 kernel") ///
label (4 "2020 kernel") size(small) position(1) ///
row(4) ring(0) subtitle(年份)) xtitle("") note("默认bwidth")
bw = 0.05
kdensity ln2017,bw(0.05)addplot((kdensity ln2018, bw(0.05)) (kdensity ln2019, bw(0.05)) (kdensity ln2020, bw(0.05))) ///
legend( label (1 "2017 kernel") label (2 "2018 kernel") label (3 "2019 kernel") ///
label (4 "2020 kernel") size(small) position(1) ///
row(4) ring(0) subtitle(年份)) xtitle("") note("bwidth=0.05")
bw = 0.5
kdensity ln2017,bw(0.5)addplot((kdensity ln2018, bw(0.5)) (kdensity ln2019, bw(0.5)) (kdensity ln2020, bw(0.5))) ///
legend( label (1 "2017 kernel") label (2 "2018 kernel") label (3 "2019 kernel") ///
label (4 "2020 kernel") size(small) position(1) ///
row(4) ring(0) subtitle(年份)) xtitle("") note("bwidth=0.5")
bw = 1
kdensity ln2017,bw(1)addplot((kdensity ln2018, bw(1)) (kdensity ln2019, bw(1)) (kdensity ln2020, bw(1))) ///
legend( label (1 "2017 kernel") label (2 "2018 kernel") label (3 "2019 kernel") ///
label (4 "2020 kernel") size(small) position(1) ///
row(4) ring(0) subtitle(年份)) xtitle("") note("bwidth=1")
3、其他
关于Stata绘制核密度的更多调整可参考Stata绘图 | Stata绘制多期（年）核密度图或使用help kdensity命令查看帮助文档。关于如何利用R语言绘制核密度图可参考R语言绘图 | 多期核密度图绘制。
如有帮助请多多点赞哦！
参考资料
国家统计局: https://data.stats.gov.cn/index.htm
文章转载自日常分享的小懒猫，如果涉嫌侵权，请发送邮件至：contact@modb.pro进行举报，并提供相关证据，一经查实，墨天轮将立刻删除相关内容。
