I wrote a function in R to do hierarchical clustering.
Klasteryzacja_hierarchiczna = function(ExprSet,a,b,c,d) {
Klasteryzacja_h =exprs(ExprSet)'
clusters = hclust(dist(Klasteryzacja_h[a:b,c:d]))'
plot(clusters,xlab="Clusters")'
}
but it gave me this:
You need to set the sub
parameter.
plot(clusters, xlab="Clusters", sub="")