我的热图功能和我自己的群集功能有问题......
我试着这样说但是它不起作用:
heatmap(data, hclustfun=function(x) myclust(x))
其中myclust是一个简单的函数
myclust= function(data) {
#D=dist(data)
hc=hclust(D)
#mclust= as.mclust(hc$merge)
#m=reorder.mclust(data)
mc=hc
mc$order=leaf
return(mc=mc)
}
它在heatmap(data, hclustfun = function(x) myclust(x))
中返回错误:
列树形图排序给出了错误长度的索引
实际上它正在使用其他群集功能,例如diana()
或pam()
,但不能使用我自己的功能......