Description Usage Arguments Details Value Author(s) References See Also Examples
Indicator power calculation of Halme et al. (2009) or the congruence between indicator and target species.
1 |
x |
Community data frame or matrix. |
type |
The type of statistic to be returned. See Details for explanation. |
Halme et al. (2009) described an index of indicator power defined as
IP_I = sqrt(a*b), where a = S / O_I and
b = 1 - (O_T - S) / (N - O_I).
N is the number of sites,
S is the number of shared occurrences of the indicator (I)
and the target (T) species. O_I and O_T are number
of occurrences of the indicator and target species. The type
argument in the function call enables to choose which statistic to
return. type = 0
returns IP_I, type = 1
returns
a, type = 2
returns b.
Total indicator power (TIP) of an indicator species is the column mean
(without its own value, see examples).
Halme et al. (2009) explain how to calculate confidence
intervals for these statistics, see Examples.
A matrix with indicator species as rows and target species as columns (this is indicated by the first letters of the row/column names).
Peter Solymos
Halme, P., Mönkkönen, M., Kotiaho, J. S, Ylisirniö, A-L. 2009. Quantifying the indicator power of an indicator species. Conservation Biology 23: 1008–1016.
indval
(package labdsv) for the
indicator species analysis of Dufrêne & Legendre.
Function beals
estimates individual cell probabilities
of species occurrences.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 | data(dune)
## IP values
ip <- indpower(dune)
## and TIP values
diag(ip) <- NA
(TIP <- rowMeans(ip, na.rm=TRUE))
## p value calculation for a species
## from Halme et al. 2009
## i is ID for the species
i <- 1
fun <- function(x, i) indpower(x)[i,-i]
## 'c0' randomizes species occurrences
os <- oecosimu(dune, fun, "c0", i=i, nsimul=99)
## get z values from oecosimu output
z <- os$oecosimu$z
## p-value
(p <- sum(z) / sqrt(length(z)))
## 'heterogeneity' measure
(chi2 <- sum((z - mean(z))^2))
pchisq(chi2, df=length(z)-1)
## Halme et al.'s suggested output
out <- c(TIP=TIP[i],
significance=p,
heterogeneity=chi2,
minIP=min(fun(dune, i=i)),
varIP=sd(fun(dune, i=i)^2))
out
|
Loading required package: permute
Loading required package: lattice
This is vegan 2.4-3
i.Achimill i.Agrostol i.Airaprae i.Alopgeni i.Anthodor i.Bellpere i.Bromhord
0.3186250 0.3342800 0.2168133 0.3416198 0.3567884 0.3432281 0.3665632
i.Chenalbu i.Cirsarve i.Comapalu i.Eleopalu i.Elymrepe i.Empenigr i.Hyporadi
0.2095044 0.2781640 0.1713273 0.2414787 0.3263516 0.2016196 0.2378197
i.Juncarti i.Juncbufo i.Lolipere i.Planlanc i.Poaprat i.Poatriv i.Ranuflam
0.2915850 0.3331330 0.3998442 0.3426064 0.4094319 0.3929520 0.2663080
i.Rumeacet i.Sagiproc i.Salirepe i.Scorautu i.Trifprat i.Trifrepe i.Vicilath
0.3484684 0.3788905 0.2898512 0.4362493 0.3145854 0.4503764 0.2605349
i.Bracruta i.Callcusp
0.4252676 0.2070766
Warning message:
In oecosimu(dune, fun, "c0", i = i, nsimul = 99) :
nullmodel transformed 'comm' to binary data
[1] -1.833791
[1] 86.72236
[1] 0.9999999
TIP.i.Achimill significance heterogeneity minIP varIP
0.3186250 -1.8337915 86.7223575 0.0000000 0.2142097
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.