Форумы GIS-Lab.info

Код: Выделить всё

x <- read.table('quake.dat', sep=' ', header= FALSE) # читаем данные 
x<-scale(x,center = TRUE, scale = TRUE) # центрируем, масштабируем

library (StatMatch)
d2 <- mahalanobis.dist(data.x=x, data.y=NULL, vc=NULL) # расчет дистанции
hcd <- hclust(as.dist(d2), 'ward') # кластеризация

gc(reset=TRUE) # оптимизация памяти
library (MASS) # пакет для записи матриц в файл
write.matrix(hcd$height, file='cl_mah.dat')
write.matrix(hcd$merge, file='cl_mah.txt')
plot(hcd)

rm(hcd); rm(d2);gc(reset=TRUE)# удаляем ненужные объекты и оптимизируем память

d <-  dist(x, method = "euclidean", diag = FALSE, upper = FALSE) # расчет дистанции
hc <- hclust(d^2, 'ward') # кластеризация
write.matrix(hc$height, file='cl_quake.dat')
write.matrix(hc$merge, file='cl_quake.txt')
plot(hcd)