Komorbiditet
Komorbiditet som covariat: enkelte diagnoser eller en samlet score
Warning
Under udvikling. Strukturel skitse - udbygges med konkrete eksempler.
Komorbiditet ved baseline er en hyppig covariat. Hvordan du håndterer den hører til din analyseplan og styres typisk af din DAG (se Fase 1 - Vælg dine kovariater med en DAG). Der er grundlæggende to tilgange:
- Justér for enkelte komorbiditeter som separate covariater, fx specifikke diagnoser du i din DAG har identificeret som confoundere. Dem udtrækker du som almindelige diagnose-variable (samme LPR-mønster som Udfald).
- Justér for den samlede komorbiditetsbyrde med ét sammenfattende indeks, fx NMI (Nordic Multimorbidity Index), en valideret, vægtet multimorbiditets-score. Brug et færdigt, valideret indeks frem for at kode det fra bunden, se NMI og Algoritmer & specialpakker.
Hvad er bedst: enkelte komorbiditeter eller en samlet score?
Det afhænger af din DAG og dit outcome:
- Enkelte komorbiditeter giver mest kontrol og gennemsigtighed, men kræver at du på forhånd ved præcis hvilke diagnoser der er confoundere. Brug det, når din DAG peger på nogle få, specifikke komorbiditeter.
- En samlet score (fx NMI) fanger den generelle comorbiditetsbyrde i én variabel og er praktisk, når mange komorbiditeter er potentielle confoundere, eller når du blot vil beskrive byrden i Table 1.
- Pas på overjustering: en score som NMI er udviklet til at forudsige morbiditets- og mortalitetsoutcomes, og dit eget outcome må ikke indgå i den score du justerer for. Hvis dit outcome (eller en stærk prædiktor for det) er en af de diagnoser scoren vægter, justerer du delvist for outcomet selv. NMI håndterer fx demens-studier ved at udelade demens-prædiktoren (se NMI).
Bruger du et samlet indeks, er mønstret: tag dine LPR-diagnoser (filtreret til kohorten og til tiden før index), kør dem gennem algoritmen, og få en score per pnr, som du gemmer som .rds og kobler på i Fase 12.
kohort_pnrs <- unique(readRDS("sti/til/full_cohort.rds")$pnr)
diagnoser_foer_index <- open_dataset("sti/til/lpr_diag/") %>%
rename_with(tolower) %>%
semi_join(tibble(pnr = kohort_pnrs), by = "pnr") %>% # KUN din kohorte
# ... hent kontaktdato, behold diagnoser før index_date ...
collect()
# → giv diagnoserne til NMI-algoritmen for at få en score per pnr