作者:
locka (locka)
2019-03-27 23:19:58感謝 celestialgod 版主大大提點:
以前以為 *apply 家族的函數就已經是向量化(vectorized)的寫法了
查了資料才發現其實底層背後還是有 for 迴圈 (覺得震撼啊...)
試試看這樣的寫法
theta <- seq(0,1,len=100)
df <- rep(19,len=100)
n <- rep(20,len=100)
vrt <- Vectorize(rt)
x <- vrt(n=n, df=df, ncp=1/theta)
於是 x[,1] ... x[,100] 就是100個 n 等於20 然後對應各自 delta 值的 t 分配樣本了
(但是不知道 df, n 的預先定義有沒有意義?)
請版上各位高手再指點~ 謝謝大家
======
補充:
但還是有查到 *apply function 的好處:
1. 程式易讀性
2. 會 pre-allocate 向量的記憶體空間
2. 只影響區域變數不會改變全域變數
ref: https://www.r-bloggers.com/vectorization-in-r-why/
※ 引述《ntpuisbest (阿龍)》之銘言:
: n <- 20
: theta=seq(0,1,len=100)
: rt(n ,1/theta )
: 如題
: 我想要生100組 ,每組都是n=20的t分配樣本
: 只是這100組的theta都不一樣
: 我像上面那樣打 只會回傳20個樣本
: 並不是我想要的 2000個樣本 請問要如何打才能要我要的結果
: 想避免for loop
: 用loop的話 我知道怎麼做
作者:
VIATOR (阿布拉卡達不拉)
2019-04-11 23:57:00*apply各種變形各有些微不同,之後可能還要用do.call,rbind而且如果想用平行運算,可能還要改回foreach我雖然會用*apply了,但常常懷疑*apply的好處preallocate 大概是我目前唯一覺得的好處但懂得自己preallocate 的話,其實for迴圈也沒什麼不好