Durga Lal Shrestha

Sorry for previous rating, it was problem due to the FLEX website. I clicked back button on browser not "Submit" button.

I was searching the function that finds the kth largest or smallest values and locations of 2d array and thought that this is the function. But it appears it produces different results that I expected

list = 1:16;
list = reshape(list,4,4)
[res loc]= mink(list,5)

list =

1 5 9 13
2 6 10 14
3 7 11 15
4 8 12 16

res =

1 5 9 13
2 6 10 14
3 7 11 15
4 8 12 16

loc =

1 1 1 1
2 2 2 2
3 3 3 3
4 4 4 4

I thought it produces the max or min of the whole array. So I have written my own program minmax which produces the results

[res,loci,locj] = minmax(list,5)

out =

1
2
3
4
5

loci =

1
2
3
4
1

locj =

1
1
1
1
2

In terms of speed, I do not think it is efficient, see below

list = 1:1000000;
list = reshape(list,1000,1000);
tic
[res loc]= mink(list,5);
toc

%%
tic
[out,loci,locj] = minmax(list,5);
toc

Elapsed time is 0.046343 seconds.
Elapsed time is 0.039491 seconds.

Does not return handles of subplot in SubHdl as these are deleted in line 160.
Handles should be of new object created as:

SubHdl(k) = copyobj(OldGca, FigHdl);
set(SubHdl(k), 'Position', NewPos);

Also it does not produce legends, otherwise it is useful.

It does not seems to compatible to subplot.

For those looking for a weighted median: wmedian = @(X) wprctile(X,50);

How should I adjust it so that I can show gradients between labels?