We introduce a discrete‐time anthrax disease model in human and herbivore populations, where the herbivore recruitment function is the classic Beverton–Holt model while that of humans is a constant recruitment function. We use the next generation method to compute the basic reproduction numbers for the submodels consisting of only herbivores
(
scriptR
0
A
) and only humans
(
scriptR
0
H
). In addition, we compute it for the full human–herbivore model
(
scriptR
0
A
H
). When
R
0
A
<
1 or
R
0
H
<
1 or
R
0
A
H
<
1, the number of anthrax infections decreases and the disease does not invade the corresponding population. Whereas, the number of anthrax infections increases and the disease invades the corresponding population when
R
0
A
>
1 or
R
0
H
>
1 or
R
0
A
H
>
1. Our simulations show that it is possible for the anthrax disease to invade a human–herbivore population where the disease does not invade either herbivores or humans in isolation. Recommendations for resource managers
To help guide anthrax disease prevention and control strategies, we provide basic reproduction numbers for the herbivore only population
(
scriptR
0
A
), human‐only population
(
scriptR
0
A
) and human–herbivore population
(
scriptR
0
A
H
).
R
0
<
1 implies the number of anthrax infections deceases and anthrax does not invade the population, while
R
0
>
1 implies the number of anthrax cases increases and anthrax invasion occurs.
It is possible for anthrax to invade a population consisting of herbivores and humans where the disease does not invade either the herbivores or humans in isolation. Thus, more anthrax surveillance is required in mixed populations of herbivores and humans.
Our human–herbivore model can be used to predict the number of anthrax‐infected humans and animals given the number of known anthrax infections in populations.