如何解决您可以在R沸腾事件模拟中检查到达属性吗?
我是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
) %>%
...
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。