R中的快速採樣

Question

有沒有比base::sample函數更快速的方法來獲取隨機子樣本（無需替換）？

Answer 1

沒有

我可以在3毫秒10,000個樣本在我的筆記本電腦更換。沒有更換，我可以在5毫秒內得到它們。從500個分佈中多次繪製，需要66 ms。你需要多快？

Answer 2

缺乏使用CUDA/gputools，我不認爲你會找到更快的方式。

你應該有興趣在走這樣一條路：

http://brainarray.mbni.med.umich.edu/brainarray/rgpgpu/

http://cran.r-project.org/web/packages/gputools/index.html

Answer 3

您可以通過消除基地::樣本函數調用得到一個加速的一點點：

> x<- rnorm(10000) 
> system.time(for(i in 1:100000) x[.Internal(sample(10000L, 10L, FALSE, NULL))]) 
    user system elapsed 
    2.873 0.017 2.851 
> system.time(for(i in 1:100000) sample(x,10)) 
    user system elapsed 
    3.420 0.025 3.258 

根據您的問題，可能還有其他更加巧妙的方法來加速您的代碼。想想辦法來取代許多小電話，以一個大的樣本。