我使用glm()
和rxGlm()
计算了几个模型(第二个来自Microsoft R
)。很遗憾,rxGlm()
未存储stargazer
所需的所有信息。因此,在尝试创建摘要表时(即使通过RxGlm
调整as.glm()
数据后),我收到以下错误消息:
Error in qr.lm(object) : lm object does not have a proper 'qr' component.
Rank zero or should not have used lm(.., qr=FALSE).
我已经逐步读出t统计数据和p值,并将它们分别反馈给stargazer
。但是,stargazer仍然需要将模型输出对象存储在工作区中,否则会发送错误消息。
这是我从模型输出中提取统计信息的方法:
obj1.t <- summary(obj1)$coef[ , "z value"]
obj1.p <- summary(obj1)$coef[ , "Pr(>|z|)"]
这是我stargazer
命令的简化形式,其中se =
和p =
用于反馈以前提取的统计信息。
stargazer(list(obj1, obj2),
type = "html", table.layout = "cd=!t-s-!a=!n", star.cutoffs=c(0.05,0.01,0.001), no.space = TRUE,
omit = c(1989:2015), font.size = "normalsize",
out = "Test.html",
df= FALSE,
column.labels = c("(1)", "(2)"),
add.lines = list(c("fixed effects", "No", "Yes")),
dep.var.labels = c("Dummy"),
title = "GLM PROBIT MODEL",
se= list(obj1.t, obj2.t),
p = list(obj1.p, obj2.p),
notes = "t statistics shown in parentheses")
现在我的问题是:有没有办法用stargazer创建回归输出表而不必提供模型输出对象?那么基本上将所有需要的数据存储在不同的向量中,然后将它们反馈给观星者? RxGlm
摘要提供了手动填充回归结果表所需的所有信息。但是,我正在寻找一种自动完成的方法。