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

看板R_Language作者 (dirk)時間11年前 (2013/04/08 21:57), 編輯推噓1(105)
留言6則, 5人參與, 最新討論串4/8 (看更多)
: : [問題敘述]: : : 想生成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

04/08 22:00, , 1F
04/08 22:00, 1F
謝謝喔^^ ※ 編輯: goddirk 來自: 122.117.36.192 (04/08 23:12)

04/09 01:12, , 2F
之前都聽人家說要少用R的迴圈
04/09 01:12, 2F

04/09 02:14, , 3F
能用向量就用向量. 等原po資料夠大就會明白了.
04/09 02:14, 3F

04/09 08:30, , 4F
loop真的要少用 資料不用太大 5萬筆就感覺得出差別了
04/09 08:30, 4F

04/09 21:41, , 5F
jit後應該還好吧
04/09 21:41, 5F
是RUN的時間會變很長嗎? ※ 編輯: goddirk 來自: 122.117.36.192 (04/10 00:22)

04/10 12:18, , 6F
jit? 願聞其詳
04/10 12:18, 6F
文章代碼(AID): #1HOip6H5 (R_Language)
文章代碼(AID): #1HOip6H5 (R_Language)