问题描述
我是R的新手,而离散事件仿真的新手。 r.simmer似乎是一个很好的解决方案,因此尝试使用它为医院病床入院策略建模,对我来说,轨迹分支选项能够同时检查到达的资源状态和套利属性对我来说似乎很自然。
library(simmer)
daypatient<-trajectory() %>%
log_("sent home") %>%
timeout(300) %>%
seize("admin",1) %>%
release("admin",1) %>%
seize("bed",1)
patient <-trajectory() %>%
log_("start pat") %>%
branch(
function() rnorm(1)>0.5,# actually want to inspect arrival priority and queues to decide which trajectory path to follow
continue= TRUE,daypatient
) %>%
seize("bed",1) %>%
seize("theatre",1) %>%
timeout(rnorm(1,5)) %>%
release("theatre",1) %>%
timeout(50) %>%
release("bed",1)
# some patients need treatment pdq...
serIoUs<-data.frame(time=rexp(10,1),priority=2)
# ... other less so
wounded<-data.frame(time=rexp(10,priority=1)
hos <- simmer() %>%
add_resource("bed",2) %>%
add_resource("admin",3) %>%
add_resource("theatre",1) %>%
add_dataframe("serIoUs",patient,serIoUs) %>%
add_dataframe("wounded",wounded)
hos %>% run(until = 500)
有什么办法做到这一点,还是我从根本上误解了r.simmer / DES框架?
解决方法
毅力发现它
...
branch(
function () {
crit<-get_attribute(hos,"crit")
avail<-get_queue_count(hos,"bed")
if( sendhomeoption && avail==0 && crit == 1 ){
return (TRUE)
} else {
return(FALSE)
}
},# actually want to inspect arrival and queues
continue= c(TRUE,FALSE),daypatient,immediate
) %>%
...