library(readxl)
d <- read_excel("percentageofcovid19casesthatarepositivefororf1abandngenes.xlsx",
sheet = "Sheet4")
d$N2<-d$Tests-d$N
str(d)
## tibble [11 × 5] (S3: tbl_df/tbl/data.frame)
## $ Date : POSIXct[1:11], format: "2020-09-30" "2020-10-07" ...
## $ Tests : num [1:11] 555 845 1104 1076 1105 ...
## $ Percent: num [1:11] 11 8 7 9 11 10 13 13 20 20 ...
## $ N : num [1:11] 61 67.6 77.3 96.8 121.5 ...
## $ N2 : num [1:11] 494 777 1027 979 983 ...
theme_set(theme_bw())
ggplot(d,aes(x=Date,y=N)) + geom_line() + geom_point() + geom_vline(xintercept = d$Date[8], col="red") + geom_line(aes(y=N2), col="blue") + geom_point(aes(y=N2), col="blue") + scale_y_continuous(trans='log2')
aqm::dt(d)
theme_set(theme_bw())
ggplot(d,aes(x=Date,y=Percent)) + geom_line() + geom_point() + geom_vline(xintercept = d$Date[8], col="red") + scale_y_continuous(trans='log2')
theme_set(theme_bw())
ggplot(d,aes(x=Date,y=Tests)) + geom_line() + geom_point() + geom_vline(xintercept = d$Date[8], col="red") + scale_y_continuous(trans='log10')
theme_set(theme_bw())
ggplot(d,aes(x=Date,y=N)) + geom_line() + geom_point() + geom_vline(xintercept = d$Date[8], col="red")