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]]