我正在尝试为观星者不支持的模型复制漂亮的观星者模型(lm(输出。
线性模型观星器的输出可以手工制作吗？由于我们可以从每个模型中创建一个数据帧，然后将创建的数据帧插入stargezer：
library(spdep)
data(afcon, package="spData")
afcon$Y = rnorm(42, 50, 20)
cns <- knearneigh(cbind(afcon$x, afcon$y), k=7, longlat=T)
scnsn <- knn2nb(cns, row.names = NULL, sym = T)
W <- nb2listw(scnsn, zero.policy = TRUE)
ols <- lm(totcon ~ Y, data = afcon)
spatial.lag <- lagsarlm(totcon ~ Y, data = afcon, W)
summary(model)
stargazer(ols, type = "text")
summary(spatial.lag)
data.frame(
spatial.lag$coefficients,
spatial.lag$rest.se
) %>%
rename(coeffs = spatial.lag.coefficients,
se = spatial.lag.rest.se) %>%
stargazer(type = "text", summary = F)
当我们做stargazer(ols)的输出非常好时，我想手工为spatial.lag复制相同的输出，有没有办法做到这一点，如何上标等等…
你的意思是^｛*｝？如果是这样的话，在《观星者》中是不可能的！！我已经试过了，所以我建议你像我在这里一样检查xtable包。
我将展示一种可以使用的方法：stargazer真的很好，即使使用尚未支持的模型对象，你也可以创建如上所述的表，例如，假设分位数回归模型不受stargazher支持(甚至思想是(：
技巧是，你需要能够获得系数和标准误差，例如作为矢量。然后为观星者提供支持的模型对象，例如lm作为模板，然后机械地指定应该使用的系数和标准误差：
library(stargazer)
library(tidyverse)
library(quantreg)
df <- mtcars
model1 <- lm(hp ~ factor(gear) + qsec + disp, data = df)
quantreg <- rq(hp ~ factor(gear) + qsec + disp, data = df)
summary_qr <- summary(quantreg, se = "boot")
# Standart Error for quant reg
se_qr = c(211.78266, 29.17307, 58.61105, 9.70908, 0.12090)
stargazer(model1, model1,
coef = list(NULL, summary_qr$coefficients),
se = list(NULL, se_qr),
type = "text")