ggbiplot - how not to use the feature vectors in t

2020-06-27 18:54发布

I have a dataset data$cell_line.sva, which has dim of 313 11875.

cc.pca <- prcomp(data$cell_line.sva, center = TRUE, scale. = TRUE, retx = TRUE) 
g <- ggbiplot(cc.pca, obs.scale = 1, var.scale = 1, groups = as.factor(cgpResponse), ellipse = TRUE, circle = FALSE)

enter image description here

How can I get rid of the feature names? (the red text)

标签: r ggbiplot
4条回答
Summer. ? 凉城
2楼-- · 2020-06-27 18:57

You need to use the varname.size argument to do that.

Using the example from the documentation:

data(wine)
wine.pca <- prcomp(wine, scale. = TRUE)
print(ggbiplot(wine.pca, obs.scale = 1, var.scale = 1, 
               groups = wine.class, ellipse = TRUE, circle = TRUE))

enter image description here

and then add the varname.size argument:

print(ggbiplot(wine.pca, obs.scale = 1, var.scale = 1, 
               groups = wine.class, ellipse = TRUE, circle = TRUE,
               varname.size=0)) #set it to zero

enter image description here

And you have what you want!

查看更多
祖国的老花朵
3楼-- · 2020-06-27 19:05

I can't just comment because I don't have the necessary reputation points yet, but removing the arrows and the names works quite easily by setting var.axes to false:

ggbiplot(wine.pca, obs.scale = 1, var.scale = 1, 
         groups = wine.class, ellipse = TRUE, 
         circle = TRUE, var.axes=FALSE)
查看更多
成全新的幸福
4楼-- · 2020-06-27 19:14

I cannot quite figure out how this would yield a useful result, but here goes anyway. The names are not something that the function lets you suppress via parameter settings, at least to my reading of the code and help page. So looking at the code it looks as though the labels for the factors are pulled from the $rotations element of the prcomp object. Trying to set those names all to a blank character created an error so I succeeded by setting to the a value of varying lengths of blanks.

data(wine)    # need a reproducible example so use the help page 
 wine.pca <- prcomp(wine, scale. = TRUE)
print(ggbiplot(wine.pca, obs.scale = 1, var.scale = 1, groups = wine.class, ellipse = TRUE, circle = TRUE))
# that was the equivalent of your plot
# Now change the input value

dimnames(wine.pca$rotation)[[1]] <- 
   Reduce(function(x,y) paste0(x,y),    # function to concatentate the lanks
          rep(" ",dim(wine.pca$rotation)[1]),   # corrrect number of blanks
           acc=TRUE)                    # save all intermediate strings
 print(ggbiplot(wine.pca, obs.scale = 1, var.scale = 1, groups = wine.class, 
        ellipse = TRUE, circle = TRUE))
 #Look, Ma! No labels

enter image description here

查看更多
Ridiculous、
5楼-- · 2020-06-27 19:21

Ok, so I do this in a very brute way that does the job.

    print(ggbiplot(wine.pca, obs.scale = 1, var.scale = 1,
           groups = wine.class, ellipse = TRUE, circle = TRUE,
           varname.size=0, varname.adjust = 20))

It simply sets the offset for variable names beyond the plot limits.

查看更多
登录 后发表回答