Траектория в рукоположении с использованием metaMDS в R

Я использую vegan пакет в R для выполнения неметрического многомерного масштабирования с использованием функции metaMDS, Для этого я использую файл equilibria0.txt, который содержит 101 вид в столбцах (первый столбец представляет пустые сайты в моей модели решетки) и 500 сообществ в строках. Я использую следующий код R для выполнения NMDS:

library(vegan)        # perform NMDS
library(cluster)      # calculate dissimilarities in data, silhouette scores, ...
library(RColorBrewer) # better colour options

comm <- read.table("equilibria0.txt")   # communities in rows, species in columns
comm <- comm[,-1]     # delete the first column (the empty sites)
comm <- decostand(comm, method = "total")*100

NMDS <- metaMDS(comm, # our community-by-species matrix
                k=2)  # the number of reduced dimensions
k <- 2                # number of clusters

# introduce some colours for the clusters
lightcolours <- c("darkolivegreen3", "cadetblue3")

# calculate Bray-Curtis distance among samples
comm.bc.dist <- vegdist(comm, method = "bray")
# cluster communities using average-linkage algorithm
comm.bc.clust <- hclust(comm.bc.dist, method = "average")

# k is the number of clusters we want
# every LC is assigned a number from 1 to k
comm.bc.cut <- cutree(comm.bc.clust, k=2)

# set up the plotting area but don't plot anything yet
mds.fig <- ordiplot(NMDS, type = "none")
# plot just the samples, colour by cluster, pch=19 means plot a circle
for (i in 1:k) {
   points(mds.fig,
          "sites",
           pch = 19,
           col = lightcolours[i],
           select = comm.bc.cut == i)
}

Графический вывод показывает мне, что система имеет два разных кластера, как я и ожидал:

Теперь у меня есть еще 30 файлов, все идентичные equilibria0.txt, за исключением первой строки. Иными словами, из 500 местных сообществ 499 не изменились, а одна (первая строка) изменилась. Это потому, что я делаю небольшое изменение в местном сообществе, чтобы посмотреть, будет ли оно вовремя эволюционировать в другой кластер после этого изменения. Поэтому я хотел бы получить всю информацию из этих 30 файлов на одном графике, в частности, показать траекторию, как на следующем графике, который можно найти в статье Gibson et al.

Конкретно, есть одно местное сообщество (точка в плоскости NMDS), представляющее интерес. Его координаты относительно других изменится, если вы перейдете из equilibria0.txt в equilibria1.txt в equilibria2.txt и т.д., которые на самом деле все файлы данных на разных временных шагах. Поэтому я хочу представить эволюцию этого конкретного местного сообщества, предпочтительно по линии траектории, подобной приведенной выше. Я предполагаю, что остальные 499 точек (сообществ) не будут значительно сдвигаться в своих координатах.

0 ответов

Другие вопросы по тегам