可视化R中posthoc Tukey的临界值/成对比较

Question

我试图从posthoc Tukey获得关键值的细粒度可视化.有一些很好的指导方针可用于可视化成对比较,但我需要更精确的东西.我的想法是,我会得到一个图,其中每个小方块代表下面矩阵的临界值,编码方式如下:

• 如果该值高于或等于5.45 - 它是一个黑色方块;
• 如果该值低于或等于-5.45 - 则为灰色方块;
• 如果该值介于-5.65和5.65之间 - 它是一个白色方块.

数据矩阵在这里.

或者你可能会更好地建议如何将这些关键值可视化？

编辑:根据@Aaron和@DWin的评论,我想为上述数据提供更多的背景,并为我的问题提供理由.我正在研究七个虚拟角色的可接受性的平均评级,每个虚拟角色在5个不同的级别上进行动画制作.所以,我有两个因素 - 角色(7个级别)和动作(5个级别).因为我发现了这两个因素之间的相互作用,所以我决定研究所有运动水平的所有角色的均值之间的差异,这导致了这个庞大的矩阵,作为posthoc Tukey的输出.现在可能太详细了,但请不要把我扔给Cross Validated,他们会活着吃我...

Answer 1

这非常简单image:

d <- as.matrix(read.table("http://dl.dropbox.com/u/2505196/postH.dat"))    
image(x=1:35, y=1:35, as.matrix(d), breaks=c(min(d), -5.45, 5.45, max(d)), 
      col=c("grey", "white", "black"))

只有一半,将一半设置为缺失d[upper.tri(d)] <- NA并添加na.rm=TRUE到 min和max功能.