library(wooldridge)
## Warning: package 'wooldridge' was built under R version 4.1.3
library(rmarkdown)
data("murder")
paged_table(murder)
library(plm)
## Warning: package 'plm' was built under R version 4.1.3
Murder2<-pdata.frame(murder, index = c("year", "state"))
pdim(Murder2)
## Balanced Panel: n = 3, T = 51, N = 153
pooling model (OLS)
pooling_modle1 <- plm(mrdrte~exec+unem, data = Murder2, model = "pooling")
summary(pooling_modle1)
## Pooling Model
##
## Call:
## plm(formula = mrdrte ~ exec + unem, data = Murder2, model = "pooling")
##
## Balanced Panel: n = 3, T = 51, N = 153
##
## Residuals:
## Min. 1st Qu. Median 3rd Qu. Max.
## -9.1747 -3.4721 -1.4156 1.1137 69.1431
##
## Coefficients:
## Estimate Std. Error t-value Pr(>|t|)
## (Intercept) 0.34812 2.68724 0.1295 0.897100
## exec 0.16502 0.19387 0.8512 0.396008
## unem 1.25891 0.43736 2.8784 0.004581 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Total Sum of Squares: 12845
## Residual Sum of Squares: 12046
## R-Squared: 0.062264
## Adj. R-Squared: 0.04976
## F-statistic: 4.97983 on 2 and 150 DF, p-value: 0.0080549
Fixed effect modeli (within)
within_modeli <- plm(mrdrte~unem+exec, data = Murder2,model = "within")
summary(within_modeli)
## Oneway (individual) effect Within Model
##
## Call:
## plm(formula = mrdrte ~ unem + exec, data = Murder2, model = "within")
##
## Balanced Panel: n = 3, T = 51, N = 153
##
## Residuals:
## Min. 1st Qu. Median 3rd Qu. Max.
## -9.1301 -3.1194 -1.2107 1.3794 67.8099
##
## Coefficients:
## Estimate Std. Error t-value Pr(>|t|)
## unem 1.39079 0.45087 3.0847 0.002432 **
## exec 0.16275 0.19393 0.8392 0.402685
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Total Sum of Squares: 12762
## Residual Sum of Squares: 11868
## R-Squared: 0.070092
## Adj. R-Squared: 0.044959
## F-statistic: 5.57773 on 2 and 148 DF, p-value: 0.0046193
stargezer
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(pooling_modle1,within_modeli, type = "text")
##
## ==========================================================
## Dependent variable:
## ---------------------------------------------
## mrdrte
## (1) (2)
## ----------------------------------------------------------
## exec 0.165 0.163
## (0.194) (0.194)
##
## unem 1.259*** 1.391***
## (0.437) (0.451)
##
## Constant 0.348
## (2.687)
##
## ----------------------------------------------------------
## Observations 153 153
## R2 0.062 0.070
## Adjusted R2 0.050 0.045
## F Statistic 4.980*** (df = 2; 150) 5.578*** (df = 2; 148)
## ==========================================================
## Note: *p<0.1; **p<0.05; ***p<0.01
Normal Regresyon
normal<- lm(mrdrte~unem+exec, data = Murder2)
summary(normal)
##
## Call:
## lm(formula = mrdrte ~ unem + exec, data = Murder2)
##
## Residuals:
## Min 1Q Median 3Q Max
## -9.175 -3.472 -1.416 1.114 69.143
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.3481 2.6872 0.130 0.89710
## unem 1.2589 0.4374 2.878 0.00458 **
## exec 0.1650 0.1939 0.851 0.39601
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 8.961 on 150 degrees of freedom
## Multiple R-squared: 0.06226, Adjusted R-squared: 0.04976
## F-statistic: 4.98 on 2 and 150 DF, p-value: 0.008055