Skip to 0 minutes and 10 seconds In this video, we will do some mathematics. But don’t worry, if you don’t understand the logic behind this at once, important is only to know that we use the following formulas and matrices to visualise growth rates of cattle and to understand how diseases can influence the population of a herd. And please feel free to pause the video at every moment to have a look at each individual slide. So let’s start with some definitions. Animal populations can be represented using
Skip to 0 minutes and 50 seconds the two aspects: age and sex. Here you see a representation of a cattle population that is divided in females and males
Skip to 1 minute and 1 second and three age classes: newborn young ones, subadult juveniles which are sexually immature, and adult ones, cows and bulls who reproduce. The arrows indicate the birth and death of animals and the transitions between age classes. The calculation of the dynamics of such a population is straightforward.
Skip to 1 minute and 30 seconds The different classes of a population can be considered as a vector called the population vector N. Here we look only at the female group. The different age classes can be represented as Y for young ones, S for subadults, and A for adults. Further, we have a so-called projection matrix that describes the demographic behaviour of the population. Such a projection matrix has a typical structure. It consists of the birth rate, the survival rate of the young animals, and the persistence of the subadults and adults. In the top right corner, you have the birth rate b. On the diagonal, you have the rate
Skip to 2 minutes and 21 seconds of persistence in the same class: for subadults ps, and for adults pa. This is the inverse of the duration in that class. On the first subdiagonal is the survival rate of the young ones, sy, and the subadult, ss.
Skip to 2 minutes and 47 seconds If we want to project a population N at a certain time t plus 1, we can take the population at a given time t and multiply it with the projection matrix P. After several iterations, typically 20 to 30, we can write N t plus 1 equal Pn times Nt. The projection matrix can be written as Pn. After several multiplications with itself, the matrix reduces itself to diagonal values of lambda called the ’dominant Eigenvalue’ of the projection matrix P. Since P has now become a diagonal matrix, we can replace it with lambda to the power of n. After several iterations– 20 to 30, for example– the population structure becomes stable. The proportions between the vector elements are constant.
Skip to 3 minutes and 54 seconds The stabilised population vector becomes the Eigenvector. We can standardise the Eigenvector to a length of 1, which is also called the ’normed Eigenvector’. Let us develop a simple matrix model for a population of female cattle, with the following numbers. We have 50,000 calves, 20,000 heifers, and 40,000 cows. The time step is one year. The annual birth rate is 60% calving per year. The survival rate of the calves is 70%, which means 30% die or are taken off per year. Heifers stay two years as heifers. Their persistence rate is 50%. The survival rate of heifers is 60%, and the persistence of cows is 65%. These figures are not real and only for demonstration purposes.
Skip to 5 minutes and 4 seconds If you don’t have a spreadsheet application at hand, you can solve the matrix as follows. To see a projection of the population of the next year t plus 1, you have to sum up the multiplication of every element of the matrix rows with an element of the vector. Let us do it together. 0 times Y is 0, plus 0 times S is 0, plus b times A results in bA, which is understandable as multiplying the annual birth rate with the number of cows. In the same way, you deal with the other roles.
Skip to 5 minutes and 54 seconds This is a basic matrix calculus, which you may also find in other sources.
Skip to 6 minutes and 3 seconds Let’s solve the first matrix together. Y is equal 24,000, S is equal 45,000, and A is 38,000. Now, solve the first seven years by using a calculator. Make sure to always multiply the vector of the following year with the same values of the birth rate, the survival rate, and so on. You can download a PDF with the solutions to compare your results. If you have a table calculator programme, you can download the spreadsheet file afterwards and continue the rows to the year 21. What do you notice in this year? And what happens after year 10? Write you insights in the comments section.
Solving the formula (part one): introducing the matrix
If you want to convince with your arguments in the field of One Health, it is indispensable that you have calculations and graphs supporting your results. In this video Jakob Zinsstag explains some of the theoretical work when dealing with One Health.
Formulae and matrices are used to describe growth rates of cattle and to indicate how diseases can influence the population of a herd. If you don’t understand the logic immediately, do not be alarmed. As you finish the video in the following step, the theoretical aspects described at the beginning of this video should be much clearer. Feel free to pause the video and to retrace its content with your own calculations.
At the end of the video we ask you to finish the calculations in the chart for the first seven years by only using a calculator. Make sure to round the digit after a decimal point up (≥ 0,5) or down (< 0,5) to get correct integers. You find a document with the solutions in the ‘downloads’ section. If you have a table calculator programme, you can download the spreadsheet file afterwards and continue to year 21.
Note: if you do not have your own spreadsheet calculator programme you might upload it to a free viewer like this. (If you calculate using this free online app, be aware that when you have a lot of formulas and calculations going on in your sheet, a change of variables may need up to 15 seconds to affect all the relevant cells).
What do you notice in year 21? What happens after 10 years? Write your insights in the ‘comments’ section below.
© University of Basel