6

我想知道是否可以为图中的列名而不是标签设置不同的形状。我已经展示了一个可重现的例子来获得情节。

#dataset
mtcars

#split the data into different blocks
a<-mtcars[,1:3]
b<-mtcars[,4:6]
c<-mtcars[,7:9]
d<-mtcars[,10:11]

### create a list of dataframes
A<-list(blockA=a, blockB=b,blockC=c,blockD=d)

### design matrix
C <- matrix(c(0, 0, 0, 0,  
              0, 0, 0, 0, 
              0, 0, 0, 0, 
              1, 1, 1, 1), 4, 4)
tau = c(1, 1, 1, 1)
library(mixOmics)

##### computing the sGCCA
wrap.result.sgcca = wrapper.sgcca(data = A, design = C, penalty = rep(1, length(A)),
  ncomp = c(2, 2, 2,2), scheme = "centroid", verbose = FALSE)

##### visualization
col <- rep(c("blue", "red", "darkgreen", "darkviolet"))

plotVar(wrap.result.sgcca, col=col, comp = c(1,2), 
  block = c(1,2,3,4), ncomp.select = c(1,2), labels = TRUE)

相关圆图

在这个情节中，我设法用独特的颜色给块上色。我还想用形状来表示列。例如，当我们考虑数据集 mtcars 的不同列名时

colnames(mtcars)
[1] "mpg"  "cyl"  "disp" "hp"   "drat" "wt"   "qsec" "vs"   "am"   "gear" "carb"

根据我拆分的块 a(columns1:3)、b(columns(4:6)、c(columns 7:9)、d(columns 10:11)，我想要第 1、4、7,10 列具有相同的形状（例如圆形），第 2,5,8,11 列（例如正方形），第 3,6,9 列（例如星形）。我尝试设置 pch 值，但未能根据列名生成形状这个情节。有人可以帮我吗？谢谢

4

0 回答 0