f.Hclust_col_pc <- function ( expr.mx, title.ch) { f.distCol.Pearson <- function (expr.mx) # correlation among columns { dist.mx <- 1 - cor (expr.mx, method = "pearson") dist.dmx <- as.dist (dist.mx) return (dist.dmx) } x.hclust <- hclust (d = f.distCol.Pearson (expr.mx)) x11 () plot ( x.hclust, main = title.ch, xlab = "1 - Pearson Correlation") }