几个有趣的R代码_生物软件圈_商圈

几个有趣的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)

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


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/

***我也来一段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

几个有趣的R代码

楼主 | 收藏 | 举报 2018-12-18 00:00 浏览: 118 回复: 3

此图代表我的心


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)

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


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)

楼主 | 收藏 | 举报 2018-12-18 00:00 浏览: 118 回复: 3

此图代表我的心


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)

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


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)

全部行业分类

生物软件

几个有趣的R代码

此图代表我的心

我心因你而有爱

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

本文作者

Ta的圈贴

最新创建圈子

最新话题

几个有趣的R代码

此图代表我的心

我心因你而有爱

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

此图代表我的心

我心因你而有爱

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

全部行业分类

生物软件

几个有趣的R代码

此图代表我的心

我心因你而有爱

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

相关推荐：

本文作者

Ta的圈贴

最新创建圈子

最新话题

几个有趣的R代码

此图代表我的心

我心因你而有爱

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

此图代表我的心

我心因你而有爱

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