All installed packages were loaded to be used in the R Program.
library(plm)
library(knitr)
library(tidyverse)## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr 1.1.1 ✔ readr 2.1.4
## ✔ forcats 1.0.0 ✔ stringr 1.5.0
## ✔ ggplot2 3.4.2 ✔ tibble 3.2.1
## ✔ lubridate 1.9.2 ✔ tidyr 1.3.0
## ✔ purrr 1.0.1
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::between() masks plm::between()
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag() masks plm::lag(), stats::lag()
## ✖ dplyr::lead() masks plm::lead()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errorslibrary(lubridate)
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=stargazertest.dat <- load("Ozone_Drought_Final.RData")
df <- combinedAir.final
regions <- read.csv("region_code.csv")
df1 <- df %>%
merge(regions, by = "State.Code")
df2 <- df1 %>%
mutate(no_drought = ifelse(USDM.categorical == "NoDrought", 1, 0),
moderate = ifelse(USDM.categorical == "ModerateDrought", 1, 0),
severe = ifelse(USDM.categorical == "SevereDrought", 1, 0))md_reg <- lm(Max.Ozone ~ moderate + severe,
data = df2)
stargazer(md_reg, type = "text", title = "All Regions", align = TRUE)##
## All Regions
## ===================================================
## Dependent variable:
## -------------------------------
## Max.Ozone
## ---------------------------------------------------
## moderate 1.924***
## (0.017)
##
## severe 3.994***
## (0.027)
##
## Constant 41.280***
## (0.007)
##
## ---------------------------------------------------
## Observations 4,799,111
## R2 0.006
## Adjusted R2 0.006
## Residual Std. Error 14.074 (df = 4799108)
## F Statistic 15,420.680*** (df = 2; 4799108)
## ===================================================
## Note: *p<0.1; **p<0.05; ***p<0.01ne <- df2 %>%
filter(noaa_region == "northeast")
ne_reg <- lm(Max.Ozone ~ moderate + severe,
data = ne)
stargazer(ne_reg, type = "text", title = "Northeast", align = TRUE)##
## Northeast
## ===============================================
## Dependent variable:
## ---------------------------
## Max.Ozone
## -----------------------------------------------
## moderate 0.328***
## (0.052)
##
## severe -1.027***
## (0.172)
##
## Constant 40.068***
## (0.016)
##
## -----------------------------------------------
## Observations 914,988
## R2 0.0001
## Adjusted R2 0.0001
## Residual Std. Error 14.598 (df = 914985)
## F Statistic 38.647*** (df = 2; 914985)
## ===============================================
## Note: *p<0.1; **p<0.05; ***p<0.01nr <- df2 %>%
filter(noaa_region == "northern_rockies")
nr_reg <- lm(Max.Ozone ~ moderate + severe,
data = nr)
stargazer(nr_reg, type = "text", title = "Northern Rockies", align = TRUE)##
## Northern Rockies
## ===============================================
## Dependent variable:
## ---------------------------
## Max.Ozone
## -----------------------------------------------
## moderate 1.264***
## (0.051)
##
## severe 3.088***
## (0.080)
##
## Constant 40.445***
## (0.027)
##
## -----------------------------------------------
## Observations 226,906
## R2 0.008
## Adjusted R2 0.008
## Residual Std. Error 10.477 (df = 226903)
## F Statistic 905.122*** (df = 2; 226903)
## ===============================================
## Note: *p<0.1; **p<0.05; ***p<0.01nw <- df2 %>%
filter(noaa_region == "northwest")
nw_reg <- lm(Max.Ozone ~ moderate + severe,
data = nw)
stargazer(nw_reg, type = "text", title = "Northwest", align = TRUE)##
## Northwest
## =================================================
## Dependent variable:
## -----------------------------
## Max.Ozone
## -------------------------------------------------
## moderate 3.534***
## (0.089)
##
## severe 6.177***
## (0.201)
##
## Constant 35.588***
## (0.045)
##
## -------------------------------------------------
## Observations 105,338
## R2 0.021
## Adjusted R2 0.021
## Residual Std. Error 12.445 (df = 105335)
## F Statistic 1,132.407*** (df = 2; 105335)
## =================================================
## Note: *p<0.1; **p<0.05; ***p<0.01ohv <- df2 %>%
filter(noaa_region == "ohio_valley")
ohv_reg <- lm(Max.Ozone ~ moderate + severe,
data = ohv)
stargazer(ohv_reg, type = "text", title = "Ohio Valley", align = TRUE)##
## Ohio Valley
## ==================================================
## Dependent variable:
## ------------------------------
## Max.Ozone
## --------------------------------------------------
## moderate 1.974***
## (0.046)
##
## severe 3.933***
## (0.108)
##
## Constant 43.448***
## (0.016)
##
## --------------------------------------------------
## Observations 1,006,480
## R2 0.003
## Adjusted R2 0.003
## Residual Std. Error 14.543 (df = 1006477)
## F Statistic 1,505.052*** (df = 2; 1006477)
## ==================================================
## Note: *p<0.1; **p<0.05; ***p<0.01sth <- df2 %>%
filter(noaa_region == "south")
sth_reg <- lm(Max.Ozone ~ moderate + severe,
data = sth)
stargazer(sth_reg, type = "text", title = "South", align = TRUE)##
## South
## =================================================
## Dependent variable:
## -----------------------------
## Max.Ozone
## -------------------------------------------------
## moderate 1.395***
## (0.036)
##
## severe 2.846***
## (0.053)
##
## Constant 39.519***
## (0.019)
##
## -------------------------------------------------
## Observations 901,864
## R2 0.004
## Adjusted R2 0.004
## Residual Std. Error 14.439 (df = 901861)
## F Statistic 1,863.638*** (df = 2; 901861)
## =================================================
## Note: *p<0.1; **p<0.05; ***p<0.01se <- df2 %>%
filter(noaa_region == "southeast")
se_reg <- lm(Max.Ozone ~ moderate + severe,
data = se)
stargazer(se_reg, type = "text", title = "Southeast", align = TRUE)##
## Southeast
## =================================================
## Dependent variable:
## -----------------------------
## Max.Ozone
## -------------------------------------------------
## moderate 4.057***
## (0.037)
##
## severe 5.616***
## (0.056)
##
## Constant 40.997***
## (0.017)
##
## -------------------------------------------------
## Observations 871,694
## R2 0.022
## Adjusted R2 0.022
## Residual Std. Error 13.697 (df = 871691)
## F Statistic 9,625.759*** (df = 2; 871691)
## =================================================
## Note: *p<0.1; **p<0.05; ***p<0.01sw <- df2 %>%
filter(noaa_region == "southwest")
sw_reg <- lm(Max.Ozone ~ moderate + severe,
data = sw)
stargazer(sw_reg, type = "text", title = "Southwest", align = TRUE)##
## Southwest
## =================================================
## Dependent variable:
## -----------------------------
## Max.Ozone
## -------------------------------------------------
## moderate 0.412***
## (0.054)
##
## severe 3.734***
## (0.071)
##
## Constant 46.152***
## (0.038)
##
## -------------------------------------------------
## Observations 241,725
## R2 0.012
## Adjusted R2 0.012
## Residual Std. Error 12.165 (df = 241722)
## F Statistic 1,453.451*** (df = 2; 241722)
## =================================================
## Note: *p<0.1; **p<0.05; ***p<0.01uw <- df2 %>%
filter(noaa_region == "upper_midwest")
uw_reg <- lm(Max.Ozone ~ moderate + severe,
data = uw)
stargazer(uw_reg, type = "text", title = "Upper Midwest", align = TRUE)##
## Upper Midwest
## ===============================================
## Dependent variable:
## ---------------------------
## Max.Ozone
## -----------------------------------------------
## moderate -0.745***
## (0.063)
##
## severe -1.035***
## (0.212)
##
## Constant 41.718***
## (0.023)
##
## -----------------------------------------------
## Observations 382,355
## R2 0.0004
## Adjusted R2 0.0004
## Residual Std. Error 13.058 (df = 382352)
## F Statistic 79.860*** (df = 2; 382352)
## ===============================================
## Note: *p<0.1; **p<0.05; ***p<0.01wst <- df2 %>%
filter(noaa_region == "west")
wst_reg <- lm(Max.Ozone ~ moderate + severe,
data = wst)
stargazer(wst_reg, type = "text", title = "West", align = TRUE)##
## West
## ===============================================
## Dependent variable:
## ---------------------------
## Max.Ozone
## -----------------------------------------------
## moderate -0.341***
## (0.084)
##
## severe 0.348***
## (0.091)
##
## Constant 45.955***
## (0.060)
##
## -----------------------------------------------
## Observations 147,761
## R2 0.0004
## Adjusted R2 0.0004
## Residual Std. Error 13.793 (df = 147758)
## F Statistic 29.007*** (df = 2; 147758)
## ===============================================
## Note: *p<0.1; **p<0.05; ***p<0.01