查看R包meboot中meboot函数如下：

<br />
function (x, reps = 999, trim = 0.1, reachbnd = TRUE, expand.sd = TRUE,<br />
    force.clt = TRUE, elaps = FALSE, colsubj, coldata, coltimes,<br />
    ...)<br />
{<br />
    if ("pdata.frame" %in% class(x)) {<br />
        res <- meboot.pdata.frame(x, reps, trim, reachbnd, expand.sd,<br />
            force.clt, elaps, colsubj, coldata, coltimes, ...)<br />
        return(res)<br />
    }<br />
    ptm1 <- proc.time()<br />
    n <- length(x)<br />
    xx <- sort(x)<br />
    ordxx <- order(x)<br />
    z <- rowMeans(embed(xx, 2))<br />
    dv <- abs(diff(as.numeric(x)))<br />
    dvtrim <- mean(dv, trim = trim)<br />
    xmin <- xx[1] - dvtrim<br />
    xmax <- xx[n] + dvtrim<br />
    aux <- colSums(t(embed(xx, 3)) * c(0.25, 0.5, 0.25))<br />
    desintxb <- c(0.75 * xx[1] + 0.25 * xx[2], aux, 0.25 * xx[n -<br />
        1] + 0.75 * xx[n])<br />
    ensemble <- matrix(x, nrow = n, ncol = reps)<br />
    ensemble <- apply(ensemble, 2, meboot.part, n, z, xmin, xmax,<br />
        desintxb, reachbnd)<br />
    qseq <- apply(ensemble, 2, sort)<br />
    ensemble[ordxx, ] <- qseq<br />
    if (expand.sd)<br />
        ensemble <- expand.sd(x = x, ensemble = ensemble, ...)<br />
    if (force.clt)<br />
        ensemble <- force.clt(x = x, ensemble = ensemble)<br />
    if (is.ts(x)) {<br />
        ensemble <- ts(ensemble, frequency = frequency(x), start = start(x))<br />
        dimnames(ensemble)[[2]] <- paste("Series", 1:reps)<br />
    }<br />
    ptm2 <- proc.time()<br />
    elapsr <- elapsedtime(ptm1, ptm2)<br />
    if (elaps)<br />
        cat("\n  Elapsed time:", elapsr$elaps, paste(elapsr$units,<br />
            ".", sep = ""), "\n")<br />
    list(x = x, ensemble = ensemble, xx = xx, z = z, dv = dv,<br />
        dvtrim = dvtrim, xmin = xmin, xmax = xmax, desintxb = desintxb,<br />
        ordxx = ordxx, elaps = elapsr)<br />
}<br />

为了一步步的看这个算法是怎么实现的，我将头部文件改为：

meboottest <- function (x, reps = 10, trim = 0.1, reachbnd = TRUE, expand.sd = TRUE,<br />
    force.clt = TRUE, elaps = FALSE, colsubj, coldata, coltimes,<br />
    ...)

然后我输如测试数据：

<br />
library(meboot)<br />
x <- c(4,12,36,20,8)<br />
meboottest(x)<br />

然后出现错误：

<br />
Error in match.fun(FUN) : object 'meboot.part' not found<br />

当然了在

<br />
library(meboot)<br />
meboot(x)<br />

则可以运行成功。不知道怎么回事？
</p>