我有很多模特,而且我已经厌倦了手动获取每一个结果。首先,我为解释的方差创建一个函数。
exp.dev <- function(x) {
1-(x$deviance/x$null.deviance)
}
但我无法进行此循环。例如,假设我有三个glm模型&#34; mod1,mod2,mod3&#34;我的目的是写另一个函数,包括我的exp.dev函数,并将这三个模型一次循环。
我试过这个
for (i in 1:3) exp.dev(modi)
但我得到了
Error in exp.dev(modi) : object 'modi' not found
我尝试了一些不同的版本,但我认为这并不像我想的那么简单。