library(tidyverse)
## Warning: package 'tidyverse' was built under R version 4.0.2
## -- Attaching packages ----------------------------------------------------------------------------------------------- tidyverse 1.3.0 --
## v ggplot2 3.3.2 v purrr 0.3.4
## v tibble 3.0.3 v dplyr 1.0.2
## v tidyr 1.1.2 v stringr 1.4.0
## v readr 1.3.1 v forcats 0.5.0
## Warning: package 'tibble' was built under R version 4.0.2
## Warning: package 'tidyr' was built under R version 4.0.2
## Warning: package 'dplyr' was built under R version 4.0.2
## -- Conflicts -------------------------------------------------------------------------------------------------- tidyverse_conflicts() --
## x dplyr::filter() masks stats::filter()
## x dplyr::lag() masks stats::lag()
library(ggplot2)
library(openintro)
## Warning: package 'openintro' was built under R version 4.0.2
## Loading required package: airports
## Warning: package 'airports' was built under R version 4.0.2
## Loading required package: cherryblossom
## Warning: package 'cherryblossom' was built under R version 4.0.2
## Loading required package: usdata
## Warning: package 'usdata' was built under R version 4.0.2
library(dplyr)
lax_flights <- read.csv("C:/Users/Jerome/Documents/Math_217/lax_flights.csv")
sfo_feb_flights <- read.csv ("C:/Users/Jerome/Documents/Math_217/sfo_feb_flights.csv")
head(sfo_feb_flights)
## X year month day dep_time dep_delay arr_time arr_delay carrier tailnum flight
## 1 1 2013 2 18 1527 57 1903 48 DL N711ZX 1322
## 2 2 2013 2 3 613 14 1008 38 UA N502UA 691
## 3 3 2013 2 15 955 -5 1313 -28 DL N717TW 1765
## 4 4 2013 2 18 1928 15 2239 -6 UA N24212 1214
## 5 5 2013 2 24 1340 2 1644 -21 UA N76269 1111
## 6 6 2013 2 25 1415 -10 1737 -13 UA N532UA 394
## origin dest air_time distance hour minute
## 1 JFK SFO 358 2586 15 27
## 2 JFK SFO 367 2586 6 13
## 3 JFK SFO 338 2586 9 55
## 4 EWR SFO 353 2565 19 28
## 5 EWR SFO 341 2565 13 40
## 6 JFK SFO 355 2586 14 15
ggplot (data = sfo_feb_flights, aes(x = dep_delay)) +
geom_histogram(bin = 5)
## Warning: Ignoring unknown parameters: bin
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.
summary(sfo_feb_flights$dep_delay)
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## -10.0 -5.0 -2.0 10.5 9.0 209.0
rowan <- read.csv ("C:/Users/Jerome/Documents/Math_217/rowan_data.csv")
plot (rowan$altitude, rowan$respiration, main = "Scatterplot of Rowan Altitude and Respiration", xlab = "Altitude", ylab = "Respiration")
fit <- lm(rowan$respiration ~ rowan$altitude)
abline(fit)
