Применил следующий набор команд статистического пакета R ?
Код: Выделить всё
> setwd("C:/data/R Project/earthqwake") # установил рабочую директорию
> x <- read.table('earth.csv', sep=',', header=TRUE) # прочитал данные из файла в матрицу Х
> x <- x[-1] # удалил первый столбец с номером записи
> d <- dist(x, method = "euclidean", diag = FALSE, upper = FALSE) # расчет матрицы расстояний
> hc <- hclust(d) # построение дерева кластеров
> plot(hc) # отрисовал результат
Могу ли я в этом пакете получить каким-нибудь образом матрицу объединения наподобие той, которую выдает STATISTICA? Кстати, STATISTICA с поставленной задачей не справилась: вылетала молча на этапе объединения
Код: Выделить всё
Связь Объект 1 Объект 2 Объект 3 Объект 4 ... ... ... ... ... ... ... ... Объект N
0.0001 c_1 c_2
0.0001 c_5 c_345 c_515
...........................
0.2324 c_3 c_6 c_453 c_345 c_123 c_543 c_231
...........................
2.6543 c_3 c_6 c_453 c_345 c_123 c_543 c_231 ... ... ... ... ... ... ... ... ... c_N
где с_i - порядковый номер записи в файле данных i [1, N]
Стоит "обратная" задача: полученные кластеры с определенным уровнем связи "вынести на местность" и заполнить значениями.
Задача является продолжением и расширением темы