找回密码
 注册
查看: 9261|回复: 9

R语言 浪漫的R代码

[复制链接]
发表于 2011-3-31 20:34:24 | 显示全部楼层 |阅读模式
本帖最后由 genechip 于 2011-3-31 20:53 编辑

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 代码

浪漫的R 代码
回复

使用道具 举报

发表于 2011-3-31 20:55:38 | 显示全部楼层
沙发了。哈哈哈
回复 支持 反对

使用道具 举报

发表于 2011-3-31 21:22:35 | 显示全部楼层
这个更炫,Biostatistician(2)群小马巅巅  21:08:02 提供


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)
回复 支持 反对

使用道具 举报

发表于 2012-6-19 16:12:47 | 显示全部楼层
你们太厉害了{:soso_e179:}
回复 支持 反对

使用道具 举报

发表于 2012-6-21 19:53:32 | 显示全部楼层
太酷了!!!!!!soooooo 牛!
回复 支持 反对

使用道具 举报

发表于 2012-7-5 16:10:03 | 显示全部楼层
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)
######################
z=-Y
layout(matrix(1:2, 1, 2, byrow = TRUE))
plot(X,z,col=0,xlim=c(-4,0))
text(X,z,"O",col=a)
plot(X,Y,col=0,xlim=c(0,4))
text(X,Y,"O",col=a)
回复 支持 反对

使用道具 举报

发表于 2013-3-13 16:25:20 | 显示全部楼层
的确是好东西
回复 支持 反对

使用道具 举报

发表于 2013-4-26 10:06:20 | 显示全部楼层
spectacular{:soso_e113:}
回复 支持 反对

使用道具 举报

发表于 2014-4-28 18:41:39 | 显示全部楼层
哈哈 so cool  厉害
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

手机版|小黑屋|生物统计家园 网站价格

GMT+8, 2024-11-23 20:01 , Processed in 0.035604 second(s), 18 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表