作者goddirk (dirk)
看板R_Language
标题Re: [问题] 生成随机数据矩阵
时间Mon Apr 8 21:57:23 2013
: : [问题叙述]:
: : 想生成n=5 , 500组的exp随机数据
: : 然後要求每组的mean跟sd
: x <- matrix( rexp(2500), ncol=500 )
: x_mean <- apply(x, 2, mean)
: x_sd <- apply(x, 2, sd)
上面那个方法我有用过类似的
x<-array(rexp(2500),c(5,500))
不过好像都是先生成2500个资料在弄成矩阵
我想要的是一次5个资料产生,然後弄回圈做500组
在这边先感谢几位大大为我解惑
=我尽量养成习惯弄<-
-----------------------------------------------------------
以下是今天打出来的
e5m<-numeric(500)
e5s<-numeric(500)
for(i in 1:500) {
e5=rexp(5,rate=1);e5m[i]=mean(
e5);e5s[i]=sd(
e5)} #3
e5m
e5s
meane5<-mean(e5m)
sde5<-mean(e5s)
meane5
sde5
照我这部分打出来的部分
#3里面的
e5是否为同一组数据啊?
我觉得理论上应该是一样啦.....不过还是问清楚!!!
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 59.126.30.68
1F:→ andrew43:是 04/08 22:00
谢谢喔^^
※ 编辑: goddirk 来自: 122.117.36.192 (04/08 23:12)
2F:推 Yukirin:之前都听人家说要少用R的回圈 04/09 01:12
3F:→ andrew43:能用向量就用向量. 等原po资料够大就会明白了. 04/09 02:14
4F:→ MIZUYAMA:loop真的要少用 资料不用太大 5万笔就感觉得出差别了 04/09 08:30
5F:→ obarisk:jit後应该还好吧 04/09 21:41
是RUN的时间会变很长吗?
※ 编辑: goddirk 来自: 122.117.36.192 (04/10 00:22)
6F:→ Wush978:jit? 愿闻其详 04/10 12:18