[英]How to change color palette of mosaic plot
我有以下数据和代码:
> mat <- matrix(c(120,230,84,70,130,83,13,26,18),3)
> dimnames(mat) <- list(c("good","fair","poor"),c("a","b","c"))
> mat
> mat
a b c
good 120 70 13
fair 230 130 26
poor 84 83 18
>
> chisq.test(mat)
Pearson's Chi-squared test
data: mat
X-squared = 11.411, df = 4, p-value = 0.02231
> library(vcd)
> mosaic(mat, shade=T)
卡方检验显示出显着差异,但由于残差很小,因此看不到颜色。 残差较高时,可以看到颜色:
?mosaic 显示此信息:
highlighting_fill: color vector or palette function used for a highlighted variable, if any.
如何更改上述绘图的填充颜色?
没有真正优雅的解决方案,但您可以使用 gpar 手动定义颜色
mosaic(mat, shade=T, colorize = T, gp = gpar(fill=matrix(c("red","blue","light green", "green", "black","red","brown","white","blue"), 3, 3)))
下一步是用合理的颜色填充 gpar 中的矩阵。 这里有一篇关于如何根据值分配色带的很好的帖子。 https://stackoverflow.com/a/6320940/1502898
默认情况下, vcd
中马赛克图的阴影使用截止值 2 和 4 作为与每个单元格关联的残差的绝对值,如 Friendly (1994, Journal of the American Statistical Association , 89 , 190-200) 所建议的。 这背后的想法是,如果模型拟合良好(由于残差的近似正态性),这些截止值应该与大约 5% 和 0.01% 交叉。 然而,这些临界值对应于细胞的个体显着性,而卡方检验评估整体显着性。
有两种策略可以解决这个问题:(1) 由于 Friendly (1994) shading 旨在显示偏差模式而不是直接对应显着性,因此可以简单地将截止值修改为较低的值对于给定的表(例如,此处为 1 和 1.8)。 (2) 切换到最大绝对残差(而不是它们的平方和)的测试,以便在某些显着性水平(例如,10% 和 1%)上显着的所有单元格都可以被着色。 这两种策略都在vcd
中实现,并在?shading_hcl
的示例中进行了说明(使用与mat
属性非常相似的示例)。 Zeileis, Meyer, Hornik (2007) 解释了基本思想。 “用于可视化(条件)独立性的基于残差的阴影”,计算和图形统计杂志, 16 (3),507-525。 doi:10.1198/106186007X237856
策略 (1) 可以通过以下方式应用于您的数据:
mosaic(mat, gp = shading_hcl, gp_args = list(interpolate = c(1, 1.8)))
策略 (2) 实现为:
set.seed(403)
mosaic(mat, gp = shading_max)
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.