There are 3 sets of data included in this package from the US Centers
for Disease Control and Prevention (CDC) and the US Food and Drug
Administration (FDA) United States Vaccine Adverse
Event Reporting System (VAERS) program. There is the vax
table which includes information about the vaccines. There is the
data table which has data from the VAERS-1 form. There is
the symptoms table which has the adverse events.
Please see the README, DESCRIPTION, and the individual R data files for more information.
install.load::load_package("vaers", "vaersND", "data.table", "pivottabler", "skimr",
"DT", "formatR")
# load needed packages using the load_package function from the install.load
# package (it is assumed that you have already installed these packages)
## load vaers_data
data(vaers_data)
# How many reports for each sex?
setkey(vaers_data, VAERS_ID)
vaers_data[, .N, by = list(SEX)] # Source 1## SEX N
## <char> <int>
## 1: F 428001
## 2: M 239499
## 3: U 92268# How many people died, listed by ID (only using the vaers_data]?
vaers_data[DIED == "Y", list(VAERS_ID, DIED)]## Key: <VAERS_ID>
## VAERS_ID DIED
## <num> <char>
## 1: 25003 Y
## 2: 25026 Y
## 3: 25028 Y
## 4: 25520 Y
## 5: 25526 Y
## ---
## 7110: 1157506 Y
## 7111: 1157515 Y
## 7112: 1157534 Y
## 7113: 1157683 Y
## 7114: 1157733 Y## Key: <VAERS_ID>
## VAERS_ID SEX
## <num> <char>
## 1: 25001 F
## 2: 25006 F
## 3: 25008 F
## 4: 25016 F
## 5: 25017 F
## ---
## 427997: 1157040 F
## 427998: 1157133 F
## 427999: 1157506 F
## 428000: 1157534 F
## 428001: 1157733 F## load vaers_symptoms
data(vaers_symptoms)
# How many reports of autism for SYMPTOM1?
setkey(vaers_symptoms, VAERS_ID)
vaers_symptoms[SYMPTOM1 == "Autism", .N] # Source 2## [1] 974## load vaers_vax
data(vaers_vax)
# What are the counts for each of the VAX_TYPEs?
setkey(vaers_vax, VAERS_ID)
vaers_vax[, .N, by = list(VAX_TYPE)][order(VAX_TYPE)] # Source 1## VAX_TYPE N
## <char> <int>
## 1: 6VAX-F 190
## 2: ADEN 8
## 3: ADEN_4_7 149
## 4: ANTH 9137
## 5: BCG 193
## 6: CEE 1
## 7: CHOL 172
## 8: COVID19 57048
## 9: DF 1
## 10: DPP 200
## 11: DT 2379
## 12: DTAP 56061
## 13: DTAPH 701
## 14: DTAPHEPBIP 11036
## 15: DTAPIPV 8458
## 16: DTAPIPVHIB 8247
## 17: DTIPV 29
## 18: DTOX 76
## 19: DTP 20758
## 20: DTPHEP 40
## 21: DTPHIB 6070
## 22: DTPIHI 57
## 23: DTPIPV 88
## 24: DTPPHIB 54
## 25: DTaPIPV 16
## 26: DTaPIPVHib 22
## 27: FLU(H1N1) 8471
## 28: FLU3 107570
## 29: FLU4 29476
## 30: FLUA3 2344
## 31: FLUA4 691
## 32: FLUC3 999
## 33: FLUC4 4388
## 34: FLUN(H1N1) 2892
## 35: FLUN3 5819
## 36: FLUN4 2993
## 37: FLUR3 198
## 38: FLUR4 1743
## 39: FLUX 5449
## 40: FLUX(H1N1) 2026
## 41: H5N1 5
## 42: HBHEPB 5410
## 43: HBPV 259
## 44: HEP 58146
## 45: HEPA 36045
## 46: HEPAB 2992
## 47: HEPATYP 5
## 48: HIBV 48149
## 49: HPV2 260
## 50: HPV4 36902
## 51: HPV9 13529
## 52: HPVX 1056
## 53: IPV 36349
## 54: JEV 484
## 55: JEV1 366
## 56: JEVX 37
## 57: LYME 2224
## 58: MEA 638
## 59: MEN 2692
## 60: MENB 5839
## 61: MENHIB 19
## 62: MER 161
## 63: MM 66
## 64: MMR 75552
## 65: MMRV 15623
## 66: MNC 1
## 67: MNQ 23079
## 68: MNQHIB 21
## 69: MU 155
## 70: MUR 29
## 71: OPV 23671
## 72: PER 143
## 73: PLAGUE 28
## 74: PNC 23316
## 75: PNC10 16
## 76: PNC13 26595
## 77: PPV 58896
## 78: RAB 4012
## 79: RUB 886
## 80: RV 693
## 81: RV1 2380
## 82: RV5 18489
## 83: RVX 880
## 84: SMALL 5736
## 85: SSEV 4
## 86: TBE 5
## 87: TD 17483
## 88: TDAP 37270
## 89: TDAPIPV 7
## 90: TTOX 3204
## 91: TYP 7764
## 92: UNK 6872
## 93: VARCEL 78021
## 94: VARZOS 82794
## 95: YF 3388
## VAX_TYPE N
## <char> <int>## [1] 75552# How many reports of COVID19 as the VAX_TYPE or the VAX_NAME contains COVID19?
vaers_vax[VAX_TYPE == "COVID19" | VAX_NAME %like% "COVID19", .N] # Source 2## [1] 57048# join vaers_data and vaers_vax to get access to the type and name of the
# vaccination$
data_vax <- vaers_data[vaers_vax]
setkey(data_vax, VAERS_ID)
# get the total number of VAERS_IDs in data_vax
length(data_vax$VAERS_ID)## [1] 1124902## [1] 759765# remove the duplicated rows
data_vax2 <- unique(data_vax)
setkey(data_vax2, VAERS_ID)
# How many people have died from all vaccines since 1990? Include the number
# of entries for each VAERS_ID
data_vax2[DIED == "Y", .N, list(VAERS_ID, DIED)] # Source 2## VAERS_ID DIED N
## <num> <char> <int>
## 1: 25003 Y 2
## 2: 25026 Y 3
## 3: 25028 Y 1
## 4: 25520 Y 1
## 5: 25526 Y 1
## ---
## 7110: 1157506 Y 1
## 7111: 1157515 Y 1
## 7112: 1157534 Y 1
## 7113: 1157683 Y 1
## 7114: 1157733 Y 1# How many people have died from all vaccines since 1990?
data_vax2[DIED == "Y", list(VAERS_ID, DIED)] # Source 2## Key: <VAERS_ID>
## VAERS_ID DIED
## <num> <char>
## 1: 25003 Y
## 2: 25003 Y
## 3: 25026 Y
## 4: 25026 Y
## 5: 25026 Y
## ---
## 12665: 1157506 Y
## 12666: 1157515 Y
## 12667: 1157534 Y
## 12668: 1157683 Y
## 12669: 1157733 Y# How many people have died from the COVID19 Experimental Gene Therapy
# Technology {Search for discussion by Judy Mikovits, Ph.D.]?
data_vax2[VAX_TYPE == "COVID19" & DIED == "Y", list(VAERS_ID, DIED)] # Source 2## Key: <VAERS_ID>
## VAERS_ID DIED
## <num> <char>
## 1: 909095 Y
## 2: 910363 Y
## 3: 913143 Y
## 4: 913733 Y
## 5: 914604 Y
## ---
## 2279: 1157506 Y
## 2280: 1157515 Y
## 2281: 1157534 Y
## 2282: 1157683 Y
## 2283: 1157733 Y# What is the ratio, expressed as a percentage, of the COVID19 Experimental
# Gene Therapy Technology {Search for discussion by Judy Mikovits, Ph.D.] to
# all deaths associated with vaccination$ since 1990?
(nrow(data_vax2[VAX_TYPE == "COVID19" & DIED == "Y"])/nrow(data_vax2[DIED == "Y"])) *
100 # Source 2## [1] 18.02036# What is the first date reported with COVID19 as the VAX_TYPE or the VAX_NAME
# contains COVID19?
data_vax2ed <- data_vax2[VAX_TYPE == "COVID19" | VAX_NAME %like% "COVID19"]
data_vax2ed[1, "VAX_DATE"]## VAX_DATE
## <Date>
## 1: 2020-12-15# How many days from the vaccination date to the person's death date? copy the
# data_vax2 data.table object
data_vax2er <- copy(data_vax2)
# create a new column with the value of the time difference between the date
# the person died and the vaccination date
data_vax2er[, `Days_Vaccination$_Death` := data_vax2$DATEDIED - data_vax2$VAX_DATE]
# group the number of days and sort by the greatest number of days between the
# date the person died and the vaccination date
data_vax2err <- data_vax2er[, .(.N), by = .(`Days_Vaccination$_Death`)][order(-N)] # Source 2
data_vax2err## Days_Vaccination$_Death N
## <difftime> <int>
## 1: NA days 1113385
## 2: 1 days 2264
## 3: 2 days 1226
## 4: 3 days 909
## 5: 0 days 673
## ---
## 367: -15 days 1
## 368: 367 days 1
## 369: -23 days 1
## 370: 352 days 1
## 371: -9 days 1
# Create a pivot table of the data_vax2 data based on VAX_NAME and if they DIED
qpvt(data_vax2, c("VAX_NAME"), "DIED", "n()")## Y NA Total
## ADENOVIRUS (NO BRAND NAME) 3 3
## ADENOVIRUS (TYPE 4, NO BRAND NAME) 1 1
## ADENOVIRUS (TYPE 7, NO BRAND NAME) 4 4
## ADENOVIRUS TYPES 4 & 7, LIVE, ORAL (NO BRAND NAME) 1 148 149
## ANTHRAX (BIOTHRAX) 27 6182 6209
## ANTHRAX (NO BRAND NAME) 4 2921 2925
## BCG (MYCOBAX) 1 36 37
## BCG (NO BRAND NAME) 2 87 89
## BCG (TICE) 67 67
## CHOLERA (NO BRAND NAME) 35 35
## CHOLERA (USP) 77 77
## CHOLERA (VAXCHORA) 60 60
## COVID19 (COVID19 (JANSSEN)) 54 5277 5331
## COVID19 (COVID19 (MODERNA)) 1167 22321 23488
## COVID19 (COVID19 (PFIZER-BIONTECH)) 1050 27036 28086
## COVID19 (COVID19 (UNKNOWN)) 12 71 83
## DENGUE TETRAVALENT (DENGVAXIA) 1 1
## DIPHTHERIA TOXOID + PERTUSSIS + IPV (NO BRAND NAME) 1 176 177
## DIPHTHERIA TOXOID + PERTUSSIS + IPV (QUATRO VIRELON) 23 23
## DIPHTHERIA TOXOIDS (NO BRAND NAME) 76 76
## DT + IPV (NO BRAND NAME) 29 29
## DT ADSORBED (DECAVAC) 237 237
## DT ADSORBED (DITANRIX) 1 78 79
## DT ADSORBED (NO BRAND NAME) 11 2052 2063
## DT+IPV+HIB+HEPB (NO BRAND NAME) 57 57
## DTAP (ACEL-IMUNE) 74 3270 3344
## DTAP (CERTIVA) 4 98 102
## DTAP (DAPTACEL) 127 14795 14922
## DTAP (INFANRIX) 233 18561 18794
## DTAP (NO BRAND NAME) 48 4592 4640
## DTAP (TRIPEDIA) 204 14053 14257
## DTAP + HEPB + IPV (FOREIGN) 14 14
## DTAP + HEPB + IPV (INFANRIX PENTA) 4 4
## DTAP + HEPB + IPV (NO BRAND NAME) 1 164 165
## DTAP + HEPB + IPV (PEDIARIX) 402 10448 10850
## DTAP + HIB (NO BRAND NAME) 1 152 153
## DTAP + HIB (TRIHIBIT) 2 546 548
## DTAP + IPV (FOREIGN) 12 12
## DTAP + IPV (INFANRIX TETRA) 5 5
## DTAP + IPV (KINRIX) 8 7403 7411
## DTAP + IPV (QUADRACEL) 899 899
## DTAP + IPV (UNKNOWN) 146 146
## DTAP + IPV + HIB (FOREIGN) 29 29
## DTAP + IPV + HIB (INFANRIX QUINTA) 3 3
## DTAP + IPV + HIB (NO BRAND NAME) 2 77 79
## DTAP + IPV + HIB (PENTACEL) 211 7855 8066
## DTAP + IPV + HIB (UNKNOWN) 87 87
## DTAP+IPV+HEPB+HIB (FOREIGN) 23 23
## DTAP+IPV+HEPB+HIB (HEXAVAC) 2 2
## DTAP+IPV+HEPB+HIB (HEXAVAX) 26 26
## DTAP+IPV+HEPB+HIB (INFANRIX HEXA) 15 15
## DTAP+IPV+HEPB+HIB (NO BRAND NAME) 3 121 124
## DTP (NO BRAND NAME) 359 12273 12632
## DTP (TRI-IMMUNOL) 312 7814 8126
## DTP + HEP B (NO BRAND NAME) 1 23 24
## DTP + HEP B (TRITANRIX) 16 16
## DTP + HIB (DTP + ACTHIB) 1 46 47
## DTP + HIB (NO BRAND NAME) 2 102 104
## DTP + HIB (TETRAMUNE) 218 5701 5919
## DTP + IPV (NO BRAND NAME) 87 87
## DTP + IPV + ACT-HIB (NO BRAND NAME) 54 54
## FSME-IMMUN. (NO BRAND NAME) 1 1
## HEP A (HAVRIX) 56 19187 19243
## HEP A (NO BRAND NAME) 9 2383 2392
## HEP A (VAQTA) 31 14378 14409
## HEP A + HEP B (NO BRAND NAME) 1 115 116
## HEP A + HEP B (TWINRIX) 8 2868 2876
## HEP A + TYP (FOREIGN) 3 3
## HEP A + TYP (HEPATYRIX) 2 2
## HEP B (ENGERIX-B) 369 26306 26675
## HEP B (FOREIGN) 34 1634 1668
## HEP B (GENHEVAC B) 53 53
## HEP B (HEPLISAV-B) 1 253 254
## HEP B (NO BRAND NAME) 72 3668 3740
## HEP B (RECOMBIVAX HB) 419 25324 25743
## HIB (ACTHIB) 415 16126 16541
## HIB (HIBERIX) 17 569 586
## HIB (HIBTITER) 562 13783 14345
## HIB (NO BRAND NAME) 53 2896 2949
## HIB (OMNIHIB) 59 59
## HIB (PEDVAXHIB) 278 11974 12252
## HIB (PROHIBIT) 35 1361 1396
## HIB (TETRACOQ) 8 8
## HIB + HEP B (COMVAX) 193 5204 5397
## HIB + HEP B (NO BRAND NAME) 13 13
## HIB POLYSACCHARIDE (FOREIGN) 11 11
## HIB POLYSACCHARIDE (HIBIMUNE) 9 228 237
## HIB POLYSACCHARIDE (NO BRAND NAME) 22 22
## HPV (CERVARIX) 2 258 260
## HPV (GARDASIL 9) 24 13505 13529
## HPV (GARDASIL) 191 36711 36902
## HPV (NO BRAND NAME) 26 1030 1056
## INFLUENZA (H1N1) (H1N1 (MONOVALENT) (CSL) 2 271 273
## INFLUENZA (H1N1) (H1N1 (MONOVALENT) (GSK) 53 53
## INFLUENZA (H1N1) (H1N1 (MONOVALENT) (MEDIMMUNE) 8 2882 2890
## INFLUENZA (H1N1) (H1N1 (MONOVALENT) (MEDIMMUNE)) 1 1 2
## INFLUENZA (H1N1) (H1N1 (MONOVALENT) (NOVARTIS) 19 3387 3406
## INFLUENZA (H1N1) (H1N1 (MONOVALENT) (NOVARTIS)) 6 6
## INFLUENZA (H1N1) (H1N1 (MONOVALENT) (SANOFI) 37 4688 4725
## INFLUENZA (H1N1) (H1N1 (MONOVALENT) (SANOFI)) 8 8
## INFLUENZA (H1N1) (H1N1 (MONOVALENT) (UNKNOWN) 22 1972 1994
## INFLUENZA (H1N1) (H1N1 (MONOVALENT) (UNKNOWN)) 32 32
## INFLUENZA (SEASONAL) (AFLURIA QUADRIVALENT) 8 3281 3289
## INFLUENZA (SEASONAL) (AFLURIA) 20 7168 7188
## INFLUENZA (SEASONAL) (AGRIFLU) 1 48 49
## INFLUENZA (SEASONAL) (FLU-IMUNE) 10 278 288
## INFLUENZA (SEASONAL) (FLUAD QUADRIVALENT) 1 689 690
## INFLUENZA (SEASONAL) (FLUAD) 8 2329 2337
## INFLUENZA (SEASONAL) (FLUARIX QUADRIVALENT) 24 8334 8358
## INFLUENZA (SEASONAL) (FLUARIX) 19 4612 4631
## INFLUENZA (SEASONAL) (FLUBLOK QUADRIVALENT) 3 1739 1742
## INFLUENZA (SEASONAL) (FLUBLOK) 198 198
## INFLUENZA (SEASONAL) (FLUCELVAX FOREIGN (OPTAFLU) 5 5
## INFLUENZA (SEASONAL) (FLUCELVAX QUADRIVALENT) 6 4379 4385
## INFLUENZA (SEASONAL) (FLUCELVAX) 2 990 992
## INFLUENZA (SEASONAL) (FLUENZ TETRA) 10 10
## INFLUENZA (SEASONAL) (FLUENZ) 1 1
## INFLUENZA (SEASONAL) (FLULAVAL QUADRIVALENT) 15 3369 3384
## INFLUENZA (SEASONAL) (FLULAVAL) 22 3150 3172
## INFLUENZA (SEASONAL) (FLUMIST QUADRIVALENT) 10 2973 2983
## INFLUENZA (SEASONAL) (FLUMIST) 20 5798 5818
## INFLUENZA (SEASONAL) (FLUOGEN) 44 1765 1809
## INFLUENZA (SEASONAL) (FLUSHIELD) 45 4197 4242
## INFLUENZA (SEASONAL) (FLUVAX) 3 3
## INFLUENZA (SEASONAL) (FLUVIRIN) 96 17963 18059
## INFLUENZA (SEASONAL) (FLUZONE HIGH-DOSE QUADRIVALENT) 11 1921 1932
## INFLUENZA (SEASONAL) (FLUZONE HIGH-DOSE) 59 16297 16356
## INFLUENZA (SEASONAL) (FLUZONE INTRADERMAL QUADRIVALENT) 235 235
## INFLUENZA (SEASONAL) (FLUZONE INTRADERMAL) 1 1162 1163
## INFLUENZA (SEASONAL) (FLUZONE QUADRIVALENT) 46 12151 12197
## INFLUENZA (SEASONAL) (FLUZONE) 324 35551 35875
## INFLUENZA (SEASONAL) (FOREIGN) 14 14
## INFLUENZA (SEASONAL) (NO BRAND NAME) 288 19839 20127
## INFLUENZA (SEASONAL) (OPTAFLU) 2 2
## INFLUENZA (SEASONAL) (PANDEMIC FLU VACCINE (H5N1) 5 5
## INFLUENZA (SEASONAL) (QIV DRESDEN) 9 54 63
## INFLUENZA (SEASONAL) (QIV QUEBEC) 1 1
## INFLUENZA (SEASONAL) (TIV DRESDEN) 38 38
## JAPANESE ENCEPHALITIS (IXIARO) 1 364 365
## JAPANESE ENCEPHALITIS (J-VAX) 2 187 189
## JAPANESE ENCEPHALITIS (JE-VAX) 295 295
## JAPANESE ENCEPHALITIS (NO BRAND NAME) 37 37
## LYME (LYMERIX) 9 2215 2224
## MEASLES (ATTENUVAX) 7 562 569
## MEASLES (FOREIGN) 3 3
## MEASLES (NO BRAND NAME) 1 65 66
## MEASLES + MUMPS (MM-VAX) 53 53
## MEASLES + MUMPS (NO BRAND NAME) 13 13
## MEASLES + MUMPS + RUBELLA (MMR I) 481 481
## MEASLES + MUMPS + RUBELLA (MMR II) 238 69348 69586
## MEASLES + MUMPS + RUBELLA (NO BRAND NAME) 12 3798 3810
## MEASLES + MUMPS + RUBELLA (PRIORIX) 97 97
## MEASLES + MUMPS + RUBELLA (VIRIVAC) 2 1575 1577
## MEASLES + MUMPS + RUBELLA + VARICELLA (NO BRAND NAME) 2 343 345
## MEASLES + MUMPS + RUBELLA + VARICELLA (PROQUAD) 20 15257 15277
## MEASLES + RUBELLA (MR-VAX II) 156 156
## MEASLES + RUBELLA (NO BRAND NAME) 5 5
## MENINGICOCCAL C & Y + HIB (MENHIBRIX) 8 8
## MENINGICOCCAL CONJUGATE + HIB (FOREIGN) 2 2
## MENINGICOCCAL CONJUGATE + HIB (UNKNOWN) 8 8
## MENINGOCOCCAL (MENINGITEC) 1 1
## MENINGOCOCCAL (MENOMUNE) 10 1409 1419
## MENINGOCOCCAL (MENOMUNE-A/C) 1 68 69
## MENINGOCOCCAL (NO BRAND NAME) 28 1175 1203
## MENINGOCOCCAL B (BEXSERO) 6 2979 2985
## MENINGOCOCCAL B (TRUMENBA) 2853 2853
## MENINGOCOCCAL C & Y + HIB (MENHIBRIX) 13 13
## MENINGOCOCCAL CONJUGATE (MENACTRA) 39 17285 17324
## MENINGOCOCCAL CONJUGATE (MENVEO) 5 5674 5679
## MENINGOCOCCAL CONJUGATE (NO BRAND NAME) 76 76
## MENINGOCOCCAL CONJUGATE + HIB (MENITORIX) 5 5
## MENINGOCOCCAL CONJUGATE + HIB (UNKNOWN) 4 4
## MUMPS (MUMPSVAX I) 2 44 46
## MUMPS (MUMPSVAX II) 2 90 92
## MUMPS (NO BRAND NAME) 17 17
## MUMPS + RUBELLA (FOREIGN) 27 27
## MUMPS + RUBELLA (NO BRAND NAME) 2 2
## PERTUSSIS (NO BRAND NAME) 3 140 143
## PLAGUE (NO BRAND NAME) 2 15 17
## PLAGUE (USP) 1 10 11
## PNEUMO (FOREIGN) 21 21
## PNEUMO (NO BRAND NAME) 24 4472 4496
## PNEUMO (PNEUMOVAX) 147 51294 51441
## PNEUMO (PNU-IMUNE) 28 2911 2939
## PNEUMO (PREVNAR) 686 22630 23316
## PNEUMO (PREVNAR13) 384 26207 26591
## PNEUMO (SYNFLORIX) 16 16
## POLIO VIRUS, INACT. (IPOL) 239 23372 23611
## POLIO VIRUS, INACT. (NO BRAND NAME) 86 4275 4361
## POLIO VIRUS, INACT. (POLIOVAX) 249 8128 8377
## POLIO VIRUS, ORAL (NO BRAND NAME) 51 1252 1303
## POLIO VIRUS, ORAL (ORIMUNE) 787 21581 22368
## RABIES (IMOVAX ID) 1 300 301
## RABIES (IMOVAX) 5 1923 1928
## RABIES (NO BRAND NAME) 3 344 347
## RABIES (RABAVERT) 4 1305 1309
## RABIES (RABIE-VAX) 49 49
## RABIES (RABIPUR) 2 51 53
## RABIES (RABIVAC) 25 25
## ROTAVIRUS (NO BRAND NAME) 6 871 877
## ROTAVIRUS (ROTARIX) 86 2294 2380
## ROTAVIRUS (ROTASHIELD) 10 683 693
## ROTAVIRUS (ROTATEQ) 367 18121 18488
## RUBELLA (MERUVAX I) 2 2
## RUBELLA (MERUVAX II) 809 809
## RUBELLA (NO BRAND NAME) 71 71
## RUBELLA (RUDIVAX) 4 4
## SMALLPOX (ACAM2000) 3 1876 1879
## SMALLPOX (DRYVAX) 13 3343 3356
## SMALLPOX (NO BRAND NAME) 1 499 500
## SUMMER/SPRING ENCEPH (NO BRAND NAME) 4 4
## TD ADSORBED (DITANRIX) 35 35
## TD ADSORBED (NO BRAND NAME) 26 14514 14540
## TD ADSORBED (TD-RIX) 16 16
## TD ADSORBED (TDVAX) 80 80
## TD ADSORBED (TENIVAC) 1 476 477
## TDAP (ADACEL) 36 22060 22096
## TDAP (BOOSTRIX) 21 12581 12602
## TDAP (NO BRAND NAME) 9 2561 2570
## TDAP + IPV (FOREIGN) 7 7
## TETANUS DIPHTHERIA (NO BRAND NAME) 3 2134 2137
## TETANUS DIPHTHERIA (TENIVAC) 198 198
## TETANUS TOXOID (NO BRAND NAME) 8 2195 2203
## TETANUS TOXOID (TEVAX) 31 31
## TETANUS TOXOID, ADSORBED (NO BRAND NAME) 2 968 970
## TICK-BORNE ENCEPH (NO BRAND NAME) 5 5
## TYPHOID LIVE ORAL TY21A (NO BRAND NAME) 94 94
## TYPHOID LIVE ORAL TY21A (VIVOTIF) 3 2327 2330
## TYPHOID VI POLYSACCHARIDE (ACETONE INACTIVATED DRIED) 1 166 167
## TYPHOID VI POLYSACCHARIDE (NO BRAND NAME) 4 1377 1381
## TYPHOID VI POLYSACCHARIDE (TYPHIM VI) 10 3781 3791
## VACCINE NOT SPECIFIED (FOREIGN) 9 9
## VACCINE NOT SPECIFIED (NO BRAND NAME) 69 5934 6003
## VACCINE NOT SPECIFIED (OTHER) 4 652 656
## VARICELLA (NO BRAND NAME) 4 2512 2516
## VARICELLA (VARILRIX) 1 173 174
## VARICELLA (VARIVAX) 154 75176 75330
## YELLOW FEVER (NO BRAND NAME) 4 264 268
## YELLOW FEVER (STAMARIL) 121 121
## YELLOW FEVER (YF-VAX) 19 2980 2999
## ZOSTER (NO BRAND NAME) 5 1099 1104
## ZOSTER (SHINGRIX) 47 41564 41611
## ZOSTER (ZOSTAVAX) 94 33582 33676
## ZOSTER LIVE (ZOSTAVAX) 39 6336 6375
## Total 12669 1111853 1124522
| Name | vaers_data |
| Number of rows | 759768 |
| Number of columns | 35 |
| Key | VAERS_ID |
| _______________________ | |
| Column type frequency: | |
| character | 22 |
| Date | 6 |
| numeric | 7 |
| ________________________ | |
| Group variables | None |
Variable type: character
| skim_variable | n_missing | complete_rate | min | max | empty | n_unique | whitespace |
|---|---|---|---|---|---|---|---|
| STATE | 122740 | 0.84 | 2 | 2 | 0 | 108 | 0 |
| SEX | 0 | 1.00 | 1 | 1 | 0 | 3 | 0 |
| SYMPTOM_TEXT | 308 | 1.00 | 1 | 26639 | 0 | 743266 | 0 |
| DIED | 752654 | 0.01 | 1 | 1 | 0 | 1 | 0 |
| L_THREAT | 748807 | 0.01 | 1 | 1 | 0 | 1 | 0 |
| ER_VISIT | 573298 | 0.25 | 1 | 1 | 0 | 1 | 0 |
| HOSPITAL | 717501 | 0.06 | 1 | 1 | 0 | 1 | 0 |
| X_STAY | 756575 | 0.00 | 1 | 1 | 0 | 1 | 0 |
| DISABLE | 746984 | 0.02 | 1 | 1 | 0 | 1 | 0 |
| RECOVD | 43644 | 0.94 | 1 | 1 | 0 | 3 | 0 |
| LAB_DATA | 442705 | 0.42 | 1 | 32000 | 0 | 125035 | 0 |
| V_ADMINBY | 1054 | 1.00 | 3 | 3 | 0 | 9 | 0 |
| V_FUNDBY | 203963 | 0.73 | 3 | 3 | 0 | 5 | 0 |
| OTHER_MEDS | 356307 | 0.53 | 1 | 247 | 0 | 140276 | 0 |
| CUR_ILL | 267601 | 0.65 | 1 | 4668 | 0 | 60142 | 0 |
| HISTORY | 288263 | 0.62 | 1 | 10374 | 0 | 161474 | 0 |
| PRIOR_VAX | 708934 | 0.07 | 1 | 132 | 0 | 26268 | 0 |
| SPLTTYPE | 457825 | 0.40 | 1 | 25 | 0 | 281418 | 0 |
| BIRTH_DEFECT | 759635 | 0.00 | 1 | 1 | 0 | 1 | 0 |
| OFC_VISIT | 712169 | 0.06 | 1 | 1 | 0 | 1 | 0 |
| ER_ED_VISIT | 737980 | 0.03 | 1 | 1 | 0 | 1 | 0 |
| ALLERGIES | 657041 | 0.14 | 1 | 2819 | 0 | 33203 | 0 |
Variable type: Date
| skim_variable | n_missing | complete_rate | min | max | median | n_unique |
|---|---|---|---|---|---|---|
| RECVDATE | 0 | 1.00 | 1990-07-02 | 2021-04-01 | 2013-01-28 | 9687 |
| RPT_DATE | 232626 | 0.69 | 1922-08-17 | 2021-02-13 | 2009-11-15 | 10731 |
| DATEDIED | 753786 | 0.01 | 1968-02-03 | 2021-04-01 | 2011-02-18 | 3343 |
| VAX_DATE | 93699 | 0.88 | 1901-12-28 | 2021-03-30 | 2012-05-02 | 12218 |
| ONSET_DATE | 134035 | 0.82 | 1920-11-29 | 2021-04-01 | 2012-05-08 | 12101 |
| TODAYS_DATE | 558896 | 0.26 | 2000-10-15 | 2021-04-01 | 2019-10-21 | 1416 |
Variable type: numeric
| skim_variable | n_missing | complete_rate | mean | sd | p0 | p25 | p50 | p75 | p100 | hist |
|---|---|---|---|---|---|---|---|---|---|---|
| VAERS_ID | 0 | 1.00 | 493465.15 | 280822.63 | 25001 | 254683.8 | 482767.5 | 732374.2 | 1157733 | ▇▇▇▇▂ |
| AGE_YRS | 126079 | 0.83 | 32.07 | 26.79 | 0 | 5.0 | 28.0 | 57.0 | 117 | ▇▃▅▂▁ |
| CAGE_YR | 184610 | 0.76 | 31.92 | 27.06 | -30 | 4.0 | 28.0 | 57.0 | 123 | ▂▇▅▃▁ |
| CAGE_MO | 642906 | 0.15 | 0.33 | 0.25 | 0 | 0.2 | 0.3 | 0.5 | 1 | ▇▅▃▁▁ |
| HOSPDAYS | 729830 | 0.04 | 6.17 | 17.99 | 0 | 2.0 | 3.0 | 5.0 | 999 | ▇▁▁▁▁ |
| NUMDAYS | 146512 | 0.81 | 42.46 | 402.95 | -2189 | 0.0 | 1.0 | 2.0 | 36896 | ▇▁▁▁▁ |
| FORM_VERS | 485117 | 0.36 | 1.74 | 0.44 | 1 | 1.0 | 2.0 | 2.0 | 2 | ▃▁▁▁▇ |
| Name | vaers_symptoms |
| Number of rows | 931884 |
| Number of columns | 11 |
| Key | VAERS_ID |
| _______________________ | |
| Column type frequency: | |
| character | 5 |
| numeric | 6 |
| ________________________ | |
| Group variables | None |
Variable type: character
| skim_variable | n_missing | complete_rate | min | max | empty | n_unique | whitespace |
|---|---|---|---|---|---|---|---|
| SYMPTOM1 | 22 | 1.00 | 3 | 88 | 0 | 6989 | 0 |
| SYMPTOM2 | 188393 | 0.80 | 3 | 88 | 0 | 6789 | 0 |
| SYMPTOM3 | 424937 | 0.54 | 3 | 88 | 0 | 6381 | 0 |
| SYMPTOM4 | 579215 | 0.38 | 3 | 88 | 0 | 5826 | 0 |
| SYMPTOM5 | 687630 | 0.26 | 3 | 88 | 0 | 5477 | 0 |
Variable type: numeric
| skim_variable | n_missing | complete_rate | mean | sd | p0 | p25 | p50 | p75 | p100 | hist |
|---|---|---|---|---|---|---|---|---|---|---|
| VAERS_ID | 0 | 1.00 | 500099.33 | 281227.94 | 25001.0 | 270614.8 | 480683.5 | 742128.2 | 1157733 | ▇▇▆▇▂ |
| SYMPTOMVERSION1 | 22 | 1.00 | 15.24 | 5.59 | 8.1 | 9.1 | 15.1 | 21.0 | 24 | ▇▃▃▃▆ |
| SYMPTOMVERSION2 | 188393 | 0.80 | 15.33 | 5.56 | 8.1 | 10.0 | 16.0 | 21.0 | 24 | ▇▃▃▅▆ |
| SYMPTOMVERSION3 | 424937 | 0.54 | 15.17 | 5.64 | 8.1 | 10.0 | 15.0 | 21.0 | 24 | ▇▃▃▃▆ |
| SYMPTOMVERSION4 | 579215 | 0.38 | 15.42 | 5.60 | 8.1 | 10.1 | 15.1 | 21.0 | 24 | ▇▅▃▃▇ |
| SYMPTOMVERSION5 | 687630 | 0.26 | 15.58 | 5.55 | 8.1 | 11.0 | 15.1 | 21.1 | 24 | ▇▅▃▃▇ |
| Name | vaers_vax |
| Number of rows | 1124896 |
| Number of columns | 9 |
| Key | VAERS_ID |
| _______________________ | |
| Column type frequency: | |
| character | 8 |
| numeric | 1 |
| ________________________ | |
| Group variables | None |
Variable type: character
| skim_variable | n_missing | complete_rate | min | max | empty | n_unique | whitespace |
|---|---|---|---|---|---|---|---|
| VAX_TYPE | 0 | 1.00 | 2 | 10 | 0 | 95 | 0 |
| VAX_MANU | 0 | 1.00 | 3 | 39 | 0 | 45 | 0 |
| VAX_LOT | 264874 | 0.76 | 1 | 15 | 0 | 136931 | 0 |
| VAX_DOSE | 578655 | 0.49 | 1 | 4 | 0 | 49 | 0 |
| VAX_ROUTE | 185420 | 0.84 | 1 | 3 | 0 | 16 | 0 |
| VAX_SITE | 221733 | 0.80 | 1 | 3 | 0 | 18 | 0 |
| VAX_NAME | 0 | 1.00 | 10 | 55 | 0 | 244 | 0 |
| VAX_DOSE_SERIES | 792504 | 0.30 | 1 | 3 | 0 | 8 | 0 |
Variable type: numeric
| skim_variable | n_missing | complete_rate | mean | sd | p0 | p25 | p50 | p75 | p100 | hist |
|---|---|---|---|---|---|---|---|---|---|---|
| VAERS_ID | 0 | 1 | 451117.4 | 273028.3 | 25001 | 221903.8 | 417843 | 672404.5 | 1157733 | ▇▇▆▆▁ |
# load the data sets
data(vaersND_data)
data(vaersND_symptoms)
data(vaersND_vax)
# use skimr to summarize each of the 3 Non-Domestic data sets
skim(vaersND_data)| Name | vaersND_data |
| Number of rows | 99613 |
| Number of columns | 35 |
| Key | NULL |
| _______________________ | |
| Column type frequency: | |
| character | 22 |
| Date | 6 |
| numeric | 7 |
| ________________________ | |
| Group variables | None |
Variable type: character
| skim_variable | n_missing | complete_rate | min | max | empty | n_unique | whitespace |
|---|---|---|---|---|---|---|---|
| STATE | 0 | 1.00 | 2 | 2 | 0 | 1 | 0 |
| SEX | 0 | 1.00 | 1 | 1 | 0 | 3 | 0 |
| SYMPTOM_TEXT | 1 | 1.00 | 3 | 31056 | 0 | 99258 | 0 |
| DIED | 95690 | 0.04 | 1 | 1 | 0 | 1 | 0 |
| L_THREAT | 95815 | 0.04 | 1 | 1 | 0 | 1 | 0 |
| ER_VISIT | 91745 | 0.08 | 1 | 1 | 0 | 1 | 0 |
| HOSPITAL | 59784 | 0.40 | 1 | 1 | 0 | 1 | 0 |
| X_STAY | 98210 | 0.01 | 1 | 1 | 0 | 1 | 0 |
| DISABLE | 92516 | 0.07 | 1 | 1 | 0 | 1 | 0 |
| RECOVD | 3696 | 0.96 | 1 | 1 | 0 | 3 | 0 |
| LAB_DATA | 43716 | 0.56 | 1 | 32000 | 0 | 35679 | 0 |
| V_ADMINBY | 0 | 1.00 | 3 | 3 | 0 | 8 | 0 |
| V_FUNDBY | 23198 | 0.77 | 3 | 3 | 0 | 4 | 0 |
| OTHER_MEDS | 52327 | 0.47 | 1 | 240 | 0 | 10382 | 0 |
| CUR_ILL | 52529 | 0.47 | 1 | 1847 | 0 | 8410 | 0 |
| HISTORY | 59251 | 0.41 | 1 | 3010 | 0 | 19156 | 0 |
| PRIOR_VAX | 95534 | 0.04 | 11 | 128 | 0 | 725 | 0 |
| SPLTTYPE | 693 | 0.99 | 4 | 25 | 0 | 91684 | 0 |
| BIRTH_DEFECT | 99568 | 0.00 | 1 | 1 | 0 | 1 | 0 |
| OFC_VISIT | 98385 | 0.01 | 1 | 1 | 0 | 1 | 0 |
| ER_ED_VISIT | 98655 | 0.01 | 1 | 1 | 0 | 1 | 0 |
| ALLERGIES | 99399 | 0.00 | 2 | 145 | 0 | 113 | 0 |
Variable type: Date
| skim_variable | n_missing | complete_rate | min | max | median | n_unique |
|---|---|---|---|---|---|---|
| RECVDATE | 0 | 1.00 | 1990-07-02 | 2021-02-26 | 2014-06-24 | 6448 |
| RPT_DATE | 23807 | 0.76 | 1988-10-26 | 2019-08-27 | 2012-09-26 | 6367 |
| DATEDIED | 97761 | 0.02 | 1974-10-05 | 2021-01-29 | 2012-11-15 | 1502 |
| VAX_DATE | 23926 | 0.76 | 1955-10-20 | 2021-01-28 | 2012-09-14 | 9201 |
| ONSET_DATE | 30408 | 0.69 | 1955-09-01 | 2021-01-28 | 2013-02-01 | 9398 |
| TODAYS_DATE | 76418 | 0.23 | 2017-06-30 | 2021-02-24 | 2019-04-24 | 1127 |
Variable type: numeric
| skim_variable | n_missing | complete_rate | mean | sd | p0 | p25 | p50 | p75 | p100 | hist |
|---|---|---|---|---|---|---|---|---|---|---|
| VAERS_ID | 0 | 1.00 | 536415.28 | 213347.88 | 25002 | 383778.00 | 534804.0 | 721165.0 | 1055830 | ▂▆▇▇▂ |
| AGE_YRS | 28681 | 0.71 | 18.65 | 23.56 | 0 | 0.75 | 9.0 | 29.0 | 103 | ▇▂▁▁▁ |
| CAGE_YR | 70771 | 0.29 | 16.26 | 21.82 | 0 | 0.00 | 5.0 | 27.0 | 100 | ▇▂▁▁▁ |
| CAGE_MO | 86302 | 0.13 | 0.34 | 0.22 | 0 | 0.20 | 0.3 | 0.4 | 1 | ▇▆▂▁▁ |
| HOSPDAYS | 87368 | 0.12 | 4.38 | 18.21 | 0 | 0.00 | 2.0 | 5.0 | 1096 | ▇▁▁▁▁ |
| NUMDAYS | 35233 | 0.65 | 70.30 | 382.34 | 0 | 0.00 | 1.0 | 8.0 | 16443 | ▇▁▁▁▁ |
| FORM_VERS | 0 | 1.00 | 1.23 | 0.42 | 1 | 1.00 | 1.0 | 1.0 | 2 | ▇▁▁▁▂ |
| Name | vaersND_symptoms |
| Number of rows | 161306 |
| Number of columns | 11 |
| Key | NULL |
| _______________________ | |
| Column type frequency: | |
| character | 5 |
| numeric | 6 |
| ________________________ | |
| Group variables | None |
Variable type: character
| skim_variable | n_missing | complete_rate | min | max | empty | n_unique | whitespace |
|---|---|---|---|---|---|---|---|
| SYMPTOM1 | 0 | 1.00 | 3 | 88 | 0 | 6349 | 0 |
| SYMPTOM2 | 24539 | 0.85 | 3 | 63 | 0 | 6139 | 0 |
| SYMPTOM3 | 47696 | 0.70 | 3 | 63 | 0 | 5717 | 0 |
| SYMPTOM4 | 68500 | 0.58 | 3 | 66 | 0 | 5333 | 0 |
| SYMPTOM5 | 85992 | 0.47 | 3 | 66 | 0 | 4944 | 0 |
Variable type: numeric
| skim_variable | n_missing | complete_rate | mean | sd | p0 | p25 | p50 | p75 | p100 | hist |
|---|---|---|---|---|---|---|---|---|---|---|
| VAERS_ID | 0 | 1.00 | 537920.23 | 206779.06 | 25002.0 | 389004 | 535738.5 | 707393.0 | 1055830.0 | ▂▆▇▇▂ |
| SYMPTOMVERSION1 | 0 | 1.00 | 16.56 | 4.78 | 8.1 | 13 | 16.1 | 21.0 | 24.0 | ▆▅▆▇▇ |
| SYMPTOMVERSION2 | 24539 | 0.85 | 16.55 | 4.76 | 8.1 | 13 | 16.1 | 21.0 | 23.1 | ▆▅▆▇▇ |
| SYMPTOMVERSION3 | 47696 | 0.70 | 16.55 | 4.73 | 8.1 | 13 | 16.1 | 20.1 | 23.1 | ▆▅▆▇▇ |
| SYMPTOMVERSION4 | 68500 | 0.58 | 16.62 | 4.67 | 8.1 | 13 | 16.1 | 20.1 | 23.1 | ▆▅▆▇▇ |
| SYMPTOMVERSION5 | 85992 | 0.47 | 16.62 | 4.63 | 8.1 | 13 | 16.1 | 20.1 | 23.1 | ▆▅▆▇▇ |
| Name | vaersND_vax |
| Number of rows | 144725 |
| Number of columns | 8 |
| Key | NULL |
| _______________________ | |
| Column type frequency: | |
| character | 7 |
| numeric | 1 |
| ________________________ | |
| Group variables | None |
Variable type: character
| skim_variable | n_missing | complete_rate | min | max | empty | n_unique | whitespace |
|---|---|---|---|---|---|---|---|
| VAX_TYPE | 0 | 1.00 | 2 | 10 | 0 | 90 | 0 |
| VAX_MANU | 0 | 1.00 | 3 | 33 | 0 | 30 | 0 |
| VAX_LOT | 78974 | 0.45 | 1 | 15 | 0 | 30254 | 0 |
| VAX_DOSE_SERIES | 4 | 1.00 | 1 | 3 | 0 | 8 | 0 |
| VAX_ROUTE | 17867 | 0.88 | 1 | 4 | 0 | 15 | 0 |
| VAX_SITE | 38147 | 0.74 | 2 | 3 | 0 | 12 | 0 |
| VAX_NAME | 0 | 1.00 | 10 | 55 | 0 | 215 | 0 |
Variable type: numeric
| skim_variable | n_missing | complete_rate | mean | sd | p0 | p25 | p50 | p75 | p100 | hist |
|---|---|---|---|---|---|---|---|---|---|---|
| VAERS_ID | 0 | 1 | 536190.1 | 215641.8 | 25002 | 374347 | 534336 | 726722 | 1055830 | ▂▆▇▇▂ |
US Centers for Disease Control and Prevention (CDC) and the US Food and Drug Administration (FDA) Vaccine Adverse Event Reporting System (VAERS), https://vaers.hhs.gov/.
Source 1
r - Is there an aggregate FUN option to count
occurrences? - Stack Overflow answered by mnel on Aug 2 2013. See https://stackoverflow.com/questions/9809166/is-there-an-aggregate-fun-option-to-count-occurrences.
Source 2
https://rdatatable.gitlab.io/data.table/articles/datatable-intro.html
Introduction to data.table
Questioning the Universe Publishing (QUP) {Irucka Embry]: Resources
to help us rethink the current crisis (in perception):
“Coronavirus” Vaccination$ Experimental Gene
Therapy Technology {Search for discussion by Judy Mikovits,
Ph.D.] https://www.questionuniverse.com/rethink_current_crisis.html#cov_vaccin
Questioning the Universe Publishing (QUP) {Irucka Embry]: Resources
to help us rethink, reimagine, & reFeel our world:
Vaccination$
https://www.questionuniverse.com/rethink.html#vaccine
EcoC²S [Irucka Embry): Resources for a Healthier You (Resources to
help us take responsibility back to heal ourselves and each other):
Vaccination$
https://www.ecoccs.com/healing.html#vaccines