#Load Ozone Dataset
load("Ozone Data with Region.RData")
combinedAir.final<-df4
combinedAir.final<-arrange(combinedAir.final,GEOID,Year)
combinedAir.final$USDM.categorical<-factor(combinedAir.final$USDM.categorical,levels=c("NoDrought","ModerateDrought","SevereDrought"))md2_random<-plm(Max.Ozone~USDM.categorical +elevation+ Longitude + Latitude+tdmean+tmean+ppt,
data =combinedAir.final,
index = c("GEOID","Year","month"),
model = "random")
md2_fix<-plm(Max.Ozone~USDM.categorical +elevation+ Longitude + Latitude+tdmean+tmean+ppt,
data =combinedAir.final,
index = c("GEOID","Year","month"),
model = "between")
stargazer(md2_random,md2_fix,type = "text", title = "All Regions", align = TRUE,
add.lines = list(c("Model Name", "Random Effect Model", "Fixed Effect Model")),
column.labels = c("Random Effect Model", "Fixed Effect Model"),
dep.var.caption = "Dependent Variable: Ozone Max",
dep.var.labels.include = FALSE,
omit.table.layout = "#")##
## All Regions
## ============================================================================
## Dependent Variable: Ozone Max
## --------------------------------------------
## Random Effect Model Fixed Effect Model
## ----------------------------------------------------------------------------
## USDM.categoricalModerateDrought 0.531*** 1.736
## (0.012) (1.079)
##
## USDM.categoricalSevereDrought 0.304*** 1.834
## (0.018) (1.614)
##
## elevation 0.008*** 0.007***
## (0.00003) (0.0003)
##
## Longitude 0.080*** 0.106***
## (0.008) (0.012)
##
## Latitude -0.052** 0.381***
## (0.020) (0.037)
##
## tdmean -0.827*** -0.165**
## (0.001) (0.076)
##
## tmean 1.346*** 1.067***
## (0.001) (0.070)
##
## ppt -0.136*** 0.272
## (0.001) (0.216)
##
## Constant 35.119*** 17.854***
## (1.028) (2.526)
##
## ----------------------------------------------------------------------------
## Model Name Random Effect Model Fixed Effect Model
## Observations 6,589,427 923
## R2 0.227 0.520
## Adjusted R2 0.227 0.515
## F Statistic 1,896,923.000*** 123.617*** (df = 8; 914)
## ============================================================================
## Note: *p<0.1; **p<0.05; ***p<0.01ne<- combinedAir.final%>%
filter(noaa_region=="northeast")
ne_random<-plm(Max.Ozone~USDM.categorical +elevation+ Longitude + Latitude+tdmean+tmean+ppt,
data =ne,
index = c("GEOID","Year","month"),
model = "random")
ne_fix<-plm(Max.Ozone~USDM.categorical +elevation+ Longitude + Latitude+tdmean+tmean+ppt,
data =ne,
index = c("GEOID","Year","month"),
model = "between")
stargazer(ne_random,ne_fix,type = "text", title = "Northeast", align = TRUE,
add.lines = list(c("Model Name", "Random Effect Model", "Fixed Effect Model")),
column.labels = c("Random Effect Model", "Fixed Effect Model"),
dep.var.caption = "Dependent Variable: Ozone Max",
dep.var.labels.include = FALSE,
omit.table.layout = "#")##
## Northeast
## ===========================================================================
## Dependent Variable: Ozone Max
## -------------------------------------------
## Random Effect Model Fixed Effect Model
## ---------------------------------------------------------------------------
## USDM.categoricalModerateDrought -0.404*** -5.784
## (0.043) (4.344)
##
## USDM.categoricalSevereDrought -1.368*** 39.585**
## (0.141) (18.575)
##
## elevation 0.017*** 0.009***
## (0.0002) (0.002)
##
## Longitude 1.473*** 0.272**
## (0.060) (0.119)
##
## Latitude -2.153*** -0.737***
## (0.099) (0.202)
##
## tdmean -1.325*** 1.893***
## (0.004) (0.317)
##
## tmean 1.923*** -0.891***
## (0.005) (0.322)
##
## ppt -0.062*** -0.971
## (0.002) (0.683)
##
## Constant 221.680*** 90.955***
## (7.495) (17.406)
##
## ---------------------------------------------------------------------------
## Model Name Random Effect Model Fixed Effect Model
## Observations 960,225 153
## R2 0.247 0.667
## Adjusted R2 0.247 0.649
## F Statistic 310,341.900*** 36.079*** (df = 8; 144)
## ===========================================================================
## Note: *p<0.1; **p<0.05; ***p<0.01nr<- combinedAir.final%>%
filter(noaa_region=="northern_rockies")
nr_random<-plm(Max.Ozone~USDM.categorical +elevation+ Longitude + Latitude+tdmean+tmean+ppt,
data =nr,
index = c("GEOID","Year","month"),
model = "random")
nr_fix<-plm(Max.Ozone~USDM.categorical +elevation+ Longitude + Latitude+tdmean+tmean+ppt,
data =nr,
index = c("GEOID","Year","month"),
model = "between")
stargazer(nr_random,nr_fix,type = "text", title = "Northern Rockies", align = TRUE,
add.lines = list(c("Model Name", "Random Effect Model", "Fixed Effect Model")),
column.labels = c("Random Effect Model", "Fixed Effect Model"),
dep.var.caption = "Dependent Variable: Ozone Max",
dep.var.labels.include = FALSE,
omit.table.layout = "#")##
## Northern Rockies
## ==========================================================================
## Dependent Variable: Ozone Max
## ------------------------------------------
## Random Effect Model Fixed Effect Model
## --------------------------------------------------------------------------
## USDM.categoricalModerateDrought 0.217*** 2.129
## (0.043) (2.128)
##
## USDM.categoricalSevereDrought -0.176*** 2.123
## (0.068) (4.972)
##
## elevation 0.007*** 0.010***
## (0.0002) (0.001)
##
## Longitude 0.117*** 0.521***
## (0.041) (0.117)
##
## Latitude 0.316*** 0.445*
## (0.072) (0.228)
##
## tdmean -0.475*** 0.091
## (0.005) (0.398)
##
## tmean 0.708*** 0.624**
## (0.004) (0.281)
##
## ppt -0.045*** -1.286
## (0.004) (0.987)
##
## Constant 25.545*** 59.968***
## (4.758) (10.058)
##
## --------------------------------------------------------------------------
## Model Name Random Effect Model Fixed Effect Model
## Observations 226,906 51
## R2 0.220 0.869
## Adjusted R2 0.220 0.844
## F Statistic 56,241.860*** 34.779*** (df = 8; 42)
## ==========================================================================
## Note: *p<0.1; **p<0.05; ***p<0.01nw<- combinedAir.final%>%
filter(noaa_region=="northwest")
nw_random<-plm(Max.Ozone~USDM.categorical +elevation+ Longitude + Latitude+tdmean+tmean+ppt,
data =nw,
index = c("GEOID","Year","month"),
model = "random")
nw_fix<-plm(Max.Ozone~USDM.categorical +elevation+ Longitude + Latitude+tdmean+tmean+ppt,
data =nw,
index = c("GEOID","Year","month"),
model = "between")
stargazer(nw_random,nw_fix,type = "text", title = "Northwest", align = TRUE,
add.lines = list(c("Model Name", "Random Effect Model", "Fixed Effect Model")),
column.labels = c("Random Effect Model", "Fixed Effect Model"),
dep.var.caption = "Dependent Variable: Ozone Max",
dep.var.labels.include = FALSE,
omit.table.layout = "#")##
## Northwest
## ==========================================================================
## Dependent Variable: Ozone Max
## ------------------------------------------
## Random Effect Model Fixed Effect Model
## --------------------------------------------------------------------------
## USDM.categoricalModerateDrought -0.325*** 1.924
## (0.074) (2.736)
##
## USDM.categoricalSevereDrought -1.848*** 5.122
## (0.171) (4.065)
##
## elevation 0.009*** 0.008***
## (0.0002) (0.002)
##
## Longitude -2.028*** 0.347
## (0.094) (0.239)
##
## Latitude -1.659*** 0.428
## (0.169) (0.302)
##
## tdmean -1.282*** -0.642
## (0.011) (0.437)
##
## tmean 1.358*** 1.553***
## (0.009) (0.249)
##
## ppt 0.156*** 0.512
## (0.005) (0.545)
##
## Constant -146.165*** 36.740
## (12.049) (34.239)
##
## --------------------------------------------------------------------------
## Model Name Random Effect Model Fixed Effect Model
## Observations 105,338 30
## R2 0.249 0.928
## Adjusted R2 0.249 0.901
## F Statistic 30,909.910*** 33.908*** (df = 8; 21)
## ==========================================================================
## Note: *p<0.1; **p<0.05; ***p<0.01ohv<- combinedAir.final%>%
filter(noaa_region=="ohio_valley")
ohv_random<-plm(Max.Ozone~USDM.categorical +elevation+ Longitude + Latitude+tdmean+tmean+ppt,
data =ohv,
index = c("GEOID","Year","month"),
model = "random")
ohv_fix<-plm(Max.Ozone~USDM.categorical +elevation+ Longitude + Latitude+tdmean+tmean+ppt,
data =ohv,
index = c("GEOID","Year","month"),
model = "between")
stargazer(ohv_random,ohv_fix,type = "text", title = "Ohio Valley", align = TRUE,
add.lines = list(c("Model Name", "Random Effect Model", "Fixed Effect Model")),
column.labels = c("Random Effect Model", "Fixed Effect Model"),
dep.var.caption = "Dependent Variable: Ozone Max",
dep.var.labels.include = FALSE,
omit.table.layout = "#")##
## Ohio Valley
## ===========================================================================
## Dependent Variable: Ozone Max
## -------------------------------------------
## Random Effect Model Fixed Effect Model
## ---------------------------------------------------------------------------
## USDM.categoricalModerateDrought 0.352*** 6.416***
## (0.038) (1.696)
##
## USDM.categoricalSevereDrought -0.658*** 9.057
## (0.090) (7.025)
##
## elevation 0.008*** 0.006***
## (0.0003) (0.002)
##
## Longitude 0.253*** 0.248***
## (0.034) (0.052)
##
## Latitude -0.256*** -0.072
## (0.062) (0.107)
##
## tdmean -1.578*** 2.677***
## (0.005) (0.468)
##
## tmean 2.199*** -1.255***
## (0.005) (0.414)
##
## ppt -0.145*** -1.975***
## (0.001) (0.507)
##
## Constant 54.996*** 64.189***
## (4.050) (7.800)
##
## ---------------------------------------------------------------------------
## Model Name Random Effect Model Fixed Effect Model
## Observations 1,006,480 185
## R2 0.284 0.733
## Adjusted R2 0.284 0.721
## F Statistic 379,350.300*** 60.382*** (df = 8; 176)
## ===========================================================================
## Note: *p<0.1; **p<0.05; ***p<0.01sth<- combinedAir.final%>%
filter(noaa_region=="south")
sth_random<-plm(Max.Ozone~USDM.categorical +elevation+ Longitude + Latitude+tdmean+tmean+ppt,
data =sth,
index = c("GEOID","Year","month"),
model = "random")
sth_fix<-plm(Max.Ozone~USDM.categorical +elevation+ Longitude + Latitude+tdmean+tmean+ppt,
data =sth,
index = c("GEOID","Year","month"),
model = "between")
stargazer(sth_random,sth_fix,type = "text", title = "South", align = TRUE,
add.lines = list(c("Model Name", "Random Effect Model", "Fixed Effect Model")),
column.labels = c("Random Effect Model", "Fixed Effect Model"),
dep.var.caption = "Dependent Variable: Ozone Max",
dep.var.labels.include = FALSE,
omit.table.layout = "#")##
## South
## ==========================================================================
## Dependent Variable: Ozone Max
## ------------------------------------------
## Random Effect Model Fixed Effect Model
## --------------------------------------------------------------------------
## USDM.categoricalModerateDrought -0.003 -0.863
## (0.032) (2.700)
##
## USDM.categoricalSevereDrought -0.425*** 5.270
## (0.047) (3.656)
##
## elevation 0.017*** 0.006*
## (0.001) (0.003)
##
## Longitude 1.084*** -0.037
## (0.061) (0.188)
##
## Latitude 0.650*** 1.025***
## (0.063) (0.161)
##
## tdmean -1.335*** -0.465
## (0.004) (0.329)
##
## tmean 1.765*** 1.370***
## (0.004) (0.331)
##
## ppt -0.080*** 1.751**
## (0.001) (0.802)
##
## Constant 102.535*** -23.524
## (6.134) (23.651)
##
## --------------------------------------------------------------------------
## Model Name Random Effect Model Fixed Effect Model
## Observations 951,486 128
## R2 0.219 0.398
## Adjusted R2 0.219 0.358
## F Statistic 254,749.200*** 9.849*** (df = 8; 119)
## ==========================================================================
## Note: *p<0.1; **p<0.05; ***p<0.01se<- combinedAir.final%>%
filter(noaa_region=="southeast")
se_random<-plm(Max.Ozone~USDM.categorical +elevation+ Longitude + Latitude+tdmean+tmean+ppt,
data =se,
index = c("GEOID","Year","month"),
model = "random")
se_fix<-plm(Max.Ozone~USDM.categorical +elevation+ Longitude + Latitude+tdmean+tmean+ppt,
data =se,
index = c("GEOID","Year","month"),
model = "between")
stargazer(se_random,se_fix,type = "text", title = "Southeast", align = TRUE,
add.lines = list(c("Model Name", "Random Effect Model", "Fixed Effect Model")),
column.labels = c("Random Effect Model", "Fixed Effect Model"),
dep.var.caption = "Dependent Variable: Ozone Max",
dep.var.labels.include = FALSE,
omit.table.layout = "#")##
## Southeast
## ===========================================================================
## Dependent Variable: Ozone Max
## -------------------------------------------
## Random Effect Model Fixed Effect Model
## ---------------------------------------------------------------------------
## USDM.categoricalModerateDrought 2.455*** -1.317
## (0.031) (2.726)
##
## USDM.categoricalSevereDrought 3.340*** 18.408***
## (0.046) (4.216)
##
## elevation 0.012*** 0.006***
## (0.0004) (0.001)
##
## Longitude -0.018 -0.102
## (0.052) (0.070)
##
## Latitude 0.456*** 1.124***
## (0.051) (0.083)
##
## tdmean -2.049*** 0.634*
## (0.005) (0.369)
##
## tmean 2.353*** 0.264
## (0.005) (0.372)
##
## ppt -0.080*** -0.860*
## (0.001) (0.487)
##
## Constant 6.062 -16.536**
## (5.174) (7.654)
##
## ---------------------------------------------------------------------------
## Model Name Random Effect Model Fixed Effect Model
## Observations 982,391 165
## R2 0.212 0.727
## Adjusted R2 0.212 0.713
## F Statistic 251,302.800*** 52.043*** (df = 8; 156)
## ===========================================================================
## Note: *p<0.1; **p<0.05; ***p<0.01sw<- combinedAir.final%>%
filter(noaa_region=="southwest")
sw_random<-plm(Max.Ozone~USDM.categorical +elevation+ Longitude + Latitude+tdmean+tmean+ppt,
data =sw,
index = c("GEOID","Year","month"),
model = "random")
sw_fix<-plm(Max.Ozone~USDM.categorical +elevation+ Longitude + Latitude+tdmean+tmean+ppt,
data =sw,
index = c("GEOID","Year","month"),
model = "between")
stargazer(sw_random,sw_fix,type = "text", title = "Southwest", align = TRUE,
add.lines = list(c("Model Name", "Random Effect Model", "Fixed Effect Model")),
column.labels = c("Random Effect Model", "Fixed Effect Model"),
dep.var.caption = "Dependent Variable: Ozone Max",
dep.var.labels.include = FALSE,
omit.table.layout = "#")##
## Southwest
## ==========================================================================
## Dependent Variable: Ozone Max
## ------------------------------------------
## Random Effect Model Fixed Effect Model
## --------------------------------------------------------------------------
## USDM.categoricalModerateDrought 0.184*** 4.011
## (0.028) (2.874)
##
## USDM.categoricalSevereDrought 1.184*** -6.469***
## (0.036) (2.377)
##
## elevation 0.010*** 0.007***
## (0.0001) (0.001)
##
## Longitude 1.004*** -0.397**
## (0.053) (0.175)
##
## Latitude 0.821*** -0.048
## (0.059) (0.160)
##
## tdmean -0.359*** 0.097
## (0.002) (0.263)
##
## tmean 0.961*** 0.631***
## (0.002) (0.195)
##
## ppt 0.021*** 0.911
## (0.004) (0.878)
##
## Constant 95.789*** -14.990
## (6.322) (15.261)
##
## --------------------------------------------------------------------------
## Model Name Random Effect Model Fixed Effect Model
## Observations 683,530 66
## R2 0.353 0.647
## Adjusted R2 0.353 0.598
## F Statistic 360,596.300*** 13.079*** (df = 8; 57)
## ==========================================================================
## Note: *p<0.1; **p<0.05; ***p<0.01um<- combinedAir.final%>%
filter(noaa_region=="upper_midwest")
um_random<-plm(Max.Ozone~USDM.categorical +elevation+ Longitude + Latitude+tdmean+tmean+ppt,
data =um,
index = c("GEOID","Year","month"),
model = "random")
um_fix<-plm(Max.Ozone~USDM.categorical +elevation+ Longitude + Latitude+tdmean+tmean+ppt,
data =um,
index = c("GEOID","Year","month"),
model = "between")
stargazer(um_random,um_fix,type = "text", title = "Upper Midwest", align = TRUE,
add.lines = list(c("Model Name", "Random Effect Model", "Fixed Effect Model")),
column.labels = c("Random Effect Model", "Fixed Effect Model"),
dep.var.caption = "Dependent Variable: Ozone Max",
dep.var.labels.include = FALSE,
omit.table.layout = "#")##
## Upper Midwest
## ==========================================================================
## Dependent Variable: Ozone Max
## ------------------------------------------
## Random Effect Model Fixed Effect Model
## --------------------------------------------------------------------------
## USDM.categoricalModerateDrought -0.361*** -0.731
## (0.056) (4.874)
##
## USDM.categoricalSevereDrought -1.657*** -2.016
## (0.188) (21.967)
##
## elevation 0.021*** 0.006
## (0.002) (0.004)
##
## Longitude 0.877*** 0.411***
## (0.058) (0.106)
##
## Latitude 1.881*** -0.089
## (0.105) (0.214)
##
## tdmean -1.417*** 1.562***
## (0.007) (0.550)
##
## tmean 1.752*** -0.554
## (0.007) (0.479)
##
## ppt -0.081*** -0.710
## (0.002) (1.128)
##
## Constant 19.250*** 77.189***
## (6.189) (11.747)
##
## --------------------------------------------------------------------------
## Model Name Random Effect Model Fixed Effect Model
## Observations 382,355 86
## R2 0.187 0.677
## Adjusted R2 0.187 0.643
## F Statistic 84,727.550*** 20.159*** (df = 8; 77)
## ==========================================================================
## Note: *p<0.1; **p<0.05; ***p<0.01wst<- combinedAir.final%>%
filter(noaa_region=="upper_midwest")
wst_random<-plm(Max.Ozone~USDM.categorical +elevation+ Longitude + Latitude+tdmean+tmean+ppt,
data =wst,
index = c("GEOID","Year","month"),
model = "random")
wst_fix<-plm(Max.Ozone~USDM.categorical +elevation+ Longitude + Latitude+tdmean+tmean+ppt,
data =wst,
index = c("GEOID","Year","month"),
model = "between")
stargazer(wst_random,wst_fix,type = "text", title = "West", align = TRUE,
add.lines = list(c("Model Name", "Random Effect Model", "Fixed Effect Model")),
column.labels = c("Random Effect Model", "Fixed Effect Model"),
dep.var.caption = "Dependent Variable: Ozone Max",
dep.var.labels.include = FALSE,
omit.table.layout = "#")##
## West
## ==========================================================================
## Dependent Variable: Ozone Max
## ------------------------------------------
## Random Effect Model Fixed Effect Model
## --------------------------------------------------------------------------
## USDM.categoricalModerateDrought -0.361*** -0.731
## (0.056) (4.874)
##
## USDM.categoricalSevereDrought -1.657*** -2.016
## (0.188) (21.967)
##
## elevation 0.021*** 0.006
## (0.002) (0.004)
##
## Longitude 0.877*** 0.411***
## (0.058) (0.106)
##
## Latitude 1.881*** -0.089
## (0.105) (0.214)
##
## tdmean -1.417*** 1.562***
## (0.007) (0.550)
##
## tmean 1.752*** -0.554
## (0.007) (0.479)
##
## ppt -0.081*** -0.710
## (0.002) (1.128)
##
## Constant 19.250*** 77.189***
## (6.189) (11.747)
##
## --------------------------------------------------------------------------
## Model Name Random Effect Model Fixed Effect Model
## Observations 382,355 86
## R2 0.187 0.677
## Adjusted R2 0.187 0.643
## F Statistic 84,727.550*** 20.159*** (df = 8; 77)
## ==========================================================================
## Note: *p<0.1; **p<0.05; ***p<0.01