Computes median and interquartile range of transit times for each
processing stage along the entity flow. If a reached_stageN
column is present in the entity data, computes the fraction of
entities completing each stage.
Usage
stage_throughput(data, stages = c("Stage1", "Stage2", "Stage3", "Stage4"))