Estimation function for Shannon's index. Internal use in
sbdiv for methods rpht, tsht,
asht. Sums up species counts in each columns for every
treatment group and estimates Shannon's index with bias correction on
the resulting vectors of summed up species counts.
$$\widehat{HBC}_{i} = \hat{H}_{i} + (S_i -1)/(2N_{i\bullet}) -
(1-\sum(1/\hat{p}_{i\bullet s}))/(12N_{i\bullet}^2) -
\sum((1/\hat{p}_{i\bullet s})-(1/(\hat{p}_{i\bullet
s}^2)))/(12N_{i\bullet}^3);$$
\(i=1,...,k;s=1,...,S;p_{i \bullet s}=\frac{\sum_{j=1}^{n}x_{sj}}{N_{i\bullet}}\);
$$\hat{H}_i=(-1)\sum_{s=1}^{S}(\hat{p}_{i \bullet s} log(\hat{p}_{i \bullet s}))$$
\(N_{i\bullet}= \sum_{j=1}^{n}N_{ij}\) Number of observed individuals in treatment \(i\).
estShannonf(X, f)
Estimated Shannon-Wiener index for treatment groups
Estimated variance of Shannon-Wiener index for treatment groups
\(n\) times \(p\) matrix containing species in \(p\) columns and replicates in \(n\) rows.
Factor variable containing treatment groups. Must be of length: replicates times treatment groups.