解释一下“par(mar=c(2,2,2,4)) circos.par(gap.degree=c(3,rep(2, nrow(cor1)-1)), start.degree = 180) chordDiagram(cor1, grid.col=rainbow(ncol(rt)), col=col1, transparency = 0.5, symmetric = T) par(xpd=T)”
时间: 2023-08-31 21:13:31 浏览: 294
这段代码主要用于绘制一个基于圆环的相关系数矩阵可视化图,主要包含以下几个步骤:
1. `par(mar=c(2,2,2,4))` 设置绘图的边缘空白区域大小,其中 mar 代表边缘空白区域大小(上下左右),这里分别设置为 2, 2, 2, 4。
2. `circos.par(gap.degree=c(3,rep(2, nrow(cor1)-1)), start.degree = 180)` 设定 circos 图的参数,其中 `gap.degree` 定义了每个扇形之间的间隔大小,第一个元素为扇形与画布环之间的距离,后面的元素表示每个扇形之间的距离,这里分别设置为 3 和 2,共有 nrow(cor1)-1 个扇形;`start.degree` 定义了起始角度,这里设置为 180 度。
3. `chordDiagram(cor1, grid.col=rainbow(ncol(rt)), col=col1, transparency = 0.5, symmetric = T)` 核心函数 `chordDiagram` 用于绘制基于圆环的相关系数矩阵可视化图。其中 `cor1` 表示相关系数矩阵;`grid.col=rainbow(ncol(rt))` 表示象征性地为每一列设置不同的颜色,这里使用了彩虹色;`col=col1` 表示用颜色矩阵 `col1` 对相关系数矩阵中的元素进行着色;`transparency = 0.5` 表示设置透明度为 0.5;`symmetric = T` 表示相关系数矩阵是对称的。
4. `par(xpd=T)` 用于在绘制完 circos 图后,将画布扩展到绘图边缘空白区域,以防止绘图被边缘空白区域截断。
阅读全文