我是R的新手,我正在尝试在数据表上进行一些聚类,其中行表示单个对象,列表示已为这些对象测量的特征 . 我已经完成了一些聚类教程,但我得到了一些输出,但是,我在聚类后获得的热图与根据同一数据表与另一个程序生成的热图完全不对应 . 虽然该程序的热图确实表明了对象之间标记表达的明显差异,但我的热图并没有显示出太大的差异,我无法识别热图上的任何聚类(即颜色)图案,它似乎只是一个随机混乱的集合颜色彼此接近(没有大的对比度) . 这是我正在使用的代码的一个例子,也许有人知道我可能做错了什么 .
mydata
datamat
datalog
我正在使用聚类的日志值,因为我知道其他程序也这样做
library(gplots)
hr
mycl
mycol
heatmap(datamat, Rowv=as.dendrogram(hr), Colv=NA,
col=colorpanel(40, "black","yellow","green"),
scale="column", RowSideColors=mycol)
再次,我绘制原始颜色,但使用日志集群,因为我知道这是其他程序的作用 .
我尝试使用这些方法,但我没有得到任何至少在某种程度上看起来像聚集热图的东西 . 当我取出缩放时,热图变得非常暗(我实际上确信我已经以某种方式按列扩展或规范化数据) . 我也尝试用k-means进行聚类,但同样,这没有帮助 . 我的想法是,由于两个异常值,颜色标度可能无法完全使用,但是尽管去除它们会略微增加热图上绘制的颜色范围,但这仍然没有显示出正确的聚类 .
还有什么我可以玩的吗?
是否可以使用热图更改色阶,以便在最后一个具有“一切都大于特定值”范围的仓中找到异常值?我尝试用heatmap.2(参数“break”)来做这个,但是我没有成功,而且我没有设法把我使用的行侧颜色放在热图功能上 .