生物软件
加入本圈
收藏本贴

几个有趣的R代码

楼主  收藏   举报   帖子创建时间:  2018-12-18 00:00 回复：3 关注量：118

此图代表我的心


n=50000;
r=0.7;r_e=(1-r*r)^.5;
X=rnorm(n);
Y=X*r+r_e*rnorm(n);
Y=ifelse(X>0,Y,-Y);
plot(X,Y,col="pink")

R语言脚本，有趣的图

我心因你而有爱


n=50000;
r=0.7;r_e=(1-r*r)^.5;
X=rnorm(n);
Y=X*r+r_e*rnorm(n);
Y=ifelse(X>0,Y,-Y);
a<-sample(c(2,6,7,8),50000,T)
b<-sample(c(76,79,86,69),50000,T)
plot(X,Y,col=0)
text(X,Y,"lOVE",col=a)

用R语言画心

此图太猥琐了，运行结果就不贴上了


xrange = c(-15, 15)
yrange = c(0, 16)
plot(0, xlim = xrange, ylim = yrange, type = "n")
#plot(xrange, yrange, type="n", ann=FALSE, axes=FALSE)
yr = seq(yrange[1], yrange[2], len = 50)
offsetFn = function(y) {
2 * sin(0 + y/3)
}
offset = offsetFn(yr)
leftE = function(y) {
-10 - offsetFn(y)
}
rightE = function(y) {
10 + offsetFn(y)
}
xp = c(leftE(yr), rev(rightE(yr)))
yp = c(yr, rev(yr))
polygon(xp, yp, col = "#ffeecc", border = NA)
h = 9
xt = seq(0, rightE(h), len = 100)
yt = log(1 + log(1 + log(xt + 1)))
yt = yt - min(yt)
yt = h * yt/max(yt)
x = c(leftE(h), rightE(h), rev(xt), -xt)
y = c(h, h, rev(yt), yt)
polygon(x, y, col = "red", border = NA)

来源：http://yixf.name/2011/04/02/%E6%B5%AA%E6%BC%AB%E7%9A%84%E8%A1%A8%E7%99%BD%E4%B8%8E%E7%8C%A5%E7%90%90%E7%9A%84%E7%9B%B8%E6%80%9D-%E5%87%A0%E4%B8%AA%E6%9C%89%E8%B6%A3%E7%9A%84r%E4%BB%A3%E7%A0%81/

分享到：
  2. 评论

相关推荐：

  1. NCBI在线BLAST使用方法与结果详解 2938
  2. 神经网络术语：Epoch、Batch Size和迭代 527
  3. Consed的安装与使用教程 465
  4. 陈连福的NGS生物信息学培训教材V2.1 277
  5. WGCNA分析使用教程 272
更多圈贴本圈首页
手机扫一扫，随时随地看商友圈
  • ***我也来一段matlab代码： 1970-01-01 08:00
    #1

  • ***x=[-1.65:0.00001:1.65];hold 1970-01-01 08:00
    #2

  • ***1boya888 2012-06-12 15:50
    #3

    xysheep是来踢馆的？！ :razz:matlab PK R

本文作者

Ta的圈贴

更多

最新创建圈子

最新话题

我要发帖