getDailyMaxResults: getDailyMaxResults

In babsim.hospital: Bartz & Bartz Simulation Hospital

getDailyMaxResults

Description

Combine babsim simulation results with real (field) data. Use daily max resources. Input: simmer simulation environment and field data formatted using getRealBeds. The formatted filed date has dim (nxm, 5), the output data has dimension (nxm, 15). The method get_mon_resources function is used to extract information from the babsim.hospital simulation. The function is used by modelResultHospital to prepare the calculation of the error.

Usage

getDailyMaxResults(envs, fieldEvents, conf)

Arguments

envs	simmer simulation environment. Result from babsim.hospital simulation runs, e.g., output from babsimHospital.
fieldEvents	Real values. Output from getRealBeds, i.e., a (nxm, 5)-dim data.frame with the following variables: resource chr: 'bed' 'bed' 'bed' 'bed' ... time int: 1 2 3 4 5 6 7 8 9 10 ... med int: 2 2 3 3 3 3 3 3 4 4 ... source chr: 'GA' 'GA' 'GA' 'GA' ... date POSIXct, format: '2020-03-03 01:00:00' '2020-03-04 01:00:00' '2020-03-05 01:00:00' '2020-03-06 01:00:00' ...
conf	list with the following entries (generated with babsimToolsConf): seed seed. Change the seed value to get different output for the same input parameters. Default: 123 simRepeats simmer repeats parallel simmer parallel runs. Default: FALSE perCores percentage of cores used for parallel simmer simulations. Default: 0.5 (=50 percent) ICU use ICU infection data. Default: FALSE logLevel log leved (0 or 1). Default: 0 (no output) maxCapacity max capacity of resources. Default: 1e6 dataset char name of the data set. Default: 'GA' simulationDates list with StartDate and EndDate. Period that is used for the simulation (babsim, simmer). Default: list(StartDate = '2020-03-03', EndDate = '2020-06-24') fieldDates list with StartDate and EndDate. Period when real data is available (resource usage). Default: list(StartDate = '2020-03-03', EndDate = '2020-06-24') simulationData data frame. Data used for the simulation. Default: dataCovidBeds20200624 fieldEvents data frame. Data used for the evaluation (error). Default: GABeds220200624 resource vector with resource names. Default: c('bed', 'intensiveBed', 'intensiveBedVentilation')

Details

get_mon_resources returns state changes in resources:

• 'resource': resources name

• 'time': time instant of the event that triggered the state change

• 'server': server count

• 'queue': queue count

• 'capacity': capacity

• 'queue_size': queue size

• 'system': system count (server + queue). If no queues are used, system values equal server values.

• 'system_limit': system limit (capacity + queue_size)

Value

This function returns an env data frame (tibble [nxm, 15] (S3: grouped_df/tbl_df/tbl/data.frame)) with the following entries:

resource (chr)

name of the seized resource: 'bed' 'bed' 'bed' 'bed' ...

time (num)

time step: 3 10 12 13 14 15 15 15 15 16 ...

server (int)

server: 1 2 3 2 3 4 3 4 5 6 ...

limit (num)

limit: Inf Inf Inf Inf Inf ...

replication (int)

replication: 1 1 1 1 1 1 1 1 1 1 ...

upper (int)

upper: 1 2 3 2 3 5 5 5 5 7 ...

lower (int)

lower: 1 2 3 2 3 3 3 3 3 5 ...

med (num)

med: 1 2 3 2 3 4 4 4 4 6 ...

date (POSIXct)

time, format: yyyy-mm-dd hh:mm.ss

rwdate (POSIXct)

format: '2020-03-01' '2020-03-08' '2020-03-15' '2020-03-15' ...

source (chr)

name of the simulation that was used: 'babsim' 'babsim' 'babsim' 'babsim' ...

See Also

modelResultHospital

Examples

data <- getSyntheticData()
para <- babsimHospitalPara()
conf <- babsimToolsConf()
conf <- getConfFromData(
  conf = conf,
  simData = data$simData,
  fieldData = data$fieldData
)
arrivalTimes <- getArrivalTimes(data$simData$Infected)
envs <- babsimHospital(arrivalTimes = arrivalTimes, conf = conf, para = para)
fieldEvents <- getRealBeds(
  data = data$fieldData,
  resource = c("bed", "intensiveBed", "intensiveBedVentilation")
)
res <- getDailyMaxResults(envs = envs, fieldEvents = fieldEvents, conf = conf)

babsim.hospital documentation built on May 30, 2022, 9:05 a.m.