当我试图执行下面的代码时,我得到了一个错误"Error in statistic(data, original, ...) : unused argument (original)"。具体来说,我试图为每个plot#应用一个函数并引导这些结果。我是不是错过了直截了当的东西?
library(codyn) # has sample dataset called pplots
library(boot)
str(pplots)
stability <- function(x){
mean(multivariate_change(x,
species.var = "species",
time.var = "year",
abundance.var = "relative_cover",
replicate.var = "plot")$composition_change)
}
boot_obj <- lapply(splitspplots,boot,statistic = stability,R = 20)
boot_obj <- lapply(splitspplots,myBootFun)
myBootFun <- function(x,i) {
lapply(x[i],boot, statistic=stability, R = 10)
}
myBootFun <- function(x,i) {
boot(x[i], statistic=stability, R = 10)
}
splitspplots <- split(pplots,pplots$plot)
lapply(splitspplots, print(myBootFun))发布于 2020-03-13 12:57:02
我想出了一种不同的方法来做到这一点。
我使用for循环在不同的组/站点上循环,但在这里没有显示这一点。这是引导代码。
data_sampled <- data_grouped %>%
group_by(SITE_ID) %>%
crossing(id = seq(100)) %>% # create replicate dataframes
sample_n(size = n(), replace = T) # sample dataframeshttps://stackoverflow.com/questions/60670522
复制相似问题