`isfe.fts` <- function (data, max.order = N - 3, N = 10, h = 5:10, method = c("classical", 
              "M", "rapca"), mean = TRUE, level = FALSE, fmethod = c("arima", "ets", 
			  "ets.na", "struct", "rwdrift", "rw"), lambda = 3, ...) 
{
    method <- match.arg(method)
    fmethod <- match.arg(fmethod)
    isfe.vec <- matrix(NA, max.order + 1, length(h))
    for (i in 0:max.order) {
        isfe.vec[i + 1, ] <- colMeans(aveISFE(data, 
            order = i, N = N, h = h, method = method, fmethod = fmethod, 
            lambda = lambda, mean = mean, level = level, ...))
    }
    rownames(isfe.vec) <- 0:max.order
    colnames(isfe.vec) <- h
    return(isfe.vec)
}