Я пытаюсь воспроизвести хороший вывод модели stargazer (lm) для модели, не поддерживаемой stargazer.
можно ли выводить линейную модель stargazer вручную? Поскольку мы можем создать фрейм данных из каждой модели и затем вставить созданный фрейм данных в stargazer:
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
есть ли способ, как это сделать, как верхний индекс et c ...