Nothing
test_that("$N$ patients can stay home", {
skip_if_quicktest()
set.seed(1111)
require("simmer")
require("dplyr")
N <- 2
conf <- babsimToolsConf()
ex50I <- ex1InfectedDf
ex50A <- data.frame(getArrivalTimes(ex1InfectedDf$Infected))
ex50A <- data.frame(time = ex50A[1:N, ])
para <- babsimHospitalPara()
para$FactorPatientsInfectedToHospital <- 0.0
y <- babsimHospital(
arrivalTimes = ex50A,
conf = conf,
para = para
)
# value of attribute "No Hospital Required" should be equal with number of patients
attributes <- get_mon_attributes(y)
s <- attributes[attributes$key == "No Hospital Required", ]
expect_equal(N, nrow(s))
})
########################################################################################
test_that("$N$ patients arrive. Everybody goes to the normal station and is healed", {
skip_if_quicktest()
set.seed(1112)
require("simmer")
require("dplyr")
N <- 4
conf <- babsimToolsConf()
ex51I <- ex1InfectedDf
ex51A <- data.frame(getArrivalTimes(ex1InfectedDf$Infected))
ex51A <- data.frame(time = ex51A[1:N, ])
para <- babsimHospitalPara()
#
para$FactorPatientsInfectedToHospital <- 1.0
#
para$FactorPatientsHospitalToIntensive <- 0.0
para$FactorPatientsHospitalToVentilation <- 0.0
#
para$FactorPatientsNormalToDeath <- 0.0
para$FactorPatientsNormalToVentilation <- 0.0
para$FactorPatientsNormalToIntensive <- 0.0
y <- babsimHospital(
arrivalTimes = ex51A,
conf = conf,
para = para
)
# value of attribute "Healed" should be equal with number of patients
attributes <- get_mon_attributes(y)
s <- attributes[attributes$key == "Healed", ]
expect_equal(N, nrow(s))
})
#######################################################################################
test_that("$N$ patients arrive. Everybody goes to the normal station and dies", {
skip_if_quicktest()
set.seed(1113)
require("simmer")
require("dplyr")
N <- 4
conf <- babsimToolsConf()
ex52I <- ex1InfectedDf
ex52A <- data.frame(getArrivalTimes(ex1InfectedDf$Infected))
ex52A <- data.frame(time = ex52A[1:N, ])
para <- babsimHospitalPara()
#
para$FactorPatientsInfectedToHospital <- 1.0
#
# everybody goes to normal:
para$FactorPatientsHospitalToIntensive <- 0
para$FactorPatientsHospitalToVentilation <- 0
#
para$FactorPatientsNormalToDeath <- 1.0
para$FactorPatientsNormalToVentilation <- 0
para$FactorPatientsNormalToIntensive <- 1.0
y <- babsimHospital(
arrivalTimes = ex52A,
conf = conf,
para = para
)
# value of attribute "Dead" should be equal with number of patients
attributes <- get_mon_attributes(y)
s <- attributes[attributes$key == "Dead", ]
expect_equal(N, nrow(s))
})
#######################################################################################
test_that("$N$ patients arrive. Everybody goes to the normal station.
normal station -> death: 50%
normal station -> healed: 50%
", {
skip_if_quicktest()
set.seed(1114)
require("simmer")
require("dplyr")
N <- 8
conf <- babsimToolsConf()
ex53I <- ex1InfectedDf
ex53A <- data.frame(getArrivalTimes(ex1InfectedDf$Infected))
ex53A <- data.frame(time = ex53A[1:N, ])
para <- babsimHospitalPara()
#
para$FactorPatientsInfectedToHospital <- 1.0
#
# everybody goes to normal:
para$FactorPatientsHospitalToIntensive <- 0
para$FactorPatientsHospitalToVentilation <- 0
#
# 50% healthy, 50% death
para$FactorPatientsNormalToDeath <- 0.5
para$FactorPatientsNormalToVentilation <- 0
para$FactorPatientsNormalToIntensive <- 0
y <- babsimHospital(
arrivalTimes = ex53A,
conf = conf,
para = para
)
# value of attribute "Healed" should be less than N
# value of attribute "Dead" should be lass than N
attributes <- get_mon_attributes(y)
s <- attributes[attributes$key == "Healed", ]
expect_lt(nrow(s), N)
attributes <- get_mon_attributes(y)
s <- attributes[attributes$key == "Dead", ]
expect_lt(nrow(s), N)
})
#######################################################################################
test_that("$N$ patients generated. Real simulation factors", {
skip_if_quicktest()
set.seed(1115)
require("simmer")
require("dplyr")
N <- 400
conf <- babsimToolsConf()
ex54I <- ex1InfectedDf
ex54A <- data.frame(getArrivalTimes(ex1InfectedDf$Infected))
ex54A <- data.frame(time = ex54A[1:N, ])
para <- babsimHospitalPara()
y <- babsimHospital(
arrivalTimes = ex54A,
conf = conf,
para = para
)
# sum of attributes should be equal with number of patients
attributes <- get_mon_attributes(y)
expect_equal(N, nrow(attributes))
})
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.