Funktion Risk_Pneumonie

Berechnet das Mortalitätsrisiko einer Pneumonie bei jedem Fall der Eingabe und fügt eine numerische Spalte mit dem Ergebniss hinzu.(Risikoadjustierter Indikator - D1_15_M - 14_26). Die Funktion getCoeff gibt sämtliche logit-Regressionskoeffizienten des ausgewählten Referenzjahres als log(odds) zurück. Die Koeffizienten der im Fall vorhandenen Risikofaktoren werden summiert und die Summe wird zum Risikowert konvertiert.

Risk_Pneumonie(x)

Arguments

x

ein tibble (siehe Format des Rückgabewertes von iqi.input)

Value

Mortalitätrisiko der Pneumonie

Spezifikation

rcoef <- RACoefficients %>% filter(Indic == 'D1_15') %>% select(-Indic) %>% deframe
logodds = rcoef['Intercept'] +
if_else(ageyears > 0, ageyears * rcoef['RF_alte'], 0) +
if_else(sex == '2', rcoef['RF_sex'], 0) +
if_else(any(sdx %in_table% R_ICD_hrs), rcoef['RF_hrs'], 0) +
if_else(any(sdx %in_table% R_ICD_hfa), rcoef['RF_hfa'], 0) +
if_else(any(sdx %in_table% R_ICD_khk), rcoef['RF_khk'], 0) +
if_else(any(sdx %in_table% R_ICD_hypt_orig), rcoef['RF_hypt'], 0) +
if_else(any(sdx %in_table% R_ICD_vit), rcoef['RF_vit'], 0) +
if_else(any(sdx %in_table% R_ICD_aex), rcoef['RF_aex'], 0) +
if_else(any(sdx %in_table% R_ICD_leb), rcoef['RF_leb'], 0) +
if_else(any(sdx %in_table% R_ICD_ren), rcoef['RF_ren'], 0) +
if_else(any(ddx %in_table% R_ICD_diab), rcoef['RF_diab'], 0) +
if_else(any(sdx %in_table% R_ICD_adi3), rcoef['RF_adi3'], 0) +
if_else(any(sdx %in_table% R_ICD_kex), rcoef['RF_kex'], 0) +
if_else(any(sdx %in_table% R_ICD_coa), rcoef['RF_coa'], 0) +
if_else(any(sdx %in_table% R_ICD_FSTR), rcoef['RF_fstr'], 0)
Risk_Pneumonie = exp(logodds)/(1 + exp(logodds))

Examples

if (FALSE)  iqi.input(path2file) %>% Risk_Pneumonie()  # \dontrun{}