Shannon diversity index formula

P

peter

I'm trying to create a single formula that evaluates to the Shannon
diversity index (SDI) on a single array. An explanation of SDI can be
found on http://en.wikipedia.org/wiki/Diversity_index.

SDI = -E { p(i)*log(p(i)) }

In words, that's the negative sum over the (proportion times the log
of the proportion)

It's a simple thing to do if all the values in the array are > 0, but
if one is zero, that screws up the LOG function, because LOG(0) is
undefined and so it returns an error.

This is my formula at present: =-SUMPRODUCT((B2:B20/SUM(B$2:B
$20)),LOG(B2:B20/SUM(B$2:B$20)))
And below is the sample data...

A B
1 Species TOTAL
2 Species 1 2
3 Species 2 2
4 Species 3 2
5 Species 4 2
6 Species 5 2
7 Species 6 2
8 Species 7 2
9 Species 8 2
10 Species 9 2
11 Species 10 2
12 Species 11 2
13 Species 12 2
14 Species 13 2
15 Species 14 2
16 Species 15 2
17 Species 16 2
18 Species 17 2
19 Species 18 2
20 Species 19 2

The TOTAL column is the sum of all the following columns. However,
there will always be a number of zeros in the TOTAL column, so the
function returns an error. I've tried using conditional tests as
explained on http://www.xldynamic.com/source/xld.SUMPRODUCT.html, but
I can't figure out how to get those to work inside the LOG function
which is the source of the error.

The SDI is also easy to calculate if I have a couple of columns in
which to put the proportions and the LOG(proportions), but that's not
what I'm looking for.

I have seen a few functions made by various folks around on the
internet, but 1) some have this same problem, and 2) some require an
add-in, neither of which are acceptable.

What say you all?
Peter
 
B

Barb Reinhardt

Try this

=-SUMPRODUCT(--(C2:C20>0),(C2:C20/SUM(C2:C20)),(LOG(C2:C20/SUM(C2:C20))))
 
G

Gary''s Student

p should never be zero. Because pi is the fraction of individuals belonging
to the i-th species, by definition it can't be zero. Including a zero is
like including the number of whales in the Amazon river basin - there aren't
any. Just exclude zeros, polar bears, kangaroos, etc.
 
P

peter

Well, when you've got a datasheet with a list of 400 species on it,
and the polar bear doesn't happen to show up for the survey that day,
then the polar bear gets a zero. So what you're saying is what I'm
trying to do... exclude them without having to do lots of
manipulations.

And I'm quite happy to say that there are, in fact, whales in the
Amazon River basin. The Amazon River dolphin: http://www.acsonline.org/factpack/Boto.htm.
You're correct that the polar bears and kangaroos are elsewhere.

However, I'm not asking a theoretical question here, I'm asking for a
formula. So, back on track... anyone with anything helpful?
 
G

Gary''s Student

Sorry about the whales thing... I have nothing against whales, some of my
best friend are cetaceans.


In any event in C2 thru C20 enter:

=-(B2/SUM($B$2:$B$20))*IF(B2=0,0,LOG(B2/SUM(B$2:B$20)))
and in another cell:

=sum(C2:C20)
 
P

peter

Thanks, but I'm really trying to consolidate the whole thing into a
single formula in a single cell.
 
B

Barb Reinhardt

Try this then:

=-SUMPRODUCT((C$2:C$20/SUM(C$2:C$20)),(IF(C$2:C$20>0,(LOG(C$2:C$20/SUM(C$2:C$20))))))
Activate with CTRL SHIFT ENTER
You should see {} around the formula when you are done.
 
B

Barb Reinhardt

Try this then:

=-SUMPRODUCT((C$2:C$20/SUM(C$2:C$20)),(IF(C$2:C$20>0,(LOG(C$2:C$20/SUM(C$2:C$20))))))
Activate with CTRL SHIFT ENTER
You should see {} around the formula when you are done.
 
H

Harlan Grove

peter said:
This function still evaluates to a #NUM error.

Barb's formula doesn't work because --(C2:C20>0) doesn't exclude
calculating LOG(C2:C20/SUM(C2:C20)) for zero values in C2:C20. There's
a fix for that. I've also rearranged terms to take advantage of the
distributive law.

=-SUMPRODUCT(C2:C20,LOG(C2:C20+(C2:C20=0)))/SUM(C2:C20)
+LOG(SUM(C2:C20))
 

Ask a Question

Want to reply to this thread or ask your own question?

You'll need to choose a username for the site, which only take a couple of moments. After that, you can post your question and our members will help you out.

Ask a Question

Top