library(wooldridge)
library(rmarkdown)
data("prison")
library(plm)
paged_table(prison)
pdim(prison)
## Balanced Panel: n = 51, T = 14, N = 714
model1 <- plm(black ~ unem + factor(year), data = prison , model = "within" )
summary(model1)
## Oneway (individual) effect Within Model
##
## Call:
## plm(formula = black ~ unem + factor(year), data = prison, model = "within")
##
## Balanced Panel: n = 51, T = 14, N = 714
##
## Residuals:
## Min. 1st Qu. Median 3rd Qu. Max.
## -3.1294e-02 -8.3058e-04 -7.2992e-06 8.8317e-04 3.1236e-02
##
## Coefficients:
## Estimate Std. Error t-value Pr(>|t|)
## unem 0.00954518 0.01297528 0.7356 0.4622135
## factor(year)81 0.00025240 0.00074501 0.3388 0.7348770
## factor(year)82 0.00036440 0.00080646 0.4519 0.6515240
## factor(year)83 0.00066001 0.00080775 0.8171 0.4141763
## factor(year)84 0.00114655 0.00074501 1.5390 0.1242966
## factor(year)85 0.00146273 0.00074336 1.9677 0.0495247 *
## factor(year)86 0.00177265 0.00074272 2.3867 0.0172863 *
## factor(year)87 0.00213930 0.00074591 2.8680 0.0042645 **
## factor(year)88 0.00251437 0.00076318 3.2946 0.0010393 **
## factor(year)89 0.00284322 0.00077396 3.6736 0.0002589 ***
## factor(year)90 0.00311748 0.00076545 4.0727 5.221e-05 ***
## factor(year)91 0.00332039 0.00074461 4.4592 9.692e-06 ***
## factor(year)92 0.00357968 0.00074243 4.8215 1.776e-06 ***
## factor(year)93 0.00392563 0.00074567 5.2646 1.913e-07 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Total Sum of Squares: 0.01016
## Residual Sum of Squares: 0.0091219
## R-Squared: 0.10219
## Adj. R-Squared: 0.01365
## F-statistic: 5.27621 on 14 and 649 DF, p-value: 1.5045e-09
model2 <- plm(black ~ unem + factor(year), data = prison )
summary(model2)
## Oneway (individual) effect Within Model
##
## Call:
## plm(formula = black ~ unem + factor(year), data = prison)
##
## Balanced Panel: n = 51, T = 14, N = 714
##
## Residuals:
## Min. 1st Qu. Median 3rd Qu. Max.
## -3.1294e-02 -8.3058e-04 -7.2992e-06 8.8317e-04 3.1236e-02
##
## Coefficients:
## Estimate Std. Error t-value Pr(>|t|)
## unem 0.00954518 0.01297528 0.7356 0.4622135
## factor(year)81 0.00025240 0.00074501 0.3388 0.7348770
## factor(year)82 0.00036440 0.00080646 0.4519 0.6515240
## factor(year)83 0.00066001 0.00080775 0.8171 0.4141763
## factor(year)84 0.00114655 0.00074501 1.5390 0.1242966
## factor(year)85 0.00146273 0.00074336 1.9677 0.0495247 *
## factor(year)86 0.00177265 0.00074272 2.3867 0.0172863 *
## factor(year)87 0.00213930 0.00074591 2.8680 0.0042645 **
## factor(year)88 0.00251437 0.00076318 3.2946 0.0010393 **
## factor(year)89 0.00284322 0.00077396 3.6736 0.0002589 ***
## factor(year)90 0.00311748 0.00076545 4.0727 5.221e-05 ***
## factor(year)91 0.00332039 0.00074461 4.4592 9.692e-06 ***
## factor(year)92 0.00357968 0.00074243 4.8215 1.776e-06 ***
## factor(year)93 0.00392563 0.00074567 5.2646 1.913e-07 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Total Sum of Squares: 0.01016
## Residual Sum of Squares: 0.0091219
## R-Squared: 0.10219
## Adj. R-Squared: 0.01365
## F-statistic: 5.27621 on 14 and 649 DF, p-value: 1.5045e-09
model3 <- plm(black ~ unem + factor(year)-1 , data = prison )
summary(model3)
## Oneway (individual) effect Within Model
##
## Call:
## plm(formula = black ~ unem + factor(year) - 1, data = prison)
##
## Balanced Panel: n = 51, T = 14, N = 714
##
## Residuals:
## Min. 1st Qu. Median 3rd Qu. Max.
## -3.1294e-02 -8.3058e-04 -7.2992e-06 8.8317e-04 3.1236e-02
##
## Coefficients: (1 dropped because of singularities)
## Estimate Std. Error t-value Pr(>|t|)
## unem 0.00954518 0.01297528 0.7356 0.4622135
## factor(year)80 -0.00392563 0.00074567 -5.2646 1.913e-07 ***
## factor(year)81 -0.00367323 0.00075398 -4.8718 1.391e-06 ***
## factor(year)82 -0.00356123 0.00083605 -4.2596 2.351e-05 ***
## factor(year)83 -0.00326563 0.00083757 -3.8989 0.0001067 ***
## factor(year)84 -0.00277908 0.00075398 -3.6859 0.0002469 ***
## factor(year)85 -0.00246290 0.00075007 -3.2836 0.0010801 **
## factor(year)86 -0.00215298 0.00074794 -2.8786 0.0041266 **
## factor(year)87 -0.00178633 0.00074243 -2.4061 0.0164037 *
## factor(year)88 -0.00141126 0.00075014 -1.8813 0.0603731 .
## factor(year)89 -0.00108242 0.00075726 -1.4294 0.1533744
## factor(year)90 -0.00080815 0.00075156 -1.0753 0.2826405
## factor(year)91 -0.00060525 0.00074253 -0.8151 0.4153027
## factor(year)92 -0.00034595 0.00074527 -0.4642 0.6426646
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Total Sum of Squares: 0.01016
## Residual Sum of Squares: 0.0091219
## R-Squared: 0.10219
## Adj. R-Squared: 0.01365
## F-statistic: 5.27621 on 14 and 649 DF, p-value: 1.5045e-09
library(stargazer)
##
## Please cite as:
## Hlavac, Marek (2022). stargazer: Well-Formatted Regression and Summary Statistics Tables.
## R package version 5.2.3. https://CRAN.R-project.org/package=stargazer
stargazer(list(model1 , model2 ,model3) , type = "text" )
##
## ========================================================
## Dependent variable:
## -----------------------------
## black
## (1) (2) (3)
## --------------------------------------------------------
## unem 0.010 0.010 0.010
## (0.013) (0.013) (0.013)
##
## factor(year)80 -0.004***
## (0.001)
##
## factor(year)81 0.0003 0.0003 -0.004***
## (0.001) (0.001) (0.001)
##
## factor(year)82 0.0004 0.0004 -0.004***
## (0.001) (0.001) (0.001)
##
## factor(year)83 0.001 0.001 -0.003***
## (0.001) (0.001) (0.001)
##
## factor(year)84 0.001 0.001 -0.003***
## (0.001) (0.001) (0.001)
##
## factor(year)85 0.001** 0.001** -0.002***
## (0.001) (0.001) (0.001)
##
## factor(year)86 0.002** 0.002** -0.002***
## (0.001) (0.001) (0.001)
##
## factor(year)87 0.002*** 0.002*** -0.002**
## (0.001) (0.001) (0.001)
##
## factor(year)88 0.003*** 0.003*** -0.001*
## (0.001) (0.001) (0.001)
##
## factor(year)89 0.003*** 0.003*** -0.001
## (0.001) (0.001) (0.001)
##
## factor(year)90 0.003*** 0.003*** -0.001
## (0.001) (0.001) (0.001)
##
## factor(year)91 0.003*** 0.003*** -0.001
## (0.001) (0.001) (0.001)
##
## factor(year)92 0.004*** 0.004*** -0.0003
## (0.001) (0.001) (0.001)
##
## factor(year)93 0.004*** 0.004***
## (0.001) (0.001)
##
## --------------------------------------------------------
## Observations 714 714 714
## R2 0.102 0.102 0.102
## Adjusted R2 0.014 0.014 0.014
## F Statistic (df = 14; 649) 5.276*** 5.276*** 5.276***
## ========================================================
## Note: *p<0.1; **p<0.05; ***p<0.01