facilityeq: Calculate the equilibrium of a facility transmission model

In facilityepimath: Analyze Mathematical Models of Healthcare Facility Transmission

Calculate the equilibrium of a facility transmission model

Description

Calculate the equilibrium of a facility transmission model

Usage

facilityeq(S, C, A, R, transm, init, mgf = NULL)

Arguments

S	A matrix of state transition rates between and removal from the susceptible states in the absence of colonized individuals
C	A matrix of state transition rates between and removal from the colonized states
A	A matrix describing transitions from susceptible to colonized states at acquisition
R	A matrix of recovery rates: state transition rates from colonized to susceptible states
transm	A vector of transmission rates from each colonized state
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; the vector contains the equilibrium S states followed by C states

Examples

S <- 0
C <- rbind(c(-0.38,0),c(0.08,0))
A <- rbind(1,0)
R <- cbind(0.3,0)
transm <- c(0.1,0.05)
init <- c(0.99,0.01,0)
mgf <- function(x, deriv=0) MGFgamma(x, rate=0.2, shape=3, deriv)
facilityeq(S, C, A, R, transm, init, mgf)

facilityepimath documentation built on April 4, 2025, 12:59 a.m.