[英]How to Create AIC Model Selection Table in R in LaTex format?
希望为LaTex格式的出版物创建AIC选择表,但我似乎无法获得所需的表格。 我已经用谷歌搜索将其杀死,很惊讶我找不到答案。 我在R中找到了更多晦涩问题的答案。
下面是一些代码,我制作了一些我不太喜欢的表,最下面是我想制作的一般结构,但是像注视者软件包一样采用了很好的乳胶表格式。
我试图对这两个软件包使用额外的参数来达到我想要的但没有成功。
##Create dummy variables
a<-1:10
b<-c(10:3,1,2)
c<-c(1,4,5,3,7,3,6,2,4,5)
##Create df
df<-data.frame(a,b,c)
##Build models
m1<-lm(a~b,data=df)
summary(m1)
m2<-lm(a~c,data=df)
m3<-lm(a~b+c,data=df)
m4<-lm(a~b*c,data=df)
##View list of AIC values
AIC(m1,m2,m3,m4)
########################CREATE AIC SELECTION TABLE
##Using MuMIn Package
library(MuMIn)
modelTABLE <- model.sel(m1,m2,m3,m4)
View(modelTABLE) ##No AIC values, just AICc, no R-squared, and model name (i.e, a~b) not present
##Using stargazer Package
library(stargazer)
test<-stargazer(m1,m2,m3,m4 ,
type = "text",
title="Regression Results",
align=TRUE,
style="default",
dep.var.labels.include=TRUE,
flip=FALSE
## ,out="models.htm"
)
View(test) ##More of a table depicting individual covariate attributes, bottom of table doesn't have AIC
###Would like a table similar to the following
Model ModelName df logLik AIC delta AICweight R2
m1 a ~ b 3 -6.111801 18.2 0 0.95 0.976
m3 a ~ b + c 4 -5.993613 20 1.8 0.05 0.976
m4 a ~ b * c 5 -5.784843 21.6 3.4 0.00 0.977
m2 a ~ c 3 -24.386821 54.8 36.6 0.00 0.068
`
model.sel
结果是data.frame
,这样你就可以对其进行修改(添加型号名称,圆号码等)和出口使用例如乳胶latex
从Hmisc
包。
# include R^2:
R2 <- function(x) summary(x)$r.squared
ms <- model.sel(m1, m2, m3, m4, extra = "R2")
i <- 1:4 # indices of columns with model terms
response <- "a"
res <- as.data.frame(ms)
v <- names(ms)[i]
v[v == "(Intercept)"] <- 1
# create formula-like model names:
mnames <- apply(res[, i], 1, function(x)
deparse(simplify.formula(reformulate(v[!is.na(x)], response = response))))
## OR
# mnames <- apply(res[, i], 1, function(x)
# sapply(attr(ms, "modelList"), function(x) deparse(formula(x)))
res <- cbind(model = mnames, res[, -i])
Hmisc::latex(res, file = "")
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.