R中有没有一种方法可以计算用Bray-Curtis相异度矩阵生成的分层聚类的差距统计?

我已经试过了 clusGap 和 fviz_gap_stat ;两者都需要调用 hcut 然而,该函数不允许使用Bray-Curtis作为度量来计算观测值之间的相异性。因此,我试着 vegdist 并在命令中输入Bray-Curtis矩阵,但没有成功。

用于估计簇数的间隙统计

gap_stat <- clusGap(otu_matrix,FUN=hcut,hc_func="hclust",hc_method="ward.D",isdiss=TRUE,Braymatrix,K.max = 50, B = 500)

Clustering k = 1,2,..., K.max (= 50): .. 

if(is.na(n)||n>65536L)停止(“大小不能为NA,也不能超过65536”): 缺少值,需要TRUE/FALSE`