我想问一下在做面板数据模型时:使用lm
函数的随机,固定和LSDV。当使用观星仪获得漂亮的结果表时,我们为所有横截面虚拟系数都获得了一个额外的变量系数。如下面的示例:
#Load packages
library(foreign)
library(plm)
library(stargazer)
#Load in Wooldridge data on crime
crime <- read.dta("http://fmwww.bc.edu/ec-p/data/wooldridge/crime4.dta")
#Declare our data to be a panel data set
crime.p <- pdata.frame(crime,index=c("county","year"))
#Run a panel model
#fixed effects / within
fixedeff <- plm(log(crmrte) ~ polpc + lwtuc + avgsen + wfed + d82 + d82 + d84 + d85,data=crime.p,model="within")
#Random effects
randomeff <- plm(log(crmrte) ~ polpc + lwtuc + avgsen + wfed + d82 + d82 + d84 + d85,data=crime.p,model="random")
#LSDV
LSDV <- lm(log(crmrte) ~ polpc + lwtuc + avgsen + wfed + d82 + d82 + d84 + d85 + factor(county)-1, data=crime.p)
stargazer(fixedeff, randomeff, LSDV, type = "text")
有没有一种方法来拥有一个漂亮的表并且其中没有所有的factor(county)
系数?
我知道我们可以在打印表格时“手动”完成此操作,但是我需要在函数内部执行此操作,因此我需要type = "text"
。
正如coeffeinjunky所说,它应该与omit
一起使用。将其作为选项添加到stargazer
:
stargazer(fixedeff, randomeff, LSDV, type = "text", omit=c("county"))
这是您要查找的输出吗?