OK, so from what I understand, you are clustering the genes, and not the samples, which means that you need to transpose your sample-gene matrix, and then run kmeans. If you then type in:
IDX = kmeans(genedata',20);
IDX will contain 100 elements, which contains numbers from 1 to 20, corresponding to which cluster the gene is in.
for i = 1:20
genenum(i) = sum(IDX == i);
end
and genenum will contain the number of genes in each cluster. You could alternatively take the histogram of IDX, and the resulting array will contain the same info as genenum.
If you'd want to include the genes included in each cluster as well, you can't store it in the same matrix as the numbers of elements are not the same in each 'cell' of the matrix, and hence I advise using cell variables instead. Hope this answers your questions!
Best Answer