Linear program for ecology

Recently worked with Dale Nimmo and our honours student Andrew on optimizing abundance (or the geometric mean of abundance) based on land types.  I found it a really interesting problem, as trying to using a general solver generally just didn’t work.  The optimization is essentially with respect to the objective:

∏_i=1…n  (∑w_j x_ij)

and so non-linear.  However by representing the geometric mean as the sum of logs and, in turn, the logarithmic function as the max min of linear functions then we’re able to solve it linearly.  I also looked at a bilevel approach for maximizing Shannon’s diversity.

The code is available here:


Leave a Reply

Fill in your details below or click an icon to log in: Logo

You are commenting using your account. Log Out /  Change )

Google photo

You are commenting using your Google account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )

Connecting to %s