This is an R Markdown
Notebook. When you execute code within the notebook, the results appear
beneath the code.
Try executing this chunk by clicking the Run button within
the chunk or by placing your cursor inside it and pressing
Cmd+Shift+Enter.
library(dplyr)
Attaching package: ‘dplyr’
The following objects are masked from ‘package:stats’:
filter, lag
The following objects are masked from ‘package:base’:
intersect, setdiff, setequal, union
data <- read.csv("death-rates-from-air-pollution-2.csv")
data
data2 <- read.csv("owid-co2-data.csv")
data2
library(dplyr)
library(ggplot2)
data_table <- data %>%
filter(Country %in% c("United States"))
data_table
data_table2 <- data2 %>%
filter(country %in% c("United States"))
data_table2
NA
library(ggthemes)
ggplot(data_table, aes(x = Year, y = Air.pollution..total...deaths.per.100.000., group = Country, color = Country)) +
labs(y = "Air Pollution (Total Deaths Per 100,000)", x = "Year",
title = "Change in Death Rate of Air Pollution Over Time in the United States",
color = "Country")+
geom_line(size = 1.5) +
theme_minimal()
Warning: Using `size` aesthetic for lines was deprecated in ggplot2 3.4.0.
Please use `linewidth` instead.

ggplot(data_table2) +
geom_line(aes(y = co2_per_capita, x = year, color = "Carbon Dioxide Emissions (Overall)")) +
geom_line(aes(y = gas_co2_per_capita, x = year, color = "Carbon Dioxide Emissions (Gas)")) +
labs(y = "Carbon Dioxide Emissions Per Capita", x = "Year",
title = "Change in Carbon Dioxide Emissions Per Capita Over Time in the United States",
color = "Emissions") +
scale_color_manual(labels = c("Carbon Dioxide Emissions (Overall)", "Carbon Dioxide Emissions (Gas)"), values = c("red", "blue")) +
theme_classic()

data_table3 <- data_table2
data_table3[data_table3$year >= 1990 & data_table2$year <= 2020,]
ggplot(data_table3[data_table3$year >= 1990 & data_table2$year <= 2020,], aes(x = year, y = nitrous_oxide_per_capita, group = country, color = country)) +
labs(y = "Nitrous Oxide Emissions Per Capita", x = "Year",
title = "Change in Nitrous Oxide Emissions Per Capita Over Time in the United States",
color = "Country")+
geom_line(size = 1.5) +
theme_light()

data_table4 <- data_table3
data_table4[data_table4$year %in% c('1800', '1850', '1900', '1950', '2000'), ]
year <- c("1850", "1850", "1900", "1900", "1950", "1950", "2000", "2000")
co2 <- c("Overall","Gas",
"Overall","Gas", "Overall","Gas", "Overall","Gas")
amount <- c(19.793, 0.000, 662.738, 11.542, 2541.485, 319.219, 6010.136, 1251.263)
circle <- data.frame(amount,year,co2)
ggplot(circle, aes(year, amount, fill = co2)) +
geom_bar(stat="identity", position = "dodge") +
labs(y = "Amount in Million Metric Tons", x = "Year",
title = "Amount of Carbon Dioxide Released in the United States From 1850 to 2000",
color = "Carbon Dioxide", fill = "Carbon Dioxide")

Add a new chunk by clicking the Insert Chunk button on the
toolbar or by pressing Cmd+Option+I.
When you save the notebook, an HTML file containing the code and
output will be saved alongside it (click the Preview button or
press Cmd+Shift+K to preview the HTML file).
The preview shows you a rendered HTML copy of the contents of the
editor. Consequently, unlike Knit, Preview does not
run any R code chunks. Instead, the output of the chunk when it was last
run in the editor is displayed.
LS0tCnRpdGxlOiAiUiBOb3RlYm9vayIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKVGhpcyBpcyBhbiBbUiBNYXJrZG93bl0oaHR0cDovL3JtYXJrZG93bi5yc3R1ZGlvLmNvbSkgTm90ZWJvb2suIFdoZW4geW91IGV4ZWN1dGUgY29kZSB3aXRoaW4gdGhlIG5vdGVib29rLCB0aGUgcmVzdWx0cyBhcHBlYXIgYmVuZWF0aCB0aGUgY29kZS4gCgpUcnkgZXhlY3V0aW5nIHRoaXMgY2h1bmsgYnkgY2xpY2tpbmcgdGhlICpSdW4qIGJ1dHRvbiB3aXRoaW4gdGhlIGNodW5rIG9yIGJ5IHBsYWNpbmcgeW91ciBjdXJzb3IgaW5zaWRlIGl0IGFuZCBwcmVzc2luZyAqQ21kK1NoaWZ0K0VudGVyKi4gCgpgYGB7cn0KbGlicmFyeShkcGx5cikKZGF0YSA8LSByZWFkLmNzdigiZGVhdGgtcmF0ZXMtZnJvbS1haXItcG9sbHV0aW9uLTIuY3N2IikKZGF0YQoKZGF0YTIgPC0gcmVhZC5jc3YoIm93aWQtY28yLWRhdGEuY3N2IikKZGF0YTIKYGBgCmBgYHtyfQpsaWJyYXJ5KGRwbHlyKQpsaWJyYXJ5KGdncGxvdDIpCgpkYXRhX3RhYmxlIDwtIGRhdGEgJT4lIAogIGZpbHRlcihDb3VudHJ5ICVpbiUgYygiVW5pdGVkIFN0YXRlcyIpKQpkYXRhX3RhYmxlCgpkYXRhX3RhYmxlMiA8LSBkYXRhMiAlPiUgCiAgZmlsdGVyKGNvdW50cnkgJWluJSBjKCJVbml0ZWQgU3RhdGVzIikpCmRhdGFfdGFibGUyCgpgYGAKYGBge3J9CmxpYnJhcnkoZ2d0aGVtZXMpCgpnZ3Bsb3QoZGF0YV90YWJsZSwgYWVzKHggPSBZZWFyLCB5ID0gQWlyLnBvbGx1dGlvbi4udG90YWwuLi5kZWF0aHMucGVyLjEwMC4wMDAuLCBncm91cCA9IENvdW50cnksIGNvbG9yID0gQ291bnRyeSkpICsKICBsYWJzKHkgPSAiQWlyIFBvbGx1dGlvbiAoVG90YWwgRGVhdGhzIFBlciAxMDAsMDAwKSIsIHggPSAiWWVhciIsIAogICAgICAgdGl0bGUgPSAiQ2hhbmdlIGluIERlYXRoIFJhdGUgb2YgQWlyIFBvbGx1dGlvbiBPdmVyIFRpbWUgaW4gdGhlIFVuaXRlZCBTdGF0ZXMiLAogICAgICAgY29sb3IgPSAiQ291bnRyeSIpKwogIGdlb21fbGluZShzaXplID0gMS41KSArCiAgCiAgdGhlbWVfbWluaW1hbCgpCmBgYApgYGB7cn0KZ2dwbG90KGRhdGFfdGFibGUyKSArCiAgZ2VvbV9saW5lKGFlcyh5ID0gY28yX3Blcl9jYXBpdGEsIHggPSB5ZWFyLCBjb2xvciA9ICJDYXJib24gRGlveGlkZSBFbWlzc2lvbnMgKE92ZXJhbGwpIikpICsKICBnZW9tX2xpbmUoYWVzKHkgPSBnYXNfY28yX3Blcl9jYXBpdGEsIHggPSB5ZWFyLCBjb2xvciA9ICJDYXJib24gRGlveGlkZSBFbWlzc2lvbnMgKEdhcykiKSkgKwogIGxhYnMoeSA9ICJDYXJib24gRGlveGlkZSBFbWlzc2lvbnMgUGVyIENhcGl0YSIsIHggPSAiWWVhciIsIAogICAgICAgdGl0bGUgPSAiQ2hhbmdlIGluIENhcmJvbiBEaW94aWRlIEVtaXNzaW9ucyBQZXIgQ2FwaXRhIE92ZXIgVGltZSBpbiB0aGUgVW5pdGVkIFN0YXRlcyIsCiAgICAgICBjb2xvciA9ICJFbWlzc2lvbnMiKSArCiAgc2NhbGVfY29sb3JfbWFudWFsKGxhYmVscyA9IGMoIkNhcmJvbiBEaW94aWRlIEVtaXNzaW9ucyAoT3ZlcmFsbCkiLCAiQ2FyYm9uIERpb3hpZGUgRW1pc3Npb25zIChHYXMpIiksIHZhbHVlcyA9IGMoInJlZCIsICJibHVlIikpICsKICAKICB0aGVtZV9jbGFzc2ljKCkKYGBgCmBgYHtyfQpkYXRhX3RhYmxlMyA8LSBkYXRhX3RhYmxlMgpkYXRhX3RhYmxlM1tkYXRhX3RhYmxlMyR5ZWFyID49IDE5OTAgJiBkYXRhX3RhYmxlMiR5ZWFyIDw9IDIwMjAsXQoKZ2dwbG90KGRhdGFfdGFibGUzW2RhdGFfdGFibGUzJHllYXIgPj0gMTk5MCAmIGRhdGFfdGFibGUyJHllYXIgPD0gMjAyMCxdLCBhZXMoeCA9IHllYXIsIHkgPSBuaXRyb3VzX294aWRlX3Blcl9jYXBpdGEsIGdyb3VwID0gY291bnRyeSwgY29sb3IgPSBjb3VudHJ5KSkgKwogIGxhYnMoeSA9ICJOaXRyb3VzIE94aWRlIEVtaXNzaW9ucyBQZXIgQ2FwaXRhIiwgeCA9ICJZZWFyIiwgCiAgICAgICB0aXRsZSA9ICJDaGFuZ2UgaW4gTml0cm91cyBPeGlkZSBFbWlzc2lvbnMgUGVyIENhcGl0YSBPdmVyIFRpbWUgaW4gdGhlIFVuaXRlZCBTdGF0ZXMiLAogICAgICAgY29sb3IgPSAiQ291bnRyeSIpKwogIGdlb21fbGluZShzaXplID0gMS41KSArCiAgCiAgdGhlbWVfbGlnaHQoKQpgYGAKYGBge3J9CmRhdGFfdGFibGU0IDwtIGRhdGFfdGFibGUzCmRhdGFfdGFibGU0W2RhdGFfdGFibGU0JHllYXIgJWluJSBjKCcxODAwJywgJzE4NTAnLCAnMTkwMCcsICcxOTUwJywgJzIwMDAnKSwgXQoKeWVhciA8LSBjKCIxODUwIiwgIjE4NTAiLCAiMTkwMCIsICIxOTAwIiwgIjE5NTAiLCAiMTk1MCIsICIyMDAwIiwgIjIwMDAiKQpjbzIgPC0gYygiT3ZlcmFsbCIsIkdhcyIsCiAgICAgICAgICAgICJPdmVyYWxsIiwiR2FzIiwgIk92ZXJhbGwiLCJHYXMiLCAiT3ZlcmFsbCIsIkdhcyIpCmFtb3VudCA8LSBjKDE5Ljc5MywgMC4wMDAsIDY2Mi43MzgsIDExLjU0MiwgMjU0MS40ODUsIDMxOS4yMTksIDYwMTAuMTM2LCAxMjUxLjI2MykKCmNpcmNsZSA8LSBkYXRhLmZyYW1lKGFtb3VudCx5ZWFyLGNvMikKCmdncGxvdChjaXJjbGUsIGFlcyh5ZWFyLCBhbW91bnQsIGZpbGwgPSBjbzIpKSArCiAgZ2VvbV9iYXIoc3RhdD0iaWRlbnRpdHkiLCBwb3NpdGlvbiA9ICJkb2RnZSIpICsKICBsYWJzKHkgPSAiQW1vdW50IGluIE1pbGxpb24gTWV0cmljIFRvbnMiLCB4ID0gIlllYXIiLCAKICAgICAgIHRpdGxlID0gIkFtb3VudCBvZiBDYXJib24gRGlveGlkZSBSZWxlYXNlZCBpbiB0aGUgVW5pdGVkIFN0YXRlcyBGcm9tIDE4NTAgdG8gMjAwMCIsCiAgICAgICBjb2xvciA9ICJDYXJib24gRGlveGlkZSIsIGZpbGwgPSAiQ2FyYm9uIERpb3hpZGUiKQpgYGAKQWRkIGEgbmV3IGNodW5rIGJ5IGNsaWNraW5nIHRoZSAqSW5zZXJ0IENodW5rKiBidXR0b24gb24gdGhlIHRvb2xiYXIgb3IgYnkgcHJlc3NpbmcgKkNtZCtPcHRpb24rSSouCgpXaGVuIHlvdSBzYXZlIHRoZSBub3RlYm9vaywgYW4gSFRNTCBmaWxlIGNvbnRhaW5pbmcgdGhlIGNvZGUgYW5kIG91dHB1dCB3aWxsIGJlIHNhdmVkIGFsb25nc2lkZSBpdCAoY2xpY2sgdGhlICpQcmV2aWV3KiBidXR0b24gb3IgcHJlc3MgKkNtZCtTaGlmdCtLKiB0byBwcmV2aWV3IHRoZSBIVE1MIGZpbGUpLiAKClRoZSBwcmV2aWV3IHNob3dzIHlvdSBhIHJlbmRlcmVkIEhUTUwgY29weSBvZiB0aGUgY29udGVudHMgb2YgdGhlIGVkaXRvci4gQ29uc2VxdWVudGx5LCB1bmxpa2UgKktuaXQqLCAqUHJldmlldyogZG9lcyBub3QgcnVuIGFueSBSIGNvZGUgY2h1bmtzLiBJbnN0ZWFkLCB0aGUgb3V0cHV0IG9mIHRoZSBjaHVuayB3aGVuIGl0IHdhcyBsYXN0IHJ1biBpbiB0aGUgZWRpdG9yIGlzIGRpc3BsYXllZC4KCg==