I have a dataset of flower shape characteristics (zygomorphic or actinomorphic, tubular or non-tubular) for 2 forests in 2 different regions. I want to see if the flower shape characteristics are significantly different in the 2 regions. Is using the anosim (Jaccard) the right test? The line of code below

anosim(flower_characteristics.matrix, Flower$region, permutations = 999, distance = "jaccard", strata = NULL,

parallel = getOption("mc.cores"))