Merry Xmas: Blessing from Stata

2017-12-25 李金洋 爬虫俱乐部 爬虫俱乐部


有问题,不要怕!点击推文底部“阅读原文”下载爬虫俱乐部用户问题登记表并按要求填写后发送至邮箱statatraining@163.com,我们会及时为您解答哟~

ps:(1)爬虫俱乐部将于2018年1月20日至28日在武汉举行两期Stata编程技术定制培训。详情请戳《爬虫俱乐部Stata编程技术定制培训班——2018年1月武汉专场

(2)喜大普奔~爬虫俱乐部的github主站正式上线了!我们的网站地址是:https://stata-club.github.io,粉丝们可以通过该网站访问过去的推文哟~

♪当当当……随 着平安夜的落幕,考研的童鞋们终于得到了解放!与此同时,一年一度的圣诞节终于到啦,大家有没有收到苹果呢*^_^*尽管圣诞节是西方的宗教节日,但飘落 的雪花,祝福的歌声,温暖的炉火,和蔼的圣诞老人,可爱的麋鹿,无不构成记忆中的小确幸。最最重要的是,能和好朋友一起分享美食,互送礼物,回忆相处的点 滴,还有什么比这样的时光更值得珍惜呢♫在这样美好的日子中,爬虫俱乐部全体成员衷心祝愿屏幕前的小可爱们圣诞快乐^_^。同时,送上一份来自stata的绘图礼物——圣诞树


想知道怎么用stata画圣诞树吗??贴心的小编已经为您双手奉上程序啦,感兴趣的小伙伴可以下载我们提供的do文件哦,快拿去运行送给喜欢的人吧!

第一步

我们使用stata中的矩阵函数以及循环变量赋值生成圣诞树的主干的数据,并生成散点图

clear

matrix L = (0.03, 0.85, 0.8, 0.20, -0.20, 0.25, -0.20\ 0.00 ,0.00, 0.0, -0.08, 0.08, -0.10, 0.10\ 0.00, 0.00, 0.0 , 0.15, 0.15, 0.12, 0.12\ 0.10, 0.85, 0.80, 0.22, 0.22, 0.25, 0.20)

mat list L

matrix B = ( 0, 0.0, 0.0,0.00, 0.00, 0.0, 0.0\ 0, 1.5, 1.5, 0.85, 0.85, 0.3, 0.4)

mat list B

matrix prob = (0.02, 0.6,.08, 0.07, 0.07, 0.07, 0.07)

set obs 100000

gen x1 = .

gen x2 = .

gen k = int(runiform(1,7))

replace x1 = 0 if _n==1

replace x2 = 2 if _n==1

forvalues i = 2(1)100000 {

       mat x = (x1[`i'-1],x2[`i'-1])

       mat a = (L[1,k[`i']],L[3,k[`i']]\L[2,k[`i']],L[4,k[`i']])

       mat b = a'*x'+ B[1...,k[`i']]

       replace x1 = b[1,1] if _n == `i'

       replace x2 = b[2,1] if _n == `i'

       di `i'

}

scatter x2 x1 ,color(green) jitter(2) title("Merry Christmas")

第二步

根据往期推文《绘图神器—gr_edit,妈妈再也不用担心我的画图》,用gr_edit命令为圣诞树添加装饰。

首先进入图形编辑器,点击红色按钮修改图片:

其次改变图像尺寸Y size

我们便可以将圣诞树的主干变得更加美观:

最后使用绘图工具为圣诞树添加装饰,点击红色按钮结束录制,根据往期推文《绘图神器—gr_edit,妈妈再也不用担心我的画图》将命令保存到指定文件:

打开命令文件夹如下图所示:

运用绘图神器gr_edit构建命令,通过命令绘制圣诞树的装饰命令。

gr_edit style.editstyle declared_ysize(8) editcopy // Graph size

gr_edit plotregion1.AddMarker added_markers editor -.1173865320507047 6.508580319160336

gr_edit plotregion1.added_markers_new = 1

gr_edit plotregion1.added_markers_rec = 1

gr_edit plotregion1.added_markers[1].style.editstyle marker( symbol(circle) linestyle( width(thin) color(red) pattern(solid)) fillcolor(red) size(large) angle(stdarrow) symangle(zero) backsymbol(none) backline( width(thin) color(black) pattern(solid)) backcolor(black) backsize(zero) backangle(stdarrow) backsymangle(zero)) line( width(thin) color(black) pattern(solid)) area( linestyle( width(thin) color(ltbluishgray) pattern(solid)) shadestyle( color(ltbluishgray) intensity(inten100) fill(pattern10))) label( textstyle( horizontal(center) vertical(middle) angle(default) size(small) color(black) position() margin(zero) linestyle( width(thin) color(black) pattern(solid))) position(6) textgap(half_tiny) horizontal(default) vertical(default)) dots( symbol(circle) linestyle( width(thin) color(black) pattern(solid)) fillcolor(black) size(vtiny) angle(horizontal) symangle(zero) backsymbol(none) backline( width(thin) color(black) pattern(solid)) backcolor(black) backsize(medium) backangle(horizontal) backsymangle(zero)) connect(direct) connect_missings(yes) editcopy 

运用上述的第一步和第二步内容,对圣诞树进行各种颜色的点缀,最终我们就生成漂亮的圣诞树啦。

由于添加点缀的命令太长,考虑推文篇幅,我们就不将程序附在后面啦,有兴趣的读者可以通过爬虫俱乐部官网下载(https://raw.githubusercontent.com/Stata-Club/Sharing-Center-of-Stata-Club/master/article/%E5%9C%A3%E8%AF%9E%E6%A0%91.do)。

最后,祝大家双旦快乐,么么哒。

注:此推文中的图片及封面(除操作部分的)均来源于网络!如有雷同,纯属巧合!


2017圣诞快乐

愿每一个心愿都能成真

Wish you a merry Christmas!

以上就是今天给大家分享的内容了,说得好就赏个铜板呗!有钱的捧个钱场,有人的捧个人场~。另外,我们开通了苹果手机打赏通道,只要扫描下方的二维码,就可以打赏啦!

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

文字编辑/闫续文

技术总编/刘贝贝


往期推文推荐:

1.高校学术大神:你的导师上榜了吗?

2.中国高校财经、管理与综合类期刊灌水排行榜

3.命令sum2docx输出统计量表到docx文件

4.reg2docx报告你的实证结果吧!

5.爬虫俱乐部又出新命令了——wordconvert转换你的word文件

6.putdocx+wordconvert—将实证结果输出到Word.docx)文档

7.Stata 15Markdown——没有做不到,只有想不到!

8.矩阵和宏的故事

9.shelloutopen anything

10.Putpdf--神奇的转换工具



关于我们

微信公众号“爬虫俱乐部”分享实用的stata命令,欢迎转载、打赏。爬虫俱乐部是由李春涛教授领导下的研究生及本科生组成的大数据分析和数据挖掘团队。

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

投稿邮箱:statatraining@163.com

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

欢迎关注爬虫俱乐部

微信扫一扫
关注该公众号

公开收集正在保存拖拽至此处以完成收集私密