R 在函数内粘贴对象名称

R 在函数内粘贴对象名称,r,R,这是的后续问题(请参见数据和以前的命令) 从mods中的模型列表开始,我现在能够找到AIC最少的模型(对应于最佳模型): 棘手的部分是,我如何告诉glht使用放入best.mod的模型(注意:这都是在一个循环中发生的)。我试过了 由于模型存储在列表mods中,您可以使用which.min(list.AIC)的索引访问“最佳模型”: list.AIC您是否知道可以使用[[?即mods[[best.mod]]?按名称选择列表元素?不,我没有。似乎运行得很好!一般来说,使用[[从列表中提取项目更灵活、

这是的后续问题(请参见数据和以前的命令)

mods
中的模型列表开始,我现在能够找到AIC最少的模型(对应于最佳模型):

棘手的部分是,我如何告诉
glht
使用放入
best.mod
的模型(注意:这都是在一个循环中发生的)。我试过了


由于模型存储在列表
mods
中,您可以使用
which.min(list.AIC)
的索引访问“最佳模型”:


list.AIC您是否知道可以使用
[[
?即
mods[[best.mod]]
?按名称选择列表元素?不,我没有。似乎运行得很好!一般来说,使用
[[
从列表中提取项目更灵活、更有用,尤其是在处理其他变量时。
$主要是在仅为保存键入而交互工作时更方便。
    mods <- lapply(methods, function(m) 
      update(amod.null, correlation = getFunction(m)(1, form = ~ x + y), method="ML"))
    names(mods) <- methods
    list.AIC <- lapply(mods, function(x) AIC(x))
    best.mod <- names(which.min(list.AIC))
library(multcomp)
res <- glht(amod.null, mcp(Date = "Tukey"))
res <- glht(paste("mods$", as.factor(best.mod),sep = "") , mcp(Date = "Tukey"))
   names(mods)
[1] "corExp"    "corGaus"   "corLin"    "corRatio"  "corSpher"
list.AIC <- sapply(mods, AIC)
best.mod <- mods[which.min(list.AIC)]
best.mod[[1]]