如何加速我的代码?

时间:2013-10-01 20:58:16

标签: r bootstrapping

如果可能,请帮助我加快下面的引导代码。

sampleOne   <- function(x) sample(x, replace = TRUE)
sampleMany  <- function(x, n) replicate(n, sampleOne(x), simplify = FALSE)
listMeans <- function(x, n) lapply(sampleMany(x, n), mean)
bootData <- function(x,n) do.call(rbind, listMeans(x,n))
sampleSize <- 100000
numBoots <- 1000
# Left Skewed distribution # shape1 = a and shape2 = b
set.seed(400)
popSkewLeft <- rbeta(sampleSize, shape1 = 5, shape2 = 1)
skewLeftbootData <- bootData(popSkewLeft, numBoots)
(bootSd <- sd(skewLeftbootData) * sqrt(sampleSize))

我知道我可以使用启动包,但我的目标是在基础R中磨练我的技能

0 个答案:

没有答案