Calculates the low-level diversity component necessary for calculating gamma diversity.

raw_gamma(meta)

Arguments

meta

object of class metacommunity

Value

raw_gamma returns an object of class powermean

Details

Values generated from raw_gamma() may be input into subdiv() and metadiv() to calculate subcommunity and metacommunity gamma diversity.

References

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.

Examples

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