520!!!Stata和R助力表白
发布:Captain-CUI | 分类:会计库
关于本站
人大经济论坛-经管之家:分享大学、考研、论文、会计、留学、数据、经济学、金融学、管理学、统计学、博弈论、统计年鉴、行业分析包括等相关资源。
经管之家是国内活跃的在线教育咨询平台!
经管之家新媒体交易平台
提供"微信号、微博、抖音、快手、头条、小红书、百家号、企鹅号、UC号、一点资讯"等虚拟账号交易,真正实现买卖双方的共赢。【请点击这里访问】
论文
- 毕业论文 | 写毕业论文
- 毕业论文 | 为毕业论文找思路
- 毕业论文 | 可以有时间好好写 ...
- 毕业论文 | 毕业论文如何选较 ...
- 毕业论文 | 毕业论文选题通过 ...
- 毕业论文 | 还有三人的毕业论 ...
- 毕业论文 | 毕业论文答辩过程 ...
- 毕业论文 | 本科毕业论文,wi ...
考研考博
- 考博 | 南大考博经济类资 ...
- 考博 | 考博英语10000词汇 ...
- 考博 | 如果复旦、南大这 ...
- 考博 | 有谁知道春招秋季 ...
- 考博 | 工作与考博?到底 ...
- 考博 | 考博应该如何选择 ...
- 考博 | 考博失败了
- 考博 | 考博考研英语作文 ...
留学
- 日本留学 | 在日本留学心得
- 日本留学 | 日本留学生活必需 ...
- 日本留学 | 【留学日本】2015 ...
- 日本留学 | 日本海外留学8年来 ...
- 日本留学 | 日本留学费用_日本 ...
- 日本留学 | 求在日本留学的师 ...
- 日本留学 | 日本留学的有没有 ...
- 日本留学 | 日本留学
TOP热门关键词
一切都是从一个故事开始的:1650年,斯德哥尔摩街头,一个宁静的午后,笛卡尔邂逅了18岁的瑞典公主克里斯汀。机遇巧合,一段纯粹、美好的爱情悄然萌发。然而,没过多久,他们的恋情传到了国王的耳朵里。国王大怒,将 ...
免费学术公开课,扫码加入![]() |
一切都是从一个故事开始的:
1650年,斯德哥尔摩街头,一个宁静的午后,笛卡尔邂逅了18岁的瑞典公主克里斯汀。
机遇巧合,一段纯粹、美好的爱情悄然萌发。
然而,没过多久,他们的恋情传到了国王的耳朵里。国王大怒,将笛卡尔放逐,而公主被软禁在宫中。 身体孱弱的笛卡尔不久便染上重病,在生命进入倒计时的那段日子,他日夜思念的还是街头偶遇的那张温暖的笑脸。他每天坚持给她写信,盼望着她的回音。然而,这些信都被国王拦截下来。在笛卡尔给克里斯汀寄出第十三封信后,他永远地离开了这个世界。此时,被软禁在宫中的小公主依然徘徊在皇宫的走廊里,思念着远方的情人。
这最后一封信上没有写一句话,只有一个方程:r=a(1-sinθ)。——传说中著名的心形函数。
不过,事实告诉我们,除非你的目标妹子也是一只Geeker(至少会用Mathematica或者MATLAB等软件),否则像笛卡尔这样单给一个函数的结局大概就是你推公式别人推妹子了……表白什么的还是选择更浅显易懂的方法吧。
在这个特殊的日子里,特意搜集了利用Stata绘制心形线的程序,代码如下
- clear
- range t 0 2*_pi 1000
- gen x=16*sin(t)^3
- gen y=13*cos(t)-5*cos(2*t)-2*cos(3*t)-cos(4*t)
- egen x_min=min(x)
- egen x_max=max(x)
- egen y_min=min(y)
- egen y_max=max(y)
- gen a=(x-x_min)/(x_max-x_min)
- gen b=(y-y_min)/(y_max-y_min)
- line b a
- gr_edit yaxis1.draw_view.setstyle, style(no)
- gr_edit xaxis1.draw_view.setstyle, style(no)
- gr_edit plotregion1.AddTextBox added_text editor .7055394244311991 .2810707216715078
- gr_edit plotregion1.added_text_new = 1
- gr_edit plotregion1.added_text_rec = 1
- gr_edit plotregion1.added_text[1].style.editstyleangle(default) size(medsmall) color(red) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(thin) color(red) pattern(solid)) box_alignment(east) editcopy
- gr_edit plotregion1.added_text[1].style.editstyle size(large) editcopy
- gr_edit plotregion1.added_text[1].text = {}
- gr_edit plotregion1.added_text[1].text.Arrpush" I LOVE YOU"
- graph export "C:\Desktop\520.png", as(png) wid(800)hei(600) replace
以上程序基于stata12运行,如果有stata14的话,还可以输入中文!祝大家好运!
注:如果图片中字符没有居中的话,请在前面输入空格调整!附上do文件
以上程序来自:http://mp.weixin.qq.com/s?__biz= ... =21#wechat_redirect
接下来是R绘制心形线,代码如下:
- 1)加加载程序包
- library(grid)
- 2)绘制心形函数
- heart <- function(lcolor){
- t=seq(0, 2*pi, by=0.1)
- x=16*sin(t)^3
- y=13*cos(t)-5*cos(2*t)-2*cos(3*t)-cos(4*t)
- a=(x-min(x))/(max(x)-min(x))
- b=(y-min(y))/(max(y)-min(y))
- grid.lines(a,b,gp=gpar(col=lcolor,lty = "solid",lwd = 6))
- }
- heart("hotpink")
- grid.newpage()
- 2)绘制玫瑰函数
- rose=function(){
- grid.circle(x=0.5, y=0.5, r=0.5,gp=gpar(fill="red",lwd = 3))
- vp <- viewport(.5, .5, w=.9, h=.9)
- pushViewport(vp)
- grid.polygon(x=c(0.08, .5, 0.94),y=c(.22, 1.03, .22),gp=gpar(lwd =
- vp2 <- viewport(.5, .5, w=.4, h=.4)
- pushViewport(vp2)
- grid.circle(x=0.5, y=0.5, r=0.5,gp=gpar(fill="red",lwd = 3))
- vp3 <- viewport(.5, .5, w=.9, h=.9,angle=180)
- pushViewport(vp3)
- grid.polygon(x=c(0.08, .5, 0.94),y=c(.22, 1.03, .22),gp=gpar(lwd = 3))}
- rose()
- 3)调整图形
- grid.newpage()
- pushViewport(viewport(x=0.1, y=0.1,w=0.2, h=0.2))
- grid.newpage()
- for (j in 1:30) {
- vp <- viewport(.5, .5, w=.9, h=.9)
- pushViewport(vp)
- heart("hotpink")
- }
- 4)调整图形及添加文字
- grid.newpage()
- vp1 <- viewport(.4, .5, w=.5, h=.5,angle=15)
- pushViewport(vp1)
- heart("red")
- vp2 <- viewport(0.9, .27, w=.7, h=.7,angle=-30)
- pushViewport(vp2)
- heart("hotpink")
- grid.text("执手偕老",
- x=0.2,y =1.2, just = c("center", "bottom"),
- gp = gpar(fontsize=30), vp = vp)
- vp3 <- viewport(-0.65, 1.2, w=.3, h=.3,angle=-30)
- pushViewport(vp3)
- rose()
以上程序来自http://mp.weixin.qq.com/s?__biz=MzA4NTAyMjQ2Mg==&mid=206788501&idx=1&sn=53e07723fdb395a787e57ea53bffd332&scene=21#wechat_redirect
没有贴出matlab的程序,表示略遗憾!已经利用matlab抱得美人归的朋友不妨分享一下程序
![](static/image/smiley/default/titter.gif)
最后,感谢那些做出巨大贡献的程序猿,没有他们的代码,表白怎能如此风骚!!!
「经管之家」APP:经管人学习、答疑、交友,就上经管之家!
免流量费下载资料----在经管之家app可以下载论坛上的所有资源,并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库,各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里,你可以找到任何学科方向、有共同话题的朋友。
经管之家(原人大经济论坛),跨越高校的围墙,带你走进经管知识的新世界。
扫描下方二维码下载并注册APP
免流量费下载资料----在经管之家app可以下载论坛上的所有资源,并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库,各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里,你可以找到任何学科方向、有共同话题的朋友。
经管之家(原人大经济论坛),跨越高校的围墙,带你走进经管知识的新世界。
扫描下方二维码下载并注册APP
![](static/2016/reg/jgewm.png)
您可能感兴趣的文章
本站推荐的文章
- 哲学名言 | 【独家发布】经典哲学名言
- 哲学书籍 | 求推荐一本讲人生目标的哲学书籍 ...
- 哲学书籍 | 20部必读的哲学书籍
- 哲学书籍 | 经济人,开拓你逻辑思维的哲学书 ...
- 哲学书籍 | 哲学书籍
- 哲学书籍 | 哲学书籍
- 哲学书籍 | 哲学书籍
- 哲学书籍 | 经典的哲学书籍
人气文章
1.凡人大经济论坛-经管之家转载的文章,均出自其它媒体或其他官网介绍,目的在于传递更多的信息,并不代表本站赞同其观点和其真实性负责;
2.转载的文章仅代表原创作者观点,与本站无关。其原创性以及文中陈述文字和内容未经本站证实,本站对该文以及其中全部或者部分内容、文字的真实性、完整性、及时性,不作出任何保证或承若;
3.如本站转载稿涉及版权等问题,请作者及时联系本站,我们会及时处理。
2.转载的文章仅代表原创作者观点,与本站无关。其原创性以及文中陈述文字和内容未经本站证实,本站对该文以及其中全部或者部分内容、文字的真实性、完整性、及时性,不作出任何保证或承若;
3.如本站转载稿涉及版权等问题,请作者及时联系本站,我们会及时处理。