作者:絮尘飘雪_896 | 来源:互联网 | 2022-12-02 00:51
我希望能够将命名的模型列表(merMod对象)传递给anova()并在输出中保留模型名称.在使用mclapply()以更高效并行运行一批慢速模型(如glmer)的情况下,这尤其有用.我提出的最好的是在模型列表的去命名版本上使用do.call,但这并不理想,因为我可能有名为(例如)"mod12","mod17"和"mod16"的模型并且这些模型名称在输出中被转换为"MODEL1","MODEL2"和"MODEL3".(在查看单个批次时,这看起来似乎微不足道,但在长时间的建模会话中,有几十个模型,这是一个混乱的确定方法.)
请注意,这与从列表创建和调用线性模型的问题不同,因为我不是要尝试跨列表比较模型对.它比在模型列表上使用lapply更复杂,因为我以非一元的方式使用anova().
这是一个最小的代表:
library(lme4)
formList <- c(mod12 = angle ~ recipe + temp + (1|recipe:replicate),
mod17 = angle ~ recipe + temperature + (1|recipe:replicate),
mod16 = angle ~ recipe * temperature + (1|recipe:replicate))
modList <- lapply(formList, FUN=lmer, data=cake)
# Fails because modList is named so it's interpreted as arg-name:arg pairs
do.call(anova, modList)
# Suboptimal because model names aren't preserved
do.call(anova, unname(modList))
# Fails because object isn't merMod (or some other class covered by methods("anova"))
do.call(anova, list(object=modList[1], ...=modList[-1], model.names=names(modList)))
第二个do.call返回:
Data: ..1
Models:
MODEL1: angle ~ recipe + temp + (1 | recipe:replicate)
MODEL2: angle ~ recipe + temperature + (1 | recipe:replicate)
MODEL3: angle ~ recipe * temperature + (1 | recipe:replicate)
Df AIC BIC logLik deviance Chisq Chi Df Pr(>Chisq)
MODEL1 6 1708.2 1729.8 -848.08 1696.2
MODEL2 10 1709.6 1745.6 -844.79 1689.6 6.5755 4 0.1601
MODEL3 20 1719.0 1791.0 -839.53 1679.0 10.5304 10 0.3953
理想情况下,输出看起来像这样:
Data: ..1
Models:
mod12: angle ~ recipe + temp + (1 | recipe:replicate)
mod17: angle ~ recipe + temperature + (1 | recipe:replicate)
mod16: angle ~ recipe * temperature + (1 | recipe:replicate)
Df AIC BIC logLik deviance Chisq Chi Df Pr(>Chisq)
mod12 6 1708.2 1729.8 -848.08 1696.2
mod17 10 1709.6 1745.6 -844.79 1689.6 6.5755 4 0.1601
mod16 20 1719.0 1791.0 -839.53 1679.0 10.5304 10 0.3953
我该怎么做呢?如果这意味着我得到一个更容易理解的输出,我对anova()周围的丑陋包装非常满意.
1> Roland..:
您可以传递这样的符号列表:
with(modList,
do.call(anova,
lapply(names(modList), as.name)))
#refitting model(s) with ML (instead of REML)
#Data: ..1
#Models:
#mod12: angle ~ recipe + temp + (1 | recipe:replicate)
#mod17: angle ~ recipe + temperature + (1 | recipe:replicate)
#mod16: angle ~ recipe * temperature + (1 | recipe:replicate)
# Df AIC BIC logLik deviance Chisq Chi Df Pr(>Chisq)
#mod12 6 1708.2 1729.8 -848.08 1696.2
#mod17 10 1709.6 1745.6 -844.79 1689.6 6.5755 4 0.1601
#mod16 20 1719.0 1791.0 -839.53 1679.0 10.5304 10 0.3953