Applied and Computational Mathematics
Volume 4, Issue 5, October 2015, Pages: 363-368

Sensitivity Analysis of Parameters in a Competition Model

Frank Nathan Ngoteya, Yaw Nkansah-Gyekye

Applied Mathematics Department, the Nelson Mandela African Institute of Science and Technology, Arusha, Tanzania

(F. N. Ngoteya)
(Y. Nkansah-Gyekye)

Frank Nathan Ngoteya, Yaw Nkansah-Gyekye. Sensitivity Analysis of Parameters in a Competition Model.Applied and Computational Mathematics.Vol.4, No. 5, 2015, pp. 363-368. doi: 10.11648/j.acm.20150405.15

Abstract: Competition models may be used by ecologists as basic guidelines in analyzing issues that contribute to the decline of species population. To be used for that purpose, a mathematical model must be prudently parameterized. Therefore, this paper examines the effects of interference competition by lions with human–related mortality to the population dynamics of African wild dogs. The model was carefully parameterized and validated with estimated data by employing the theory of basic reproduction number and by running the sensitivity analysis. Numerical simulation of the model was executed by using MATLAB to explore the outcome of certain key parameters when changes are applied on them.

Keywords: Interference Competition, African Wild Dog, Lion, Basic Reproduction Number, Sensitivity Analysis

1. Introduction

The African wild dog was formerly present throughout sub-Saharan Africa, although it was never locally abundant. This species is today restricted to fragment population mainly in southern and eastern Africa (Woodroffe et al., 1997). Potentially viable populations currently exist in Botswana, Kenya, Mozambique, South Africa, Tanzania, Zambia and Zimbabwe (Sillero-Zubiri et al., 2004).

The status of the African wild dog is classified as endangered (EN) (IUCN Red List, 2008). Wild dogs have declined across their range; they are now extinct in 25 of the 39 countries in which they were formerly recorded, with the total population estimated to be 5,750 individuals in not more than 1,000 packs (Woodroffe et al., 2004) cited by (Dutson and Sillero-Zubiri, 2005).

Competition is defined as a negative interaction that occurs among organisms whenever two or more organisms require the same limited resource (Cronin and Carson, 2015). While, Interference competition is defined as direct interactions between organisms that decrease use of the common resources (Park, 1954) cited by (Nakayama and Fuiman, 2010).

A study in Australia conducted by (Caut et al., 2007) on rats and mice they employed a classical Lotka – Volterra two species competition model and they named the model as the Competitor Release Effect. They assumed that a rat is globally a better competitor over the mouse either by better exploiting resources or by generally winning interference interactions. They called rat species the superior competitor and the mouse species the inferior competitor. According to Caut et al. (2007), an inferior competitor can increase in numbers if the superior is controlled, due to the competitive interactions between them.

On the other hand, Blanco et al. (2014) studied population dynamics of wolves and coyotes at Yellowstone National Park, and the study was titled as modeling interference competition with an infectious disease. The goal of the study was to understand how different factors may affect the decline in population size of a dominant predator, the wolf and a subordinate predator, the coyote in relation to human-related mortality, disease and other factors.

Caut et al. (2007) studied a two species competition model which have twelve parameters but authors do not explain how they came up with most influence parameters in their model. Moreover, Blanco et al. (2014) studied an interference competition with an infectious disease, whereby this disease component was a reason of employing basic reproduction number for their model’s sensitivity analysis.

Therefore, having noted that, and with exception that our competition model does not have disease component, we studied a competition model by evaluating and analyzing the sensitivity indices of the basic reproduction number () in order to determine the importance of each parameter in the competition model. However, the intention is to simplify the work of numerical analysis by considering parameters that bring much effect in the model.

2. Model Formulation

The mathematical model for two competing species with human-related mortality to the species is formulated in this section.

2.1. Assumptions

In order to formulate the models, the following assumptions were made:

The lions and African wild dogs species live in an ecosystem where external factors such as droughts, epidemics are stable or have similar effect on the competing species;

There is no African wild dog related mortality to lions;

Lions kill African wild dogs to eliminate number of threats and competition;

African wild dogs and lions are both victims to humans.

2.2. Mathematical Model

A nonlinear mathematical model is formulated and analyzed.

(1)

(2)

where  and are all positive constants and and  measure the competitive effect of  on  and  and  respectively and that.

This model is modified from the classical Lotka-Volterra competition model as provided by (Baigent, 2010) in his lecture notes on Introduction to Lotka-Volterra Dynamics.

Table 1. Parameters used in the Competition model (1) – (2).

 Parameter Description Value rate of mortality caused by human beings to African wild dogs 0.0005 rate of mortality caused by human beings to lions 0.0001 δ rate of mortality caused by lions to African wild dogs 0.00076 competition coefficient of lion species on wild dog species 1.00014 competition coefficient of wild dog species on lion species 0.039 K1 carrying capacity for African wild dogs 5800 K2 carrying capacity for lions 1020 population growth rate of African wild dogs 0.05 Population growth rate of lions 0.025

Table 2. Variables used in the Competition model and their initial values.

 Variables Description Initial Value N1 Population of African wild dogs 580 N2 Population of lions 48

On the other hand, we have chosen classical Lotka-Volterra competition model for the reason that, it is a standard model for interference and exploitative categories of competition phenomenon. Besides, to the best knowledge of the authors, the Lotka-Volterra competition model has not yet been employed to study competition between African wild dog and lion species. Therefore, we have used that opportunity to use the model for these two competing species.

However, we have modified the classical Lotka-Volterra competition model by adding the following components: the rate of mortality caused by human to the competing species ( and) and the rate of mortality caused by lions to African wild dogs (). We have modified the model in this way for the reasons that, the wild dog species are endangered (IUCN Red List, 2008) and major threats to their existence are human and lion species (WWF, 2014; Perry, 2014).

3. Sensitivity Analysis

3.1. The Basic Reproduction Number, R0

For the sake of conducting sensitivity analysis we have borrowed the concept of basic reproduction number, denoted by R0 from epidemiological modeling but with the same concept that, R0 should be less than unity so that the population under consideration can withstand.

The basic reproduction number R0 of the model (1) and (2) is calculated by using the next generation matrix of an ODE (van den Driessche and Watmough, 2002). The R0 is obtained by taking the dominant eigenvalue of , whereby;

and  for

Then from model (1) and (2) when  we obtain:

Therefore if we let wild dogs’ killings and competition from lion species and the rest of the components in the model (1) and (2).

Then, it implies that:

and

Hence by using linearization method, the associated matrices are given by

We get

and

with

Therefore:

The eigenvalues of the product of  are:

and

It follows that the reproductive number which is given by the dominant eigenvalue for model system (1) and (2) is:

If, it means interaction between wild dogs and lions is minimal, hence, the rate of competition is very small and lion–related mortality to wild dogs is less or there is no such an incident and the rate of human–related mortality to the competing species is very small. Therefore, under these circumstances, the wild dog population will prevail from being endangered and start to increase in number. Whereas, if, then, the status of wild dog species will change to being extinct.

3.2. Sensitivity Analysis of Model Parameters

Sensitivity analysis state how important each parameter is to the competition model. Besides, sensitivity analysis is normally used to determine the strength of model predictions to parameter values, since there can be errors in data collection and postulated parameter values. Moreover, it is used to determine parameters that have greater influence on R0 and the one that has less influence. Therefore, we need to calculate the sensitivity indices of the reproduction number R0 to each parameter in the model. The parameters and their respective indices are ordered from most sensitive to the least.

Besides, a highly sensitive parameter should be carefully estimated, because a small variation in that parameter will lead to large quantitative changes. An insensitive parameter, on the other hand, does not require as much effort to estimate, since a small variation in that parameter will not produce large changes to the quantity of interest (Mikucki, 2012).

Table 3. Sensitivity indices of  calculated at the baseline parameter values given in Table 1.

 S/N Parameter Sensitivity index 1 2 3 4 5 6 7 8 9

3.3. Definition

The normalized forward sensitivity index of  that depends differentially on a parameter  is defined as (Chitnis et al., 2008:

By using the explicit formula for, one can easily derive an analytical expression for the sensitivity of  with respect to each parameter that comprise it. The obtained values are described in Table 3, which presents the sensitivity indices for the baseline parameter values in Table 1.

3.4. Interpretation of Sensitivity Indices

Commencing on Table 3, while being arranged in an order of most to least sensitive parameters, it shows that, if there is an individual increase amongst the parameters  and   while keeping other parameters constant, they tend to increase the value of  indicating that they increase pressure on already declined wild dogs’ population as they have positive indices. However, for the parameters  and  when each increases while keeping other parameters constant, they decrease the value of  suggesting that they decrease pressure on the decline of wild dog population as they have negative indices.

3.5. Numerical Simulation Analysis

Here, we illustrate the analytical results of the study by carrying out numerical simulations of the model system (1) and (2) using the set of estimated parameter values given as shown below.

The final time was years, and computations were run in Matlab with the  routine. This function implements a Runge–Kutta method with a variable time step for efficient computation. Figure 1 shows the solution to the system (1) and (2) with the baseline parameter values given in Table 1 for wild dogs and lions respectively.

Figure 1. Show variables of the ODE system (1) – (2) with baseline parameter values and initial variable values as in Table 1 and Table 2.

Figure 2. Changes in the solution of the system (1) and (2) on increasing and decreasing the parameter value by 10 percent.

Figure 2 shows a graph that reflects the effects on the population of wild dogs when delta,  was increased and decreased by 10 percent. The effect of increase in delta,  is shown by cyan dotted line which indicates that wild dog population will go to extinct in less than 80 years. Whereas, the effect of decrease in delta is shown by green dotted line, implying that wild dog population will die out in more than 80 years but less than 100 years. On contrary, when the parameter delta,  is decreased by 95 percent, wild dog population will grow and it will take more than 100 years to reach its carrying capacity,  this is indicated by Figure 3. The obtained graphics support the sensitivity analysis prepared in section 3.2. The most positive sensitive parameter is the rate of wild dogs’ mortality caused by lions, where.

Figure 3. Changes in the solution of the system (1) – (2) on decreasing the parameter value by 95 percent.

Figure 4. Changes in the solution of the system (1) – (2) on increasing the parameter value by 96 percent.

Figure 5. Changes in the solution of the system (1) – (2) on increasing the parameter value by 96 percent.

Figure 6. Changes in the solution of the system (1) – (2) on the effect of increasing the parameter value by 50 percent.

Figure 7. Changes in the solution of the system (1) – (2) on the effect of increasing parameter values by 10 percent.

On the other hand, the rest of positive sensitivity indices of the parameters and are lesser compared to the sensitivity index of the parameter delta,  by more than 96 percent, therefore when each of the parameter value was increased by 96 percent the changes were difficult to observe. This is shown by Figure 4 and Figure 5, as examples for the rest of the parameters in the group of least positive sensitivity indices.

The parameters  and have a negative sensitivity index. The most negative sensitive parameter is the intrinsic growth rate of wild dog population,  with. If  is increased by 50 percent, then the basic reproduction number  decreases approximately by 33 percent respectively. In this situation the wild dog population increase accordingly and reach its peak in not more than 25 years and start to fall since parameters that hinder the growth of wild dog population are held constant as can be seen in Figure 6.

Figure 7 presents the comparison of wild dog population when the baseline parameters are considered and all the parameters are increased by percent, whereby, the population of wild dog seemed to die out in less than  years, this is shown by green dotted line, while the population of lion continued to raise as indicated by purple dotted line.

4. Discussion and Conclusion

4.1. Discussion

A nonlinear mathematical model has been analyzed to study the effects of lions’ interference competition on population dynamics of wild dogs. Nevertheless, with the consideration that our research data were estimated, we had to borrow the concept of basic reproduction number , which show that the population of wild dog will survive much longer, that is 80+ years when , while, when , the wild dog population will die out in less than 20 years.

Sensitivity analysis shows that, the rate of wild dog mortality caused by lion,  is the most sensitive parameter on  whereby when  is increased by 50 percent  increases by 49.44 percent appropriately, and the least sensitive is intrinsic growth rate of wild dog population,  which when increased by 50 percent  decreases approximately by 33 percent respectively.

In numerical simulation it is observed that the increase of rate of wild dog mortality caused by lion, tends to decrease the population of wild dogs. But the decrease or absence of wild dog mortality raises the population of wild dog towards its carrying capacity. On the other hand, the indices of least positive sensitive parameter,  and, and least negative sensitive parameters,  and  have small influence on  and their changes are not graphically definite. However, the increase of intrinsic growth rate,  while holding other parameters constant, tends to increase the population of wild dog for some years.

This indicate that much work should be on controlling or rather eradicating factors that hinder the growth of wild dog population.

4.2. Conclusion

A competition model for wild dog population dynamics towards interference competition by lions and the effects of human related–mortality was presented, with lion and wild dog populations being variables under consideration. Simulation shows that when the rate of wild dog mortality caused by lion is reduced or eliminated then there will be growth of wild dog population. Moreover, when the rest of parameters are controlled while increasing carrying capacity of wild dog, the intrinsic growth rate of wild dog will be favoured and the population will withstand. This is only possible when protected areas for wild dog population are limited with lion population because the two species, lion and wild dog cannot share the same territory as most of the time the lion species would eliminate or threaten its competitor in resource utilization. On the other hand, human activities that jeopardize lives of wild dogs should be prohibited, that is local people should be stopped from hunting and grazing in the protected areas as well as sport hunting tourism should be abolished and much effort should be on practising sustainable tourism.

References

1. Baigent, S., 2010. Lotka-Volterra Dynamics- An Introduction. [Online]Available at:http://www.ltcc.ac.uk/courses/BioMathematics/LTCC LV2010[Accessed 8 4 2015].
2. Blanco, K., Barley, K., and Mubayi, A., 2014. Population Dynamics of Wolves and Coyotes at Yellowstone National Park: Modeling Interference Competition with an Infectious Disease. pp. 2-5.
3. Caut, S., Casanovas, J. G., Virgos, E., Lozano, J., Witmer, G. W., and Courchamp, F., 2007. Rats dying for mice: Modelling the competitor Release Effect. Austral Ecology, Issue 32, pp. 858 - 868.
4. Chitnis, N., Hyman, J. and Cusching, J. M., 2008. Determining important parameters in the spread of malaria through the sensitivity analysis of a mathematical model. Bull. Math. Biol., 5(70), pp. 1272-1296.
5. Cronin, J. P. and Carson, W. P., 2015. Biology Encyclopedia Forum. [Online] Available at:http://www.biologyreference.com/Ce-Co/Competition.html[Accessed 18 August 2015].
6. Dutson, G. and Sillero-Zubiri, C., 2005. Forest-dwelling African wild dogs in the Bale Mountains, Ethiopia. [Online] Available at:http://www.canids.org/canisnews/8/African_wild_dogs_in_Ethiopia.pdf[Accessed 31 May 2015].
7. IUCN Red List, 2008. Endangered Species. [Online] Available at:http://www.iucnredlist.org/[Accessed 27 May 2015].
8. Lozano, J., Casanovas, J. G., Virgos, E. and Zorrila, J., 2013. The competitor Release Effect applied to carnivore species: How Red Foxes can Increase in numbers when persecuted. Animal Biodiversity and Conservation, Issue 36.1, p. 39.
9. Mikucki, M. A., 2012. Sensitivity analysis of the basic reproduction number and other quantities for infectious disease models, Colorado: Colorado State University.
10. Nakayama, S. and Fuiman, L. A., 2010. Body size and vigilance mediate asymmetric interference competition for food in fish larvae, Oxford: Oxford University Press.
11. Painted Dog Conservation Incorporated, 2014. Threats. [Online] Available at:http://www.painteddogconservation.iinet.net.au[Accessed 28 May 2015].
12. Sillero-Zubiri, C., Hoffmann, M. and Macdonald, D. W., 2004. Canids: Foxes, Wolves, Jackals and Dogs: Status and Conservation Action Plan, Gland, Switzerland and Cambridge, UK: IUCN/SSC Canid Specialist Group.
13. Van den Driessche, P. and Watmough, J., 2002. Reproduction Numbers and Sub-Threshold Endemic Equilibria for Compartmental Models of Disease Transmission. Mathematical Biosciences, Issue 180, pp. 29 - 48.
14. Woodroffe, R., Ginsberg, J., Macdonald, D. and IUCN/SSC Canid Specialists, 1997. The African Wild Dog - Status Survey and Conservation Action Plan. Gland, Switzerland: IUCN.

 Contents 1. 2. 2.1. 2.2. 3. 3.1. 3.2. 3.3. 3.4. 3.5. 4. 4.1. 4.2.
Article Tools