Skip to content

Instantly share code, notes, and snippets.

@MJacobs1985
Created May 16, 2022 11:23
Show Gist options
  • Save MJacobs1985/ecb6a924cba48a11295b5fccf1063170 to your computer and use it in GitHub Desktop.
Save MJacobs1985/ecb6a924cba48a11295b5fccf1063170 to your computer and use it in GitHub Desktop.
library(readxl)
library(dplyr)
library(ggplot2)
options(max.print=25000)
# https://www.rivm.nl/documenten/effectiviteit-van-covid-19-vaccinatie-tegen-ziekenhuis-en-intensive-care-opname-in-12#
df <- read_excel("VAX/Book1.xlsx")
View(df)
df%>%print(n=60)
# A tibble: 54 x 4
#Maand Leeftijd Status Opname
#<dttm> <chr> <chr> <dbl>
#1 2022-02-01 00:00:00 "\"12-49\"" "\"Niet gevaccineerd\"" 394
#2 2022-02-01 00:00:00 "\"12-49\"" "\"Basisserie deels afgerond\"" 29
#3 2022-02-01 00:00:00 "\"12-49\"" "\"Basisserie afgerond\"" 274
#4 2022-02-01 00:00:00 "\"12-49\"" "\"Boostervaccinatie ontvangen\"" 130
#5 2022-02-01 00:00:00 "\"12-49\"" "\"Herhaalprik ontvangen\"" 0
#6 2022-02-01 00:00:00 "\"12-49\"" "\"Onbekend\"" 21
#7 2022-02-01 00:00:00 "\"50-69\"" "\"Niet gevaccineerd\"" 498
#8 2022-02-01 00:00:00 "\"50-69\"" "\"Basisserie deels afgerond\"" 33
#9 2022-02-01 00:00:00 "\"50-69\"" "\"Basisserie afgerond\"" 389
#10 2022-02-01 00:00:00 "\"50-69\"" "\"Boostervaccinatie ontvangen\"" 236
#11 2022-02-01 00:00:00 "\"50-69\"" "\"Herhaalprik ontvangen\"" 0
#12 2022-02-01 00:00:00 "\"50-69\"" "\"Onbekend\"" 23
#13 2022-02-01 00:00:00 "\"70+\"" "\"Niet gevaccineerd\"" 702
#14 2022-02-01 00:00:00 "\"70+\"" "\"Basisserie deels afgerond\"" 31
#15 2022-02-01 00:00:00 "\"70+\"" "\"Basisserie afgerond\"" 520
#16 2022-02-01 00:00:00 "\"70+\"" "\"Boostervaccinatie ontvangen\"" 669
#17 2022-02-01 00:00:00 "\"70+\"" "\"Herhaalprik ontvangen\"" 0
#18 2022-02-01 00:00:00 "\"70+\"" "\"Onbekend\"" 14
#19 2022-03-01 00:00:00 "\"12-49\"" "\"Niet gevaccineerd\"" 284
#20 2022-03-01 00:00:00 "\"12-49\"" "\"Basisserie deels afgerond\"" 24
#21 2022-03-01 00:00:00 "\"12-49\"" "\"Basisserie afgerond\"" 208
#22 2022-03-01 00:00:00 "\"12-49\"" "\"Boostervaccinatie ontvangen\"" 162
#23 2022-03-01 00:00:00 "\"12-49\"" "\"Herhaalprik ontvangen\"" 0
#24 2022-03-01 00:00:00 "\"12-49\"" "\"Onbekend\"" 22
#25 2022-03-01 00:00:00 "\"50-69\"" "\"Niet gevaccineerd\"" 389
#26 2022-03-01 00:00:00 "\"50-69\"" "\"Basisserie deels afgerond\"" 26
#27 2022-03-01 00:00:00 "\"50-69\"" "\"Basisserie afgerond\"" 490
#28 2022-03-01 00:00:00 "\"50-69\"" "\"Boostervaccinatie ontvangen\"" 592
#29 2022-03-01 00:00:00 "\"50-69\"" "\"Herhaalprik ontvangen\"" 0
#30 2022-03-01 00:00:00 "\"50-69\"" "\"Onbekend\"" 20
#31 2022-03-01 00:00:00 "\"70+\"" "\"Niet gevaccineerd\"" 666
#32 2022-03-01 00:00:00 "\"70+\"" "\"Basisserie deels afgerond\"" 27
#33 2022-03-01 00:00:00 "\"70+\"" "\"Basisserie afgerond\"" 663
#34 2022-03-01 00:00:00 "\"70+\"" "\"Boostervaccinatie ontvangen\"" 1723
#35 2022-03-01 00:00:00 "\"70+\"" "\"Herhaalprik ontvangen\"" 35
#36 2022-03-01 00:00:00 "\"70+\"" "\"Onbekend\"" 24
#37 2022-04-01 00:00:00 "\"12-49\"" "\"Niet gevaccineerd\"" 63
#38 2022-04-01 00:00:00 "\"12-49\"" "\"Basisserie deels afgerond\"" 8
#39 2022-04-01 00:00:00 "\"12-49\"" "\"Basisserie afgerond\"" 64
#40 2022-04-01 00:00:00 "\"12-49\"" "\"Boostervaccinatie ontvangen\"" 61
#41 2022-04-01 00:00:00 "\"12-49\"" "\"Herhaalprik ontvangen\"" 0
#42 2022-04-01 00:00:00 "\"12-49\"" "\"Onbekend\"" 7
#43 2022-04-01 00:00:00 "\"50-69\"" "\"Niet gevaccineerd\"" 120
#44 2022-04-01 00:00:00 "\"50-69\"" "\"Basisserie deels afgerond\"" 5
#45 2022-04-01 00:00:00 "\"50-69\"" "\"Basisserie afgerond\"" 218
#46 2022-04-01 00:00:00 "\"50-69\"" "\"Boostervaccinatie ontvangen\"" 272
#47 2022-04-01 00:00:00 "\"50-69\"" "\"Herhaalprik ontvangen\"" 6
#48 2022-04-01 00:00:00 "\"50-69\"" "\"Onbekend\"" 20
#49 2022-04-01 00:00:00 "\"70+\"" "\"Niet gevaccineerd\"" 253
#50 2022-04-01 00:00:00 "\"70+\"" "\"Basisserie deels afgerond\"" 7
#51 2022-04-01 00:00:00 "\"70+\"" "\"Basisserie afgerond\"" 308
#52 2022-04-01 00:00:00 "\"70+\"" "\"Boostervaccinatie ontvangen\"" 692
#53 2022-04-01 00:00:00 "\"70+\"" "\"Herhaalprik ontvangen\"" 186
#54 2022-04-01 00:00:00 "\"70+\"" "\"Onbekend\"" 36
ggplot(df,aes(x=as.factor(Maand),
y=Opname,
col=as.factor(Status),
group=as.factor(Status)))+
geom_point(size=5)+
geom_line()+
facet_wrap(~Leeftijd, scales="free")+
theme_bw()+
labs(col="Status", x="Maand", y="IC-opname")
ggplot(df,aes(x=Maand,
y=Opname,
fill=as.factor(Status)))+
geom_bar(stat="identity",
position=position_dodge())+
facet_wrap(~Leeftijd, scales="free")+
theme_bw()+
labs(col="Status", x="Maand", y="IC-opname")
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment