Calculates the low-level diversity component necessary for calculating gamma diversity.
raw_gamma(meta)object of class metacommunity
raw_gamma returns an object of class powermean
Values generated from raw_gamma() may be input into subdiv() and
metadiv() to calculate subcommunity and metacommunity gamma diversity.
R. Reeve, T. Leinster, C. Cobbold, J. Thompson, N. Brummitt, S. Mitchell, and L. Matthews. 2016. How to partition diversity. arXiv 1404.6520v3:1–9.
pop <- cbind.data.frame(A = c(1,1), B = c(2,0), C = c(3,1))
row.names(pop) <- paste0("sp", 1:2)
pop <- pop/sum(pop)
meta <- metacommunity(pop)
# Calculate gamma component
g <- raw_gamma(meta)
subdiv(g, 1)
#> measure q type_level type_name partition_level partition_name diversity
#> 1 gamma 1 types subcommunity A 2.309401
#> 2 gamma 1 types subcommunity B 1.333333
#> 3 gamma 1 types subcommunity C 1.754765
#> dat_id transformation normalised k max_d
#> 1 naive NA NA NA NA
#> 2 naive NA NA NA NA
#> 3 naive NA NA NA NA
metadiv(g, 1)
#> measure q type_level type_name partition_level partition_name diversity
#> 1 gamma 1 types metacommunity 1.754765
#> dat_id transformation normalised k max_d
#> 1 naive NA NA NA NA