作者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