Мне интересно, какой будет правильный код в R для расчета расстояния между групповыми центроидами?Я использовал функцию betadisper в веганской системе на матрице различий, используя различия между Брейем-Кертисом, которая дает среднее расстояние до центроидов внутри групп, но теперь я хочу вычислить разницу между этими двумя центроидами.
Я пробовал функцию dist_between_centroids в пакете usedist
data(iris)
dist.iris<-vegdist(iris[,1:4],method = "bray") #calculate Bray-Curtis dissim
dist.iris
#avg distance to group centroids
permdisp<-betadisper(dist.iris,iris$Species,type="centroid")
permdisp
library(usedist)
dist_between_centroids(dist.iris,1:50, 51:100) #1:50 = group setosa, 51:100 = group versicolor
Таким образом, среднее расстояние до центроида для группового сетаза составляет 0,04, а для группового многоцветного также 0,04.Для расстояния между групповыми центроидами сетозы и версиколора я получаю вывод 0,22.Тогда это правильное расстояние между центром тяжести сетозой группы и центроидом версиколора или я что-то упустил?
Спасибо!