最近看到sas画图娱乐比较火，用R的去pk一下

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

基于SAS/GRAPH绘制动画雪景图 ：
https://bbs.pinggu.org/thread-2137713-1-1.html

用SAS画蝴蝶：
https://bbs.pinggu.org/thread-2137037-1-1.html
（这个已经用R复制了）

上面的那个有谁用R复制一下，验证一下R的强大。

ps： 茶余饭后可以娱乐一下，希望不要占用科研时间。

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：pinggu thread GRAPH HTML READ 画图 娱乐

支持啊

支持一下

支持一下

library(ggplot2)npoints=500npointsb=1200opt=theme(legend.position="none",panel.background = element_blank(),panel.grid = element_blank(),axis.ticks=element_blank(),axis.title=element_blank(),axis.text =element_blank())t=seq(0,10*pi,length=npointsb)butterfly=data.frame(x=sin(t)*(exp(1)^cos(t)-2*cos(4*t)-(sin(t/12))^5), y=cos(t)*(exp(1)^cos(t)-2*cos(4*t)-(sin(t/12))^5), s=runif(npointsb, min=.1, max=10), f=factor(sample(1:10,npointsb,TRUE)), a=runif(npointsb,min=.1, max=.4))points=data.frame(x=runif(npoints,-4,4), y=runif(npoints,-3,5), s=runif(npoints,min=30, max=50), f=factor(sample(1:10,npoints,TRUE)), a=runif(npoints,min=.05, max=.15))data=rbind(butterfly, points)ggplot(data, aes(x, y, colour=f))+geom_point(alpha=data$a,size=data$s)+opt

1. t=seq(-400, 4, by=0.01)
   
2. a=cos(t)
   
3. b=cos(4*t)
   
4. x=sin(t)*(exp(a)-2*b-(sin(t/12))*sin(t/12)*sin(t/12)*sin(t/12)*sin(t/12))
   
5. y = cos(t)*(exp(a)-2*b-(sin(t/12)*sin(t/12)*sin(t/12)*sin(t/12)*sin(t/12)))
   
6. x1=seq(0.20, (pi/4+0.20), by=0.01)
   
7. y1=3*sin(2*(x1-0.20))+2.5;
   
8. x2=-1*x1;
   
9. 
   
10. plot(x, y, xlim=c(-4,4), ylim=c(-3,6),col=4, lty=2, pch=3, cex=0.7)
    
11. points(x1,y1, col=2, pch="*", cex=1.5)
    
12. points(x2,y1, col=2, pch="*", cex=1.5)

复制代码