K means for multidimensional data
显示 更早的评论
Hi everyone. I am trying to perform Raman spectral analysis using K-means clustering . I have 100 spectrums over 534 variables(in a matrix of 100 x 534).
Now I want to cluster 100 objects .How can I do so?
I am trying with this code, K= 12 found out by iteration. Now I have to find a plot of this for my data . Please help .
K=[ ];
sa=[ ];
for k=1:20
[idx c sumd]= kmeans(matrix,k);
sa= [sa sum(sumd)];
K= [K k];
end
plot(K,sa);// to find appropriate k
idx = kmeans(matrix,12);
gscatter(scoress(:,1),scoress(:,2),scoress(:,3),idx);//
now here I need to plot the data for all the columns rather than just 2 columns. How can I do so?
1 个评论
Image Analyst
2020-6-13
编辑：Image Analyst
2020-6-13
So you have 100 observations for each absorbance (wavenumber). The absorbance at each wavenumber are the features. And now you want 12 clusters which will classify each spectrum into one of 12 possible classes? Can you attach your matrix so we can try it?
采纳的回答
Image Analyst
2020-6-14
Well this is what I got so far
clc; % Clear the command window.
fprintf('Beginning to run %s.m ...\n', mfilename);
close all; % Close all figures (except those of imtool.)
clear; % Erase all existing variables. Or clearvars if you want.
workspace; % Make sure the workspace panel is showing.
format short g;
format compact;
fontSize = 15;
[numbers, strings, raw] = xlsread('data matrix.xlsx');
[rows, columns] = size(numbers)
wavenumbers = numbers(:, 1);
for col = 2 : columns
thisSpectrum = numbers(:, col);
plot(wavenumbers, thisSpectrum, '-');
grid on;
hold on;
end
title('All Raman Spectra', 'FontSize', 20);
xlabel('Wavenumber', 'FontSize', 20);
ylabel('Absorbance', 'FontSize', 20);
[classNumber, classCentroid] = kmeans(numbers(:, 2:end)', 12)
% Plot each clas separately
hFig = figure;
for col = 2 : columns
thisClass = classNumber(col - 1);
thisSpectrum = numbers(:, col);
subplot(3, 4, thisClass);
plot(wavenumbers, thisSpectrum, '-');
grid on;
hold on;
caption = sprintf('Class #%d', thisClass);
title(caption, 'FontSize', 20, 'Interpreter', 'none');
xlabel('Wavenumber', 'FontSize', 20);
ylabel('Absorbance', 'FontSize', 20);
end
hFig.WindowState = 'maximized'
but I'm not really sure kmeans is what you want to do, as you can see from the spectra plotted for each class. I might talk to my spectroscopists tomorrow and see if they have any ideas. They are really world class. What do you want me to ask him or her?
6 个评论
ananya mittal
2020-6-14
编辑：ananya mittal
2020-6-14
Image Analyst
2020-6-14
The means (cluster centers) are returned as classCentroid.
ananya mittal
2020-6-15
Image Analyst
2020-6-15
My spectroscopist say that kmeans is one way to group data. There may be better ways but they wouldn't know without more context. Why do you think there are 12 classes? Why not some other number? Would you rather compare these against some known reference spectra to see which reference spectra each one matches best?
I ran it again and this time I told kmeans to run 15 times to get the best guess at what the classes are. The new m-file is attached. Plus I told it to show how many spectra are in each class. You can see that some classes have only 1 spectra in them. Be aware that because of the random seeding nature of kmeans, the numbering of the classes could be different on each run, so class 1 from run 1 might not have the same spectra as class 1 on run 2. They might be in class 7 on run 2 instead of in class 1 again.
ananya mittal
2020-6-15
Image Analyst
2020-6-15
I don't know what the elbow method is. But there are ways to have kmeans decide what the best value of k is. Some other function I think - I don't remember what it is off the top of my head. Maybe it's 12 but maybe it's not.
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 k-Means and k-Medoids Clustering 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
