我正在使用pca3d
包来绘制基因表达数据的主要组成部分1-5。我很难标准化3d绘图,所以相反我绘制了PC1与PC2,PC2与PC3等。我的代码如下
par(mfrow = c(2,2))
for(n in 1:4){
pca2d(pcaSix$x[,n:(n+1)], # Mclust Grouping - module Genes
title = "Six Genes From Abel, et al. 2011",
shape = 16,
group = clustersMod1six,
radius = 1,
legend = NULL,
show.centroids = FALSE,
show.ellipses = TRUE,
col = ClusterColorsSix,
axe.titles = c(paste("PCA",
n,
sep = ""),
paste("PCA",
n+1,
sep = "")
),
xlim = c(-1,1),
ylim = c(-1,1)
)
}
我遇到的问题是xlim
和ylim
论证似乎没有对情节有任何影响,我发现this post建议使用asp
参数,但设置asp
不允许我改变轴。
有谁知道我怎么能强迫pca2d
改变轴限制?
在不修改功能的情况下,无法控制轴。见line 32:
plot(NULL, type= "n",
xlim= prange$x,
ylim= prange$y,
xlab= xlab,
ylab= ylab,
bty= "none"
)