# using the example sepsis data set
phoenix_cardiovascular(
   vasoactives = dobutamine + dopamine + epinephrine + milrinone + norepinephrine + vasopressin,
   lactate = lactate,
   age = age,
   map = dbp + (sbp - dbp)/3,
   data = sepsis
)
# example data set to get all the possible scores
DF <-
  expand.grid(vasos = c(NA, 0:6),
              lactate = c(NA, 3.2, 5, 7.8, 11, 14), # units of mmol/L
              age = c(NA, 0.4, 1, 3, 12, 18, 24, 45, 60, 61, 144, 145), # months
              map = c(NA, 16:52)) # mmHg
DF$card <- phoenix_cardiovascular(vasos, lactate, age, map, DF)
head(DF)
# what if lactate is unknown for all records? - set the value either in the
# data object or the arguement value to NA
DF2 <-
  expand.grid(vasos = c(NA, 0:6),
              age = c(NA, 0.4, 1, 3, 12, 18, 24, 45, 60, 61, 144, 145), # months
              map = c(NA, 16:52)) # mmHg
DF2$card <- phoenix_cardiovascular(vasos, lactate = NA, age, map, DF2)
DF3 <-
  expand.grid(vasos = c(NA, 0:6),
              lactate = NA, # mmol/L
              age = c(NA, 0.4, 1, 3, 12, 18, 24, 45, 60, 61, 144, 145), # months
              map = c(NA, 16:52)) # mmHg
DF3$card <- phoenix_cardiovascular(vasos, lactate, age, map, DF3)
identical(DF2$card, DF3$card)
Run the code above in your browser using DataLab