Я пытаюсь визуализировать результаты PCoA{ape}
путем создания биполя в R. Теперь оси получают метки меток по умолчанию 1 и ось 2, но я хочу отредактировать это.Изменение меток осей для biplot() в R
Это код, который я пробовал:
biplot(pcoa.ntK, Y=NULL, plot.axes=c(1,2), rn=ntnames,
xlabs="PC1 (%)", ylabs="PC2 (%)")
Но метки не изменяются. Может кто-нибудь сказать мне, что я делаю неправильно здесь? И я также хотел бы отредактировать название, все советы для этого?
Мои данные:
ntK <- matrix(
c(0.00000, 0.01500, 0.01832, 0.02061, 0.01902, 0.01270, 0.02111, 0.01655, 0.01520, 0.01691,
0.01667, 0.00000, 0.01175, 0.01911, 0.01759, 0.01127, 0.01854, 0.01041, 0.00741, 0.02007,
0.02432, 0.01404, 0.00000, 0.02551, 0.01972, 0.01838, 0.02505, 0.01484, 0.01391, 0.02687,
0.01501, 0.01252, 0.01399, 0.00000, 0.01442, 0.01294, 0.01402, 0.01132, 0., 0.01455,
0.02343, 0.01951, 0.01830, 0.02440, 0.00000, 0.01727, 0.02470, 0.02021, 0.01699, 0.02482,
0.01320, 0.01054, 0.01439, 0.01847, 0.01457, 0.00000, 0.01818, 0.01366, 0.00977, 0.01394,
0.02468, 0.01950, 0.02206, 0.02251, 0.02343, 0.02040, 0.00000, 0.02028, 0.01875, 0.02558,
0.02254, 0.01276, 0.01522, 0.02117, 0.02234, 0.01790, 0.02363, 0.00000, 0.01152, 0.02557,
0.01804, 0.00792, 0.01244, 0.02019, 0.01637, 0.01116, 0.01904, 0.01004, 0.00000, 0.02099,
0.01862, 0.01988, 0.02227, 0.02200, 0.02218, 0.01476, 0.02408, 0.02066, 0.01947, 0.00000),
nrow=10,
ncol=10)
library(ape)
ntnames <- c("A","B","C","D","E","F","G","H","I","J")
pcoa.ntK <- pcoa(ntK)
как насчет данных 'rednt'? – 989
Кажется, что 'rednt' должен быть' pcoa.ntK', а 'namesaa' должен быть' ntnames' на вызове 'biplot'. – Molx
Как было сказано ниже, заголовок оси жестко закодирован, где они берут свои имена из названий столбцов 'pcoa.ntK $ vectors'. Поэтому быстрое решение - переименовать это. Так что если вы планируете первые два компьютера. используйте 'colnames (векторы pcoa.ntK $) [1: 2] <- c (" Mytitle1 "," Mytitle2 ")', затем зарисуйте, как раньше. – user20650