load("Ozone Data with Region.RData")
##Sampling 30%
df4$USDM.categorical<-factor(df4$USDM.categorical,levels=c("NoDrought","ModerateDrought","SevereDrought"))
df4$modest<-ifelse(df4$USDM.categorical=="ModerateDrought",1,0)
df4$severe<-ifelse(df4$USDM.categorical=="SevereDrought",1,0)
df5<-df4%>%
group_by(USDM.categorical,Year,month)%>%
sample_frac(0.3)
#felm Overall Model
ov05<-felm(Max.Ozone~modest+severe+elevation+Longitude+Latitude |GEOID + Year | 0 |GEOID,
data=df5)
ran05<-plm(Max.Ozone~modest+severe+elevation+Longitude+Latitude,
data=df5,
index=c("GEOID","Year","month"),
model="random")
stargazer(ov05,ran05,type="text",title = "All Region",align = TRUE,
column.labels = c("FELM Model","Random Effect Model"),
add.lines = list(c("Model Name","FELM Model","Random Effect Model")),
dep.var.caption = "Dependent Variable: Max Ozone",
dep.var.labels.include = FALSE,
omit.table.layout = "#")##
## All Region
## =============================================================
## Dependent Variable: Max Ozone
## -----------------------------------------
## felm panel
## linear
## FELM Model Random Effect Model
## -------------------------------------------------------------
## modest 1.424*** 1.310***
## (0.161) (0.025)
##
## severe 3.081*** 2.610***
## (0.243) (0.037)
##
## elevation 0.007*** 0.006***
## (0.001) (0.0001)
##
## Longitude 0.287 0.126***
## (0.836) (0.010)
##
## Latitude 0.470 -0.134***
## (1.315) (0.027)
##
## Constant 56.221***
## (1.342)
##
## -------------------------------------------------------------
## Model Name FELM Model Random Effect Model
## Observations 1,976,817 1,976,817
## R2 0.127 0.018
## Adjusted R2 0.126 0.018
## Residual Std. Error 13.609 (df = 1975871)
## F Statistic 15,382.490***
## =============================================================
## Note: *p<0.1; **p<0.05; ***p<0.01region<-df5$noaa_region%>%unique()%>%sort()
# Region Lopping
for(i in 1:length(region)){
df4_region<-df4%>%
filter(noaa_region==region[i])
felm_region<-felm(Max.Ozone~modest+severe+elevation+Longitude+Latitude |GEOID + Year | 0 |GEOID,
data=df4_region)
ran_model<-plm(Max.Ozone~modest+severe+elevation+Longitude+Latitude,
data=df4_region,
index=c("GEOID","Year","month"),
model="random")
stargazer(felm_region,ran_model,type="text",title = paste0(region[i]," ","Region"),align = TRUE,
column.labels = c("FELM Model","Random Effect Model"),
add.lines = list(c("Model Name","FELM Model","Random Effect Model")),
dep.var.caption = "Dependent Variable: Max Ozone",
dep.var.labels.include = FALSE,omit.table.layout = "#")
}##
## northeast Region
## ============================================================
## Dependent Variable: Max Ozone
## ----------------------------------------
## felm panel
## linear
## FELM Model Random Effect Model
## ------------------------------------------------------------
## modest -0.495* 0.497***
## (0.289) (0.050)
##
## severe -4.222*** -1.826***
## (0.540) (0.162)
##
## elevation 0.013*** 0.011***
## (0.001) (0.0002)
##
## Longitude 4.292** 1.855***
## (1.746) (0.087)
##
## Latitude -4.644** -3.442***
## (2.245) (0.141)
##
## Constant 319.751***
## (10.608)
##
## ------------------------------------------------------------
## Model Name FELM Model Random Effect Model
## Observations 960,225 960,225
## R2 0.089 0.005
## Adjusted R2 0.089 0.005
## Residual Std. Error 13.985 (df = 960049)
## F Statistic 3,234.418***
## ============================================================
## Note: *p<0.1; **p<0.05; ***p<0.01
##
## northernRegion
## ===========================================================
## Dependent Variable: Max Ozone
## ---------------------------------------
## felm panel
## linear
## FELM Model Random Effect Model
## -----------------------------------------------------------
## modest 0.180 0.512***
## (0.412) (0.048)
##
## severe -0.018 0.803***
## (0.623) (0.076)
##
## elevation 0.005*** 0.005***
## (0.001) (0.0002)
##
## Longitude -1.073 -0.111**
## (0.810) (0.045)
##
## Latitude 0.705 -0.270***
## (0.828) (0.080)
##
## Constant 35.327***
## (5.288)
##
## -----------------------------------------------------------
## Model Name FELM Model Random Effect Model
## Observations 226,906 226,906
## R2 0.180 0.032
## Adjusted R2 0.180 0.032
## Residual Std. Error 9.528 (df = 226832)
## F Statistic 1,181.555***
## ===========================================================
## Note: *p<0.1; **p<0.05; ***p<0.01
##
## northwest Region
## ============================================================
## Dependent Variable: Max Ozone
## ----------------------------------------
## felm panel
## linear
## FELM Model Random Effect Model
## ------------------------------------------------------------
## modest 0.228 0.471***
## (0.430) (0.082)
##
## severe -2.230*** -0.993***
## (0.708) (0.190)
##
## elevation 0.012*** 0.011***
## (0.002) (0.0002)
##
## Longitude -4.407 -2.543***
## (2.901) (0.137)
##
## Latitude -3.378 -2.535***
## (5.914) (0.277)
##
## Constant -157.097***
## (18.929)
##
## ------------------------------------------------------------
## Model Name FELM Model Random Effect Model
## Observations 105,338 105,338
## R2 0.268 0.046
## Adjusted R2 0.268 0.046
## Residual Std. Error 10.764 (df = 105285)
## F Statistic 4,064.467***
## ============================================================
## Note: *p<0.1; **p<0.05; ***p<0.01
##
## ohioRegion
## =============================================================
## Dependent Variable: Max Ozone
## -----------------------------------------
## felm panel
## linear
## FELM Model Random Effect Model
## -------------------------------------------------------------
## modest 0.153 1.872***
## (0.279) (0.045)
##
## severe 0.627 3.628***
## (0.854) (0.105)
##
## elevation 0.004** 0.005***
## (0.002) (0.0003)
##
## Longitude 3.304** 0.627***
## (1.436) (0.069)
##
## Latitude -1.994 -0.917***
## (1.698) (0.117)
##
## Constant 132.649***
## (8.026)
##
## -------------------------------------------------------------
## Model Name FELM Model Random Effect Model
## Observations 1,006,480 1,006,480
## R2 0.112 0.008
## Adjusted R2 0.111 0.008
## Residual Std. Error 13.729 (df = 1006272)
## F Statistic 3,183.263***
## =============================================================
## Note: *p<0.1; **p<0.05; ***p<0.01
##
## south Region
## ============================================================
## Dependent Variable: Max Ozone
## ----------------------------------------
## felm panel
## linear
## FELM Model Random Effect Model
## ------------------------------------------------------------
## modest 0.961*** 1.363***
## (0.251) (0.035)
##
## severe 2.183*** 2.938***
## (0.358) (0.052)
##
## elevation 0.014*** 0.014***
## (0.004) (0.001)
##
## Longitude -0.014 0.450***
## (0.808) (0.069)
##
## Latitude 3.482*** 0.500***
## (1.272) (0.070)
##
## Constant 63.687***
## (6.890)
##
## ------------------------------------------------------------
## Model Name FELM Model Random Effect Model
## Observations 951,486 951,486
## R2 0.062 0.015
## Adjusted R2 0.062 0.015
## Residual Std. Error 13.963 (df = 951335)
## F Statistic 4,518.735***
## ============================================================
## Note: *p<0.1; **p<0.05; ***p<0.01
##
## southeast Region
## ============================================================
## Dependent Variable: Max Ozone
## ----------------------------------------
## felm panel
## linear
## FELM Model Random Effect Model
## ------------------------------------------------------------
## modest 2.831*** 3.868***
## (0.194) (0.034)
##
## severe 3.592*** 5.240***
## (0.215) (0.051)
##
## elevation 0.009*** 0.006***
## (0.003) (0.0004)
##
## Longitude -0.026 -0.043
## (1.657) (0.060)
##
## Latitude -0.398 0.651***
## (2.475) (0.060)
##
## Constant 15.664***
## (6.029)
##
## ------------------------------------------------------------
## Model Name FELM Model Random Effect Model
## Observations 982,391 982,391
## R2 0.129 0.030
## Adjusted R2 0.128 0.030
## Residual Std. Error 12.919 (df = 982203)
## F Statistic 20,731.880***
## ============================================================
## Note: *p<0.1; **p<0.05; ***p<0.01
##
## southwest Region
## ============================================================
## Dependent Variable: Max Ozone
## ----------------------------------------
## felm panel
## linear
## FELM Model Random Effect Model
## ------------------------------------------------------------
## modest 0.089 -0.444***
## (0.363) (0.034)
##
## severe 4.295*** 3.715***
## (0.729) (0.043)
##
## elevation 0.008*** 0.008***
## (0.001) (0.0001)
##
## Longitude 1.148 0.255***
## (1.111) (0.063)
##
## Latitude 1.668 0.510***
## (1.151) (0.069)
##
## Constant 42.653***
## (7.484)
##
## ------------------------------------------------------------
## Model Name FELM Model Random Effect Model
## Observations 683,530 683,530
## R2 0.075 0.045
## Adjusted R2 0.075 0.045
## Residual Std. Error 12.014 (df = 683441)
## F Statistic 22,580.700***
## ============================================================
## Note: *p<0.1; **p<0.05; ***p<0.01
##
## upperRegion
## ============================================================
## Dependent Variable: Max Ozone
## ----------------------------------------
## felm panel
## linear
## FELM Model Random Effect Model
## ------------------------------------------------------------
## modest -1.140*** 0.149**
## (0.263) (0.062)
##
## severe -2.467*** 0.283
## (0.843) (0.207)
##
## elevation 0.001 0.029***
## (0.019) (0.002)
##
## Longitude 11.215** 1.245***
## (4.279) (0.080)
##
## Latitude 6.775*** 1.309***
## (2.192) (0.130)
##
## Constant 87.542***
## (8.051)
##
## ------------------------------------------------------------
## Model Name FELM Model Random Effect Model
## Observations 382,355 382,355
## R2 0.092 0.006
## Adjusted R2 0.092 0.006
## Residual Std. Error 12.449 (df = 382246)
## F Statistic 646.255***
## ============================================================
## Note: *p<0.1; **p<0.05; ***p<0.01
##
## west Region
## =============================================================
## Dependent Variable: Max Ozone
## -----------------------------------------
## felm panel
## linear
## FELM Model Random Effect Model
## -------------------------------------------------------------
## modest 1.171** 1.056***
## (0.501) (0.031)
##
## severe 1.449** 0.993***
## (0.616) (0.036)
##
## elevation 0.006*** 0.006***
## (0.002) (0.0001)
##
## Longitude 0.283 0.305***
## (1.346) (0.044)
##
## Latitude 0.109 0.036
## (2.385) (0.069)
##
## Constant 74.249***
## (5.543)
##
## -------------------------------------------------------------
## Model Name FELM Model Random Effect Model
## Observations 1,290,716 1,290,716
## R2 0.160 0.013
## Adjusted R2 0.160 0.013
## Residual Std. Error 15.003 (df = 1290634)
## F Statistic 15,613.030***
## =============================================================
## Note: *p<0.1; **p<0.05; ***p<0.01