i使用stargazer
软件包制作面板型号的乳胶表(plm
软件包)。如果第一个plm
模型是单变量,则stargazer
会掉落该模型的回归器,但仍显示模型。
如果我交换模型顺序,以使第一个模型是多变量的,则所有模型显示。我没有lm
型号的问题。有工作吗?lm
包括截距,而plm
不包括?
谢谢!这是一个例子。
dat <- data.frame(y=runif(1000),
x1=runif(1000),
x2=runif(1000),
id=rep(1:50, each=20),
time=rep(1:20, times=50)
)
pdat <- pdata.frame(dat, index=c("id", "time"))
plms <- list()
plms[[1]] <- plm(y ~ x1, data=pdat)
plms[[2]] <- plm(y ~ x1 + x2, data=pdat)
plms[[3]] <- plm(y ~ x2, data=pdat)
# drops first regressor if first model is univariate
stargazer(plms[1:3], type="text")
# but not if first model is multivariate
stargazer(plms[c(2, 1, 3)], type="text")
这产生以下。
> # drops first regressor if first model is univariate
> stargazer(plms[1:3], type="text")
========================================================================
Dependent variable:
-----------------------------------------------------------
y
(1) (2) (3)
------------------------------------------------------------------------
x1 0.004
(0.032)
x2 -0.008 -0.008
(0.032) (0.032)
------------------------------------------------------------------------
Observations 1,000 1,000 1,000
R2 0.00002 0.0001 0.0001
Adjusted R2 0.00002 0.0001 0.0001
F Statistic 0.016 (df = 1; 949) 0.038 (df = 2; 948) 0.059 (df = 1; 949)
========================================================================
Note: *p<0.1; **p<0.05; ***p<0.01
> # but not if first model is multivariate
> stargazer(plms[c(2, 1, 3)], type="text")
========================================================================
Dependent variable:
-----------------------------------------------------------
y
(1) (2) (3)
------------------------------------------------------------------------
x1 0.004 0.004
(0.032) (0.032)
x2 -0.008 -0.008
(0.032) (0.032)
------------------------------------------------------------------------
Observations 1,000 1,000 1,000
R2 0.0001 0.00002 0.0001
Adjusted R2 0.0001 0.00002 0.0001
F Statistic 0.038 (df = 2; 948) 0.016 (df = 1; 949) 0.059 (df = 1; 949)
========================================================================
Note: *p<0.1; **p<0.05; ***p<0.01
此问题已在stargazer
(版本5.0)的最新版本中解决,现在可在Cran上获得。