-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathcarditis.youth.r
53 lines (37 loc) · 1.16 KB
/
carditis.youth.r
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
library(data.table)
library(lubridate)
a1 <- mergeDVS[VAX_TYPE == "COVID19" &
grepl("carditis",All_symptoms) &
AGE_YRS <= 20,
.(VAERS_ID,
VAX_TYPE,
VAX_MANU,
AGE_YRS,
VAX_DATE,
ONSET_DATE,
ONSETminusVAX=mdy(ONSET_DATE) - mdy(VAX_DATE),
All_symptoms,
abbrevSYMPTOM_TEXT=substr(SYMPTOM_TEXT,0,225))][
VAX_DATE !="" & mdy(VAX_DATE) > mdy("11/01/2020"),][order(mdy(VAX_DATE),mdy(ONSET_DATE))]
a2 <- mergeDVS[VAX_TYPE == "COVID19" &
grepl("carditis",SYMPTOM_TEXT) &
AGE_YRS <= 20,
.(VAERS_ID,
VAX_TYPE,
VAX_MANU,
AGE_YRS,
VAX_DATE,
ONSET_DATE,
ONSETminusVAX=mdy(ONSET_DATE) - mdy(VAX_DATE),
All_symptoms,
abbrevSYMPTOM_TEXT=substr(SYMPTOM_TEXT,0,225))][
VAX_DATE !="" & mdy(VAX_DATE) > mdy("11/01/2020"),][order(mdy(VAX_DATE),mdy(ONSET_DATE))]
# binding by row a1 & a2
rbind(a1,a2)[!duplicated(VAERS_ID),]
#searching for all 'carditis' in fields relevant then removing duplicates
# this returns widest possible range
all.carditis <- mergeDVS[VAX_TYPE == "COVID19" &
(grepl("carditis",All_symptoms) |
grepl("carditis",SYMPTOM_TEXT) |
grepl("carditis",LAB_DATA)),]
all.carditis[!duplicated(VAERS_ID)][between(AGE_YRS,10,30),.N]