程序员的浪漫--R语言画小心心

Fig 1.

n=50000
r=0.6
r2=(1-r*r)^0.5
X=rnorm(n)
Y=X*r+r2*rnorm(n)
Y=ifelse(X>0,Y,-Y)
plot(X,Y,col="red")

1.JPG

Fig 2

x <- seq(-1,1,length = 5000)
y1 <- (1-x²⁾(1/2)+(x²⁾(1/3)
y2 <- -(1-x²⁾(1/2)+(x²⁾(1/3)
plot(c(rev(x),x),c(rev(y2),y1), type = 'l' ,xlim = c(-1,1), xlab = "",ylab = "", axes = F,col=2)

好像有点单调对吧？？？

那我们来参考一下加州伯克利大佬们的浪漫

library("RColorBrewer")
col<-brewer.pal(3,'Set1')[1]

RColorBrewer这个包可以玩好久

t<-seq(0,60,len=100)
x<--0.01(-t^2+40t+1200)sin(pit/180)
y<- 0.01(-t^2+40t+1200)cos(pit/180)
plot(c(-12,14),c(0,20),type = 'n',axes = T,xlab = '',ylab = '')

一定要先画箭头

arrows(-11,8,14,8,lwd=8,col='gold')
arrows(14,8,-11,8,lwd=8,col='gold',angle = 120)
polygon(x,y,col=col,border = NA)
polygon(-x,y,col=col,border = NA)
lines(x,y,lwd=1)
lines(-x,y,lwd=1)

Fig3

参考文献：
https://blog.csdn.net/kMD8d5R/article/details/79327603
https://www.r-bloggers.com/make-a-valentines-heart-with-r/
​

2.JPG 3.JPG

最后附上一个大佬的教程，他开发了出名的clusterProfiler

https://mp.weixin.qq.com/s/Iqcqo5YTy1qmXdsuIrlrYQ

据说还有分手函数图。。真的是。。。

分手函数：17 x^2 – 16|x|y + 17 y^2 + 150/|5 x + sin(5 y)| < 225

http://www.matrix67.com/blog/archives/4139</pre>

PS：
据说笛卡尔最后还是没能等到 自己等了很久的那个人
但是这不重要，重要的是。。。群里好多有心型图啊
不知道作者们愿不愿意放上来分享一下（哈哈哈哈）
各位端午节快乐，没脱单的早日脱单
以及，下次再有人问你R语言能干什么，告诉他可以拿去表白
和大神们比画图我还是嫩了点（顶锅盖逃走）