【发布时间】:2019-02-25 04:11:52
【问题描述】:
我想编写一个函数并使用lapply 自动运行一些具有不同协变量的模型。但我无法更改每个模型需要使用的变量。
bootstrap_vars <- list(vars1=c(1,2,3), vars2=c(2,3,4), vars3=c(1,3,4))
full_mod <- lm(Y ~ X[,full_vars])
getFstat <- function(Vars, X, Y, full_mod) {
boot_mod <- lm(Y ~ X[,Vars], data=data)
ano <- anova(full_mod, boot_mod, test="F")
# cat("Variables:", Vars, ". F-stat: ", ano$F[2], "\n")
ano$F[2] # Return F-statistic between full model and one bootstrap model
}
all_F <- unlist(lapply(bootstrap_vars, FUN = getFstat, X=X, Y=Y, full_mod=full_mod))
无法识别lapply 中的X 并显示anova 具有不同的boot_mod 和full_mod 数据集?我怎样才能让它工作?
【问题讨论】: