# Data and Feature Reduction in Fuzzy Modeling through Particle Swarm Optimization

Authors: S. Sakinah S. Ahmad; Witold Pedrycz
Journal: Applied Computational Intelligence and Soft Computing
 (2012)

---

## Abstract

The study is concerned with data and feature reduction in fuzzy modeling. As these reduction activities are advantageous to fuzzy models in terms of both the effectiveness of their construction and the interpretation of the resulting models, their realization deserves particular attention. The formation of a subset of meaningful features and a subset of essential instances is discussed in the context of fuzzy-rule-based models. In contrast to the existing studies, which are focused predominantly on feature selection (namely, a reduction of the input space), a position advocated here is that a reduction has to involve both data and features to become efficient to the design of fuzzy model. The reduction problem is combinatorial in its nature and, as such, calls for the use of advanced optimization techniques. In this study, we use a technique of particle swarm optimization (PSO) as an optimization vehicle of forming a subset of features and data (instances) to design a fuzzy model. Given the dimensionality of the problem (as the search space involves both features and instances), we discuss a cooperative version of the PSO along with a clustering mechanism of forming a partition of the overall search space. Finally, a series of numeric experiments using several machine learning data sets is presented.


---

## Body

## 1. Introduction

In fuzzy modeling, the two main approaches for generating the rules rely on knowledge acquisition from human experts and knowledge discovery from data [1, 2]. In recent years, knowledge discovery from data or data-driven fuzzy modeling has become more important [2–4]. In many cases, the ability to develop models efficiently is hampered by the dimensionality of the input space as well as the number of data. If we are concerned with rule-based models, the high-dimensionality of the feature space along with the topology of the rules gives rise to the curse of dimensionality [1, 4]. The number of rules increases exponentially and is equal to Pn, where n is the number of features and P stands for the number of fuzzy sets defined for each feature.The factors that contribute most to the accuracy of the data-driven fuzzy modeling are associated with the size of the input space and the decomposition of the input data. A Large number of data points or instances in a continuous input-output domain exhibit a significant impact on fuzzy models. It is well known that more training data will not always lead to a better performance for data-driven models. Large amount of training data have important implications on the modeling capabilities. Since the number of fuzzy sets determines the family of realizable approximation functions, larger datasets present the possibility of over-fitting the training data [1, 4]. Thus, the effectiveness of the fuzzy models relies on the quality of the training data. In addition, the main drawback is the fuzzy models’ relative inefficiency as the size of the data increases, regarding both the number of data points in the data set and the number of features. Moreover, one of the most widely used approaches in fuzzy modeling is the fuzzy C-means (FCM) algorithm for constructing the antecedents of the rules associated with the curse of dimensionality [5, 6].The dimensionality problem can be addressed by reducing the constructed fuzzy rules. The reduction method plays two important roles: it increases the effectiveness of the learning algorithm, since the learning algorithm will concentrate only on the most useful subset of data, and it also improves the computational efficiency as the learning algorithm involves only a subset of data smaller than the original dataset [7]. This reduction can be realized by removing the redundant fuzzy rules by exploiting a concept of fuzzy similarity [3, 7, 8]. Evolutionary algorithms have also been used for building compact fuzzy rules [9–12]. An evolutionary algorithm is used to tune the structure and the rules’ parameter of the fuzzy systems [13, 14]. However, in numerous cases, some variables are not crucial to the realization of the fuzzy model. A suitable way to overcome this problem is to implement feature selection before constructing the fuzzy models. Therefore, during the last decade, feature selection methods in conjunction with constructing fuzzy models for reducing the curse of dimensionality were developed [15–22]. This process reduces the fuzzy rule search space and increases the accuracy of the model.As mentioned above, forming the best input data as the training set to construct the fuzzy modeling is also important. However, as far as we know there is no research that has been done to simultaneously select the best subset of features and input data for constructing the fuzzy model. Most of the research is focused on reducing the fuzzy rules, and the process of simplifying the system is done once the design has been completed. Here we propose a method that reduces the complexity of the system starting from the design stage. however, the process of constructing the antecedent and the consequent parts of the fuzzy model is realized using the best subset of input data.In this paper, a comprehensive framework is proposed to construct fuzzy models from the subset of numerical input-output data. First, we develop a data-driven fuzzy modeling framework for a high-dimensional large dataset, which is capable of generating a rule-based automatically from numerical data. Second, we integrate the concept of feature selection and data selection together in the unified form to further refine (reduce) the fuzzy models. In this regard, the PSO technique is applied in order to search for the best subset of data. In order to increase the effectiveness of the PSO techniques, we introduce a new cooperative PSO method based on the information granulation approach. Third, we develop a flexible setup to cope with the optimization of variables and data to be used in the design of the fuzzy model. The proposed approach allows the user to choose the predetermined fraction of variables and data that can be used to construct the fuzzy models.This paper is organized as follows. We briefly elaborate on the selected approaches to data and feature space reduction in Section2, and then in Section 3, we recall the main algorithmic features of PSO and its cooperative version, CPSO, which is of interest in problems of high-dimensionality. The proposed fuzzy modeling framework along with its main algorithmic developments is presented in Section 4. Experimental studies are presented in Section 5, and conclusions are provided in Section 6.


## 2. Selected Approaches to Data and Space Reduction

In general, reduction processes involve feature selection (FS), instances (data) selection (IS), and a combination of these two reduction processes: feature and Instances selection (FIS). Feature selection is a subject of the main reduction pursuits. The goal of FS, which is commonly encountered in problems of system modeling and pattern recognition, is to select the best subset of features so that the model formed in this new feature (input) space exhibits the highest accuracy (classification rate) being simultaneously associated with the increased transparency of the resulting construct [23]. The process aims to discard irrelevant and/or redundant features [24]. In general, the FS algorithms can be classified into three main categories: filters, wrappers, and embedded methods. The filter method selection criterion is independent of the learning algorithm. In contrast to the wrapper method, the selection criterion is dependent on the learning algorithm and uses its performance index as the evaluation criterion. The embedded method incorporates feature selection as part of the training process. The reader can refer to [23–25] for more details.Instances selection (IS), another category of reduction approaches, is concerned with the selection of the relevant data (instances) reflective of the knowledge pertinent to the problem at hand [26, 27]. The three main functions forming the essence of IS include enabling, focusing and cleaning [26].In this study, as stated earlier, instead of approaching feature selection and instances selection separately, we focus on the integration of feature selection and instances selection in the construction of the fuzzy models. Both processes are applied simultaneously to the initial dataset, in order to obtain a suitable subset of feature and data to construct the parameters for the fuzzy model. In the literature, some methods for integrating feature and instances selection are more focused on a class of classification problems [28, 29].The ideas of feature and data reduction as well as hybrid approaches have been discussed in the realm of fuzzy modeling. Table1 offers a snapshot at the diversity of the existing approaches and the advantages gained by completing the reduction processes.Table 1
A summary of selected studies in data and feature reduction in fuzzy modeling.
ReferenceFeature reduction techniqueDataset, fuzzy model and dataOriginal data used in modelingNumber of selected featuresNumber of resulting rulesNumber of instancesNumber of featuresGaweda et al.  [15]The use of sensitivity analysis Determination of essential featuresBox-Jenkins gas furnace2961032Hadjili and Wertz [16]Deviation criterion (DC): to measure the change in fuzzy partition. Removal of features that do not significantly change the fuzzy partitionNonlinear systems in noisy environment250314Nonlinear dynamical system excited by a sinusoidal signal8001068Run-out cooling table in a hot strip mill100017512Zarandi et al. [18]Heuristic method to select featuresNonlinear System used in [3]50424Supplier chance management dataset300955Du and Zhang [19]Evolutionary optimizationBox-Jenkins gas furnace2961034MR damper identification500011610Ghazavi and Liao [20](1) Mutual correlation methods,(2) gene selection criteria(3) the relief algorithmWisconsin breast cancerPIMA Indian diabetesWelding flaw identification56976839930825333250 (3)125 (3) —Zhang et al. [21]Iterative search margin based algorithm (Simba)Wisconsin breast cancer699953Wine1781345Iris150433Ionosphere35134104


## 3. Particle Swarm Optimization and Its Cooperative Version

Population-based algorithms provide interesting solutions since any constructive method can be used to generate the initial population, and any local search technique can be used to improve each solution in the population [30]. In addition, population-based methods have the advantage of being able to combine good solutions in order to obtain potentially better ones. Most of the population-based algorithm approaches in FS and IS are based on GAs. Some recent studies [28, 29, 31] have employed population-based optimization techniques to carry out search for the best subset of variables and data for solving the application problems, but all of them were carried out to solve the classification problem. Therefore, in this study, we use population-based technique for selecting the best subset of feature and data for the regression problem. Here, we implement particle swarm optimization (PSO) techniques to intelligently search for the best subset of features and data (instances).PSO, developed by Kennedy and Eberhart, inspired by the collective behavior of birds or fish [32], is a population-based algorithm where each individual, referred to as a particle, represents a candidate solution. Each particle proceeds through the search space at a given velocity v that is dynamically modified according to its own experience and results in its local best (lb) performance. It is also affected by others particles flying experience resulting in the best value, global best (gb). The underlying expression for the update of the velocity in successive generations reads as follows:
(1)vi(t+1)=w⋅vi(t)+c1⋅r1,i(t)[lbi(t)-xi(t)]+c2⋅r2,i(t)[gbi(t)-xi(t)],xi(t+1)=xi(t)+vi(t+1),
where i=1,2,…,N+n (the search space is equal to the sum of the dimensionalities of the feature space and the size of the data). The inertia weight (w) is confined to the range [0, 1]; its values can decrease over time. The cognitive factor c1 and social factor c2 determine the relative impact coming from the particle’s own experience and the local best and global best. r1 and r2 are numbers drawn from a uniform distribution over the unit interval that brings some component of randomness to the search process.In this research, we employed the PSO-based method to handle two optimization tasks, namely, (1) selection of the optimal subset of features and (2) selection of the optimal subset of instances based on the concept of information granularity. In order to reduce the computational complexity of using the standard PSO, we employed cooperative PSO method to simultaneously solve the two optimization tasks. The motivation behind the use of cooperative PSO, as advocated in [33], is to deal effectively with the dimensionality of the search space, which becomes a serious concern when a large number of data with a large dimensionality are involved. This curse of dimensionality is a significant impediment negatively impacting the effectiveness of standard PSO. The essence of the cooperative version of PSO is essentially a parallel search for optimal subset of features and its optimal subset of instances. The cooperative strategy is achieved by dividing the candidate solution vector into components, called subswarm, where each subswarm represents a small part of the overall optimization processes. By doing this, we implement the concept of divide and conquer to solve the optimization problem, so that the process will become more efficient and fast.The mechanism of information sharing of CPSO is shown in Figure1. The cooperative search between one subswarm and other is achieved by sharing the information of the global best position (PGB) across all subswarm. Here the algorithm has the advantage of taking two steps forward because the candidate solution comes from the best position for all subswarm except only for the current subswarms being evaluated. Therefore, the algorithm will not spend too much time optimizing the features or instances that have little effect to the overall solution. The rate at which each swarm converges to the solution is significantly higher than the rate of convergence of the standard PSO.Figure 1
The schematic diagram of information sharing in CPSO.The essence of the cooperative version of PSO is to split the data into several groups so that each group is handled by a separate PSO. The main design question involves splitting the variables into groups. A sound guideline is to keep the related (associated) variables within the same group. Obviously, such relationships are not known in advance. Several possible methods are available for addressing this issue in more detail in the context of the problem at hand.(a)
As we are concerned with a collection of features and data (instances), a natural way to split the variables would be to form two groups (K=2), one for the features (n) and another one for the instances (N). This split would be legitimate if the dimensionality of both subsets was quite similar.(b)
In some situations, one of the subsets (either the data or the features) might be significantly larger than the other one. We often encounter a large number of data, but in some situations, a large number of features might be present (for instance, in microarray data analysis). This particular collection of data or features is then split intoK groups. Clustering such items is a viable algorithmic approach. Running K-means or fuzzy C-means produces clusters (group) of variables that are used in the individual PSO.(c)
In case both subsets are large, the clustering is realized both for the features and data, and the resulting, structure (partition) is used to run cooperative PSO.Algorithm1 presents the Cooperative PSO pseudocode implementing the optimization process [33]. Firstly, the PSO is divided into m subspaces, called subswarms. In our case the first subswarm represents the features search space and the rest are for instances search space. Pj(xi) refers to the position of particle i of subswarms j. The global best for each subswarm defined as Pj(GB), and the local best is defined as Pj(LBi). The cooperation between the subswarms employed in the function C(j,k), which returns m-dimensional vector formed by concatenating all the global best vector across all subswarms, except for the current position j. Here the jth component is called k and represent the position of any particle from subswarm Pj.Algorithm 1:Pseudocode for cooperative PSO.
Initializem one-dimensional PSO: Pj, j∈[1,…,m]CreateC(j,k)=[P1(BG),P2(BG),…,Pj-1(BG),k,Pj+1(BG),…,Pm(BG)]Whilestop criteria not met doforeach subswarm j∈[1,…,m]  dofor each particle i∈[1,…,s]  doif fitness(C(j,Pj(xi))) > fitness(C(j,Pj(LBi)))thenPj(LBi)=Pj(xi)iffitness(C(j,Pj(LBi))) > fitness(C(j,Pj(BG)))then  Pj(BG)=Pj(lbi)end forfor each Pj dovi,j(t+1)=w·vi,j(t)+c1·r1,i(t)[PLBi,j(t)-xi,j(t)]+c2·r2,i(t)[PBGi,j(t)-xi,j(t)]xi,j(t+1)=xi,j(t)+vi,j(t+1)


## 4. PSO-Integrated Feature and Data Reduction in Fuzzy-Rule-Based Models

As the problem of feature data reduction is inherently combinatorial nature, PSO provides an interesting and computationally viable optimization alternative. In the following subsections, we start with a general optimization setting and then discuss the PSO realization of the search process (here, a crucial design phase is a formation of the search space with a suitable encoding mechanism). Although the proposed methodology is of general nature, we concentrate on rule-based models, which are commonly present in fuzzy modeling, to help offer a detailed view of the overall design process.


### 4.1. An Overall Reduction Process

As is usual in system modeling, we consider a supervised learning scenario ion which we encounter in a finite set of training data(xk,tk), k=1,2,…,N. By stressing the nature of the data and their dimensionality, the data space along with n-dimensional feature vectors can be viewed as a Cartesian product of the data and features D×F. The essence of the reduction is to arrive at the Cartesian product of the reduced data and feature spaces, D'×F', where, D'∈D and F'∈F. The cardinality of the reduced spaces is equal to N′ and n′, where N′<N and n′<n.The overall scheme of the reduction process outlining a role of the PSO-guided reduction is illustrated in Figure2. The scheme can be divided into two important parts and can be described as follows. (a)
Reduction process via PSO: a reduction process tackles both feature reduction and data reduction simultaneously. PSO algorithm is used to search for the best feature and data for constructing the fuzzy model. Here, the size of the selected features (n′) and data (N′) is provided in advance by the user. After the PSO meets the maximum generation, the process is stopped, and the last best subset of features and data is the best subset of data for constructing the fuzzy model.(b)
Evaluation process: the fuzzy C-means algorithm is used to convert the numerical data into the information granules. Here, the information granularity process deals only with the subset of the data and features (D'×F'). Next, the consequent parameter a constructed from the fuzzy models is use to evaluate the performance of the selected data and features. At this stage we access the performance of the constructed fuzzy model in terms of their capability to fit the model by using the all instances in the original data set.Figure 2
The scheme of the proposed data and reduction for fuzzy modeling.As it becomes apparent, the original spaceD×F is reduced, and in this Cartesian product a fuzzy model, denoted by FM, is designed in the usual way (we elaborate on the form of the fuzzy model in the subsequent section). Its design is guided by a certain objective function Q expressed over all elements of original instances. The quality of the reduced space is assessed by quantifying the performance of the fuzzy model operating over the original, non-reduced space. The same performance index as used in the construction of the fuzzy model in the reduced space is used to describe the quality of the fuzzy model:
(2)Q=1N∑xk∈D×F′(FM(xk)-tk)2.
Note that the summation shown above is taken over all the elements forming the data space D. By taking another look at the overall reduction scheme, it is worth noting that the reduction is realized as in the wrapper mode, in which we use a fuzzy model to evaluate the quality of the reduction mechanism.


### 4.2. The PSO-Based Representation of the Search Space

The reduction of the data and feature spaces involves a selection of a subset of the data and a subset of the features. Therefore, the problem is combinatorial in its nature. PSO is used here to form a subset of integers which are indexes of the data or features to be used in the formation ofD'×F'. For instance, D' is represented as a set of indexes {i1,i2,…,iN′} being a subset of integers {1,2,…,N}. From the perspective of the PSO, the particle is formed as a string of real numbers in [0,1] of the length of N+n; effectively, the search space is a hypercube [0,1]N+n. The first substring of length N represents the data; the second one (having n entries) is used to optimize the subset of features. The particle is decoded as follows. Each substring is processed (decoded) separately. The real number entries are ranked. The result is a list of integers viewed as the indexes of the data. The first N′ entries out of the N-position substring are selected to form D'. The same process is applied to the substring representing the set of features. An overall decoding scheme is illustrated in Figure 3.Figure 3
From a particle in[0,1]N+n search space to a subset of instances and features.The information given by the PSO is used to represent the subset of features and data to construct the data-driven fuzzy models. Then, the numerical data are represented in terms of a collection of information granules (a fuzzy sets) produced through some clustering (fuzzy clustering). The information about the granules (clusters) is then used to construct the fuzzy models.In the cooperative PSO, the formation of the search space is realized in a more sophisticated way. The cooperative facet involves mainly exchanging information about the best positions found by the different subswarms. Here, we present a new cooperative PSO (CPSO) algorithm for the data and feature reduction process. The selection of the number of cooperating swarms is important because it will affect the performance of the cooperative PSO model. Sub-swarm 1 represents the features’ column and subswarm 2 represents the instances’ row of the particular data set. Figure4 illustrates the main difference between standard PSO and cooperative PSO. The standard PSO contains one swarm with a large dimension of search space. In contrast, for the cooperative PSO, we divide the search space into two subswarms: subswarm 1 for feature representation and subswarm 2 for instances representation. All the subswarms share the same basic particles definition illustrated in Figure 4.The particle scheme of the “standard” PSO (a) and cooperative PSO (b).
(a)(b)In general, the dimensionality for the data (instances) selection is higher than that of the feature selection. In order to reduce the impact of the curse of dimensionality, we decompose the data into several groups by using the information granulation approach. In this paper, we used the fuzzy C-means (FCMs) to construct the information granules. Therefore, the number of decomposition groups is actually the number of the clusters (c) used in the FCM. For example, if we want to decompose the data into three groups, we use the number of clusters equal to three. As a result, instead of having only two subswarms, we introduce more subswarms that represent different groups of data.Figure5 presents the process of constructing the subswarms for cooperative PSO by decomposing the instances into several subswarms. As mentioned earlier, we apply the concept of information granulation to decompose the data group. In order to identify the selected data in each decomposed group, we use the information granules (membership degrees) values to identify the index of the instances in each group. Here, we employ a winner-takes-all scheme to determine a single group for each granule, that is, the index of the instances in each of the decomposition group related to the information granule that gets the highest degrees of activation. We denote the set of data associated with the ith granules by X(i0):
(3)X(i0)={xk∈X∣Ui0k=max⁡i⁡Uik}for1≤k≤M,1≤i≤c,
where X(i0) is the decomposition groups, Uik is the information granules for each data, xk is the data (instances), M and c are the number of data and the level of information granulation, respectively.Figure 5
The particle scheme of cooperative PSO with more subswarms.


### 4.3. A Category of Fuzzy-Rule-Based Models

To make an overall presentation more focused, we consider a class of fuzzy-rule-based models governed by the collection of “c” rules:
(4)ifxisAitheny=fi(x,ai),
where i=1,2,…,c (c is the number of clusters), Ai are the information granules formed in the input space, and fi is a local linear function with some parameters ai associated with the corresponding information granule. The information granules Ai are constructed by means of fuzzy clustering, namely, fuzzy C-means (FCMs). The corresponding membership functions Ai are thus described as
(5)Ai(x)=1∑j=1c(‖x-vi‖/‖x-vj‖)2(m-1),
where vi, i=1,2,…,c are the prototypes formed through clustering, and m, m>1 is a fuzzification coefficient.


## 4.1. An Overall Reduction Process

As is usual in system modeling, we consider a supervised learning scenario ion which we encounter in a finite set of training data(xk,tk), k=1,2,…,N. By stressing the nature of the data and their dimensionality, the data space along with n-dimensional feature vectors can be viewed as a Cartesian product of the data and features D×F. The essence of the reduction is to arrive at the Cartesian product of the reduced data and feature spaces, D'×F', where, D'∈D and F'∈F. The cardinality of the reduced spaces is equal to N′ and n′, where N′<N and n′<n.The overall scheme of the reduction process outlining a role of the PSO-guided reduction is illustrated in Figure2. The scheme can be divided into two important parts and can be described as follows. (a)
Reduction process via PSO: a reduction process tackles both feature reduction and data reduction simultaneously. PSO algorithm is used to search for the best feature and data for constructing the fuzzy model. Here, the size of the selected features (n′) and data (N′) is provided in advance by the user. After the PSO meets the maximum generation, the process is stopped, and the last best subset of features and data is the best subset of data for constructing the fuzzy model.(b)
Evaluation process: the fuzzy C-means algorithm is used to convert the numerical data into the information granules. Here, the information granularity process deals only with the subset of the data and features (D'×F'). Next, the consequent parameter a constructed from the fuzzy models is use to evaluate the performance of the selected data and features. At this stage we access the performance of the constructed fuzzy model in terms of their capability to fit the model by using the all instances in the original data set.Figure 2
The scheme of the proposed data and reduction for fuzzy modeling.As it becomes apparent, the original spaceD×F is reduced, and in this Cartesian product a fuzzy model, denoted by FM, is designed in the usual way (we elaborate on the form of the fuzzy model in the subsequent section). Its design is guided by a certain objective function Q expressed over all elements of original instances. The quality of the reduced space is assessed by quantifying the performance of the fuzzy model operating over the original, non-reduced space. The same performance index as used in the construction of the fuzzy model in the reduced space is used to describe the quality of the fuzzy model:
(2)Q=1N∑xk∈D×F′(FM(xk)-tk)2.
Note that the summation shown above is taken over all the elements forming the data space D. By taking another look at the overall reduction scheme, it is worth noting that the reduction is realized as in the wrapper mode, in which we use a fuzzy model to evaluate the quality of the reduction mechanism.


## 4.2. The PSO-Based Representation of the Search Space

The reduction of the data and feature spaces involves a selection of a subset of the data and a subset of the features. Therefore, the problem is combinatorial in its nature. PSO is used here to form a subset of integers which are indexes of the data or features to be used in the formation ofD'×F'. For instance, D' is represented as a set of indexes {i1,i2,…,iN′} being a subset of integers {1,2,…,N}. From the perspective of the PSO, the particle is formed as a string of real numbers in [0,1] of the length of N+n; effectively, the search space is a hypercube [0,1]N+n. The first substring of length N represents the data; the second one (having n entries) is used to optimize the subset of features. The particle is decoded as follows. Each substring is processed (decoded) separately. The real number entries are ranked. The result is a list of integers viewed as the indexes of the data. The first N′ entries out of the N-position substring are selected to form D'. The same process is applied to the substring representing the set of features. An overall decoding scheme is illustrated in Figure 3.Figure 3
From a particle in[0,1]N+n search space to a subset of instances and features.The information given by the PSO is used to represent the subset of features and data to construct the data-driven fuzzy models. Then, the numerical data are represented in terms of a collection of information granules (a fuzzy sets) produced through some clustering (fuzzy clustering). The information about the granules (clusters) is then used to construct the fuzzy models.In the cooperative PSO, the formation of the search space is realized in a more sophisticated way. The cooperative facet involves mainly exchanging information about the best positions found by the different subswarms. Here, we present a new cooperative PSO (CPSO) algorithm for the data and feature reduction process. The selection of the number of cooperating swarms is important because it will affect the performance of the cooperative PSO model. Sub-swarm 1 represents the features’ column and subswarm 2 represents the instances’ row of the particular data set. Figure4 illustrates the main difference between standard PSO and cooperative PSO. The standard PSO contains one swarm with a large dimension of search space. In contrast, for the cooperative PSO, we divide the search space into two subswarms: subswarm 1 for feature representation and subswarm 2 for instances representation. All the subswarms share the same basic particles definition illustrated in Figure 4.The particle scheme of the “standard” PSO (a) and cooperative PSO (b).
(a)(b)In general, the dimensionality for the data (instances) selection is higher than that of the feature selection. In order to reduce the impact of the curse of dimensionality, we decompose the data into several groups by using the information granulation approach. In this paper, we used the fuzzy C-means (FCMs) to construct the information granules. Therefore, the number of decomposition groups is actually the number of the clusters (c) used in the FCM. For example, if we want to decompose the data into three groups, we use the number of clusters equal to three. As a result, instead of having only two subswarms, we introduce more subswarms that represent different groups of data.Figure5 presents the process of constructing the subswarms for cooperative PSO by decomposing the instances into several subswarms. As mentioned earlier, we apply the concept of information granulation to decompose the data group. In order to identify the selected data in each decomposed group, we use the information granules (membership degrees) values to identify the index of the instances in each group. Here, we employ a winner-takes-all scheme to determine a single group for each granule, that is, the index of the instances in each of the decomposition group related to the information granule that gets the highest degrees of activation. We denote the set of data associated with the ith granules by X(i0):
(3)X(i0)={xk∈X∣Ui0k=max⁡i⁡Uik}for1≤k≤M,1≤i≤c,
where X(i0) is the decomposition groups, Uik is the information granules for each data, xk is the data (instances), M and c are the number of data and the level of information granulation, respectively.Figure 5
The particle scheme of cooperative PSO with more subswarms.


## 4.3. A Category of Fuzzy-Rule-Based Models

To make an overall presentation more focused, we consider a class of fuzzy-rule-based models governed by the collection of “c” rules:
(4)ifxisAitheny=fi(x,ai),
where i=1,2,…,c (c is the number of clusters), Ai are the information granules formed in the input space, and fi is a local linear function with some parameters ai associated with the corresponding information granule. The information granules Ai are constructed by means of fuzzy clustering, namely, fuzzy C-means (FCMs). The corresponding membership functions Ai are thus described as
(5)Ai(x)=1∑j=1c(‖x-vi‖/‖x-vj‖)2(m-1),
where vi, i=1,2,…,c are the prototypes formed through clustering, and m, m>1 is a fuzzification coefficient.


## 5. Experimental Studies

In this section, we report our results from a set of experiments, using several machine learning data sets (seehttp://www.ics.uci.edu/~mlearn/MLRepository.html and http://lib.stat.cmu.edu/datasets/). The main objective of these experiments is to show the abilities of the proposed approach, quantify the performance of the selected subsets of features and instances, and arrive at some general conclusions. A concise summary of the data sets used in the experiment is presented in Table 2. All the data concern continuous output.Table 2
Description of data used in the experiments.
Data setAbbreviationNumber of featuresNumber of dataSparsity ration,κAir pollution PM10PM10750071.43Boston housingHousing1350638.92Body fatBody fat1425218.00Parkinson’s telemonitoringParkinson175875345.59Computer activityComputer218192390.09


### 5.1. Parameter Setup

The values of the PSO and CPSO parameters were set using the standard form as follows. The values of the inertia weight,w, were linearly from 1 to 0 over the course of optimization. The values of the cognitive factor, c1, and social factor c2, were set to 0.5 and 1.5, respectively. In Table 3, we also list the numeric values of the parameters of the PSO and CPSO environment. As to the size of the population and the number of generations, we used a larger population and a larger number of generations in the generic version of the PSO than in the CPSO because of the larger search space this algorithm operates in.Table 3
The values of the parameters used in the experiments; CPSO1: swarms located in the feature space. CPSO2: swarms located in the instance (data) space.
Optimization methodSubswarmsGenerationParticlesPSO150300CPSO1250100CPSO243050The number of subswarms used in the optimization method is also shown in Table3. The PSO method comprises only a single swarm whose individuals concatenate features and instances. In contrast, for the CPSO, we divided the search space into several subswarms that can cooperate with each other and where the individuals in the subswarms are used to represent a portion of the search space. The CPSO1 contains two subswarms that cover the data and features, respectively. In CPSO2, we used three subswarms to represent data point; in the data used here, the number of data is larger than the number of features, so a better balance of the dimensionality of the spaces is achieved. The data (instances) search space is divided into three subswarms, and the decomposition process is realized by running fuzzy clustering (each cluster forms a subswarm). In the table we used a smaller size of generation compared to particles size. This is because in [34] Shi and Eberhart mentioned that the population size does not exhibit any significant impact on the performance of the PSO method. However, the size of particles is high given the size of the search space. Here we require more particles to capture the large search space of instances selection for using the standard PSO. As a result we can find the best solution faster than using a smaller particles size. On the other hand, the number of particle is decreased when we implement the CPSO method. This is because the original large search space is divided into several groups, and the processes of searching the best subset are done in parallel.


### 5.2. Results of the Experiments

In the experiments, we looked at the performance—an average root mean squared error (RMSE)—obtained for the selected combinations of the number of features and data (instances). The results obtained for the Housing data, PM10 data, and Parkinson’s data forc=4 and c=3 clusters are summarized in Tables 4, 5, and 6, respectively. The experiments were repeated 10 times, and the reported results are the average RMSE values. We also report the values of the standard deviation of the performance index to offer a better insight into the variability of the performance. It is noticeable that the standard deviation is reduced with the increase of the data involved and the decrease of the dimensionality of the feature space.Table 4
Results for housing data; the number of clusters is set to 4,c=4; κ is the ratio of the number of selected data versus the number of selected features.
FeatureData = 10%Data = 20%Data = 30%Data = 40%Data = 50%Data = 60%Data = 70%Data = 80%Data = 90%Data = 100%(# of data = 51)(# of data = 101)(# of data = 152)(# of data = 202)(# of data = 253)(# of data = 304)(# of data = 354)(# of data = 405)(# of data = 455)(# of data = 506)κRMSEκRMSEκRMSEκRMSEκRMSEκRMSEκRMSEκRMSEκRMSEκRMSE10%51.06.341±0.253101.06.262±0.171152.05.777±0.207202.06.227±0.351253.06.389±0.0633046.483±0.214354.06.610±0.211405.06.387±0.026455.06.655  ±  0.0615067.437  ±  020%17.06.664±0.23333.75.389±0.25350.75.191±0.18567.34.884±0.11884.34.805±0.0471014.882±0.080118.04.906±0.040135.05.011±0.077151.75.149  ±  0.0711695.039  ±  0.09630%12.86.127±0.24525.35.468±0.29038.04.853±0.18350.54.574±0.07863.34.619±0.294764.45±0.09288.54.398±0.035101.34.56±0.053113.84.536  ±  0.0261274.578  ±  0.0140%10.26.321±0.51820.25.122±0.24530.44.626±0.07540.44.362±0.19050.64.172±0.12360.84.06±0.09970.84.126±0.10581.04.18±0.17791.04.343  ±  0.0261014.233  ±  0.07950%7.37.126±0.83514.45.046±0.31221.74.574±0.20628.94.018±0.10936.13.916±0.22443.43.927±0.08950.64.009±0.09357.94.085±0.13265.04.077  ±  0.09272.34.005  ±  0.08260%6.48.133±0.78212.65.120±0.18919.04.504±0.20725.34.052±0.19631.63.912±0.120383.935±0.12944.33.934±0.03550.63.923±0.11756.93.931  ±  0.10263.33.803  ±  0.08870%5.79.379±0.98411.25.003±0.23216.94.345±0.13422.43.949±0.12528.13.721±0.07133.83.722±0.06539.33.722±0.06645.03.787±0.04650.63.799  ±  0.04156.23.668  ±  0.04380%5.110.57±2.25110.15.107±0.26215.24.232±0.17320.23.67±0.09325.33.617±0.10830.43.659±0.12835.43.568±0.06440.53.567±0.08645.53.645  ±  0.06550.63.526  ±  0.04990%4.324.05±7.6818.45.324±0.20712.74.173±0.18116.83.809±0.08021.13.652±0.05025.33.569±0.02829.53.555±0.01733.83.533±0.01637.93.541±0.01542.23.550±0100%3.944.39±17.657.85.409±0.20111.74.082±0.04715.53.781±0.05519.53.687±0.03823.43.654±0.01527.23.631±0.02131.23.615±0.01135.03.605±0.01538.94.023±0Table 5
Results for PM10 dataset;c=3.
FeatureData = 10%Data = 20%Data = 30%Data = 40%Data = 50%Data = 60%Data = 70%Data = 80%Data = 90%Data = 100%(# of data = 50)(# of data = 100)(# of data = 150)(# of data = 200)(# of data = 250)(# of data = 300)(# of data = 350)(# of data = 400)(# of data = 450)(# of data = 500)κRMSEκRMSEκRMSEκRMSEκRMSEκRMSEκRMSEκRMSEκRMSEκRMSE10%50.00.931  ±  0.036100.00.979  ±  0.018150.00.983  ±  0.006200.00.985  ±  0.010250.00.998  ±  0.021300.01.036  ±  0.022350.01.071  ±  0.016400.01.088  ±  0.002450.01.099  ±  0.003500.01.116  ±  020%50.00.896  ±  0.034100.00.98  ±  0.013150.00.987  ±  0.009200.00.994  ±  0.008250.01.004  ±  0.018300.01.025  ±  0.032350.01.075  ±  0.009400.01.090  ±  0.003450.01.098  ±  0.003500.01.116  ±  030%25.00.825  ±  0.08750.00.902  ±  0.0475.00.918  ±  0.007100.00.916  ±  0.003125.00.918  ±  0.010150.00.920  ±  0.005175.00.922  ±  0.003200.00.937  ±  0.001225.00.948  ±  0.005250.00.964  ±  040%16.70.829  ±  0.02333.30.877  ±  0.00750.00.877  ±  0.01066.70.862  ±  0.00983.30.865  ±  0.004100.00.869  ±  0.008116.70.887  ±  0.004133.30.892  ±  0.006150.00.902  ±  0.002166.70.915  ±  050%12.50.802  ±  0.02925.00.816  ±  0.00837.50.822  ±  0.01250.00.826  ±  0.02862.50.843  ±  0.02475.00.870  ±  0.02387.50.891  ±  0.022100.00.906  ±  0.004112.50.905  ±  0.002125.00.925  ±  060%12.50.804  ±  0.02725.00.818  ±  0.01337.50.825  ±  0.01450.00.834  ±  0.01662.50.841  ±  0.00675.00.879  ±  0.01287.50.898  ±  0.006100.00.897  ±  0.006112.50.907  ±  0.002125.00.925  ±  070%10.00.783  ±  0.03020.00.781  ±  0.03130.00.804  ±  0.01740.00.782  ±  0.03950.00.806  ±  0.01460.00.832  ±  0.01670.00.851  ±  0.01080.00.856  ±  0.00490.00.864  ±  0.003100.00.900  ±  080%8.30.768  ±  0.02416.70.769  ±  0.00725.00.776  ±  0.01733.30.768  ±  0.02441.70.796  ±  0.01050.00.805  ±  0.01658.30.826  ±  0.01366.70.839  ±  0.00775.00.847  ±  0.00183.30.878  ±  090%8.30.774  ±  0.02616.70.771  ±  0.01725.00.767  ±  0.01433.30.796  ±  0.01041.70.777  ±  0.01950.00.815  ±  0.01758.30.820  ±  0.01466.70.843  ±  0.00175.00.851  ±  0.00383.30.878  ±  0100%7.10.786  ±  0.01214.30.758  ±  0.01721.40.764  ±  0.00728.60.765  ±  0.01535.70.772  ±  0.01642.90.795  ±  0.00450.00.808  ±  0.01057.10.818  ±  0.00564.30.824  ±  0.00371.40.883  ±  0Table 6
Results for Parkinson’s data;c=3.
FeatureData = 10%Data = 20%Data = 30%Data = 40%Data = 50%Data = 60%Data = 70%Data = 80%(# of data = 346)(# of data = 691)(# of data = 1037)(# of data = 1382)(# of data = 1728)(# of data = 2074)(# of data =2419)(# of data = 2765)κRMSEκRMSEκRMSEκRMSEκRMSEκRMSEκRMSEκRMSE10%3466.388  ±  0.1826916.221  ±  0.06410376.393  ±  0.14013826.495  ±  0.10617286.644  ±  0.13720746.644  ±  0.13724196.515  ±  0.00627656.406  ±  0.00520%1736.183  ±  0.1103465.857  ±  0.0315195.932  ±  0.0236915.972  ±  0.0088646.247  ±  0.13710376.247  ±  0.13712106.066  ±  0.01813825.970  ±  0.00230%1156.152  ±  0.1512305.799  ±  0.0613465.703  ±  0.0674615.708  ±  0.0455766.077  ±  0.0486916.077  ±  0.0488065.964  ±  0.2769225.740  ±  0.08940%866.328  ±  0.4561735.900  ±  0.0872595.886  ±  0.3423466.175  ±  0.3584326.461  ±  0.2615186.461  ±  0.2616056.299  ±  0.5636916.308  ±  0.67750%696.991  ±  0.5251387.585  ±  0.7552076.448  ±  0.5802767.493  ±  0.3033468.057  ±  0.0374158.057  ±  0.0374848.160  ±  0.0185538.110  ±  0.01260%588.357  ±  0.1641158.088  ±  0.0281738.069  ±  0.0782307.960  ±  0.0212888.104  ±  0.0213468.104  ±  0.0214038.147  ±  0.0174618.101  ±  0.00270%498.401  ±  0.068998.074  ±  0.0641488.087  ±  0.0151978.008  ±  0.0112478.123  ±  0.0092968.123  ±  0.0093468.158  ±  0.0113958.107  ±  0.00780%438.419  ±  0.091868.257  ±  0.0481308.187  ±  0.0171738.092  ±  0.0202168.177  ±  0.0082598.177  ±  0.0083028.194  ±  0.0083468.138  ±  0.00590%388.500  ±  0.106778.258  ±  0.0241158.199  ±  0.0171548.139  ±  0.0131928.200  ±  0.0092308.200  ±  0.0092698.217  ±  0.0043078.157  ±  0.006100%358.560  ±  0.026698.249  ±  0.0331048.262  ±  0.0131388.223  ±  0.0071738.239  ±  0.0062078.239  ±  0.0062428.235  ±  0.0022768.223  ±  0.001The visualization of the results in the form of a series of heat maps, see Figures6, 7, and 8, helps us arrive at a number of qualitative observations as well as to look at some quantitative relationships. In most cases, the performance index remains relatively low in some regions of the heat map. This finding demonstrates that the available data come with some evident redundancy, which exhibits a negative impact on the designed model. For the PM10 data, there is a significantly reduced performance of the model when, for a low percentage of data, the number of features starts growing. This effect is present for different numbers of clusters. The same tendency is noticeable for the other data sets. There is a sound explanation to this phenomenon: simply, the structure formed by fuzzy clustering does not fully reflect the dependencies in the data (due to the effect of the sparsity of the data), and this problem, in turn, results in the deteriorating performance of the fuzzy model. In this case, one would be better off to consider a suitable reduced set of features. In all cases experimented with, we noted an optimal combination of features and data that led to the best performance of the model. Table 7 summarizes the optimal combinations of features and data.Table 7
The optimal % of features and data for different clusters.
Data set and number of clusters% of features% of dataPima withC=37040Pima withC=410020Pima withC=510020Pima withC=610040Pima withC=710080Housing withC=38070Housing withC=48050Housing withC=580100Housing withC=68080Housing withC=790100Body Fat withC=33030Body Fat withC=410070Body Fat withC=59070Body Fat withC=69090Parkinson withC=33030Heat map for PM 10 data forc varying in-between 3 to 6.
(a)(b)(c)(d)Heat map for Body fat data for varying in-between 3 to 6.
(a)(b)(c)(d)Heat map for housing data forc=3, 4, 5, and 7.
(a)(b)(c)(d)The relationships between the percentage of data used and the resulting RMSE values are displayed in Figures9 and 10. Some interesting tendencies are worth noting. A critical number of data are required to form a fuzzy model. Increasing the number of data does not produce any improvement as the curves plotted on Figures 9(a), 9(c), and 9(a) achieve a plateau or even some increase of the RMSE is noticeable.The values of RMSE versus the percentage of data for selected number clusters: (a) housing data, (b) PM10 data, and (c) body fat data.
(a)(b)(c)Plots of RMSE versus the percentage of features for selected number clusters: (a) housing data, (b) PM10 data, and (c) body Fat data.
(a)(b)(c)Considering a fixed percentage of the data used, we look at the nature of the feature sets. Tables8, 9 and 10 displays the best feature for PM10 data, Body fat data, and Housing data, respectively.  Overall, the selected subsets of features are almost the same for different numbers of the clusters being used. Furthermore, we observe that in most cases, the reduced feature spaces exhibit an interesting “nesting” property, meaning that the extended feature space constructed subsumes the one formed previously. For example, for the Housing data, we obtain the following subsets of features:
(6){feature6}⊂{feature6,feature9,feature13}⊂{feature6,feature9,feature10,feature13}.
Here, the corresponding features are as follows: 6: average number of rooms per dwelling, 9: index of accessibility to radial highways, 13: percentage of lower status population, and 10: full-value property-tax rate per $10,000. This combination is quite convincing.Table 8
Best subsets of features for PM10 data.
FSubset of featuresD = 30%D = 40%D = 50%D = 60%D = 70%D = 80%D = 90%10%111111120%111111130%1, 71, 71, 71, 71, 61, 61, 640%1, 6, 71, 6, 71, 6, 71, 6, 71, 6, 71, 6, 71, 6, 750%1, 2, 6, 71, 2, 6, 71, 2, 6, 71, 2, 6, 71, 2, 6, 71, 2, 6, 71, 2, 6, 760%1, 2, 6, 71, 2, 6, 71, 2, 6, 71, 2, 6, 71, 2, 6, 71, 2, 6, 71, 2, 6, 770%1, 2, 3, 6, 71, 2, 3, 6, 71, 2, 3, 6, 71, 2, 3, 6, 71, 2, 3, 6, 71, 2, 4, 6, 71, 2, 4, 6, 780%1, 2, 3, 4, 6, 71, 2, 3, 5, 6, 71, 2, 3, 4, 6, 71, 2, 3, 4, 6, 71, 2, 3, 4, 6, 71, 2, 3, 4, 6, 71, 2, 3, 4, 6, 790%1, 2, 3, 4, 6, 71, 2, 3, 4, 6, 71, 2, 3, 4, 6, 71, 2, 3, 5, 4, 71, 2, 3, 4, 5, 71, 2, 3, 4, 6, 71, 2, 3, 5, 4, 7Table 9
Best subsets of features for body fat data.
FSubset of featuresD = 50%D = 60%D = 70%D = 80%10%111120%1, 6, 71, 3, 71, 3, 71, 3, 730%1, 6, 7, 91, 7, 8, 91, 3, 7, 91, 3, 7, 940%1, 6, 7, 8, 9, 121, 3, 4, 7, 8, 91, 3, 6, 7, 8, 91, 7, 8, 9, 11, 1250%1, 2, 6, 7, 8, 12, 141, 4, 8, 9, 11, 12, 141, 3, 7, 8, 9, 12, 141, 3, 4, 5, 7, 8, 1260%1, 2, 6, 7, 8, 11, 12, 141, 3, 4, 7, 8, 9, 11, 141, 3, 5, 7, 9, 11, 12, 141, 3, 5, 7, 8, 11, 12, 1470%1, 3, 4, 5, 6, 8, 9, 11, 12, 141, 3, 4, 5, 7, 8, 9, 11, 12, 141, 3, 4, 5, 7, 8, 9, 11, 12, 141, 3, 4, 5, 7, 8, 9, 11, 12, 1480%1, 3, 4, 5, 6, 7, 8, 9, 10, 12, 141, 3, 4, 5, 6, 7, 8, 9, 10, 12, 141, 3, 4, 5, 6, 7, 8, 9, 11, 12, 141, 3, 4, 5, 6, 7, 8, 9, 11, 12, 1490%1, 2, 3, 4, 5, 6, 8, 9, 10, 11, 12, 13, 141, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 141, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 141, 23, 4, 5, 6, 7, 8, 9, 10, 11, 12, 14Table 10
Best subsets of features for housing data.
FSubset of featuresD = 10%D = 30%D = 50%D = 70%D = 90%10%1266121220%5, 12, 136, 9, 136, 9, 106, 9, 106, 9, 1030%6, 7, 9, 135, 6, 9, 114, 6, 9, 131, 6, 9, 102, 3, 12, 1340%1, 3, 6, 10, 134, 6, 9, 10, 136, 9, 10, 12, 133, 5, 6, 9, 101, 4, 6, 9, 1050%3, 6, 7, 8, 9, 10, 111, 2, 6, 9, 10, 11, 131, 3, 5, 6, 8, 9, 111, 5, 6, 9, 10, 11, 121, 6, 9, 10, 11, 12, 1360%3, 5, 6, 7, 8, 9, 12, 133, 5, 6, 7, 8, 9, 10, 111, 3, 5, 6, 7, 9, 10, 111, 3, 5, 6, 7, 9, 10, 111, 3, 6, 7, 8, 9, 11, 1370%3, 5, 6, 7, 8, 9, 10, 12, 131, 3, 4, 5, 6, 9, 10, 11, 131, 3, 5, 6, 7, 9, 10, 11, 131, 3, 5, 6, 7, 9, 10, 11, 131, 3, 6, 7, 8, 9, 10, 11, 1380%1, 3, 4, 5, 6, 8, 9, 10, 11, 131, 3, 5, 6, 8, 9, 10, 11, 12, 131, 3, 5, 6, 7, 8, 9, 10, 11, 131, 3, 5, 6, 7, 8, 9, 10, 11,  131, 3, 5, 6, 7, 8, 9, 10, 11, 1390%1, 2, 3, 4, 5, 6, 7, 8, 10, 11, 12, 131, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 131, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 131, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 131, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13For the PM10 data, we arrive at a series of nested collections of features:(7){feature1}⊂{feature1,feature7}⊂{feature1,feature6,feature7}⊂{feature1,feature2,feature6,feature7},
where the corresponding features include: 1: the concentration of PM10 (particles), 7: hour of experiment per day, 6: wind direction, and 2: the number of cars per hour.Turning to the comparative analysis of performance of the swarm optimization methods, we summarize the obtained results in Figure11. For all data, the CPSO performed better than the standard PSO. Although both algorithms show the same tendency when the percentage of feature is 100% however, the RMSE produced by the CPSO is lower than the one obtained when running the PSO. Furthermore, the CPSO algorithm is more stable than the standard PSO. In most cases, the standard deviations of error produced by the CPSO are smaller than the results obtained for the standard PSO (see Table 11).Table 11
Standard deviations for PSO and CPSO (housing and PM10 data sets).
Housing (D = 50%)PM10 (D = 50%)PSOCPSOPSOCPSO0.0660.0720.0210.0070.1920.0150.0180.0090.1990.0390.0100.0070.110.0930.0040.0070.1150.0790.0240.0080.0910.0710.0060.0090.0580.0940.0140.0100.0440.0640.0100.0070.0530.0420.0190.0090.0210.0420.0160.009Values of RMSE versus the percentage of features selected when running PSO and CPSO—the use of the housing dataset: (a) 20% of selected data, (b) 30% of selected data, (c) 50% of selected data, and (d) 70% of selected data.
(a)(b)(c)(d)Figure12 shows the subsets of the features selected for different percentages of the features used in construction of the fuzzy model. The CPSO algorithm is more consistent while selecting the increasing number of features. For example, features 6 and 13 were selected when using both 30% and 70% of data. In contrast to the selection made with the PSO algorithm, the subset of the features selected here is not as stable, especially when using only 30% of data.Comparison of sets of features being selected by using PSO and CPSO2 for Housing data dataset: (a) PSO method with 30% of selected data, (b) CPSO2 method with 30% of selected data, (c) PSO method with 70% of selected data, and (d) CPSO2 method with 70% of selected data.
(a)(b)(c)(d)Table12 presents the percentage of the improvement when using the CPSO algorithm compared to the PSO algorithm. Note that in this percentage we included all different combinations of the features’ percentages and the data percentages being used. The percentage of the improvement is higher when dealing with a smaller percentage of features and data. For example, the percentage of improvement is 34% for 10% of the instances and 50% of the features selected while the percentage of improvement is less than 10% for 60% of instances and features used. These results occurred because the PSO method has to deal with a large search space for selecting a small subset of features and instances. In contrast to the search space for CPSO, the large search space is decomposed into multiple subswarms that reduce the dimensionality of the original search space.Table 12
Percentage of improvement of the RMSE obtained when using CPSO over the results formed by the PSO; Housing data set.
FD = 10%D = 20%D = 30%D = 40%D = 50%D = 60%D = 70%D = 80%D = 90%10%13152026191810191120%31312024251624231730%33261920141524182140%28211914171717181950%342112799912860%23211374789770%19171449648880%33191242338990%22144521221100%1744731111Tables13, 14, 15, and 16 show the comparison of RMSE when using the proposed method and the standard fuzzy modeling method. Here the standard fuzzy model is constructed without using any feature and instances selection and the holdout method is used to select the data based on the percentage given. The experiment for using the standard fuzzy modeling is repeated for 50 times. If we analyze the tables, we can observe that our proposed method outperforms the standard method of constructing the fuzzy model from the dataset. This can be seen clearly when using the CPSO method to search for the best subset of feature and instances. For example, in Table 13 if we use the CPSO method, the RMSE for using 70% of data is 3.413, whereas the RMSE for the standard method is 8.312. The same tendency occurs for all datasets used here.Table 13
The comparison of RMSE obtained when using standard PSO, CPSO, and standard fuzzy model with holdout method for housing data withC=3.
% of dataStandard PSOCooperative PSO1Holdout method% of featureMSE% of featureMSE% of featureMSE30904.015403.47310017.59340803.699703.46410010.80350803.573703.4141009.90760803.556703.4351008.50770803.527603.4131008.31280803.654903.4491008.16490803.679903.6151007.641Table 14
The comparison of RMSE obtained when using standard PSO, CPSO, and standard fuzzy model with holdout method for body fat data withC=3.
% of dataStandard PSOCooperative PSO1Holdout method% of featureMSE% of featureMSE% of featureMSE30304.677304.684710011.58640304.717304.54091008.29150404.617304.51361007.54860504.636404.42891007.07370404.548404.42341006.65880404.553404.42331006.23990404.582404.37711006.102Table 15
The comparison of RMSE obtained when using standard PSO, CPSO, and standard fuzzy model with holdout method for PM10 data withC=3.
% of dataStandard PSOCooperative PSO1Holdout method% of featureMSE% of featureMSE% of featureMSE301000.764800.73381002.100401000.765800.74321002.01850900.777900.77901002.03060800.805800.77691001.98670900.820800.80521002.00180800.839900.82061001.98390700.847900.84171001.976Table 16
The comparison of RMSE obtained when using CPSO and standard fuzzy model with holdout method for computer data withC=3.
% of dataCooperative PSOHoldout method% of featureMSE% of featureMSE304016.44610017.453403014.71210017.524503014.35010017.680604014.93510017.837704016.23710017.918804015.12210018.351Figure13 shows the comparison plot between the proposed method and the “standard” fuzzy modeling. In most of the cases, the proposed method showed better performance.Comparison of RMSE by using proposed method (straight line) and standard fuzzy model (dotted line): (a) Housing dataset withc=4, (b) body fat dataset with c=5, (c) Parkinson’s dataset with c=3, and (d) computer dataset with c=3.
(a)(b)(c)(d)It becomes clear that one is able to reduce the input data in terms of the number features and instances. Moreover, the flexibility of choosing the reduction level helps the user focus on the most essential subsets of data and features (variables). The knowledge acquired about the best subset of data can be used for future data collection. In addition, the user can put more effort analyzing only the best subset of data that give more impact to the overall prediction.


## 5.1. Parameter Setup

The values of the PSO and CPSO parameters were set using the standard form as follows. The values of the inertia weight,w, were linearly from 1 to 0 over the course of optimization. The values of the cognitive factor, c1, and social factor c2, were set to 0.5 and 1.5, respectively. In Table 3, we also list the numeric values of the parameters of the PSO and CPSO environment. As to the size of the population and the number of generations, we used a larger population and a larger number of generations in the generic version of the PSO than in the CPSO because of the larger search space this algorithm operates in.Table 3
The values of the parameters used in the experiments; CPSO1: swarms located in the feature space. CPSO2: swarms located in the instance (data) space.
Optimization methodSubswarmsGenerationParticlesPSO150300CPSO1250100CPSO243050The number of subswarms used in the optimization method is also shown in Table3. The PSO method comprises only a single swarm whose individuals concatenate features and instances. In contrast, for the CPSO, we divided the search space into several subswarms that can cooperate with each other and where the individuals in the subswarms are used to represent a portion of the search space. The CPSO1 contains two subswarms that cover the data and features, respectively. In CPSO2, we used three subswarms to represent data point; in the data used here, the number of data is larger than the number of features, so a better balance of the dimensionality of the spaces is achieved. The data (instances) search space is divided into three subswarms, and the decomposition process is realized by running fuzzy clustering (each cluster forms a subswarm). In the table we used a smaller size of generation compared to particles size. This is because in [34] Shi and Eberhart mentioned that the population size does not exhibit any significant impact on the performance of the PSO method. However, the size of particles is high given the size of the search space. Here we require more particles to capture the large search space of instances selection for using the standard PSO. As a result we can find the best solution faster than using a smaller particles size. On the other hand, the number of particle is decreased when we implement the CPSO method. This is because the original large search space is divided into several groups, and the processes of searching the best subset are done in parallel.


## 5.2. Results of the Experiments

In the experiments, we looked at the performance—an average root mean squared error (RMSE)—obtained for the selected combinations of the number of features and data (instances). The results obtained for the Housing data, PM10 data, and Parkinson’s data forc=4 and c=3 clusters are summarized in Tables 4, 5, and 6, respectively. The experiments were repeated 10 times, and the reported results are the average RMSE values. We also report the values of the standard deviation of the performance index to offer a better insight into the variability of the performance. It is noticeable that the standard deviation is reduced with the increase of the data involved and the decrease of the dimensionality of the feature space.Table 4
Results for housing data; the number of clusters is set to 4,c=4; κ is the ratio of the number of selected data versus the number of selected features.
FeatureData = 10%Data = 20%Data = 30%Data = 40%Data = 50%Data = 60%Data = 70%Data = 80%Data = 90%Data = 100%(# of data = 51)(# of data = 101)(# of data = 152)(# of data = 202)(# of data = 253)(# of data = 304)(# of data = 354)(# of data = 405)(# of data = 455)(# of data = 506)κRMSEκRMSEκRMSEκRMSEκRMSEκRMSEκRMSEκRMSEκRMSEκRMSE10%51.06.341±0.253101.06.262±0.171152.05.777±0.207202.06.227±0.351253.06.389±0.0633046.483±0.214354.06.610±0.211405.06.387±0.026455.06.655  ±  0.0615067.437  ±  020%17.06.664±0.23333.75.389±0.25350.75.191±0.18567.34.884±0.11884.34.805±0.0471014.882±0.080118.04.906±0.040135.05.011±0.077151.75.149  ±  0.0711695.039  ±  0.09630%12.86.127±0.24525.35.468±0.29038.04.853±0.18350.54.574±0.07863.34.619±0.294764.45±0.09288.54.398±0.035101.34.56±0.053113.84.536  ±  0.0261274.578  ±  0.0140%10.26.321±0.51820.25.122±0.24530.44.626±0.07540.44.362±0.19050.64.172±0.12360.84.06±0.09970.84.126±0.10581.04.18±0.17791.04.343  ±  0.0261014.233  ±  0.07950%7.37.126±0.83514.45.046±0.31221.74.574±0.20628.94.018±0.10936.13.916±0.22443.43.927±0.08950.64.009±0.09357.94.085±0.13265.04.077  ±  0.09272.34.005  ±  0.08260%6.48.133±0.78212.65.120±0.18919.04.504±0.20725.34.052±0.19631.63.912±0.120383.935±0.12944.33.934±0.03550.63.923±0.11756.93.931  ±  0.10263.33.803  ±  0.08870%5.79.379±0.98411.25.003±0.23216.94.345±0.13422.43.949±0.12528.13.721±0.07133.83.722±0.06539.33.722±0.06645.03.787±0.04650.63.799  ±  0.04156.23.668  ±  0.04380%5.110.57±2.25110.15.107±0.26215.24.232±0.17320.23.67±0.09325.33.617±0.10830.43.659±0.12835.43.568±0.06440.53.567±0.08645.53.645  ±  0.06550.63.526  ±  0.04990%4.324.05±7.6818.45.324±0.20712.74.173±0.18116.83.809±0.08021.13.652±0.05025.33.569±0.02829.53.555±0.01733.83.533±0.01637.93.541±0.01542.23.550±0100%3.944.39±17.657.85.409±0.20111.74.082±0.04715.53.781±0.05519.53.687±0.03823.43.654±0.01527.23.631±0.02131.23.615±0.01135.03.605±0.01538.94.023±0Table 5
Results for PM10 dataset;c=3.
FeatureData = 10%Data = 20%Data = 30%Data = 40%Data = 50%Data = 60%Data = 70%Data = 80%Data = 90%Data = 100%(# of data = 50)(# of data = 100)(# of data = 150)(# of data = 200)(# of data = 250)(# of data = 300)(# of data = 350)(# of data = 400)(# of data = 450)(# of data = 500)κRMSEκRMSEκRMSEκRMSEκRMSEκRMSEκRMSEκRMSEκRMSEκRMSE10%50.00.931  ±  0.036100.00.979  ±  0.018150.00.983  ±  0.006200.00.985  ±  0.010250.00.998  ±  0.021300.01.036  ±  0.022350.01.071  ±  0.016400.01.088  ±  0.002450.01.099  ±  0.003500.01.116  ±  020%50.00.896  ±  0.034100.00.98  ±  0.013150.00.987  ±  0.009200.00.994  ±  0.008250.01.004  ±  0.018300.01.025  ±  0.032350.01.075  ±  0.009400.01.090  ±  0.003450.01.098  ±  0.003500.01.116  ±  030%25.00.825  ±  0.08750.00.902  ±  0.0475.00.918  ±  0.007100.00.916  ±  0.003125.00.918  ±  0.010150.00.920  ±  0.005175.00.922  ±  0.003200.00.937  ±  0.001225.00.948  ±  0.005250.00.964  ±  040%16.70.829  ±  0.02333.30.877  ±  0.00750.00.877  ±  0.01066.70.862  ±  0.00983.30.865  ±  0.004100.00.869  ±  0.008116.70.887  ±  0.004133.30.892  ±  0.006150.00.902  ±  0.002166.70.915  ±  050%12.50.802  ±  0.02925.00.816  ±  0.00837.50.822  ±  0.01250.00.826  ±  0.02862.50.843  ±  0.02475.00.870  ±  0.02387.50.891  ±  0.022100.00.906  ±  0.004112.50.905  ±  0.002125.00.925  ±  060%12.50.804  ±  0.02725.00.818  ±  0.01337.50.825  ±  0.01450.00.834  ±  0.01662.50.841  ±  0.00675.00.879  ±  0.01287.50.898  ±  0.006100.00.897  ±  0.006112.50.907  ±  0.002125.00.925  ±  070%10.00.783  ±  0.03020.00.781  ±  0.03130.00.804  ±  0.01740.00.782  ±  0.03950.00.806  ±  0.01460.00.832  ±  0.01670.00.851  ±  0.01080.00.856  ±  0.00490.00.864  ±  0.003100.00.900  ±  080%8.30.768  ±  0.02416.70.769  ±  0.00725.00.776  ±  0.01733.30.768  ±  0.02441.70.796  ±  0.01050.00.805  ±  0.01658.30.826  ±  0.01366.70.839  ±  0.00775.00.847  ±  0.00183.30.878  ±  090%8.30.774  ±  0.02616.70.771  ±  0.01725.00.767  ±  0.01433.30.796  ±  0.01041.70.777  ±  0.01950.00.815  ±  0.01758.30.820  ±  0.01466.70.843  ±  0.00175.00.851  ±  0.00383.30.878  ±  0100%7.10.786  ±  0.01214.30.758  ±  0.01721.40.764  ±  0.00728.60.765  ±  0.01535.70.772  ±  0.01642.90.795  ±  0.00450.00.808  ±  0.01057.10.818  ±  0.00564.30.824  ±  0.00371.40.883  ±  0Table 6
Results for Parkinson’s data;c=3.
FeatureData = 10%Data = 20%Data = 30%Data = 40%Data = 50%Data = 60%Data = 70%Data = 80%(# of data = 346)(# of data = 691)(# of data = 1037)(# of data = 1382)(# of data = 1728)(# of data = 2074)(# of data =2419)(# of data = 2765)κRMSEκRMSEκRMSEκRMSEκRMSEκRMSEκRMSEκRMSE10%3466.388  ±  0.1826916.221  ±  0.06410376.393  ±  0.14013826.495  ±  0.10617286.644  ±  0.13720746.644  ±  0.13724196.515  ±  0.00627656.406  ±  0.00520%1736.183  ±  0.1103465.857  ±  0.0315195.932  ±  0.0236915.972  ±  0.0088646.247  ±  0.13710376.247  ±  0.13712106.066  ±  0.01813825.970  ±  0.00230%1156.152  ±  0.1512305.799  ±  0.0613465.703  ±  0.0674615.708  ±  0.0455766.077  ±  0.0486916.077  ±  0.0488065.964  ±  0.2769225.740  ±  0.08940%866.328  ±  0.4561735.900  ±  0.0872595.886  ±  0.3423466.175  ±  0.3584326.461  ±  0.2615186.461  ±  0.2616056.299  ±  0.5636916.308  ±  0.67750%696.991  ±  0.5251387.585  ±  0.7552076.448  ±  0.5802767.493  ±  0.3033468.057  ±  0.0374158.057  ±  0.0374848.160  ±  0.0185538.110  ±  0.01260%588.357  ±  0.1641158.088  ±  0.0281738.069  ±  0.0782307.960  ±  0.0212888.104  ±  0.0213468.104  ±  0.0214038.147  ±  0.0174618.101  ±  0.00270%498.401  ±  0.068998.074  ±  0.0641488.087  ±  0.0151978.008  ±  0.0112478.123  ±  0.0092968.123  ±  0.0093468.158  ±  0.0113958.107  ±  0.00780%438.419  ±  0.091868.257  ±  0.0481308.187  ±  0.0171738.092  ±  0.0202168.177  ±  0.0082598.177  ±  0.0083028.194  ±  0.0083468.138  ±  0.00590%388.500  ±  0.106778.258  ±  0.0241158.199  ±  0.0171548.139  ±  0.0131928.200  ±  0.0092308.200  ±  0.0092698.217  ±  0.0043078.157  ±  0.006100%358.560  ±  0.026698.249  ±  0.0331048.262  ±  0.0131388.223  ±  0.0071738.239  ±  0.0062078.239  ±  0.0062428.235  ±  0.0022768.223  ±  0.001The visualization of the results in the form of a series of heat maps, see Figures6, 7, and 8, helps us arrive at a number of qualitative observations as well as to look at some quantitative relationships. In most cases, the performance index remains relatively low in some regions of the heat map. This finding demonstrates that the available data come with some evident redundancy, which exhibits a negative impact on the designed model. For the PM10 data, there is a significantly reduced performance of the model when, for a low percentage of data, the number of features starts growing. This effect is present for different numbers of clusters. The same tendency is noticeable for the other data sets. There is a sound explanation to this phenomenon: simply, the structure formed by fuzzy clustering does not fully reflect the dependencies in the data (due to the effect of the sparsity of the data), and this problem, in turn, results in the deteriorating performance of the fuzzy model. In this case, one would be better off to consider a suitable reduced set of features. In all cases experimented with, we noted an optimal combination of features and data that led to the best performance of the model. Table 7 summarizes the optimal combinations of features and data.Table 7
The optimal % of features and data for different clusters.
Data set and number of clusters% of features% of dataPima withC=37040Pima withC=410020Pima withC=510020Pima withC=610040Pima withC=710080Housing withC=38070Housing withC=48050Housing withC=580100Housing withC=68080Housing withC=790100Body Fat withC=33030Body Fat withC=410070Body Fat withC=59070Body Fat withC=69090Parkinson withC=33030Heat map for PM 10 data forc varying in-between 3 to 6.
(a)(b)(c)(d)Heat map for Body fat data for varying in-between 3 to 6.
(a)(b)(c)(d)Heat map for housing data forc=3, 4, 5, and 7.
(a)(b)(c)(d)The relationships between the percentage of data used and the resulting RMSE values are displayed in Figures9 and 10. Some interesting tendencies are worth noting. A critical number of data are required to form a fuzzy model. Increasing the number of data does not produce any improvement as the curves plotted on Figures 9(a), 9(c), and 9(a) achieve a plateau or even some increase of the RMSE is noticeable.The values of RMSE versus the percentage of data for selected number clusters: (a) housing data, (b) PM10 data, and (c) body fat data.
(a)(b)(c)Plots of RMSE versus the percentage of features for selected number clusters: (a) housing data, (b) PM10 data, and (c) body Fat data.
(a)(b)(c)Considering a fixed percentage of the data used, we look at the nature of the feature sets. Tables8, 9 and 10 displays the best feature for PM10 data, Body fat data, and Housing data, respectively.  Overall, the selected subsets of features are almost the same for different numbers of the clusters being used. Furthermore, we observe that in most cases, the reduced feature spaces exhibit an interesting “nesting” property, meaning that the extended feature space constructed subsumes the one formed previously. For example, for the Housing data, we obtain the following subsets of features:
(6){feature6}⊂{feature6,feature9,feature13}⊂{feature6,feature9,feature10,feature13}.
Here, the corresponding features are as follows: 6: average number of rooms per dwelling, 9: index of accessibility to radial highways, 13: percentage of lower status population, and 10: full-value property-tax rate per $10,000. This combination is quite convincing.Table 8
Best subsets of features for PM10 data.
FSubset of featuresD = 30%D = 40%D = 50%D = 60%D = 70%D = 80%D = 90%10%111111120%111111130%1, 71, 71, 71, 71, 61, 61, 640%1, 6, 71, 6, 71, 6, 71, 6, 71, 6, 71, 6, 71, 6, 750%1, 2, 6, 71, 2, 6, 71, 2, 6, 71, 2, 6, 71, 2, 6, 71, 2, 6, 71, 2, 6, 760%1, 2, 6, 71, 2, 6, 71, 2, 6, 71, 2, 6, 71, 2, 6, 71, 2, 6, 71, 2, 6, 770%1, 2, 3, 6, 71, 2, 3, 6, 71, 2, 3, 6, 71, 2, 3, 6, 71, 2, 3, 6, 71, 2, 4, 6, 71, 2, 4, 6, 780%1, 2, 3, 4, 6, 71, 2, 3, 5, 6, 71, 2, 3, 4, 6, 71, 2, 3, 4, 6, 71, 2, 3, 4, 6, 71, 2, 3, 4, 6, 71, 2, 3, 4, 6, 790%1, 2, 3, 4, 6, 71, 2, 3, 4, 6, 71, 2, 3, 4, 6, 71, 2, 3, 5, 4, 71, 2, 3, 4, 5, 71, 2, 3, 4, 6, 71, 2, 3, 5, 4, 7Table 9
Best subsets of features for body fat data.
FSubset of featuresD = 50%D = 60%D = 70%D = 80%10%111120%1, 6, 71, 3, 71, 3, 71, 3, 730%1, 6, 7, 91, 7, 8, 91, 3, 7, 91, 3, 7, 940%1, 6, 7, 8, 9, 121, 3, 4, 7, 8, 91, 3, 6, 7, 8, 91, 7, 8, 9, 11, 1250%1, 2, 6, 7, 8, 12, 141, 4, 8, 9, 11, 12, 141, 3, 7, 8, 9, 12, 141, 3, 4, 5, 7, 8, 1260%1, 2, 6, 7, 8, 11, 12, 141, 3, 4, 7, 8, 9, 11, 141, 3, 5, 7, 9, 11, 12, 141, 3, 5, 7, 8, 11, 12, 1470%1, 3, 4, 5, 6, 8, 9, 11, 12, 141, 3, 4, 5, 7, 8, 9, 11, 12, 141, 3, 4, 5, 7, 8, 9, 11, 12, 141, 3, 4, 5, 7, 8, 9, 11, 12, 1480%1, 3, 4, 5, 6, 7, 8, 9, 10, 12, 141, 3, 4, 5, 6, 7, 8, 9, 10, 12, 141, 3, 4, 5, 6, 7, 8, 9, 11, 12, 141, 3, 4, 5, 6, 7, 8, 9, 11, 12, 1490%1, 2, 3, 4, 5, 6, 8, 9, 10, 11, 12, 13, 141, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 141, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 141, 23, 4, 5, 6, 7, 8, 9, 10, 11, 12, 14Table 10
Best subsets of features for housing data.
FSubset of featuresD = 10%D = 30%D = 50%D = 70%D = 90%10%1266121220%5, 12, 136, 9, 136, 9, 106, 9, 106, 9, 1030%6, 7, 9, 135, 6, 9, 114, 6, 9, 131, 6, 9, 102, 3, 12, 1340%1, 3, 6, 10, 134, 6, 9, 10, 136, 9, 10, 12, 133, 5, 6, 9, 101, 4, 6, 9, 1050%3, 6, 7, 8, 9, 10, 111, 2, 6, 9, 10, 11, 131, 3, 5, 6, 8, 9, 111, 5, 6, 9, 10, 11, 121, 6, 9, 10, 11, 12, 1360%3, 5, 6, 7, 8, 9, 12, 133, 5, 6, 7, 8, 9, 10, 111, 3, 5, 6, 7, 9, 10, 111, 3, 5, 6, 7, 9, 10, 111, 3, 6, 7, 8, 9, 11, 1370%3, 5, 6, 7, 8, 9, 10, 12, 131, 3, 4, 5, 6, 9, 10, 11, 131, 3, 5, 6, 7, 9, 10, 11, 131, 3, 5, 6, 7, 9, 10, 11, 131, 3, 6, 7, 8, 9, 10, 11, 1380%1, 3, 4, 5, 6, 8, 9, 10, 11, 131, 3, 5, 6, 8, 9, 10, 11, 12, 131, 3, 5, 6, 7, 8, 9, 10, 11, 131, 3, 5, 6, 7, 8, 9, 10, 11,  131, 3, 5, 6, 7, 8, 9, 10, 11, 1390%1, 2, 3, 4, 5, 6, 7, 8, 10, 11, 12, 131, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 131, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 131, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 131, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13For the PM10 data, we arrive at a series of nested collections of features:(7){feature1}⊂{feature1,feature7}⊂{feature1,feature6,feature7}⊂{feature1,feature2,feature6,feature7},
where the corresponding features include: 1: the concentration of PM10 (particles), 7: hour of experiment per day, 6: wind direction, and 2: the number of cars per hour.Turning to the comparative analysis of performance of the swarm optimization methods, we summarize the obtained results in Figure11. For all data, the CPSO performed better than the standard PSO. Although both algorithms show the same tendency when the percentage of feature is 100% however, the RMSE produced by the CPSO is lower than the one obtained when running the PSO. Furthermore, the CPSO algorithm is more stable than the standard PSO. In most cases, the standard deviations of error produced by the CPSO are smaller than the results obtained for the standard PSO (see Table 11).Table 11
Standard deviations for PSO and CPSO (housing and PM10 data sets).
Housing (D = 50%)PM10 (D = 50%)PSOCPSOPSOCPSO0.0660.0720.0210.0070.1920.0150.0180.0090.1990.0390.0100.0070.110.0930.0040.0070.1150.0790.0240.0080.0910.0710.0060.0090.0580.0940.0140.0100.0440.0640.0100.0070.0530.0420.0190.0090.0210.0420.0160.009Values of RMSE versus the percentage of features selected when running PSO and CPSO—the use of the housing dataset: (a) 20% of selected data, (b) 30% of selected data, (c) 50% of selected data, and (d) 70% of selected data.
(a)(b)(c)(d)Figure12 shows the subsets of the features selected for different percentages of the features used in construction of the fuzzy model. The CPSO algorithm is more consistent while selecting the increasing number of features. For example, features 6 and 13 were selected when using both 30% and 70% of data. In contrast to the selection made with the PSO algorithm, the subset of the features selected here is not as stable, especially when using only 30% of data.Comparison of sets of features being selected by using PSO and CPSO2 for Housing data dataset: (a) PSO method with 30% of selected data, (b) CPSO2 method with 30% of selected data, (c) PSO method with 70% of selected data, and (d) CPSO2 method with 70% of selected data.
(a)(b)(c)(d)Table12 presents the percentage of the improvement when using the CPSO algorithm compared to the PSO algorithm. Note that in this percentage we included all different combinations of the features’ percentages and the data percentages being used. The percentage of the improvement is higher when dealing with a smaller percentage of features and data. For example, the percentage of improvement is 34% for 10% of the instances and 50% of the features selected while the percentage of improvement is less than 10% for 60% of instances and features used. These results occurred because the PSO method has to deal with a large search space for selecting a small subset of features and instances. In contrast to the search space for CPSO, the large search space is decomposed into multiple subswarms that reduce the dimensionality of the original search space.Table 12
Percentage of improvement of the RMSE obtained when using CPSO over the results formed by the PSO; Housing data set.
FD = 10%D = 20%D = 30%D = 40%D = 50%D = 60%D = 70%D = 80%D = 90%10%13152026191810191120%31312024251624231730%33261920141524182140%28211914171717181950%342112799912860%23211374789770%19171449648880%33191242338990%22144521221100%1744731111Tables13, 14, 15, and 16 show the comparison of RMSE when using the proposed method and the standard fuzzy modeling method. Here the standard fuzzy model is constructed without using any feature and instances selection and the holdout method is used to select the data based on the percentage given. The experiment for using the standard fuzzy modeling is repeated for 50 times. If we analyze the tables, we can observe that our proposed method outperforms the standard method of constructing the fuzzy model from the dataset. This can be seen clearly when using the CPSO method to search for the best subset of feature and instances. For example, in Table 13 if we use the CPSO method, the RMSE for using 70% of data is 3.413, whereas the RMSE for the standard method is 8.312. The same tendency occurs for all datasets used here.Table 13
The comparison of RMSE obtained when using standard PSO, CPSO, and standard fuzzy model with holdout method for housing data withC=3.
% of dataStandard PSOCooperative PSO1Holdout method% of featureMSE% of featureMSE% of featureMSE30904.015403.47310017.59340803.699703.46410010.80350803.573703.4141009.90760803.556703.4351008.50770803.527603.4131008.31280803.654903.4491008.16490803.679903.6151007.641Table 14
The comparison of RMSE obtained when using standard PSO, CPSO, and standard fuzzy model with holdout method for body fat data withC=3.
% of dataStandard PSOCooperative PSO1Holdout method% of featureMSE% of featureMSE% of featureMSE30304.677304.684710011.58640304.717304.54091008.29150404.617304.51361007.54860504.636404.42891007.07370404.548404.42341006.65880404.553404.42331006.23990404.582404.37711006.102Table 15
The comparison of RMSE obtained when using standard PSO, CPSO, and standard fuzzy model with holdout method for PM10 data withC=3.
% of dataStandard PSOCooperative PSO1Holdout method% of featureMSE% of featureMSE% of featureMSE301000.764800.73381002.100401000.765800.74321002.01850900.777900.77901002.03060800.805800.77691001.98670900.820800.80521002.00180800.839900.82061001.98390700.847900.84171001.976Table 16
The comparison of RMSE obtained when using CPSO and standard fuzzy model with holdout method for computer data withC=3.
% of dataCooperative PSOHoldout method% of featureMSE% of featureMSE304016.44610017.453403014.71210017.524503014.35010017.680604014.93510017.837704016.23710017.918804015.12210018.351Figure13 shows the comparison plot between the proposed method and the “standard” fuzzy modeling. In most of the cases, the proposed method showed better performance.Comparison of RMSE by using proposed method (straight line) and standard fuzzy model (dotted line): (a) Housing dataset withc=4, (b) body fat dataset with c=5, (c) Parkinson’s dataset with c=3, and (d) computer dataset with c=3.
(a)(b)(c)(d)It becomes clear that one is able to reduce the input data in terms of the number features and instances. Moreover, the flexibility of choosing the reduction level helps the user focus on the most essential subsets of data and features (variables). The knowledge acquired about the best subset of data can be used for future data collection. In addition, the user can put more effort analyzing only the best subset of data that give more impact to the overall prediction.


## 6. Conclusions

In this paper, we proposed a simple framework for constructing fuzzy modeling from high-dimensional and large data. This framework has several advantages that make it better suited than other frameworks for sharing various real-life problems. Firstly, the simultaneous feature and instances selection is easily adapted to construct the structure of the fuzzy model. Secondly, the best selected subset of data obtained with this framework is capable of representing the original large data set. Thirdly, we construct an optimal (or suboptimal) collection of features and data based on the PSO. In addition, a cooperative PSO is developed in order to overcome the limitation of using standard PSO when dealing with a high-dimensional search space. The size of the selected features and data used to construct the fuzzy model can be adjusted based upon the feedback provided in terms of the performance of the model constructed for the currently accepted.The effectiveness of the framework was validated by using four well-known regression data sets. The experiment results showed that the proposed fuzzy modeling framework is able to handle high dimensionality and a large data set simultaneously. Moreover, the curse of dimensionality problem in fuzzy modeling was substantially reduced.In the future work one could concentrate on improving the cooperative PSO by fine-tuning the parameters of the method such as, for example, the cognitive and social parameter.



---

*Source: 347157-2012-04-09.xml*