I performed cluster analysis in R (using hclust under vegan package) to find relationships between habitats based on species. A dendrogram/tree was derived from Bray-Curtis sim. matrix. Cluster methods are Ward, Ward.D2, Single, Complete, Average etc.
However, when I perform an ANOVA with post-test, the significant differences between pairs of habitats do not represent the differences observed in the dendrogram, when using a variety of methods. E.g. ANOVA computed that A was different from C, but the cluster analysis showed that A is similar to C.
ANOVA was performed using species abundances and cluster analysis using Bray-Curtis which converts to binary/presence-absence.
Any ideas?