Example

###################################################
#
data <- exprs (GCCL.37$eset)# eset of GCCL.37
#1. repeated centering
data.centering <- medpolish(data)$residual

#2. repeated standardization
std.data <- data.centering
for (i in 1:10)
{
#By row (gene)
std.data <- t(scale(t(std.data), scale=T))
#By column (array)
std.data <- scale(std.data, scale=T)
}

HC <- hclust(dist(t(std.data)))
plot (HC) # Hclust plot of 2 large natural clusters of 37 GCCL


cl.37 <- cutree(HC, k=2) # 2 large clusters 16 G1 and 21 G2
setwd("E:\\Projects\\8.ComBAT\\ComBat201\\G1G2")
write.table(cl.37, file="G1G2.cl.37.medpolish.txt", sep="\t")


cl.37 <- cutree(HC, k=3) # 2 large clusters 16 G1 and 21 G2
setwd("E:\\Projects\\8.ComBAT\\ComBat201\\G1G2")
write.table(cl.37, file="G1G2.cl.37.medpolish3.txt", sep="\t")