Я использую пакет pca3d
для построения основных компонентов 1-5 данных экспрессии генов. Мне было трудно стандартизировать трехмерный график, поэтому вместо этого я нарисовал ПК1 против ПК2, ПК2 против ПК3 и т. Д. Мой код для этого ниже
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
не влияют на сюжет. Я обнаружил в этом посте , в котором предлагалось использовать аргумент asp
, но установка asp
не позволила мне изменить оси.
Кто-нибудь знает, как я мог заставить pca2d
изменить пределы осей?