Re: [問題] 生成隨機數據矩陣

[goddirk (dirk)]

: : [問題敘述]:
: : 想生成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
→ andrew43:是
謝謝喔^^
※ 編輯: goddirk         來自: 122.117.36.192       (04/08 23:12)
推 Yukirin:之前都聽人家說要少用R的迴圈
→ andrew43:能用向量就用向量. 等原po資料夠大就會明白了.
→ MIZUYAMA:loop真的要少用 資料不用太大 5萬筆就感覺得出差別了
→ obarisk:jit後應該還好吧
是RUN的時間會變很長嗎?
※ 編輯: goddirk         來自: 122.117.36.192       (04/10 00:22)
→ Wush978:jit? 願聞其詳