Shannon-Weanerの多様度指数H'を返すR関数

ShannonH <- function(v, base="e") {
	v <- subset(v, v>0)
	if (base=="e") base <- exp(1)
	s <- sum(v)
	h <- -sum(v/s * logb(v/s, base))
	return(h)
}

base引数で、対数の底を指定できるようにした。