Abstract
We consider an optimization problem in ecology where our objective is to maximize biodiversity with respect to different land-use allocations. As it turns out, the main problem can be framed as learning the weights of a weighted arithmetic mean where the objective is the geometric mean of its outputs. We propose methods for approximating solutions to this and similar problems, which are non-linear by nature, using linear and bilevel techniques.
Access this chapter
Tax calculation will be finalised at checkout
Purchases are for personal use only
Notes
- 1.
- 2.
However since the constraints are linear, quadratic programming formulations would also be fine.
- 3.
As a side note the common starling is an introduced species and was not included in the previous example of 21 native bird species.
- 4.
References
Tuomisto, H.: An updated consumer’s guide to evenness and related indices. Oikos 121, 1203–1218 (2012)
Kelly, L.T., Bennett, A.F., Clarke, M.F., McCarthy, M.A.: Optimal fire histories for biodiversity conservation. Conserv. Biol. 29, 473–481 (2015)
Beliakov, G., Bustince, H., Calvo, T.: A Practical Guide to Averaging Functions. Springer, New York (2015)
Beliakov, G., Pradera, A., Calvo, T.: Aggregation Functions: A Guide for Practitioners. Springer, Heidelberg (2007)
Grabisch, M., Marichal, J.-L., Mesiar, R., Pap, E.: Aggregation Functions. Cambridge University Press, Cambridge (2009)
Torra, Y., Narukawa, V.: Modeling Decisions. Information Fusion and Aggregation Operators. Springer, Heidelberg (2007)
Grimm, N.B., Faeth, S.H., Golubiewski, N.E., Redman, C.L., Wu, J., Bai, X., Briggs, J.M.: Global change and the ecology of cities. Science 319, 756–760 (2008)
Lin, B.B., Fuller, R.A.: FORUM: sharing or sparing? how should we grow the world’s cities? J. Appl. Ecol. 50, 1161–1168 (2013)
Beliakov, G., James, S.: Citation-based journal ranks: the use of fuzzy measures. Fuzzy Sets Syst. 167, 101–119 (2011)
Beliakov, G., James, S.: Using linear programming for weights identification of generalized Bonferroni means in R. In: Torra, V., Narukawa, Y., López, B., Villaret, M. (eds.) MDAI 2012. LNCS, vol. 7647, pp. 35–44. Springer, Heidelberg (2012)
Beliakov, G., James, S., Gómez, D., RodrÃguez, J.T., Montero, J.: Learning stable weights for data of varying dimension. In: Proceedings of the 8th International Summer School on Aggregation Operators, University of Silesia, Katowice, Poland (2015)
Beliakov, G., James, S., Gómez, D., RodrÃguez, J.T., Montero, J.: Approaches to learning strictly-stable weights for data with missing values. Fuzzy Sets and Systems (2016, submitted)
Beliakov, G., James, S., Nimmo, D.: Learning aggregation weights from 3-tuple comparison sets. In: Proceedings of IFSA 2014, Edmonton, Canada. pp. 1–6 (2013)
Beliakov, G., James, S., Nimmo, D.: Using aggregation functions to model human judgements of species diversity. Inf. Sci. 306, 21–33 (2015)
Author information
Authors and Affiliations
Corresponding author
Editor information
Editors and Affiliations
Appendix: Implementations
Appendix: Implementations
We have implemented all three approaches to optimization as functions in an R library available at our websiteFootnote 4.
Optimizing total abundance: eco.opti().
Description of inputs
-
species.data - matrix of species abundances per unit of land area, i.e. with \(x_{ij}\) denoting the i-th species and its abundance for land type j;
-
densities - vector of densities per unit of land area, i.e. with \(d_j\) denoting the human population density for a given land-type j;
-
tot.pop - the total population required to be fit into the given land area;
-
tot.land - the total area over which we need to distribute the population.
Additional optional constraints
-
w.min / w.max - vectors denoting minimum or maximum bounds on the land-types, e.g. if we want to ensure that at least 20Â % of the land is populated at minimum density we incorporate the constraint \(w_1 \ge 0.2\) (assuming \(w_1\) is the land-type with minimum density), or alternatively we may wish to limit high density housing to at most 40Â % of the land area etc.;
-
spec.min / spec.max - vector placing minimum or maximum bounds for a particular species, for example, if we want to make sure that a rare species is above a given threshold \(\gamma _i\), the linear constraint \(\sum \limits _{j=1}^n w_j x_{ij} \ge \gamma _i\) is added for that species.
The function also gives as output a number of ecological indices such as the individual species abundances, and the Simpson and Shannon diversity indices.
Maximizing the geometric mean of species abundances: eco.opti.gm()
Description of inputs
In addition to all inputs and constraints used with eco.opti(), this function has two additional optional parameters to control the precision.
-
fprec - a positive integer giving the number of tangent functions to be defined. The default setting is 100 linear segments, and so gains in accuracy can be achieved with settings of 500, 1000 etc., however obviously at the cost of computation time;
-
max.x - a real number giving the maximum value for the domain over which the tangent functions are calculated, the default setting is 10000, and so depending on the scale given it could be necessary to increase this value (or decrease it for finer accuracy) or the optimization will be the same as it would be for maximizing abundance.
For the number of tangent lines K, optimizing over 5 species with 100 linear segments will require \(5\times 100 = 500\) additional constraints, use of 1000 linear segments will require 5000 additional constraints and so on. We need to be careful when reducing the precision, since the log function’s gradient changes more drastically for values closer to zero than it does for large values. Rather than taking equal step sizes in calculating our tangent lines, they were distributed using \(t_k = \exp (-k\cdot \max (s)/K)\) where \(\max (s)\) is the max.x parameter above.
We note also that by default the setting for spec.min will be 1 for all species. It could be adjusted to a fractional value if desired.
Maximizing Shannon’s diversity index: eco.opti.sh()
This function uses the same inputs as with the previous two. The program first solves a maximum and minimum problem using eco.opti() in order to find the feasible bounds to search for M. Note that \(-t \ln t\) is concave for \(t \in (0,1]\).
Another biodiversity index used as an objective and included in the code made available online is Simpson’s diversity index \(1/(\sum _{i=1}^m p_i^2)\). This is performed in a similar manner, however now we are minimizing for a convex function \(y = t^2\) rather than maximizing for a concave function and so we need to make the appropriate changes when using the linear framework above.
Rights and permissions
Copyright information
© 2016 Springer International Publishing Switzerland
About this paper
Cite this paper
Beliakov, G., Geschke, A., James, S., Nimmo, D. (2016). Linear Optimization for Ecological Indices Based on Aggregation Functions. In: Carvalho, J., Lesot, MJ., Kaymak, U., Vieira, S., Bouchon-Meunier, B., Yager, R. (eds) Information Processing and Management of Uncertainty in Knowledge-Based Systems. IPMU 2016. Communications in Computer and Information Science, vol 611. Springer, Cham. https://doi.org/10.1007/978-3-319-40581-0_33
Download citation
DOI: https://doi.org/10.1007/978-3-319-40581-0_33
Published:
Publisher Name: Springer, Cham
Print ISBN: 978-3-319-40580-3
Online ISBN: 978-3-319-40581-0
eBook Packages: Computer ScienceComputer Science (R0)