##Supplementary file 15_Comparison of total sapling density (n) among treatment groups per species (Figure 4)

##load data from "sapling density per plot" file (Supplementary file 8)

##Non-parametric test
library(dunn.test)

##Quercusacuta
dunn.test(sapling_density_data_per_plot$Quercusacuta,sapling_density_data_per_plot$Treatment,method = "bonferroni")

##Abies firma
dunn.test(sapling_density_data_per_plot$Abiesfirma,sapling_density_data_per_plot$Treatment,method = "bonferroni")

##Cinnamomum tenuifolium
dunn.test(sapling_density_data_per_plot$Cinnamomumtenuifolium,sapling_density_data_per_plot$Treatment,method = "bonferroni")

##Neolitsea sericea
dunn.test(sapling_density_data_per_plot$Cinnamomumtenuifolium,sapling_density_data_per_plot$Treatment,method = "bonferroni")

##Castanopsis sieboldii
dunn.test(sapling_density_data_per_plot$Castanopsissieboldii,sapling_density_data_per_plot$Treatment,method = "bonferroni")

##Eurya japonica
dunn.test(sapling_density_data_per_plot$Euryajaponica,sapling_density_data_per_plot$Treatment,method = "bonferroni")