# A Greedy Clustering Algorithm Based on Interval Pattern Concepts and the Problem of Optimal Box Positioning.

1. Introduction

We consider the problem of clustering, that is, splitting a finite set X [subset] [R.sup.d] into disjoint subsets (called clusters) in such a way that points from the same cluster are similar (with respect to some criterion) and points from different clusters are dissimilar (see, e.g., ). It is convenient to present the input data in the form of a numerical context (table) whose rows correspond to objects and columns correspond to attributes of the objects.

Formal concept analysis (FCA) is a data analysis method based on applied lattice theory and order theory. The object-attribute binary relation is visualized with the use of the line diagram of the concept lattice. Within the framework of this theory a formal concept is defined as a pair (extent, intent) obeying a Galois connection (for exact definitions see the monograph  by Ganter and Wille).

There exist several generalizations of FCA to fuzzy and numerical contexts. One of them is known as the theory of pattern structures introduced by Ganter and Kuznetsov in . An important particular case of pattern concepts, which are the key object in the theory of pattern structures, is interval pattern concepts with the operation of interval intersection. Interval pattern concepts allow one to apply cluster analysis to rows of formal numerical contexts. In this case the criterion of similarity consists in belonging of all the differences between the values of the corresponding attributes to given intervals.

It can be easily seen that the problem of finding an interval pattern concept of maximum extent size (i.e., cardinality) can be reformulated as the problem of the optimal positioning of a d-dimensional box with given edge lengths for the given set X, that is, finding a position of the box that maximizes the number of points of the set X enclosed by the box (the details are given below in Section 2.2).

The existing algorithms that solve the problem of finding the optimal position of a box do not allow one to obtain an exact or at least approximate solution for high-dimensional data within a reasonable time (see a detailed survey in Section 2.2). The main goal of this paper is to propose a greedy algorithm which gives an approximate solution to this problem and a clustering algorithm based on the optimal positioning problem. We propose a clustering algorithm with

O((dn log (n) + [d.sup.3][n.sup.1-1/d]/[s.sub.min] f(n, d))n/[c.sub.min]) (1)

worst-case time and O(dn) space complexity, where f(n, d) denotes the number of iterations of the main stage of the algorithm, and parameters [s.sub.min] and [c.sub.min] regulate the duration of each iteration. Greater number of iterations and greater duration of each iteration provide better approximation.

The rest of the paper is organized as follows. In Section 2 we introduce the main definitions and formalize the statement of the problem. In Sections 3 and 4 we formulate our algorithms. In Sections 5 and 6 we describe the validation results and make some concluding remarks.

2. Main Definitions and Statement of the Problem

In this section we start with the main definitions from the theory of formal concepts and then present a geometrical reformulation of the problem of finding the interval pattern concept of maximum extent size (we call it simply the maximum interval pattern concept).

2.1. Formal Concepts. Let us recall the main definitions which we need to formalize our clustering method based on interval pattern concepts. Additional details can be found in [2, 3].

Definition 1. An upper (lower) semilattice is a partially ordered set (M, [less than or equal to]) such that for any elements x, y [member of] M there exists a unique least upper bound (greatest lower bound, resp.).

Definition 2. A semilattice operation on the set M is a binary operation [??]: MxM that features the following properties for a certain e [member of] M and any elements x, y, z [member of] M:

(i) x [??] x = x (idempotency).

(ii) x [??] y = y [??] x (commutativity).

(iii) (x [??] y) [??] z = x [??] (y [??] z) (associativity).

(iv) e [??] x = e.

Definition 3. A lattice is an ordered set (L, [less than or equal to]) which is at the same time an upper and a lower semilattice.

Definition 4. Let (P, [[less than or equal to].sub.P]) and (Q, [[less than or equal to].sub.Q]) be partially ordered sets. 0A Galois connection between these sets is a pair of maps [phi]: P [right arrow] Q and [psi]: Q [right arrow] P (each of them is referred to as a Galois operator) such that the following relations hold for any [p.sub.1], [p.sub.2] [member of] P and [q.sub.1], [q.sub.2] [member of] Q:

(i) [p.sub.1] [[less than or equal to].sub.P] [p.sub.2] [right arrow] [phi]([p.sub.1]) [[greater than or equal to].sub.Q] [phi]([p.sub.2]) (anti-isotone property).

(ii) [q.sub.1] [[less than or equal to].sub.Q] [q.sub.2] [right arrow] [psi]([q.sub.1]) [[greater than or equal to].sub.P] [psi]([q.sub.2]) (anti-isotone property).

(iii) [p.sub.1] [[less than or equal to].sub.P] [psi]([phi]([p.sub.1])) and [q.sub.1] [[less than or equal to].sub.Q] [phi]([psi]([q.sub.1])) (isotone property).

Applying the Galois operator twice, namely, [psi]([phi](p)) and [phi]([psi](q)), defines a closure operator.

Definition 5. A closure operator [bar.(*)] on M is a map that assigns a closure [bar.X] [??] M to each subset X [??] M under the following conditions:

(i) X [less than or equal to] Y [right arrow] [bar.X] [less than or equal to] [bar.y] (monotony).

(ii) X [right arrow] [bar.X] (extensity).

(iii) [??] = [bar.X] (idempotency).

Definition 6. A pattern structure is a triple (G, (D, [??]), [delta]), where G is a set of objects, (D, [??]) is a meet-semilattice of potential object descriptions, and [delta]: G [right arrow] D is a function that associates descriptions with objects.

The Galois connection between the subsets of the set of objects and the set of descriptions for the pattern structure (G, (D, [??]), [delta]) is defined as follows:

[mathematical expression not reproducible]. (2)

Definition 7. A pattern concept of the pattern structure (G, (D, [??]), [delta]) is a pair (A, d), where A [??] G is a subset of the set of objects and d [member of] D is one of the descriptions in the semilattice, such that [A.sup.[??]] = d and [d.sup.[??]] = A; A is called the pattern extent of the concept and d is the pattern intent.

A particular case of a pattern concept is the interval pattern concept. The set D consists of the rows of a numerical context, which are treated as tuples of intervals of zero length. An interval pattern concept is a pair (A, d), where A is a subset of the set of objects and d is a tuple of intervals with ends determined by the smallest and the largest values of the corresponding component in the descriptions of all objects in A.

Interval pattern concepts are convenient to use in the analysis of numerical contexts, when there is a need to divide all data into clusters that comprise objects in which the numerical data is similarly "distributed" in the rows.

For each component of an interval pattern concept we introduce the width [sigma]: the difference between the largest and the smallest values of the component. Then a clustering procedure can be defined using a standard greedy approach. Specifically, at each step the maximum interval pattern concept is identified, that is, an interval pattern concept with the maximum number of objects, whose width with respect to each component does not exceed a predefined [sigma]. The objects of the identified interval pattern concept are combined into a cluster and excluded from the set of objects analyzed at subsequent steps.

In Example 1 presented in Table 1 the objects are pupils and the numerical data of the context consist of the grades they got at exams in various disciplines.

We need to divide the set of pupils into clusters in such a way that the grades of pupils in the same cluster differ by at most 1 for each of the disciplines. Such a setting corresponds to [sigma] = 1; in this case we obtain 6 clusters (interval pattern concepts whose width is not greater than 1), each containing one pupil. In the case [sigma] = 2 we arrive at the same 6 clusters.

When [sigma] = 3 we have five clusters [[A, D].sup.[??]] = {[8, 9], [9, 9], [9,10], [6,9]}, {B}, {C}, {E}, {F}, and in the case [sigma] = 4 we obtain three clusters [{A,C,D}.sup.[??]] = {[6,9], [5,9], [9,10], [6,9]}, [{B,E}.sup.[??]] = {[8, 8], [2,4], [6,6], [5,9]}, {F}.

Example 2. In the conditions of the previous example let us set [[sigma].sub.1] = 1, [[sigma].sub.2] = 1, [[sigma].sub.3] = 10, [[sigma].sub.4] = 3. Then the set of pupils can be divided into four clusters {A,D}, {C,F}, {B}, {E}:

[mathematical expression not reproducible]. (3)

Clustering methods based on interval pattern concepts find applications in the analysis of experimental data. For instance, applications of such methods to gene expression analysis were discussed in [4, 5].

2.2. Geometry. Let P be a set of n points in [R.sup.d] (d [member of] N) and [[delta].sub.1], [[delta].sub.2], ..., [[delta].sub.d] be a set of positive real numbers.

Definition 8. A d-orthotope (also called a box) with center x = ([x.sub.1], ..., [x.sub.d]) [member of] [R.sup.d] and edge lengths [[delta].sub.1], [[delta].sub.2], ..., [[delta].sub.d] is the Cartesian product of the intervals

[[x.sub.1] - [[delta].sub.1]/2, [x.sub.1] + [[delta].sub.1]/2] x *** x [[x.sub.d] - [[delta].sub.d]/2, [x.sub.d] + [[delta].sub.d]/2]. (4)

It can be easily seen that the problem of identification of maximum interval pattern concept can be reformulated in terms of finding the optimal position of the box with the edge lengths [[delta].sub.1], [[delta].sub.2],..., [[delta].sub.d], that is, maximizing the number of points of set P enclosed by the box. This formulation can be generalized to the problem of finding the optimal position of a ball in an arbitrary metric space, since any box can be treated as a ball in the stretched [L.sub.[infinity]] metric in which the distance [rho](x, y) between the points x = ([x.sub.1],..., [x.sub.d]) and y = ([y.sub.1], ..., [y.sub.d]) is defined as

[mathematical expression not reproducible]. (5)

The problem of optimal positioning has been well studied for d = 2: some lower and sharp upper bounds on complexity are known (see, e.g., [6, 7]). However, to the best of our knowledge for the case of an arbitrary dimension d no lower bounds and no efficient exact algorithms are available so far. de Figueiredo and da Fonseca noted  that the problem can be solved exactly in roughly 0([n.sup.d+1/d]) time by projecting the points onto a (d + 1)-dimensional paraboloid and using half-space range searching data structures . In the same paper for the case of weighted points under certain additional restrictions they also obtained a lower bound [OMEGA]([n.sup.d]) for exact algorithms and indicated that existing algorithms for the unweighted version of the problem do not beat this lower bound in the worst case. Eckstein et al. showed that a generalization of the problem of optimal positioning whose input also includes a set of prohibited points is NP-hard .

Known approximate algorithms for optimal positioning also have time complexity which depends on d exponentially. For example, de Figueiredo and da Fonseca suggested an approximate algorithm  which solves the problem in worst-case time 0([3.sup.d]n/[[epsilon].sup.d-1]), where 0 < [epsilon] < 1 is a given approximation parameter. Due to exponential dependence on d these approximate algorithms are also practically inapplicable in the case of high dimension, and there is a need to develop an algorithm which can produce an approximate solution in reasonable time.

3. A Greedy Algorithm for Finding an Approximately Optimal Position of a Box

In this section we present a greedy algorithm for finding an approximately optimal position of a box with edge lengths [[delta].sub.1], [[delta].sub.1], ..., [[delta].sub.d] for a set P = [{[P.sub.i]}.sup.n.sub.i=1] [subset] [R.sup.d] (the order in which points are listed in P is insignificant). This algorithm is auxiliary for the clustering method described in Section 4.

The algorithm has several input parameters: positive real numbers s, [s.sub.min], [lambda] < 1, and a function f: N x N [right arrow] N. The parameters s, [s.sub.min], and [lambda] regulate the duration of one iteration. The function f takes the values n and d as inputs and returns the number of iterations at the main stage of the algorithm. Greater number of iterations and greater duration of each iteration provide better approximation.

The algorithm includes two basic stages: the preprocessing stage and the main stage.

3.1. Preprocessing

(1) At the first stage of our algorithm the box with the edge lengths [[delta].sub.1], [[delta].sub.2], ..., [[delta].sub.d] is transformed into the unit cube (we call it simply the cube) by means of dividing the ith coordinate of each point by [[delta].sub.i], i = 1, ..., d. This stage can be performed in O(dn) operations.

(2) We consider the integer lattice with edges of length 1, compute the number of points of P in each cell, and denote the cell that contains the maximum number of points by [C.sub.0]. The cell [C.sub.0] is called the base cube. Let [y.sub.0] [member of] [R.sup.d] denote the center of [C.sub.0]. This stage requires O(dn) operations as well.

(3) At the final step of the preprocessing stage we build a k-d tree data structure (which is used at the main stage to organize the fast range search) in 0(dn log(n)) operations with the space complexity of o(dn) (see [11, 12]).

3.2. The Main Stage. Let [mathematical expression not reproducible] denote the function which counts number of points of the set P in an arbitrary subset of [R.sup.d]. The main idea of our algorithm consists in constructing a finite sequence of cubes that starts from a random point y in the base cube and satisfies the condition that the next cube contains more points than the previous one. Let [D.sup.y.sub.1], ..., [D.sup.y.sub.k(y)], denote these cubes with centers [z.sup.y.sub.1], ..., [z.sup.y.sub.k(y)], respectively. In our notation we have [z.sup.y.sub.1] = y and q([D.sup.y.sub.i]) < q([D.sup.y.sub.i+1]) for all i [member of] {1, ..., k(y) - 1}. After f(n, d) iterations the algorithm returns a locally optimal cube C.

Definition 9. The t-neighborhood of a cube D with center x = ([x.sub.1],..., [x.sub.d]) is the set consisting of all cubes with centers at points of the form ([x.sub.1],..., [x.sub.i-1], [x.sub.i] [+ or -] t, [x.sub.i+1],..., [x.sub.d]) for all i [member of] {1,...,d}, that is, all cubes obtained through translation of D along one of the axes by the distance [+ or -] t.

Now we describe the procedure of constructing the sequence of cubes. Let y be an arbitrary point in the base cube [C.sub.0] and [z.sup.y.sub.1] = y, [D.sup.y.sub.1] be the cube with center at [z.sup.y.sub.1], [s.sub.1] = s. In order to get a definite estimate on the precision of the algorithm (see Theorem 11) we initialize the first iteration deterministically by taking the center of [C.sub.0] as y. Other iterations are initialized randomly.

Suppose that the cubes [D.sup.y.sub.1], ..., [D.sup.y.sub.m] with centers [z.sup.y.sub.1],..., [z.sup.y.sub.m], respectively, and the numbers [s.sub.1],..., [s.sub.m] have been already constructed. There are two possible cases.

(1) If there exists a cube D in the sm-neighborhood of [D.sup.y.sub.m] such that q(D) > q([D.sup.y.sub.m]), then we set [D.sup.y.sub.m+1] = D, take the center of D as [z.sup.y.sub.m+1], and take [s.sub.m+1] = [s.sub.m]. In other words, if there exists a cube in the [s.sub.m]-neighborhood of the current cube which contains more points of P, then we move the current cube to this position.

(2) If there are no such cubes (i.e., all cubes in the [s.sub.m]-neighborhood of the current cube contain at most the same number of points), then we set [D.sup.y.sub.m+1] = [D.sup.y.sub.m], [z.sup.y.sub.m+1] = [z.sup.y.sub.m], and [s.sub.m+1] = [s.sub.min] (i.e., decrease the current step size). If [s.sub.m+1] < [s.sub.min] (the step size threshold is reached), then the procedure is ended and [D.sup.y.sub.m] is returned as the procedure result.

In order to obtain acceptable time complexity we impose additional restrictions on the selection of the next cube. These assumptions are necessary to avoid the situation where the length of the sequence grows exponentially with d. Validation on experimental data confirmed that these restrictions do not essentially affect the clustering results.

Restriction 1. All cubes in the sequence must have common points with the base cube [C.sub.0].

In Figure 1 we present an example of a set P for which this requirement causes a significant difference between the exact solution and the solution produced by the algorithm. However, this difference is essentially reduced at further steps of the clustering algorithm as generally it affects only the order in which clusters are constructed.

Restriction 2. For each individual coordinate it is not allowed to translate the cube in the opposite directions at different steps of the procedure described above.

The above restrictions lead to the following lemma.

Lemma 10. The main stage of the algorithm has

O ([d.sup.3] [n.sup.1-1/d]/[s.sub.min] f(n,d)) (6)

worst-case time complexity.

Proof. First we get an upper estimate for the length k(y) of the sequence

of cubes (for an arbitrary y [member of] [C.sub.0]). Due to Restrictions 1 and 2 we have

k(y) [less than or equal to] O(d/[s.sub.min]). (7)

Thus, Restrictions 1 and 2 avoid the situation where the length of the sequence grows exponentially with d. Each step of the procedure of constructing the sequence of cubes requires 2d evaluations of the function q for the cubes (i.e., 2d range searches). With the use of a k-d tree the range search can be performed with 0(d[n.sup.1-1/d]) worst-case time complexity (see ). The procedure of constructing the sequence of cubes involves f(n, d) iterations, so the above complexity bound holds.

Note that we also have a trivial estimate k(y) [less than or equal to] n, as q([D.sup.y.sub.m]) grows, and hence k(y) [less than or equal to] q([D.sup.y.sub.k(y)])) [less than or equal to] n. Thus, without the imposed restrictions the worst-case complexity estimate

0([d.sup.2] [n.sup.2-1/d] f(n,d)) (8)

holds, and hence the restrictions can be omitted without violation of practical feasibility in case if the number of objects n has the same order as the dimensionality d.

3.3. Precision and Complexity of the Algorithm

Theorem 11. Let [y.sub.0] be a center of [mathematical expression not reproducible] be a cube produced by an algorithm iteration which was initialized with [y.sub.0] (and so for this iteration [mathematical expression not reproducible], and [D.sub.opt] be an optimal cube (i.e., [D.sub.opt] [member of] arg max q(D), where maximum is taken over all unit cubes in [R.sup.d]). Then

1/[2.sup.d] [less than or equal to] q ([D.sub.alg]/[D.sub.opt]) [less than or equal to] 1, (9)

and this estimate is sharp.

Proof. The upper estimate is trivial. The lower estimate follows from the fact that [D.sub.opt] is covered by at most [2.sup.d] cells of the integer lattice with edges of length 1, and hence

q([D.sub.opt]) [less than or equal to] [2.sup.-d] q ([C.sub.0]) [less than or equal to] [2.sup.-d] q ([D.sub.alg]). (10)

An example that shows that the estimate is sharp is similar to the example from Figure 1. For example, we can locate the center of [D.sub.opt] at the integer lattice node and put [2.sup.d] points in [D.sub.opt] in such a way that each cell of the integer lattice contains at most one of these points. Then, we select an arbitrary cell of the integer lattice that is distant from [D.sub.opt] and put one point to this cell, which becomes [C.sub.0].

Theorem 12. The algorithm for finding an approximately optimal position of the box has O (dn log (n) + [d.sup.3][n.sup.1-1/d]/[s.sub.min] f(n, d))) (11)

worst-case time complexity and O(dn) space complexity.

Proof. Combining the estimates for the time and space complexity of the preprocessing stage and the main stage of the algorithm gives the bounds mentioned above.

Note that omitting Restrictions 1 and 2 results in the worst-case time complexity estimate

O([d.sup.2] [n.sup.2-1/d] f(n, d)). (12)

4. Clustering Algorithm

Now let us consider the clustering problem, that is, the problem of splitting the given set P = [{[p.sub.i]}.sup.n.sub.i=1] [subset] [R.sup.d] into mutually disjoint subsets [C.sub.1],..., [C.sub.k]. Following interval pattern concept approach, we construct clusters with controlled interval pattern concept width. We propose a clustering algorithm based on the greedy approach and the procedure for finding an approximately optimal position of a box described in Section 3. The algorithm is not sensitive to the order in which points P are given. The parameters of the algorithm include positive real numbers [[delta].sub.1], [[delta].sub.2], ..., [[delta].sub.d] and all parameters of the positioning algorithm, namely, s, [s.sub.min], [lambda], and f(n, d).

First, we put [P.sub.1] = P and find an approximately optimal position [D.sub.1] of the box with the edge lengths [[delta].sub.1],..., [[delta].sub.d] for the set [P.sub.1]. Now suppose that the sets [D.sub.1],..., [D.sub.i] and [P.sub.1],..., [P.sub.i] have been already constructed and let [P.sub.i+1] = [P.sub.i]\[D.sub.i]. If [P.sub.i+1] = 0 then the procedure is ended. Else we find an approximately optimal position [D.sub.i+1] of the box for the set [P.sub.i+1]. The output of this procedure is a set of clusters [C.sub.i] = [P.sub.i] [intersection] [D.sub.i].

In order to avoid producing a lot of small clusters consisting of outliers we impose one more restriction.

Restriction 3. The resulting clusters must include at least [c.sub.min] objects.

With this restriction if the size of [P.sub.i+1] [intersection] [D.sub.i+1] is less than [c.sub.min] then the procedure ends (and points belonging to P \ ([C.sub.1] [union] *** [union] [C.sub.i]) are considered unclustered and referred to as outliers).

Restriction 3 together with Theorem 12 immediately leads to the following theorem.

Theorem 13. The clustering algorithm has

O ((dn log (n) + [d.sup.3][n.sup.1-1/d]/[s.sub.min] f(n,d)) x n/[c.sub.min]) (13)

worst-case time complexity and O(dn) space complexity.

If Restrictions 1-3 are omitted, the worst-case time complexity estimate is

0((dn log (n) + [d.sup.2][n.sup.1-1/d]/[s.sub.min] f(n,d)). (14)

5. Validation

Validation of the clustering algorithm developed in this study was performed on a dataset of tactile images registered by the Medical Tactile Endosurgical Complex (MTEC) during examination of artificial samples. MTEC allows intraoperative mechanoreceptor tactile examination of tissues and is already used in endoscopic surgery [14-16]. As methods for automated analysis of medical tactile images are still insufficient, validation results in particular and the developed clustering algorithm in general provide new opportunities for the medical domain applications.

The key component of MTEC is a tactile mechanoreceptor [17, Fig. 1]. Its operating head is equipped with 19 pressure sensors that perform synchronous measurements 100 times per second. Each measurement result (called "a tactile frame" and consisting of 19 values) is wirelessly transmitted to a computer that performs preprocessing and visualization. The sensors are located at the operating head surface which is a circle with diameter 20 mm.

In order to create a dataset of tactile images we utilized MTEC for tactile examinations of three types of artificial samples. The samples were similar to the L-samples utilized in the study --they were made using a soft silicone (Ecoflex 00-10, Shore hardness 00-10A) according to manufacturer's instructions and had a shape of a rectangular block with length, width, and height of 40 mm, 35 mm, and 11 mm, respectively. The difference was in sizes and shapes of hard inclusions enclosed in the samples. For the first sample type (ST1) the inclusion had a form of a spherical cap with base diameter 8 mm and height 2.4 mm oriented for palpation from the convex side. For the second sample type (ST2) the inclusion had a form of a spherical cap with a base diameter 4.7 mm and height 1.7 mm also oriented for palpation from the convex side. For the third sample type (ST3) the inclusions were the same as for ST2, but they were oriented for palpation from the flat side. For all sample types the inclusions were located in the center at height of approximately 3 mm. Thus, sample types were similar with a difference in either size or convexity of the inclusion. These samples simulated tissue with malignant neoplasms.

Totally 55 tactile examinations of the described samples were performed using MTEC. The contact angle was kept approximately equal to 90[degrees], and inclusions were located close to the center of the operating head surface. We performed twenty-two, seventeen, and sixteen examinations for samples of ST1, ST2, and ST3 types, respectively. For each examination one tactile frame was selected, namely, the one with the largest standard deviation (SD) of values, and other tactile frames were disregarded. Visualization of tactile frames for each sample type is presented in Figures 2(a)-2(c).

Thus, each examination was associated with a point in R19, and the total number of points was 55. This set of points was clustered using the developed clustering algorithm, and the results were compared with the results of k-means clustering (k = 3, Euclidean distance; see, e.g., ), which was used as a reference. Scikit-learn implementation  of k-means algorithm was utilized. Adjusted and raw Rand indexes (clustering result versus original classes; see, e.g., ) were used as compared characteristics of the clustering quality. Note that both clustering algorithms use random initialization, so multiple runs were performed for clustering quality estimation (specifically, 100 runs were performed to estimate Rand index for each algorithm with given parameters).

The results produced by both the proposed algorithm and by the k-means algorithm were unsatisfactory. However, the poor quality of the resulting clustering was predictable as examining of a single sample can result in tactile frames that are essentially different with respect to representation by a point in [R.sup.19] due to rotation and slight shifts of a tactile mechanoreceptor.

To get better results we mapped the data to the new 9dimensional space of attributes. The new attributes included

(i) SD of all values in a tactile frame;

(ii) mean and SD of the values corresponding to 7 middle sensors;

(iii) mean and SD of the values corresponding to 12 outer sensors;

(iv) mean and SD of the values corresponding to sensors that belong to the main diagonals (3 diagonals each consisting of 5 sensors, 13 sensors in total; see Figure 2(d) for details);

(v) mean and SD of the values corresponding to sensors that belong to the secondary diagonals (6 diagonals each consisting of 4 sensors, 12 sensors in total; see Figure 2(d) for details).

These attributes are robust to rotations proportional to 60[degrees]. The values of mean and SD were computed after scaling the values to [0, 1] segment.

Transition to the new attribute space essentially improved the clustering quality, but our algorithm left 10-14 points as outliers ([c.sub.min] was set equal to 8; the values of s, [s.sub.min], and A were set equal to 0.9, 0.3, and 0.8 respectively, and [sigma] was set equal to 0.27 for all attributes). A representative result of one run is presented in Table 2. Then we placed outliers points to the obtained clusters by the fc-nearest neighbors algorithm (fc = 8, unweighted; see, e.g., ). A representative result for one run is presented in Table 3.

Table 4 contains mean values and SDs for Rand indices and timing information.

As one can see, the proposed algorithm has an acceptable running time, and both our and k-means algorithm reach mean quality plateau already at 20 iterations.

The advantage of the proposed algorithm over the k-means algorithm with respect to the clustering quality was statistically significant. For example, for 20 iterations and adjusted Rand index the comparison of our algorithm with outliers and the k-means on 100 runs resulted in Mann-Whitney U-test two-tailed p value equal to 1.0 x [10.sup.-10]. As outliers are the points that are the most difficult for clustering, the advantage of our algorithm complemented by kNN-attributing of outliers to clusters over the k-means was lower but still firmly significant with Mann-Whitney U-test two-tailed p value equal to 9.5 x [10.sup.-4].

Interestingly, the transition to the new attribute space improved the quality of our algorithm more than the quality of the k-means clustering. For example, for 20 iterations, adjusted Rand index, and 100 runs, the comparison of the clustering quality for the initial attribute space and the new attribute space resulted in Mann-Whitney U-test two-tailed p values not exceeding [10.sup.-12] for both "with outliers" and "no outliers" versions of our algorithm, while for k-means the p value was 0.43.

6. Conclusions

In this paper we proposed a greedy clustering algorithm based on interval pattern concepts. The obtained theoretical estimate on algorithm complexity proved computational feasibility for high-dimensional spaces, and the validation on experimental data demonstrated high quality of the resulting clustering in comparison with conventional clustering algorithms such as k-means.

Particular results obtained during validation, such as a new attribute space for tactile frames registered by the Medical Tactile Endosurgical Complex, have individual significance as they provide new opportunities for the medical domain applications aimed at automated analysis of tactile images.

Data Access

Dataset of tactile frames used for the validation and the Python script that implements the developed clustering algorithm are available upon request from the authors.

https://doi.org/10.1155/2017/4323590

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The authors thank Dr. Alexey V. Galatenko and Dr. Vladimir V. Galatenko for valuable comments and discussions. The research was supported by the Russian Science Foundation (Project 16-11-00058 "The Development of Methods and Algorithms for Automated Analysis of Medical Tactile Information and Classification of Tactile Images").

References

 B. S. Everitt, S. Landau, M. Leese, and D. Stahl, Cluster analysis, Wiley Series in Probability and Statistics, John Wiley & Sons, Chichester, UK, Fifth edition, 2011.

 B. Ganter and R. Wille, Formal Concept Analysis: Mathematical Foundations, Springer, Berlin, Germany, 1999.

 B. Ganter and S. O. Kuznetsov, "Pattern Structures and Their Projections," in Conceptual Structures: Broadening the Base, vol. 2120 of Lecture Notes in Computer Science, pp. 129-142, Springer Berlin Heidelberg, Berlin, Heidelberg, 2001.

 M. Kaytoue, S. Duplessis, S. O. Kuznetsov, and A. Napoli, "Two fca-based methods for mining gene expression data," Lecture Notes in Computer Science (including subseries Lecture Notes in Artificial Intelligence and Lecture Notes in Bioinformatics), vol. 5548, pp. 251-266, 2009.

 M. Kaytoue, S. O. Kuznetsov, A. Napoli, and S. Duplessis, "Mining gene expression data with pattern structures in formal concept analysis," Information Sciences. An International Journal, vol. 181, no. 10, pp. 1989-2001, 2011.

 B. Aronov and S. Har-Peled, "On approximating the depth and related problems," SIAM Journal on Computing, vol. 38, no. 3, pp. 899-921, 2008.

 B. M. Chazelle and D. T. Lee, "On a circle placement problem," Computing. Archives for Scientific Computing, vol. 36, no. 1-2, pp. 1-16, 1986.

 C. M. de Figueiredo and G. D. da Fonseca, "Enclosing weighted points with an almost-unit ball," Information Processing Letters, vol. 109, no. 21-22, pp. 1216-1221, 2009.

 J. Matousek, "Range searching with efficient hierarchical cuttings," Discrete & Computational Geometry, vol. 10, no. 1, pp. 157-182, 1993.

 J. Eckstein, P. L. Hammer, Y. Liu, M. Nediak, and B. Simeone, "The maximum box problem and its application to data analysis," Computational Optimization and Applications. An International Journal, vol. 23, no. 3, pp. 285-298, 2002.

 J. L. Bentley, "Multidimensional binary search trees used for associative searching," Communications of the ACM, vol. 18, no. 9, pp. 509-517, 1975.

 J. H. Freidman, L. B. Jon, and A. F. Raphael, "An algorithm for finding best matches in logarithmic expected time," ACM Transactions on Mathematical Software, vol. 3, no. 3, pp. 209-226, 1977.

 D. T. Lee and C. K. Wong, "Worst-case analysis for region and partial region searches in multidimensional binary search trees and balanced quad trees," Acta Informatica, vol. 9, no. 1, pp. 23-29, 1977.

 V. Sadovnichy, R. Gabidullina, M. Sokolov, V. Galatenko, V. Budanov, and E. Nakashidze, "Haptic device in endoscopy," in Proceedings of the 21st Medicine Meets Virtual Reality Conference, NextMed/MMVR 2014, pp. 365-368, USA, February 2014.

 V. Barmin, V. Sadovnichy, M. Sokolov, O. Pikin, and A. Amiraliev, "An original device for intraoperative detection of small indeterminate nodules," European Journal of Cardio-thoracic Surgery, vol. 46, no. 6, pp. 1027-1031, 2014.

 R. F. Solodova, V. V. Galatenko, E. R. Nakashidze et al., "Instrumental tactile diagnostics in robot-assisted surgery," Medical Devices: Evidence and Research, vol. 9, pp. 377-382, 2016.

 V. M. Staroverov, V. V. Galatenko, T. V. Zykova, Y. I. Rakhmatulin, D. V. Rukhovich, and V. E. Podol'skii, "Automated real-time correction of intraoperative medical tactile images: sensitivity adjustment and suppression of contact angle artifact," Applied Mathematical Sciences, vol. 10, pp. 2831-2842, 2016.

 F. Pedregosa, G. Varoquaux, and A. Gramfort, "Scikit-learn: machine learning in Python," Journal of Machine Learning Research, vol. 12, pp. 2825-2830, 2011.

 T. Mitchell, Machine Learning, McGraw Hill, New York, NY, USA, 1997.

Stepan A. Nersisyan, Vera V. Pankratieva, Vladimir M. Staroverov, and Vladimir E. Podolskii

Faculty of Mechanics and Mathematics, Lomonosov Moscow State University, Leninskie Gory 1, Moscow 119991, Russia

Correspondence should be addressed to Stepan A. Nersisyan; s.a.nersisyan@gmail.com

Received 19 June 2017; Accepted 15 August 2017; Published 25 September 2017

Caption: FIGURE 1: The base cube is colored red; the global optimum is blue. There is no way to move from the red cube to the blue one without losing touch with the base cube.

Caption: FIGURE 2: (a-c) Examples of tactile frames for examinations of ST1 (a), ST2 (b), and ST3 (c) type samples. Pressure values are scaled to [0,255] segment and color-coded. (d) Correspondence between sensors and attributes from the new attribute space. Each hexagon represents one sensor. Middle sensors are colored in light-gray, outer sensors are colored in dark-gray. The main diagonals are shown by orange lines; the secondary diagonals are shown by blue lines. Centers of the hexagons that represent sensors belonging to both main and secondary diagonals are colored in red; belonging only to main diagonals, in orange; belonging only to secondary diagonals, in blue.
```Table 1: A fuzzy formal context, where the objects are pupils
and the attributes are disciplines.

Arts   Mathematics   Computer   Sports
science

A    9          9           10        9
B    8          2           6         5
C    6          5           10        7
D    8          9           9         6
E    8          4           6         9
F    6          5           2         10

Table 2: Correspondence between the original classes and the
clusters constructed by the proposed algorithm (with outliers).

1st cluster   2nd cluster   3rd cluster   Unclustered
9 points      13 points     22 points     11 points

ST1 22      9 points      1 points      5 points      7 points
points
ST2 17      0 points      12 points     3 points      2 points
points
ST3 16      0 points      0 points      14 points     2 points
points

Table 3: Correspondence between the original classes and the
clusters constructed by the proposed algorithm (no outliers).

1st cluster   2nd cluster   3rd cluster
11 points     17 points     27 points

ST1 22      11 points     3 points      8 points
points
ST2 17      0 points      14 points     3 points
points
ST3 16      0 points      0 points      16 points
points

Table 4: Dependency of Rand index values and the running time for
our and k-means clustering methods on number of iterations
performed (100 program runs for each value). Values of Rand index
are presented in terms of medians and interquartile ranges (IQR).

Number of      Clustering      Rand index    Rand index    Average
iterations       method          median         IQR        running
/raw)         /raw)      seconds)

20             Our method       0.43/0.73     0.12/0.05      0.8
(with outliers)
Our method       0.39/0.73     0.10/0.04      0.8
(no outliers)
k-means        0.32/0.70     0.21/0.09      0.02
50             Our method       0.43/0.73     0.08/0.05      2.4
(with outliers)
Our method       0.39/0.73     0.06/0.03      2.5
(no outliers)
k-means        0.27/0.68     0.20/0.09      0.05

100            Our method       0.42/0.74     0.08/0.04      4.2
(with outliers)
Our method       0.39/0.73     0.06/0.03      4.3
(no outliers)
k-means        0.31/0.70     0.20/0.09      0.09
```
Title Annotation: Printer friendly Cite/link Email Feedback Research Article Nersisyan, Stepan A.; Pankratieva, Vera V.; Staroverov, Vladimir M.; Podolskii, Vladimir E. Journal of Applied Mathematics Report Jan 1, 2017 6529 Analysis of Economic Burden of Seasonal Influenza: An Actuarial Based Conceptual Model. Sparse Optimization of Vibration Signal by ADMM. Algorithms Clustering (Computers) Convergence (Mathematics) Information management Mathematical research