library(wooldridge)
library(rmarkdown)
data("crime4")
library(plm)
paged_table(crime4)
pdim(crime4)
## Balanced Panel: n = 90, T = 7, N = 630
model1 <- plm(crmrte ~ prbarr + factor(year), data = crime4 , model = "within" )
summary(model1)
## Oneway (individual) effect Within Model
## 
## Call:
## plm(formula = crmrte ~ prbarr + factor(year), data = crime4, 
##     model = "within")
## 
## Balanced Panel: n = 90, T = 7, N = 630
## 
## Residuals:
##        Min.     1st Qu.      Median     3rd Qu.        Max. 
## -0.04395851 -0.00206428 -0.00011912  0.00192094  0.09339363 
## 
## Coefficients:
##                   Estimate  Std. Error t-value Pr(>|t|)   
## prbarr         -1.1164e-03  2.6306e-03 -0.4244 0.671450   
## factor(year)82 -8.8019e-05  1.0282e-03 -0.0856 0.931815   
## factor(year)83 -2.0101e-03  1.0306e-03 -1.9505 0.051642 . 
## factor(year)84 -3.2776e-03  1.0289e-03 -3.1855 0.001529 **
## factor(year)85 -2.9853e-03  1.0278e-03 -2.9046 0.003830 **
## factor(year)86 -4.6785e-04  1.0277e-03 -0.4552 0.649136   
## factor(year)87  7.5544e-04  1.0278e-03  0.7350 0.462640   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Total Sum of Squares:    0.026722
## Residual Sum of Squares: 0.025333
## R-Squared:      0.051981
## Adj. R-Squared: -0.11877
## F-statistic: 4.17503 on 7 and 533 DF, p-value: 0.00017221
model2 <- plm(crmrte ~ prbarr + factor(year), data = crime4 )
summary(model2)
## Oneway (individual) effect Within Model
## 
## Call:
## plm(formula = crmrte ~ prbarr + factor(year), data = crime4)
## 
## Balanced Panel: n = 90, T = 7, N = 630
## 
## Residuals:
##        Min.     1st Qu.      Median     3rd Qu.        Max. 
## -0.04395851 -0.00206428 -0.00011912  0.00192094  0.09339363 
## 
## Coefficients:
##                   Estimate  Std. Error t-value Pr(>|t|)   
## prbarr         -1.1164e-03  2.6306e-03 -0.4244 0.671450   
## factor(year)82 -8.8019e-05  1.0282e-03 -0.0856 0.931815   
## factor(year)83 -2.0101e-03  1.0306e-03 -1.9505 0.051642 . 
## factor(year)84 -3.2776e-03  1.0289e-03 -3.1855 0.001529 **
## factor(year)85 -2.9853e-03  1.0278e-03 -2.9046 0.003830 **
## factor(year)86 -4.6785e-04  1.0277e-03 -0.4552 0.649136   
## factor(year)87  7.5544e-04  1.0278e-03  0.7350 0.462640   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Total Sum of Squares:    0.026722
## Residual Sum of Squares: 0.025333
## R-Squared:      0.051981
## Adj. R-Squared: -0.11877
## F-statistic: 4.17503 on 7 and 533 DF, p-value: 0.00017221
model3 <- plm(crmrte ~ prbarr + factor(year)-1 , data = crime4 )
summary(model3)
## Oneway (individual) effect Within Model
## 
## Call:
## plm(formula = crmrte ~ prbarr + factor(year) - 1, data = crime4)
## 
## Balanced Panel: n = 90, T = 7, N = 630
## 
## Residuals:
##        Min.     1st Qu.      Median     3rd Qu.        Max. 
## -0.04395851 -0.00206428 -0.00011912  0.00192094  0.09339363 
## 
## Coefficients: (1 dropped because of singularities)
##                   Estimate  Std. Error t-value  Pr(>|t|)    
## prbarr         -0.00111643  0.00263064 -0.4244 0.6714501    
## factor(year)81 -0.00075544  0.00102776 -0.7350 0.4626403    
## factor(year)82 -0.00084346  0.00102861 -0.8200 0.4125819    
## factor(year)83 -0.00276553  0.00103135 -2.6815 0.0075573 ** 
## factor(year)84 -0.00403303  0.00102942 -3.9178 0.0001010 ***
## factor(year)85 -0.00374077  0.00102797 -3.6390 0.0003004 ***
## factor(year)86 -0.00122329  0.00102771 -1.1903 0.2344554    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Total Sum of Squares:    0.026722
## Residual Sum of Squares: 0.025333
## R-Squared:      0.051981
## Adj. R-Squared: -0.11877
## F-statistic: 4.17503 on 7 and 533 DF, p-value: 0.00017221
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:     
##                           -----------------------------
##                                      crmrte            
##                              (1)       (2)       (3)   
## -------------------------------------------------------
## prbarr                     -0.001    -0.001    -0.001  
##                            (0.003)   (0.003)   (0.003) 
##                                                        
## factor(year)81                                 -0.001  
##                                                (0.001) 
##                                                        
## factor(year)82             -0.0001   -0.0001   -0.001  
##                            (0.001)   (0.001)   (0.001) 
##                                                        
## factor(year)83             -0.002*   -0.002*  -0.003***
##                            (0.001)   (0.001)   (0.001) 
##                                                        
## factor(year)84            -0.003*** -0.003*** -0.004***
##                            (0.001)   (0.001)   (0.001) 
##                                                        
## factor(year)85            -0.003*** -0.003*** -0.004***
##                            (0.001)   (0.001)   (0.001) 
##                                                        
## factor(year)86             -0.0005   -0.0005   -0.001  
##                            (0.001)   (0.001)   (0.001) 
##                                                        
## factor(year)87              0.001     0.001            
##                            (0.001)   (0.001)           
##                                                        
## -------------------------------------------------------
## Observations                 630       630       630   
## R2                          0.052     0.052     0.052  
## Adjusted R2                -0.119    -0.119    -0.119  
## F Statistic (df = 7; 533) 4.175***  4.175***  4.175*** 
## =======================================================
## Note:                       *p<0.1; **p<0.05; ***p<0.01