Calculate the equilibrium of a linear facility model
equilib(M, init, mgf = NULL)
Arguments
- M
A matrix of state transition rates between facility patient states
- init
A vector of admission state probabilities to each state
- mgf
The moment generating function characterizing a time-of-stay-dependent removal hazard
Value
A vector with the proportion of patients in each state at equilibrium
Examples
M <- rbind(c(-0.06,0.03,0),c(0.06,-0.08,0),c(0,0.05,0))
init <- c(0.95,0.05,0)
mgf <- function(x, deriv=0) MGFgamma(x, rate = 0.05, shape = 2.5, deriv)
equilib(M, init, mgf)
#> [1] 0.3558121 0.2291196 0.4150683