library(tidyverse)
## ── Attaching packages ──────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────── tidyverse 1.3.0 ──
## ✓ ggplot2 3.3.0 ✓ purrr 0.3.3
## ✓ tibble 2.1.3 ✓ dplyr 0.8.5
## ✓ tidyr 1.0.2 ✓ stringr 1.4.0
## ✓ readr 1.3.1 ✓ forcats 0.5.0
## ── Conflicts ─────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────────── tidyverse_conflicts() ──
## x dplyr::filter() masks stats::filter()
## x dplyr::lag() masks stats::lag()
##Data pulled from google sheets
deaths_317<-7122
recovered_317<-79387
world_pop<-
death_rate_317<-deaths_317/(recovered_317+deaths_317); death_rate_317
## [1] 0.08232669
deaths_318<-8054
recovered_318<-81917
death_rate_318<-deaths_318/(recovered_318+deaths_318); death_rate_318
## [1] 0.08951773
deaths_320<-9780
recovered_320<-85590
death_rate_320<-deaths_320/(recovered_320+deaths_320); death_rate_320
## [1] 0.102548
deaths_326<-23719
recovered_326<-119843
death_rate_326<-deaths_326/(recovered_326+deaths_326); death_rate_326
## [1] 0.1652178
deaths_330<-34873
recovered_330<-157463
death_rate_330<-deaths_330/(recovered_330+deaths_330); death_rate_330
## [1] 0.1813129
deaths_331<-37947
recovered_331<-166918
death_rate_331<-deaths_331/(recovered_331+deaths_331); death_rate_331
## [1] 0.1852293
deaths_401<-42196
recovered_401<-175911
death_rate_401<-deaths_401/(recovered_401+deaths_401); death_rate_401
## [1] 0.1934647
deaths_402<-48111
recovered_402<-199106
death_rate_402<-deaths_402/(recovered_402+deaths_402); death_rate_402
## [1] 0.1946104
deaths_403<-55822
recovered_403<-213639
death_rate_403<-deaths_403/(recovered_403+deaths_403); death_rate_403
## [1] 0.2071617
deaths_404<-60359
recovered_404<-239385
death_rate_404<-deaths_404/(recovered_404+deaths_404); death_rate_404
## [1] 0.2013685
'Today World'
## [1] "Today World"
deaths_405<-68935
recovered_405<-260212
death_rate_405<-deaths_405/(recovered_405+deaths_405); death_rate_405
## [1] 0.2094353
confirmed_us_330<-143724
deaths_us_330<-2592
recovered_us_330<-4050
death_rate_us_330<-deaths_us_330/(recovered_us_330+deaths_us_330); death_rate_us_330
## [1] 0.3902439
confirmed_us_331<-164359
deaths_us_331<-3173
recovered_us_331<-4784
death_rate_us_331<-deaths_us_331/(recovered_us_331+deaths_us_331); death_rate_us_331
## [1] 0.3987684
confirmed_us_401<-188280
deaths_us_401<-3976
recovered_us_401<-5020
death_rate_us_401<-deaths_us_401/(recovered_us_401+deaths_us_401); death_rate_us_401
## [1] 0.4419742
confirmed_us_402<-215608
deaths_us_402<-5110
recovered_us_402<-7760
death_rate_us_402<-deaths_us_402/(recovered_us_402+deaths_us_402); death_rate_us_402
## [1] 0.3970474
confirmed_us_403<-257379
deaths_us_403<-6558
recovered_us_403<-9359
death_rate_us_403<-deaths_us_403/(recovered_us_403+deaths_us_403); death_rate_us_403
## [1] 0.4120123
confirmed_us_404<-277621
deaths_us_404<-7403
recovered_us_404<-10946
death_rate_us_404<-deaths_us_404/(recovered_us_404+deaths_us_404); death_rate_us_404
## [1] 0.4034552
'Today US'
## [1] "Today US"
confirmed_us_405<-330263
deaths_us_405<-9444
recovered_us_405<-12713
death_rate_us_405<-deaths_us_405/(recovered_us_405+deaths_us_405); death_rate_us_405
## [1] 0.426231
death_rate_world<-c(0.08232669, 0.08951773, 0.102548, 0.1652178, 0.1813129, 0.1852293, 0.1934647, 0.1946104, 0.2071617, 0.2013685,0.2094353)
deaths_world<-c(deaths_317, deaths_318, deaths_320, deaths_326, deaths_330, deaths_331, deaths_401, deaths_402, deaths_403, deaths_404, deaths_405)
days_world<-c(1, 2, 4, 10, 14, 15, 16, 17, 18,19,20)
recover_world<-c(recovered_317, recovered_318, recovered_320, recovered_326, recovered_330, recovered_331, recovered_401, recovered_402, recovered_403, recovered_404, recovered_405)
plot(days_world, recover_world)
plot(days_world, death_rate_world)
linearModdeathrateworld<-lm(death_rate_world ~ days_world); linearModdeathrateworld
##
## Call:
## lm(formula = death_rate_world ~ days_world)
##
## Coefficients:
## (Intercept) days_world
## 0.079488 0.006896
linearModrecover<-lm(log(recover_world) ~ days_world); linearModrecover
##
## Call:
## lm(formula = log(recover_world) ~ days_world)
##
## Coefficients:
## (Intercept) days_world
## 11.15009 0.06172
plot(days_world, deaths_world)
days_us<-c(1,2,3,4, 5, 6, 7)
death_rate_us<- c(0.39024390, .3987684, 0.4419742, 0.3970474, 0.4120123, 0.4034552, 0.426231)
plot(days_us, death_rate_us)
linearModdeathrate_us<-lm(death_rate_us~ days_us); linearModdeathrate_us
##
## Call:
## lm(formula = death_rate_us ~ days_us)
##
## Coefficients:
## (Intercept) days_us
## 0.39748 0.00312
'Geometric Mean- Seems pretty useless- Not sure how to use with exponential growth'
## [1] "Geometric Mean- Seems pretty useless- Not sure how to use with exponential growth"
n<-length(death_rate_world)
death_rate_world_mean_geo<-(prod(death_rate_world))^(1/n); death_rate_world_mean_geo
## [1] 0.1565011
deaths_hour_last6days<-(deaths_326-deaths_320)/144;deaths_hour_last6days
## [1] 96.79861
deaths_hour_317<-(deaths_318-deaths_317)/24;deaths_hour_317
## [1] 38.83333
deaths_hour_320<-(deaths_320-deaths_318)/48;deaths_hour_320
## [1] 35.95833