modelResultHospital: modelResultHospital

Description

Simulate one parameter configuration from babsimHospitalPara. The simulation is by default deterministic, because conf$seed is used for set.seed.

Usage

modelResultHospital(para, conf, data)

Arguments

para

Simulation model parameters. The function babsimHospitalPara can be used to generate these parameters.

conf

list with the following entries:

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')

data

list with simData and fieldData

Value

This function returns an env data frame (tibble [560 <U+00D7> 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 ...
queue (int): 1 1 2 3 4 3 4 1 0 2 ...
capacity (num): 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 ...
queue_size (num): Inf Inf Inf Inf Inf ...
system (int): 1 1 2 3 4 3 4 1 0 2 ...
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' ...

Examples

Run this code

# NOT RUN {
# First example: OBK data
data <- getObkData()
para <- babsimHospitalPara()
para$GammaShapeParameter <- 0.8
# turn off parallelized simulation:
conf <- babsimToolsConf()
conf <- getConfFromData(
  conf = conf,
  simData = data$simData,
  fieldData = data$fieldData
)
# no logging (default)
conf$logLevel <- 0
res <- modelResultHospital(
  para = para,
  conf = conf,
  data = data
)

# Second example: synthetic data
data <- getSyntheticData()
para <- babsimHospitalPara()
conf <- babsimToolsConf()
conf <- getConfFromData(
  conf = conf,
  simData = data$simData,
  fieldData = data$fieldData
)
res <- modelResultHospital(para = para, conf = conf, data = data)
# }

Run the code above in your browser using DataLab