The following line returns an observed shannon diversity of 12.06 for girdled treatment in the spider dataset : iNEXT(spider, q=0, datatype="abundance").
Which doesn't match the observed shannon diversity of 2.49 returned when running ChaoShannon(spider$Girdled, datatype="abundance")
Doing it by hand, the correct shannon diversity for girdled is 2.49. Where does the 12.06 value of iNEXT come from?
Hi there,
The following line returns an observed shannon diversity of 12.06 for girdled treatment in the spider dataset : iNEXT(spider, q=0, datatype="abundance").
Which doesn't match the observed shannon diversity of 2.49 returned when running ChaoShannon(spider$Girdled, datatype="abundance")
Doing it by hand, the correct shannon diversity for girdled is 2.49. Where does the 12.06 value of iNEXT come from?
p<-spider$Girdled/sum(spider$Girdled) -sum(p*log(p))
Thanks!