Initially I combined all the datasets. I calculated the times between positive COVID test in healthcare worker (earliest positive test) and each of this healthcare worker intubation procedure. I kept those pair intubation - positive COVID test with the shortest time between them for each positive healthcare worker and this is what I got:
| Intubation to positive COVID test pairs | ||||||
|---|---|---|---|---|---|---|
| With lowest number of days from intubation to positive COVID test for each healthcare worker that tested positive | ||||||
| no | intubation_id | account_id | intubation_date | date_of_outcome | days_intub_to_pos_covid | patient_s_covid_19_status |
| 1 | 8388 | 5098 | 2020-11-08 | 2020-11-09 | 1 | Confirmed |
| 2 | 8282 | 5127 | 2020-10-31 | 2020-11-04 | 4 | Confirmed |
| 3 | 8263 | 5116 | 2020-10-21 | 2020-10-27 | 6 | Confirmed |
| 4 | 8195 | 5082 | 2020-10-21 | 2020-10-28 | 7 | Suspected |
| 5 | 8304 | 5068 | 2020-11-04 | 2020-11-11 | 7 | Confirmed |
| 6 | 8989 | 3729 | 2020-12-11 | 2020-12-18 | 7 | Confirmed |
| 7 | 8839 | 5132 | 2020-11-20 | 2020-11-30 | 10 | Confirmed |
| 8 | 8858 | 4439 | 2020-12-02 | 2020-12-12 | 10 | Confirmed |
| 9 | 8834 | 5175 | 2020-11-29 | 2020-12-11 | 12 | Confirmed |
| 10 | 8578 | 5130 | 2020-11-06 | 2020-11-29 | 23 | Confirmed |
| 11 | 8257 | 5114 | 2020-10-30 | 2020-11-27 | 28 | Confirmed |
| 12 | 8917 | 5259 | 2020-10-09 | 2020-11-07 | 29 | Confirmed |
| 13 | 8166 | 5065 | 2020-10-22 | 2020-12-18 | 57 | Confirmed |
| 14 | 3503 | 3733 | 2020-04-17 | 2020-11-23 | 220 | Suspected |
Should we remove those with patient’s COVID_19 status suspected?
57 and 220 days from exposure to positive test doesn’t make too much sense so I filter out those records and we are left with 12 healthcare workers that have positive COVID after intubation. I will be using these 12 workers for the rest of analysis.
Also, on the top of removing 2 healthcare workers that had + test > 30 days from intubation, I removed all healthcare workers that did not have any follow up (we don’t know if they are COVID+ or COVID-, so cannot put them in any of these categories) and we are left with total of 62
Here is our healthcare worker data divided into those with outcome and no outcome:
| Characteristic | Overall, N = 621 | No Outcome, N = 501 | Outcome, N = 121 | p-value2 |
|---|---|---|---|---|
| age | 38 (8) | 38 (8) | 39 (9) | 0.6 |
| gender | 0.5 | |||
| Female | 26 (42%) | 20 (40%) | 6 (50%) | |
| Male | 36 (58%) | 30 (60%) | 6 (50%) | |
| specialty | 0.4 | |||
| Anaesthetics | 42 (68%) | 34 (68%) | 8 (67%) | |
| Emergency Medicine | 1 (1.6%) | 1 (2.0%) | 0 (0%) | |
| Intensive Care Medicine | 18 (29%) | 15 (30%) | 3 (25%) | |
| Pre-Hospital Emergency Medicine | 1 (1.6%) | 0 (0%) | 1 (8.3%) | |
| grade | 0.5 | |||
| Advanced Critical Care Practitioner (ACCP) | 1 (1.6%) | 1 (2.0%) | 0 (0%) | |
| Cons / SAS / Attending | 33 (53%) | 26 (52%) | 7 (58%) | |
| Other | 1 (1.6%) | 1 (2.0%) | 0 (0%) | |
| Paramedic | 1 (1.6%) | 0 (0%) | 1 (8.3%) | |
| Registered Nurse | 1 (1.6%) | 1 (2.0%) | 0 (0%) | |
| Trainee grades | 25 (40%) | 21 (42%) | 4 (33%) | |
| intubator | 58 (94%) | 47 (94%) | 11 (92%) | >0.9 |
|
1
Mean (SD); n (%)
2
Wilcoxon rank sum test; Pearson's Chi-squared test; Fisher's exact test
|
||||
Will do Kaplan-Mayers soon, altough I don’t think we will get a lot of significance in Cox regression because we have so few + outcomes.
How would you count “Procedures with outcome” as part of your table 1, as we do not really know which intubation caused the outcome? You had n=14 out of 292 all intubations in your table
Please double check my analysis with the raw data as I there might be mistakes!
And let me know what changes you would suggest and what you would like me to do next.
| Most common PPD combinations | ||||||||||
|---|---|---|---|---|---|---|---|---|---|---|
| For more or equal than 5 intubation procedures | ||||||||||
| eyewear | hat | gown | apron | gloves | surgical_mask | ffp2 | ffp3 | papr | plastic_drape_box | count |
| 1 | 1 | 1 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 90 |
| 1 | 1 | 1 | 1 | 1 | 0 | 0 | 1 | 0 | 0 | 36 |
| 1 | 1 | 1 | 0 | 1 | 0 | 1 | 0 | 0 | 0 | 29 |
| 1 | 1 | 1 | 1 | 1 | 0 | 1 | 0 | 0 | 0 | 18 |
| 1 | 0 | 1 | 0 | 1 | 1 | 0 | 1 | 0 | 0 | 15 |
| 1 | 1 | 0 | 1 | 1 | 0 | 0 | 1 | 0 | 0 | 13 |
| 1 | 1 | 1 | 0 | 1 | 1 | 1 | 0 | 0 | 0 | 11 |
| 1 | 0 | 1 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 7 |
| 1 | 1 | 1 | 0 | 1 | 1 | 0 | 1 | 0 | 0 | 7 |
| 1 | 0 | 1 | 1 | 1 | 0 | 0 | 1 | 0 | 0 | 6 |
| 1 | 1 | 0 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 6 |
| 1 | 1 | 1 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 5 |
Or do we want to do combinations with less PPD elements than we have now (we have 10)?
| Symptoms in healthcare workers with + follow ups | |||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| no | account_id | fever | cough | sore_throat | sob | headache | photophobia | myalgia | fatigue | nausea_vomiting | diarrhoea | abdominal_pain | loss_of_smell_taste |
| 1 | 3729 | 0 | 0 | 1 | 0 | 0 | 0 | 1 | 1 | 0 | 0 | 0 | 0 |
| 2 | 3729 | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 1 | 0 | 0 | 0 | 0 |
| 3 | 3729 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 0 | 0 |
| 4 | 3733 | 1 | 1 | 0 | 1 | 1 | 0 | 0 | 1 | 0 | 0 | 0 | 0 |
| 5 | 4439 | 0 | 0 | 1 | 0 | 1 | 0 | 1 | 1 | 0 | 1 | 0 | 0 |
| 6 | 5065 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 1 |
| 7 | 5068 | 0 | 1 | 1 | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 1 |
| 8 | 5082 | 1 | 1 | 0 | 0 | 1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
| 9 | 5098 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
| 10 | 5114 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
| 11 | 5114 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
| 12 | 5116 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 0 | 0 |
| 13 | 5116 | 0 | 0 | 0 | 0 | 1 | 0 | 0 | 1 | 0 | 0 | 0 | 0 |
| 14 | 5127 | 1 | 0 | 0 | 0 | 1 | 0 | 1 | 1 | 0 | 0 | 0 | 0 |
| 15 | 5130 | 1 | 0 | 0 | 0 | 1 | 0 | 1 | 1 | 0 | 0 | 0 | 1 |
| 16 | 5132 | 1 | 1 | 0 | 0 | 1 | 0 | 1 | 1 | 0 | 0 | 0 | 1 |
| 17 | 5175 | 0 | 1 | 0 | 0 | 1 | 0 | 1 | 1 | 0 | 0 | 0 | 1 |
| 18 | 5259 | 1 | 1 | 0 | 1 | 0 | 1 | 0 | 1 | 0 | 1 | 0 | 1 |
Some healthcare have different symptoms on different follow up records (ex account_id=3729), so I guess we have to mark the symptom as “1” if that person had it in at least 1 follow up?