我有两个使用lm
创建的线性模型,我想与stargazer
包中的表进行比较。在大多数情况下,我喜欢我得到的结果。但Akaike信息标准没有显示出来。 The docs我可以在"aic"
参数中传递keep.stat
以包含它。但它并不存在。没有错误消息。
stargazer(model1, model2, type="text", report="vc", header=FALSE,
title="Linear Models Predicting Forest Land",
keep.stat=c("aic", "rsq", "n"), omit.table.layout="n")
Linear Models Predicting Forest Land
==========================================
Dependent variable:
--------------------
forest
(1) (2)
------------------------------------------
log.MS.MIL.XPND.GD.ZS -11.948 -12.557
log.TX.VAL.AGRI.ZS.UN 2.310 2.299
log.NY.GDP.MKTP.CD 0.505
Constant 40.857 28.365
------------------------------------------
Observations 183 183
R2 0.142 0.146
==========================================
我没有看到为什么它无法包含它的任何理由。在这些模型上调用全局AIC
函数可以正常工作。
> AIC(model1)
[1] 1586.17
> AIC(model2)
[1] 1587.208
答案 0 :(得分:4)
问题由.AIC
中定义的stargazer:::.stargazer.wrap
函数给出
可以看出,此函数不计算lm
模型的AIC:
.AIC <- function(object.name) {
model.name <- .get.model.name(object.name)
if (model.name %in% c("coeftest")) {
return(NA)
}
if (model.name %in% c("lmer", "lme", "nlme", "glmer",
"nlmer", "ergm", "gls", "Gls", "lagsarlm", "errorsarlm",
"", "Arima")) {
return(as.vector(AIC(object.name)))
}
if (model.name %in% c("censReg")) {
return(as.vector(AIC(object.name)[1]))
}
if (model.name %in% c("fGARCH")) {
return(object.name@fit$ics["AIC"])
}
if (model.name %in% c("maBina")) {
return(as.vector(object.name$w$aic))
}
if (model.name %in% c("arima")) {
return(as.vector(object.name$aic))
}
else if (!is.null(.summary.object$aic)) {
return(as.vector(.summary.object$aic))
}
else if (!is.null(object.name$AIC)) {
return(as.vector(object.name$AIC))
}
return(NA)
}
.get.model.name
中的.AIC
功能调用.model.identify
。如果模型的组件call
为lm()
,则.model.identify
会返回ls
:
if (object.name$call[1] == "lm()") {
return("ls")
}
解决方案1 :使用add.lines
。
set.seed(12345)
n <- 100
df <- data.frame(y=rnorm(n), x1=rnorm(n), x2=rnorm(n))
model1 <- lm(y ~ x1, data=df)
model2 <- lm(y ~ x2, data=df)
library(stargazer)
stargazer(model1, model2, type="text", report="vc", header=FALSE,
title="Linear Models Predicting Forest Land",
keep.stat=c("rsq", "n"), omit.table.layout="n",
add.lines=list(c("AIC", round(AIC(model1),1), round(AIC(model2),1))))
,输出为:
Linear Models Predicting Forest Land
=================================
Dependent variable:
--------------------
y
(1) (2)
---------------------------------
x1 0.115
x2 -0.052
Constant 0.240 0.243
---------------------------------
AIC 309.4 310.3
Observations 100 100
R2 0.011 0.002
=================================
解决方案2 :将组件AIC
添加到模型对象。
model1 <- lm(y ~ x1, data=df)
model2 <- lm(y ~ x2, data=df)
model1$AIC <- AIC(model1)
model2$AIC <- AIC(model2)
stargazer(model1, model2, type="text", report="vc", header=FALSE,
title="Linear Models Predicting Forest Land",
keep.stat=c("aic", "rsq", "n"), omit.table.layout="n")
,输出
Linear Models Predicting Forest Land
======================================
Dependent variable:
--------------------
y
(1) (2)
--------------------------------------
x1 0.115
x2 -0.052
Constant 0.240 0.243
--------------------------------------
Observations 100 100
R2 0.011 0.002
Akaike Inf. Crit. 309.413 310.318
======================================