Plague is an enzootic disease in the western United States, even though long-term persistent infections do not seem to occur. Enzootic persistence may occur as a function of dynamic interactions between flea vectors and transiently infected hosts, but the specific levels of vector competence, host competence, and transmission and recovery rates that would promote persistence and emergence among wild hosts and vectors are not known. We developed a mathematical model of enzootic plague in the western United States and implemented the model with the following objectives: 1) to use matrix manipulation within a classic susceptible-->infective-->resistant-->susceptible (SIRS) model framework to describe transmission of the plague bacterium Yersinia pestis among rodents and fleas in California, 2) to perform sensitivity analysis with model parameters and variables to indicate which values tended to dominate model output, and 3) to determine whether enzootic maintenance would be predicted with realistic parameter values obtained from the literature for Y. pestis in California rodents and fleas. The model PlagueSIRS was implemented in discrete time as a computer simulation incorporating environmental stochasticity and seasonality, by using matrix functions in the computer language R, allowing any number of rodent and flea species to interact through parasitism and disease transmission. Sensitivity analysis indicated that the model was sensitive to flea attack rate, host recovery rate, and rodent host carrying capacity but relatively insensitive to changes in the duration of latent infection in the flea, host and vector competence, flea recovery from infection, and host mortality attributable to plague. Realistic parameters and variable values did allow for the model to predict enzootic plague in some combinations, specifically when rodent species that were susceptible to infection but resistant to morbidity were parasitized by multiple poorly competent flea species, including some that were present year-round. This model could be extended to similar vectorborne disease systems and could be used iteratively with data collection in sylvatic plague studies to better understand plague persistence and emergence in nature.