A Comparative Study on TIBA ImputationMethods in FCMdd-Based Linear Clustering with Relational Data

Relational fuzzy clustering has been developed for extracting intrinsic cluster structures of relational data and was extended to a linear fuzzy clustering model based on Fuzzy c-Medoids (FCMdd) concept, in which Fuzzy c-Means-(FCM-) like iterative algorithm was performed by defining linear cluster prototypes using two representative medoids for each line prototype. In this paper, the FCMdd-type linear clustering model is further modified in order to handle incomplete data including missing values, and the applicability of several imputation methods is compared. In several numerical experiments, it is demonstrated that some pre-imputation strategies contribute to properly selecting representative medoids of each cluster.


Introduction
Relational fuzzy clustering is a relational extension of fuzzy clustering for revealing cluster structures buried in relational data.Relational Fuzzy c-Means (RFCM) [1] extended the Fuzzy c-Means (FCM) [2] clustering criterion with mutual dissimilarity measures instead of object-type observation in FCM.Although FCM and other variants of k-Means [3] use the clustering criterion of the distance between a data point and a cluster prototype, RFCM defines the clustering criterion by using mutual dissimilarities only.When the dissimilarities among objects are measured by squared Euclidean distances, the RFCM criterion is equivalent to the centroid-less formulation of the FCM criterion.Using other dissimilarity measures; however, the RFCM criterion has no clear connection with distances between data points and prototypes.In k-Medoids [4], cluster prototypes are selected from data points, and the clustering criterion coincides with one of mutual dissimilar degree among objects.So, k-Medoids can be directly extended to relational data analysis even if taking an average cannot be done in non-Euclidean space.Fuzzy c-Medoids (FCMdd) [5] is a fuzzy extension of k-Medoids and can deal with various dissimilarity measures.
Linear fuzzy clustering models [6,7] extract linear substructures by modifying the point prototypes of FCM into lines, planes, and linear varieties.Because the subspace learning model in each cluster can be identified with fuzzy principal component analysis (fuzzy PCA) [8], they are often regarded as a kind of local principal component analysis (local PCA) [9].This paper studies the FCMdd-based linear clustering model [10], which can reveal local linear substructures buried in relational data.In [10], Haga et al. defined each prototypical line by using two representative medoids and demonstrated that the clustering modal can be applied to Euclidean relational data.The FCMdd-type linear clustering model was further modified for dealing with non-Euclidean relational data [11,12], in which data transformation, called β-spread transformation, was performed before applying the clustering algorithm in a similar manner to Non-Euclideantype Relational Fuzzy (NERF) c-Means [13].
In this paper, a comparative study on the applicability of β-spread transformation is performed in FCMdd-based linear clustering of incomplete relational data.Hathaway and Bezdek [14] proposed several methods for imputing (predicting and substituting) missing elements of incomplete relational data and showed that imputation errors can be revised by β-spread transformation in NERF c-Means.This paper demonstrates that the performance of FCMdd-type linear fuzzy clustering for incomplete relational data can also 2 Advances in Fuzzy Systems be improved by β-spread transformation through several comparative experiments including an example of document clustering.
The remaining part of this paper is organized as follows.In Section 2, linear clustering and relational clustering are briefly reviewed.Section 3 introduces FCMdd-type linear clustering model and applies several imputation methods called TIBA.Comparative results are shown in Section 4, and conclusions are given in Section 5.

Linear Clustering and
Relational Clustering u ci is the fuzzy membership degree of pattern i to cluster c, and θ is the fuzzification parameter.The larger the θ, the fuzzier the membership assignment.D ci is the clustering criterion which measures the deviation between pattern i and the prototype of cluster c.In the original FCM clustering [2], cluster prototypes are given by the centroid vectors b c , and D ci is the squared Euclidean distance as follows: The FCM model is reduced to the hard (nonfuzzy) k-Means model [3] when θ = 1, in which cluster memberships are given by the nearest prototype principle.Besides point-type prototypes b c in FCM, Fuzzy c-Lines (FCL) [6] for extracting linear clusters used linear prototypes defined as where a c is the basis vector of the principal subspace, and b c is the centroid, which the linear prototype passes through.The clustering criterion is calculated as The updating rules for membership u ci and the cluster center b c are derived as The basis vectors a c are the principal eigenvectors of the generalized fuzzy scatter matrices: This linear clustering model has close relation with local PCA [9].Indeed, when we consider only a single cluster (C = 1), the FCL clustering model is equivalent to the conventional PCA and the basis vector a c is reduced to the principal component vector.In this sense, FCL is a type of local PCA, which simultaneously performs membership estimation (local fuzzy group extraction) and fuzzy PCA [8] in each local fuzzy group considering the fuzzy membership degree of u θ ci .The prototypical line Line c can be identified with principal subspace spanned by fuzzy principal component vector a c from the local PCA view point.
When θ = 1, the FCL model is also reduced to the hard (nonfuzzy) local PCA model [15,16], in which cluster memberships are given by the nearest prototype principle.[1] is the relational extension of FCM.When we have relational data composed of mutual relations among patterns D = {d 2 i j }, the FCM-type objective function is redefined as

FCM-Type Relational Clustering. RFCM
d i j can be any type of dissimilarity between patterns i and j but is assumed to be Euclidean-like one in RFCM.Indeed, this model is equivalent to FCM only when d 2 i j is the squared Euclidean distance, and the clustering model derives only poor results if the relational information is highly non-Euclidean.
In order to modify RFCM for handling non-Euclidean distance metrics, Hathaway and Bezdek [13] considered NERF, which includes the following β-spread transformation: where β is added to off-diagonal elements of non-Euclidean relational data D. I is a unit matrix, and β is a suitably chosen scalar.M is a matrix whose elements are all 1.Hathaway and Bezdek discussed that D β is Euclidean if PD β P with P = I − (1/n)M is negative semidefinite; that is, β is greater than or equal to the largest eigenvalue of PDP.By the way, the basic RFCM iteration can be continued when clustering criteria are all nonnegative.In NERF, β is gradually increased from 0 to a certain value by considering the negative elements of clustering criteria.

FCMdd-Type Linear Clustering and TIBA Imputation
3.1.FCMdd-Type Linear Clustering.Assume that d i j is the mutual Euclidean distance such that FCMdd [5] is a fuzzy extension of k-medoids [4], which performs an FCM-like clustering by selecting b c from patterns x i , i = 1, . . ., n.The representative objects are called "medoids" and are given by solving combinatorial optimization problems.Haga et al. [10] applied the idea to linear fuzzy clustering, in which each linear prototype is spanned by two representative medoids x c1 and x c2 as The squared Euclidean distance between object i and the prototypical line Line c is given as With fixed fuzzy memberships u ci , the optimal medoids are derived by the following combinatorial optimization problem: The optimal medoid set of (c 1 , c 2 ) is searched by enumerating all pairs of objects.In order to reduce the computational cost, a simplified medoid search process was also proposed, in which medoids are selected from a subset X c of objects: where This linear fuzzy clustering model was also extended to the 2D prototype case by spanning 2D prototypical planes using three medoids [10].
Although non-Euclidean relational data may bring negative values for the clustering criteria of (12), from the practical view point, we have no trouble in operating the conventional FCMdd-type linear clustering algorithm if all clustering criteria are not negative.
Yamamoto et al. [11] proposed a procedure for β-spread transformation so as to avoid negative criterion values in FCMdd-type linear clustering.Because a negative criterion value implies a non-Euclidean situation, relational data should be revised so that the criterion value is always nonnegative.In the previous research [12], it was shown that the clustering criterion ) is satisfied.Then, β-spread transformation should be performed so that the following triangle inequality is satisfied for all objects: A plausible value of Δβ in an iteration step is obtained as Here, Δβ is positive when some D ci are negative, while Δβ is zero when all D ci are nonnegative.Then, β is monotonically increasing.
A sample procedure including the automated β-spread transformation can be summarized as follows: Step 1. Set β = 0. Randomly initialize the prototypical medoids (two representative objects) of each cluster.
Step 2. Calculate the clustering criteria D ci by (12).
Step 5. Search medoids in each cluster.
Step 6. Repeat Steps 2-5 until a certain stopping criterion is satisfied.
In Step 6, such a stopping criterion as min |u (new Although the proposed model is in the fuzzy clustering category, it is easily seen that a hard (non-fuzzy) version can be covered when θ = 1, in which cluster memberships are given by the nearest prototype principle.

Missing Value Imputation by TIBA. Hathaway and
Bezdek [14] demonstrated that the β-spread transformation is also useful for handling missing elements in relational data matrices.Although preimputation of missing elements may cause imputation errors and bring illegal effects in clustering process, β-spread transformation can decrease the illegal effects.
This paper considers the applicability of several imputation techniques in FCMdd-type linear clustering.
Hathaway and Bezdek [14] used three imputation techniques based on triangle inequality-based approximation (TIBA).The triangle inequality, which Euclidean relational data always satisfy, is represented as follows: Advances in Fuzzy Systems Assume that an element of relational matrix d i j is missing and is to be preimputed before applying the clustering algorithm.Let K i j be the corresponding index set as For each k ∈ K i j , the triangle inequality ( 17) is given as the upper bound of d i j .Missing elements are replaced with the minimum upper bound of d i j : which is called minimax TIBA.By the way, d i j is imputed by zero value if K i j is empty.
Advances in Fuzzy Systems The triangle inequality is also represented as follows: and brings the following inequalities: (21) So, the lower bound of d i j is given as Missing elements are replaced with the maximum lower bound of d i j : which is called maximin TIBA.
It is also possible to combine the previous two imputation values for predicting a reasonable estimation of missing values.The average values of minimax TIBA and maximin TIBA are used for imputing missing elements.This TIBA is called average TIBA.
These imputation techniques based on triangle inequalities can be easily applied to relational clustering problems.In the next section, these three imputation approaches are compared in FCMdd-type linear clustering tasks in conjunction with β-spread transformation.

Numerical Experiments
Two experimental results are shown in order to consider the applicability of the three TIBA imputation techniques in FCMdd-type linear clustering with β-spread transformation.
In previous researches, it has been shown that "soft" clustering models outperformed "hard" ones in local PCA tasks [15][16][17], and "fuzzy" models can be more useful than probabilistic ones [9].Therefore, in this paper, the characteristics of the fuzzy version are investigated.

Artificial Data
Set.An artificial relational data set composed of 60 patterns was generated from a 2D data set shown in Figure 1, in which patterns form two lineshaped clusters.It is obvious that the local linear structures cannot be extracted by the conventional point-prototype models such as FCM-like models and FCMdd.We made two relational data matrices.The first relational data matrix was generated by Euclidean norm, and the second one was generated by L 1 norm, which is non-Euclidean measure.The iterative algorithm was performed until the medoids became unchanged, and the model parameters were set as (C, θ) = (2,2).In order to demonstrate the characteristics of the algorithm, the initial memberships were given in a supervised manner; that is, (u 1i , u 2i ) = (0.9, 0.1) for the first visual cluster and (u 1i , u 2i ) = (0.1, 0.9) for the second one.
In the previous research [12], it was demonstrated that the two linear substructures can be successfully revealed by the FCMdd-based linear clustering algorithm without βspread transformation for Euclidean relational data while it can be done only with β-spread transformation for L 1 norm.
First, Euclidean incomplete relational data matrices were generated by removing a part of off-diagonal elements where K i j was not empty.In order to protect tridiagonal parts of relational data, the maximum number of missing elements was set as n 2 − 3n + 2.
Clustering results are compared with those without βspread transformation in Figure 2. Objects were partitioned into two clusters of circles and times, and smaller times mean that the patterns were shared almost equally by the two clusters.Medoids and prototypical lines are indicated by black circles and lines, respectively.
Each approximation method with β-spread transformation could estimate cluster medoids for capturing the two visual linear prototypes until the numbers of missing elements are less than about 30% although patterns having ambiguous memberships increased more than complete relational data.β-spread transformation performed on each approximation, minimax TIBA: β = 0.04867, maximin TIBA: β = 0.044286, average TIBA: β = 0.051706.Here, the maximum eigenvalues of PDP after imputation were, minimax TIBA: 0.053808, maximin TIBA: 0.1711612, average TIBA: 0.083694.So, the TIBA imputation brought a slightly non-Euclidean situation, and β-spread transformation successfully modified the data set.
On the other hand, without β-spread transformation, only average TIBA made it possible to extract linear substructures while minimax TIBA and maximin TIBA brought inappropriate results where some patterns depicted by black diamonds in Figure 2 had negative clustering criterion values.
These results imply that the FCMdd-type linear clustering can successfully extract linear substructures of incomplete Euclidean relational data using β-spread transformation although the three imputation techniques cause non-Euclidean relational matrices.
Second, FCMdd-type linear fuzzy clustering was applied to non-Euclidean relational data.Incomplete relational data matrices were generated in the same manner with the Euclidean case.Clustering results are depicted in Figure 3.
Without β-spread transformation; however, all the three TIBAs brought inappropriate partitions because many patterns had negative clustering criterion values.
In this way, β-spread transformation also works well in incomplete situations.

Document Clustering.
In the second experiment, TIBA imputation methods are compared in a document classification task.A relational data set was generated using a famous Japanese novel "Kokoro" by Soseki Natsume.The novel is composed of 3 chapters (Sensei and I, My Parents and I, Sensei and His Testament), and the chapters include 36, 18, 56 sections, respectively.The text data (Japanese language) can be downloaded from Aozora Bunko (http://www.aozora.gr.jp/).The sections were used as individual text documents (n = 110), which should be partitioned without the chapter information.The text documents were preprocessed using "Chasen" morphological analysis system software (http://chasen.naist.jp/hiki/ChaSen/),which segments the Japanese text string into morphemes.Wada et al. [18] performed a PCA-based structural analysis with the 83 most frequently used substantives and verbs with their tfidf weights and revealed that the chapter structure can be emphasized by using 10 meaningful keywords as is shown in Figure 4, which is 2D biplots of principal components.Chapters 2 and 3 form two linear clusters in 10D data space, and chapter 1 exists on their intersection.In this experiment, parameters were set as C = 2, θ = 2.0 with the goal of revealing the two linear substructures.Two relational data matrices were generated considering co-occurrence information of the 10 keywords.Jaccard coefficient and Dice coefficient are the similarity measures for  Jaccard's coefficient is the similarity represented as

Advances in Fuzzy Systems
Dice's coefficient is also the similarity represented as Because the linear clustering model uses distance (dissimilarity) measures, the similarity measures s i j were transformed into dissimilarity ones d i j .
Before applying the FCMdd-based linear fuzzy clustering, randomly selected elements were withheld from the relational matrix with 11,772 elements and were imputed by the three TIBA methods.Then, the cluster partitions for Jaccard's index were derived as shown in Figure 5. Two clusters are depicted by circles and times, and small times mean ambiguous assignment.Documents were properly partitioned into two clusters considering linear substructures.
Clustering results for Dice coefficient are depicted in Figure 6.Our approach also extracted linear substructure from incomplete relational data of Dice coefficient.Minimax TIBA allowed with 48% missing values or fewer.Average TIBA tolerated 55% missing values or fewer.Maximin TIBA resulted in a good partition with 63% missing values or fewer.The parameters β in β-spread transformation were given as, minimax TIBA: β = 0.760, maximin TIBA: β = 0.383, average TIBA: β = 0.608.The derived β values are still smaller than the maximum eigenvalues of PDP, minimax TIBA: 6.169, maximin TIBA: 4.288, average TIBA: 2.219.
In the experiments, it was demonstrated that the TIBA imputation methods work well for incomplete non-Euclidean relational data in conjunction with β-spread transformation.
Finally, comparison with other methods is discussed.Although we have already many clustering algorithms, some of which are used in document clustering tasks [20], most of them are designed for finding groups composed of similar pattern from the view point of "point prototype" or "hierarchical aggregation".For example, Fuzzy c-Medoids (FCMdd) [5], which is a representative method of pointprototype models, can be applied to the relational data set of this subsection.Figure 7 shows the clustering results of finding three chapter structures of circles, times, and triangles.Small times mean ambiguous assignment as well.The conventional clustering methods are useful for finding such document groups considering mutual similarity among documents (or sometime keyword groups).
On the other hand, the proposed method is designed for a different purpose of finding "local linear structures" from the view point of local PCA, which is useful for cluster-wise information summarization such as local feature map construction.In this sense, the proposed method has different future application area from the conventional clustering tools.

Conclusion
This paper compared the applicability of TIBA imputation methods and β-spread transformation for handling incomplete relational data in FCMdd-type linear clustering.In numerical experiments, three imputation techniques of minimax TIBA, maximin TIBA, and average TIBA were compared using two data sets.The experimental results indicated that β-spread transformation still works well for incomplete data in conjunction with β-spread transformation.All the three TIBAs are useful for imputing incomplete non-Euclidean relational data.
From the view point of local PCA concept, the proposed method can be used for local information summarization or local feature map construction where data structures are visually summarized in low-dimensional space in conjunction with data clustering.The application is remained in future works.Another potential future work is an extension to the case of multidimensional prototype models, which is useful for constructing 2D feature map.

Figure 5 :
Figure 5: Comparison of cluster partition from incomplete "Kokoro" text data derived with Jaccard's coefficient imputed by three TIBAs.

Figure 6 :
Figure 6: Comparison of cluster partition from incomplete "Kokoro" text data derived with Dice's coefficient imputed by three TIBAs.

Figure 7 :
Figure 7: Comparison of cluster partition of Fuzzy c-Medoids for incomplete "Kokoro" text data derived with Jaccard's and Dice's coefficients imputed by maximin TIBA.
2.1.FCM-Type Linear Clustering.Assume that we have m-dimensional observations of n patterns x i = (x i1 , . . ., x im ) , i = 1, . . ., n.With the goal of partitioning the n patterns into C clusters, the objective function for FCM-type clustering is defined as