forming a cell array from matrix

Question

0 个投票

I want to include the genes included in each cluster .

i have genedata =[1:1:100]

the code given below is for number of genes in cluster. Now i want to include the genes included

genedata=[1:1:100]
IDX = kmeans(genedata',20)
for i = 1:20
  genenum(i) = sum(IDX == i);
end

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

Answer 1

Walter Roberson 2011-11-19

0 个投票

First off, the loop you have is okay, but you could also replace it with the more efficient

genenum = accumarray(IDX(:));

As for including the genes included: you are not clear as to what you mean by that, so I will guess:

genes = cell(20,1);
for K = 1 : 20
  genes{K} = genedata(IDX==K);
end

FIR 2011-11-21

thanks andrei

FIR 2011-11-21

andrei if i do yhe following operation i get error

Error using ==> horzcat

CAT arguments dimensions are not consistent.

Error in ==> v at 10

fix=[j genenum genes]

where j=[1:20]'

please help