clear

set obs 100

gen x = cond(_n > 50,_n*runiform(2,3), _n*runiform(4,5)) //其中第50是数据的断点位置

gen y= x^2+x+10+runiform(1,40000)

global b=x[50]

首先使用rd自带的绘图命令

rd y x, z0($b) mbw(100) graph

紧接着来观察cmogram命令绘图与rd绘图的不同

cmogram y x if x > $b-20 & x <$b+20, title(断点回归) cut($b) scatter lineat($b) lfit histopts(bin(15))

如果需要加上置信区间只需要将lift改为lfitci即可，如：

cmogram y x if x > $b-20 & x <$b+20, title(断点回归) cut($b) scatter lineat($b)  lfitci histopts(bin(15))

通过命令

cmogram y x if x > $b-20 & x <$b+20, title(断点回归) cut($b) scatter lineat($b) qfit histopts(bin(15))

同样如果需要加上置信区间只需要将qfit改为qfitci即可，如：

cmogram y x if x > $b-20 & x <$b+20, title(断点回归) cut($b) scatter lineat($b) qfitci histopts(bin(15))

以上图形均是定义窗宽为40，其中cut()是断点位置，scatter lineat()是在断点位置加上位置线，相当于scatter绘图中的xline()函数的作用，而histopts(bin(15))在之前推文中也提过是构造十五个箱体，希望通过这几个图的展示我们可以看出cmogram绘图对断点回归显著的展示作用，同时也希望该推文能够让你们对cmograph命令有更深刻的了解。如果大家对cmogram感兴趣可以自行help cmogram,里面也有很多有趣的examples,让你们熟练掌握该命令。

接下来是空气质量报告

全国空气质量如下

甘肃的朋友们请保重

以上就是今天给大家分享的内容了，说得好就赏个铜板呗！有钱的捧个钱场，有人的捧个人场~，点赞打赏随您心意，么么哒~

应广大粉丝要求，爬虫俱乐部的推文公众号打赏功能可以开发票啦，累计打赏超过1000元我们即可给您开具发票，发票类别为“咨询费”。用心做事，只为做您更贴心的小爬虫。第一批发票已经寄到各位小主的手中，大家快来给小爬虫打赏呀~

文字编辑：梅洁瓷傲

技术总编：刘贝贝

往期推文推荐：

1.合并输出回归结果和其他检验结果——esttab和estadd

2.关于RTF你不知道的命令

3.关于RTF你不知道的命令之番外篇

4.免费事件研究，一片片从邮局寄来

5.免费的股价同步性，一片片从邮局寄来

6.Stata叫你回家听音乐了！

7.一言不合就用stata写邮件（Outlook/Foxmail）

8.玩转stata之调用浏览器

9.I have a Stata, I have a python

10.I have a Stata, I have a Python之二——pdf转word

此外，欢迎大家踊跃投稿，介绍一些关于stata的数据处理和分析技巧。

投稿邮箱：

xueyuan19920310@163.com
投稿要求：
1）必须原创，禁止抄袭；
2）必须准确，详细，有例子，有截图；
注意事项：
1）所有投稿都会经过本公众号运营团队成员的审核，审核通过才可录用，一经录用，会在该推文里为作者署名，并有赏金分成。
2）邮件请注明投稿，邮件名称为“投稿”+“推文名称”。
3）应广大读者要求，现开通有偿问答服务，如果大家遇到关于stata分析数据的问题，可以在公众号中提出，只需支付少量赏金，我们会在后期的推文里给予解答。

欢迎关注爬虫俱乐部