# ----------------------
# 1. tree level mutation

# 1.1 diameter classes
# returns the diameter class
# default to 5 cm classes
# include second argument for e.g. 10 cm classes

diameter.class <- function(x, multiple = 5) {
  multiple * floor(x / multiple + 0.5)
}

# 1.2 diameter increment

# ----------------------
# 2. summarise into plot

# 2.1 q value

# 2.2 Gini coefficient

# 2.3 Shannon Diversity Index

# 2.4 Stand Diversity Index

# 2.5 Transition period

# 2.6 Annual recruitment

# 2.7 Topheight

# ----------------------
# 3. plot level mutation

# 3.1 form class

# 3.2 form height