RE: programming problem about elements taken

*To*: mathgroup at smc.vnet.net*Subject*: [mg72480] RE: [mg72470] programming problem about elements taken*From*: "David Park" <djmp at earthlink.net>*Date*: Sun, 31 Dec 2006 05:31:06 -0500 (EST)

Barrow, It's not completely clear what your grouping criterion is. After all, the skull bone is connected to the toe bone but they are quite far apart. I'm going to use a criterion of sorting all the numbers and then splitting them into groups whenever there is a gap greater than or equal to epsilon. Here is a test set. testset = Table[Random[], {20}] {0.854031, 0.1915, 0.671011, 0.875499, 0.0538857, 0.420417, 0.0711807, 0.240185, 0.672597, 0.0146134, 0.640515, 0.790636, 0.742558, 0.540169, 0.136822, 0.0476439, 0.725082, 0.482359, 0.414704, 0.411541} We sort the testset. step1 = Sort[testset] {0.0146134, 0.0476439, 0.0538857, 0.0711807, 0.136822, 0.1915, 0.240185, 0.411541, 0.414704, 0.420417, 0.482359, 0.540169, 0.640515, 0.671011, 0.672597, 0.725082, 0.742558, 0.790636, 0.854031, 0.875499} Pick an epsilon of 1/4 the maximum gap. epsilon = Max[step2]/4 0.075833 Now use the Split command to group numbers whenever the gap is less than epsilon. Split[step1, #2 - #1 < epsilon &] Length[%] {{0.0146134, 0.0476439, 0.0538857, 0.0711807, 0.136822, 0.1915, 0.240185}, {0.411541, 0.414704, 0.420417, 0.482359, 0.540169}, {0.640515, 0.671011, 0.672597, 0.725082, 0.742558, 0.790636, 0.854031, 0.875499}} 3 If that is what you want, then write and document a routine for your custom use. groupValues::usage = "groupValues[epsilon][valuelist] will sort and group the numbers in \ valuelist whenever the gap between two successive values is greater than or \ equal to epsilon."; groupValues[epsilon_?Positive][valuelist : {__?NumericQ}] := Split[Sort[valuelist], #2 - #1 < epsilon &] testset // groupValues[epsilon] Length[%] {{0.0146134, 0.0476439, 0.0538857, 0.0711807, 0.136822, 0.1915, 0.240185}, {0.411541, 0.414704, 0.420417, 0.482359, 0.540169}, {0.640515, 0.671011, 0.672597, 0.725082, 0.742558, 0.790636, 0.854031, 0.875499}} 3 It will almost always be worthwhile to write routines to extend Mathematica to conveniently do the things you want. Sometimes they might be quite extensive, and sometimes, as in this example, they might be simple combinations that are easier to use and save typing. If you build up sets of such routines you will soon find that Mathematica is much more powerful and convenient. David Park djmp at earthlink.net http://home.earthlink.net/~djmp/ From: Barrow [mailto:GRseminar at gmail.com] Dear all, I have a list of numbers (A),in fact, they are numbers distributed over [0,1]. Given parameter \epsilon, I have to choose elements from A such that their distances are larger than \epsilon, so the elements are "distringuishable". My goal is to find the maximum number of elements from A to meet the above "distinct criterion". How to do it with the functions build-in in mathematica ?? Thanks in advence. Sincerely Barrow