I decided to remake my dendrograms using a year old code which worked just fine before:
fit <- pvclust(t(matrix), method.hclust = "ward.D2", meth