kmeans matlab code feed собственный источник данных

#matlab #cluster-analysis #k-means

#matlab #кластерный анализ #k-означает

Вопрос:

Я хочу попробовать этот код кластеризации K-means в моем собственном файле, как мне изменить его, чтобы он не создавал случайную информацию, а считывал ее из моего собственного источника данных?

 %% generate sample data
K = 3;
numObservarations = 100;
dimensions = 3;
data = rand([numObservarations dimensions]);

%% cluster
opts = statset('MaxIter', 500, 'Display', 'iter');
[clustIDX, clusters, interClustSum, Dist] = kmeans(data, K, 'options',opts, ...
    'distance','sqEuclidean', 'EmptyAction','singleton', 'replicates',3);

%% plot data clusters
figure, hold on
scatter3(data(:,1),data(:,2),data(:,3), 50, clustIDX, 'filled')
scatter3(clusters(:,1),clusters(:,2),clusters(:,3), 200, (1:K)', 'filled')
hold off, xlabel('x'), ylabel('y'), zlabel('z')

%% plot clusters quality
figure
[silh,h] = silhouette(data, clustIDX);
avrgScore = mean(silh);


%% Assign data to clusters
% calculate distance (squared) of all instances to each cluster centroid
D = zeros(numObservarations, K);     % init distances
for k=1:K
    %d = sum((x-y).^2).^0.5
    D(:,k) = sum( ((data - repmat(clusters(k,:),numObservarations,1)).^2), 2);
end

% find  for all instances the cluster closet to it
[minDists, clusterIndices] = min(D, [], 2);

% compare it with what you expect it to be
sum(clusterIndices == clustIDX)
  

Комментарии:

1. Можете ли вы четко определить, какая часть опубликованного вами кода генерирует случайные данные? Если да, то это переменная, в которой вы должны хранить свой набор данных.

2. Пожалуйста, расскажите нам больше о вашем источнике данных. Находится ли он в базе данных, в файле (какого формата?), …

3. @Laurent, прочитайте другие вопросы OP.

4. @Nzbuu Не нужно читать историю вопросов OP, чтобы выяснить детали. Каждый вопрос должен быть автономным, и если OP не может приложить минимальные усилия к своему собственному вопросу (на который ему предположительно нужен ответ) и прояснить его, почему другие должны?

5. @yoda, я полностью согласен. Просто ссылаясь на объем комментариев OP по этому вопросу.

Ответ №1:

Строка, которая создает случайные данные, является:

 data = rand([numObservarations dimensions]);
  

Просто замените эту строку кодом, который считывает ваши данные (возможно, с помощью команды matlab, такой как textscan ) в переменную с именем data .