# Ant Colony Optimization with Three Stages for Independent Test Cost Attribute Reduction

Authors: Zilong Xu; Hong Zhao; Fan Min; William Zhu
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

Minimal test cost attribute reduction is an important problem in cost-sensitive learning. Recently, heuristic algorithms including the information gain-based algorithm and the genetic algorithm have been designed for this problem. However, in many cases these algorithms cannot find the optimal solution. In this paper, we develop an ant colony optimization algorithm to tackle this problem. The attribute set is represented as a graph with each vertex corresponding to an attribute and weight of each edge to pheromone. Our algorithm contains three stages, namely, the addition stage, the deletion stage, and the filtration stage. In the addition stage, each ant starts from the initial position and traverses edges probabilistically until the stopping criterion is satisfied.
The pheromone of the traveled path is also updated in this process. In the deletion stage, each ant deletes redundant attributes. Two strategies, called the centralized deletion strategy and the distributed deletion strategy, are proposed. Finally, the ant with minimal test cost is selected to construct the reduct in the filtration stage. Experimental results on UCI datasets indicate that the algorithm is significantly better than the information gain-based one. It also outperforms the genetic algorithm on medium-sized dataset Mushroom.


---

## Body

## 1. Introduction

Cost-sensitive attribute reduction has gained much research interest in rough sets. People have considered three types of costs, namely, test cost [1], misclassification cost [2], and delay cost [3, 4]. Test cost is the money and/or time that is required to obtain attribute values. In real applications, only part of tests is needed to maintain enough information for classification. We would like to choose an attribute reduct [5] that minimizes the total test cost. This issue is called the minimal test cost reduct (MTR) problem [1].We explain the MTR problem through a simple example. There is a medical decision system with five attributes, Mcv, Alkphos, Sgpt, Sgot, and Gammagt. The test costs of these five attributes are $55, $10, $15, $20, and $25, respectively. The decision system has three reducts,{Mcv,Alkphos}, {Alkphos,Sgpt,Gammagt}, and {Alkphos,Sgot,Gammagt}. Alkphos is the core attribute. The minimal reduct is {Mcv,Alkphos} including only two attributes. The minimal test cost reduct is {Alkphos,Sgpt,Gammagt} with only $50.The minimal test cost attribute reduction (MTR) is proposed for saving resources, and the problem is meaningful in applications. The MTR problem is more general than the minimal reduct problem which is NP-hard; hence the MTR problem is at least NP-hard. Consequently, heuristic algorithms are needed to deal with such problem. Heuristic algorithms including information gain-basedλ-weighted reduction algorithm [1] and genetic algorithm [6] have been designed to deal with MTR problem. Unfortunately, they often do not find the optimal solution for medium-sized datasets. Although the competition approach [1] helps to improve the performance through constructing a population of reducts, the results are still unsatisfactory. There is still room to obtain more sophisticated algorithms through other techniques.In this paper, we develop an algorithm based on ant colony optimization for the MTR problem. The attribute set is represented as a complete graph with each vertex corresponding to one attribute. Then batches of ants are generated for attribute subset selection. Our algorithm contains three stages, namely, the addition stage, the deletion stage, and the filtration stage. First, in the addition stage, core vertexes are compressed into one vertex as the initial position of all ants. From the initial position, each ant selects next vertex according to test costs of each adjacent vertex and pheromone of each adjacent edge. If the attribute vertexes an ant has traveled satisfy the positive region condition, the ant stops, otherwise continues to add attributes. Second, in the deletion stage, redundant attributes are deleted from the obtained attribute subsets, and the pheromone of each edge is updated. Two strategies including centralized deletion and distributed deletion are designed for this stage. Third, the ant with the least test cost is selected, and the attribute subset corresponding to its path is output as the result in the filtration stage.To evaluate the performance of algorithms, we adopt three measures, namely, finding optimal factor (FOF), maximal exceeding factor (MEF), and average exceeding factor (AEF) [1]. Experimental results indicate that our algorithm outperforms the information gain-based algorithm in most datasets with different test cost distributions. It can obtain better results than the genetic algorithm except some small datasets. One possible reason is that the ant colony optimization algorithm produces more diverse solutions than the existing ones. The distributed deletion strategy is superior to the centralized one on medium-sized dataset Mushroom.The rest of the paper is organized as follows. Section2 reviews the basic concepts in rough sets and decision system. Section 3 proposes the ant colony optimization to tackle the minimal test cost reduction problem. In Section 4, we present our experiment schemes and show the results. We also give a simple analysis of our experimental results. Finally, Section 5 presents the conclusion.


## 2. Preliminaries

This section reviews basic knowledge: test-cost-independent decision systems, relative reduct, minimal test cost reduct problem, genetic algorithm, and ant colony optimization.


### 2.1. Test-Cost-Independent Decision Systems

Most supervised learning approaches are based on decision systems. A decision system is often denoted byS=(U,C,D,{Va∣a∈C∪D},{Ia∣a∈C∪D}), where U is a finite set of objects called the universe, C is the set of conditional attributes, also called the set of tests, D is the set of decision attributes, also called the decision, Va is the set of values for each a∈C∪D, and Ia:U→Va is an information function for each a∈C∪D. We often denote {Va∣a∈C∪D} and {Ia∣a∈C∪D} by V and I, respectively. A decision system is often represented by a decision table, as shown in Table 1.Table 1
An exemplary decision table.
Patient
Headache
Muscle
Temperature
Snivel
Flu
x
1
No
Yes
Normal
No
No
x
2
Yes
Yes
High
Yes
Yes
x
3
Yes
Yes
Very high
No
Yes
x
4
No
Yes
Normal
No
No
x
5
No
No
High
No
No
x
6
No
Yes
Very high
Yes
YesWe consider the simplest case though most widely used type of cost-sensitive decision systems as follows.Definition 1 (see [7]).
A test-cost-independent decision system (TCI-DS)S is the 6-tuple
(1)S=(U,C,D,V,I,c),
where U,C,D,V, and I have the same meanings as in a decision system and c:C→ℝ+∪{0} is the test cost function. Test costs are independent of one another; that is, c(B)=∑a∈B‍c(a) for any B⊆C.We usually use a vectorc=[c(a1),c(a2),…,c(a|C|)] to represent the cost function. An exemplary cost vector is shown in Table 2. In fact, cost-sensitive decision systems are more general than decision systems. If all elements in c are 0, a TCI-DS coincides with a DS. For simplicity, free tests are not considered in this work. This consideration is reasonable since we always need some cost to obtain data.Table 2
The test cost vector.
a
Headache
Muscle pain
Temperature
Snivel
c
(
a
)
50
25
40
20


### 2.2. The Relative Reduct

The relative reduct is a crucial concept in rough sets, and there are many different definitions, such as positive region reducts [5], maximum distribution reducts [8], fuzzy reducts [9], and β-reduct [10]. These definitions are equivalent if the decision table is consistent. Furthermore, there are some extended concepts such as dynamic reducts [11], parallel reducts [12, 13], and M-relative reducts [14]. In some extensions of rough sets, that is, covering-based rough sets [15], there are other definitions of a reduct (see, e.g., [16–18]).Most existing reduct problems aim at finding the minimal description of the data. And the objectives include finding attribute subsets with the minimal size [5, 19], the minimal space [20], or a covering with the minimal number of subsets [16]. Since the test cost issue is the focus of this paper, we are interested in reducts with the minimal test cost. This type of reducts is defined as follows.Definition 2 (see [1]).
LetS be a TCI-DS and Red(S) the set of all reducts of S. Any R∈Red(S), where c(R)=min{c(R′)∣R′∈Red(S)}, is called a minimal test cost reduct.The set of all minimal test cost reducts is denoted byMTR(S). And the problem of constructing MTR(S) is called the minimal test cost reduct (MTR) problem. As indicated in [1], the time complexity of computing MTR(S) is the same as Red(S).


### 2.3. The Minimal Test Cost Reduct Problem

Attribute reduction is a key issue of the rough sets research. The classical [5], covering-based [16, 21], decision-theoretical [3], variable-precision [10], dominance-based [22], and neighborhood [23] rough sets models address the reduction problem from different perspectives. A number of definitions of relative reducts exist [5, 19, 23, 24] for different rough sets models. This paper employs the definition based on the positive region.


### 2.4. The Genetic Algorithm

In the computer science field of artificial intelligence, a genetic algorithm (GA) is a search heuristic that mimics the process of natural evolution [25]. This heuristic (also sometimes called a metaheuristic) is routinely used to generate useful solutions to optimization and search problems. Genetic algorithms belong to the larger class of evolutionary algorithms (EA), which generate solutions to optimization problems using techniques inspired by natural evolution, such as inheritance, mutation, selection, and crossover [25]. In [26], the genetic algorithm is employed to evolve the cost-sensitive decision trees. Recently, the genetic algorithm has been employed to tackle the minimal test cost reduction problem [6] and attribute reduction with test cost constraint [27].


### 2.5. The Ant Colony Optimization

Swarm intelligence is a relatively new approach to problem solving that takes inspiration from the social behaviors of insects and other animals [28]. The ant colony optimization (ACO) algorithm is a probabilistic technique for solving computational problems, which can be reduced to finding good paths through graphs. ACO algorithms are state-of-the-art for the sequential ordering problem [29], the vehicle routing problem with time window constraints [30], the quadratic assignment problem [31], the arc-weighted l-cardinality tree problem [32], and the shortest common supersequence problem [33]. In rough sets, the classical attribute reduction problem has been tackled by the ACO [34, 35].


### 2.6. Evaluation Measures

For evaluating the experiment results, we adopt evaluation measures proposed in [1]. The new algorithm is compared with the information gain-based heuristic algorithm on four UCI datasets [36].We need a measure to evaluate the quality of one particular reduction. Since an algorithm can run on many datasets or one dataset with different test cost settings, we adopt three metrics from a statistical viewpoint. They are finding optimal factor (FOF), maximal exceeding factor (MEF), and average exceeding factor (AEF) [1].


#### 2.6.1. Finding Optimal Factor

Let the number of experiments beK and the number of successful searches of an optimal reduct k. The finding optimal factor is defined as
(2)op=kK.


#### 2.6.2. Exceeding Factor

For a dataset with a particular test cost setting, letR′ be an optimal reduct. The exceeding factor of a reduct R is
(3)ef(R)=c*(R)-c*(R′)c*(R′).The exceeding factor provides a quantitative metric to evaluate the performance of a reduct. It indicates the badness of a reduct when it is not optimal. Naturally, ifR is an optimal reduct, the exceeding factor is 0. To demonstrate the performance of an algorithm, statistical metrics are needed. Let the number of experiments be K. In the ith experiment (1≤i≤K), the reduct computed by the algorithm is denoted by Ri. The maximal exceeding factor is defined as
(4)max1≤i≤Kef(Ri).This shows the worst case of the algorithm given some dataset. Although it relates to the performance of one particular reduct, it should be viewed as a statistical rather than an individual metric. The average exceeding factor is defined as(5)∑i=1Kef(Ri)K.
Since it is averaged on K different test-cost-sensitive decision systems, it shows the overall performance of the algorithm solely from a statistical perspective.


## 2.1. Test-Cost-Independent Decision Systems

Most supervised learning approaches are based on decision systems. A decision system is often denoted byS=(U,C,D,{Va∣a∈C∪D},{Ia∣a∈C∪D}), where U is a finite set of objects called the universe, C is the set of conditional attributes, also called the set of tests, D is the set of decision attributes, also called the decision, Va is the set of values for each a∈C∪D, and Ia:U→Va is an information function for each a∈C∪D. We often denote {Va∣a∈C∪D} and {Ia∣a∈C∪D} by V and I, respectively. A decision system is often represented by a decision table, as shown in Table 1.Table 1
An exemplary decision table.
Patient
Headache
Muscle
Temperature
Snivel
Flu
x
1
No
Yes
Normal
No
No
x
2
Yes
Yes
High
Yes
Yes
x
3
Yes
Yes
Very high
No
Yes
x
4
No
Yes
Normal
No
No
x
5
No
No
High
No
No
x
6
No
Yes
Very high
Yes
YesWe consider the simplest case though most widely used type of cost-sensitive decision systems as follows.Definition 1 (see [7]).
A test-cost-independent decision system (TCI-DS)S is the 6-tuple
(1)S=(U,C,D,V,I,c),
where U,C,D,V, and I have the same meanings as in a decision system and c:C→ℝ+∪{0} is the test cost function. Test costs are independent of one another; that is, c(B)=∑a∈B‍c(a) for any B⊆C.We usually use a vectorc=[c(a1),c(a2),…,c(a|C|)] to represent the cost function. An exemplary cost vector is shown in Table 2. In fact, cost-sensitive decision systems are more general than decision systems. If all elements in c are 0, a TCI-DS coincides with a DS. For simplicity, free tests are not considered in this work. This consideration is reasonable since we always need some cost to obtain data.Table 2
The test cost vector.
a
Headache
Muscle pain
Temperature
Snivel
c
(
a
)
50
25
40
20


## 2.2. The Relative Reduct

The relative reduct is a crucial concept in rough sets, and there are many different definitions, such as positive region reducts [5], maximum distribution reducts [8], fuzzy reducts [9], and β-reduct [10]. These definitions are equivalent if the decision table is consistent. Furthermore, there are some extended concepts such as dynamic reducts [11], parallel reducts [12, 13], and M-relative reducts [14]. In some extensions of rough sets, that is, covering-based rough sets [15], there are other definitions of a reduct (see, e.g., [16–18]).Most existing reduct problems aim at finding the minimal description of the data. And the objectives include finding attribute subsets with the minimal size [5, 19], the minimal space [20], or a covering with the minimal number of subsets [16]. Since the test cost issue is the focus of this paper, we are interested in reducts with the minimal test cost. This type of reducts is defined as follows.Definition 2 (see [1]).
LetS be a TCI-DS and Red(S) the set of all reducts of S. Any R∈Red(S), where c(R)=min{c(R′)∣R′∈Red(S)}, is called a minimal test cost reduct.The set of all minimal test cost reducts is denoted byMTR(S). And the problem of constructing MTR(S) is called the minimal test cost reduct (MTR) problem. As indicated in [1], the time complexity of computing MTR(S) is the same as Red(S).


## 2.3. The Minimal Test Cost Reduct Problem

Attribute reduction is a key issue of the rough sets research. The classical [5], covering-based [16, 21], decision-theoretical [3], variable-precision [10], dominance-based [22], and neighborhood [23] rough sets models address the reduction problem from different perspectives. A number of definitions of relative reducts exist [5, 19, 23, 24] for different rough sets models. This paper employs the definition based on the positive region.


## 2.4. The Genetic Algorithm

In the computer science field of artificial intelligence, a genetic algorithm (GA) is a search heuristic that mimics the process of natural evolution [25]. This heuristic (also sometimes called a metaheuristic) is routinely used to generate useful solutions to optimization and search problems. Genetic algorithms belong to the larger class of evolutionary algorithms (EA), which generate solutions to optimization problems using techniques inspired by natural evolution, such as inheritance, mutation, selection, and crossover [25]. In [26], the genetic algorithm is employed to evolve the cost-sensitive decision trees. Recently, the genetic algorithm has been employed to tackle the minimal test cost reduction problem [6] and attribute reduction with test cost constraint [27].


## 2.5. The Ant Colony Optimization

Swarm intelligence is a relatively new approach to problem solving that takes inspiration from the social behaviors of insects and other animals [28]. The ant colony optimization (ACO) algorithm is a probabilistic technique for solving computational problems, which can be reduced to finding good paths through graphs. ACO algorithms are state-of-the-art for the sequential ordering problem [29], the vehicle routing problem with time window constraints [30], the quadratic assignment problem [31], the arc-weighted l-cardinality tree problem [32], and the shortest common supersequence problem [33]. In rough sets, the classical attribute reduction problem has been tackled by the ACO [34, 35].


## 2.6. Evaluation Measures

For evaluating the experiment results, we adopt evaluation measures proposed in [1]. The new algorithm is compared with the information gain-based heuristic algorithm on four UCI datasets [36].We need a measure to evaluate the quality of one particular reduction. Since an algorithm can run on many datasets or one dataset with different test cost settings, we adopt three metrics from a statistical viewpoint. They are finding optimal factor (FOF), maximal exceeding factor (MEF), and average exceeding factor (AEF) [1].


### 2.6.1. Finding Optimal Factor

Let the number of experiments beK and the number of successful searches of an optimal reduct k. The finding optimal factor is defined as
(2)op=kK.


### 2.6.2. Exceeding Factor

For a dataset with a particular test cost setting, letR′ be an optimal reduct. The exceeding factor of a reduct R is
(3)ef(R)=c*(R)-c*(R′)c*(R′).The exceeding factor provides a quantitative metric to evaluate the performance of a reduct. It indicates the badness of a reduct when it is not optimal. Naturally, ifR is an optimal reduct, the exceeding factor is 0. To demonstrate the performance of an algorithm, statistical metrics are needed. Let the number of experiments be K. In the ith experiment (1≤i≤K), the reduct computed by the algorithm is denoted by Ri. The maximal exceeding factor is defined as
(4)max1≤i≤Kef(Ri).This shows the worst case of the algorithm given some dataset. Although it relates to the performance of one particular reduct, it should be viewed as a statistical rather than an individual metric. The average exceeding factor is defined as(5)∑i=1Kef(Ri)K.
Since it is averaged on K different test-cost-sensitive decision systems, it shows the overall performance of the algorithm solely from a statistical perspective.


## 2.6.1. Finding Optimal Factor

Let the number of experiments beK and the number of successful searches of an optimal reduct k. The finding optimal factor is defined as
(2)op=kK.


## 2.6.2. Exceeding Factor

For a dataset with a particular test cost setting, letR′ be an optimal reduct. The exceeding factor of a reduct R is
(3)ef(R)=c*(R)-c*(R′)c*(R′).The exceeding factor provides a quantitative metric to evaluate the performance of a reduct. It indicates the badness of a reduct when it is not optimal. Naturally, ifR is an optimal reduct, the exceeding factor is 0. To demonstrate the performance of an algorithm, statistical metrics are needed. Let the number of experiments be K. In the ith experiment (1≤i≤K), the reduct computed by the algorithm is denoted by Ri. The maximal exceeding factor is defined as
(4)max1≤i≤Kef(Ri).This shows the worst case of the algorithm given some dataset. Although it relates to the performance of one particular reduct, it should be viewed as a statistical rather than an individual metric. The average exceeding factor is defined as(5)∑i=1Kef(Ri)K.
Since it is averaged on K different test-cost-sensitive decision systems, it shows the overall performance of the algorithm solely from a statistical perspective.


## 3. The Algorithm

In this section, one hand, we revise the classical ant colony optimization. On the other hand, we present our ant colony optimization with different techniques to tackle the minimal test cost reduct problem.


### 3.1. The Problem Representation and Algorithm Framework

In order to apply ant colony optimization to tackle the minimal test cost reduction problem, we adopt the following model.Graph. The decision system is represented as a graph.Vertex. An attribute is represented as a vertex. Each feature vertex has information about test cost.Edge. There is an edge between any two vertexes. Each edge has information on pheromone density.Adjacent Matrix. The values of the matrix represent pheromone density of each edge.Because our objective is attribute reduction, not classification which produces rule sets, we do not adopt the tree structure such as the ant colony decision tree. In this section, we employ the first, simplified ant colony optimization—AS. We represent the general algorithm framework as follows.Stage 1 (addition stage).Compute the core of the dataset. K batches of ants are created with each batch containing k ants, therefore giving a total of K×k ants. Each ant takes the core as the starting position. From the initial positions, each ant traverses edges probabilistically until the stopping criterion is satisfied.Stage 2 (deletion stage).Delete redundant attributes and update the pheromone of the traveled path.Stage 3 (filtration stage).Gather the attribute subsets obtained by ants, and compute their test cost. Choose the attribute subset with minimal test cost, and output it as the result.Throughout the paper, the stopping criterion is the positive region condition. The selecting probability depends on the test cost of each adjacent attribute vertex and the pheromone of each adjacent edge. The probabilistic transition rule is(6)pijk=[τj]α[ηij]β∑l∈𝒩ik[τl]α[ηil]βj∈𝒩ik,
where k is the number of the ant, τj=1000/tc(aj) as the heuristic information, ηij means the pheromone density of the edge (i,j), α is the exponent of τj, β is the exponent of ηij, and 𝒩ik is the set of all unvisited adjacent vertexes of the attribute i.The difference between the centralized deletion strategy and the distributed strategy is the time to delete attributes. The deletion of redundant attributes follows the finish of the journey of all ants in the case of the centralized deletion strategy. When using the distributed deletion strategy, the algorithm deletes redundant attributes after each ant travels a path, and then it adjusts the path. In this situation, the pheromone of edges adjoining to redundant attribute vertexes will not be updated. The adjusting method is illustrated in Figure2. Monte Carlo methods are a broad class of computational algorithms that rely on repeated random sampling to obtain numerical results that is by running simulations many times over in order to calculate those same probabilities heuristically just like actually playing and recording your results in a real casino situation: hence the name. We employ Monte Carlo method to simulate the process that an artificial ant selects the next attribute vertex probabilistically.


### 3.2. The ACO with Distributed Deletion

We represent the substantial algorithm of the ant colony optimization with distributed deletion strategy as Algorithm1. The algorithm with centralized deletion strategy is similar.Algorithm 1:ACO with distributed deletion to MTR problem.
Input: S=(U,C,D,V,I,c)
Output: M, a reduct of S
Method: acomtr
(1)M=∅;//the best reduct obtained by ants
(2)for (i=0; i<K; i ++) do
(3)for (j=0; j<k; j++)  do
(4)    create ant[i×k+j];
(5)Bi×k+j=core(S);//takes the set of core attributes
(6)end for
//Each ant selects attributes
(7)for (each m∈[i×k,(i+1)×k] where ant[m] has not stopped) do
(8)   MonteCarlo(tc,ph);//select the next vertex
(9)if (POSBm(D)==POSC(D)) then
(10)    ant[m] stops; //Delete redundant vertexes.
(11)mtc=10000;//minimal test cost
(12)for  (j=|Bi|-1; j>=|Core(S)|; j- -) do
(13)if (POS(Bi-{aj})(D)==POSBi(D)) then
(14)Bi=Bi-{aj};
(15)end if
(16)end for
(17)    Adjust the traveled path
(18)    Update the pheromone of edges ant[m] traveled
(19)end if
(20)end for
(21)end for
//Choose the minimal test cost reduction from the lastn ants
(22)for (i=K×k-n; i<=K×k-1; i ++) do
(23)if  tc(Bi)≤mtc  then
(24)mtc=tc(Bi);
(25)M=Bi;
(26)end if
(27)end for
(28) returnM;We explain the algorithm in details as follows. Lines 1 through 8 correspond to Stage 1 in the algorithm framework. In lines 2 through 4,K batches of ants are created, with each batch containing k ants, therefore giving a total of K×k ants. Each ant takes the core, as indicated by line 5. Lines 9 through 21 represent Stage 2. Select vertexes until all ants in the batch meet the positive region condition. Whenever an ant stops, it removes the redundant attributes and adjusts the traveled path. In line 18, the ants release pheromone to the adjusted paths. After all ants finish their journey, select the best reduct as lines 22 through 27 have shown. Obviously, these lines relate to Stage 3. If we adopt the centralized deletion strategy, the algorithm deletes redundant attributes after all ants finish their journey.


### 3.3. A Running Example

The TCI-DS is given by Tables1 and 2. We illustrate the process of an ant in Figure 1, where attributes Headache, Muscle, Temperature, and  Snivel are represented by 0, 1, 2, and 3, respectively.(a) Start from the attribute Headache. (b) Select second attribute. (c) Select third attribute. (d) Stop.
(a)
(b)
(c)
(d)(a) The path which an ant traveled, (b) Delete the redundant attribute, (c) Reconstruct the path.
(a)
(b)
(c)After an ant travels a path, the redundant attributes are removed. After deletion, the algorithm reconstructs the path. This adjusting process is represented in Figure2.Each ant is placed on an attribute randomly. For simplicity, assumek=1; namely, there is one artificial ant in one batch; Consider Core(S)=∅; hence any ant initially takes no vertex.Stage 1 (addition stage)Stage 1.1 (artificial ants start). Assume ant[0] starts from the attribute[0].Stage 1.2 (artificial ants select attributes)Iteration 1.Ant[0] is at the attribute Headache. Ant[0] has three attributes to select, and they are Muscle pain, Temperature,and Snivel.According to (6), the denominator of selecting probability is (1000/25)2·12+(1000/40)2·12+(1000/20)2·12=1600+625+2500=4725.We compute the probability of three selections:(7)p01=1600/4725=0.339,p02=6254725=0.132,p03=25004725=0.529.Becausep03>p01>p02, the ant will select this as next feature with high probability. That is to say, the ant does not select the attribute necessarily. In this case, we assume the ant choose the feature Snivel. At the same time, the ant will add the selected attribute to the visited attribute set. We find the selected attribute set the ant has visited does not satisfy the positive region condition, so it will continue to select next attribute.Iteration 2.Ant[0] is at the attribute Snivel currently.The candidate features areMuscle pain and Temperature. According to the probability function, the denominator of selecting probability is (1000/25)2·12+(1000/40)2·12=1600+625=2225.Consider(8)p31=16002225=0.719,p32=6252225=0.281.More overp31>p32; then the artificial ant[0] almost selects the attribute Muscle pain.Now, the attribute set obtained by ant[0] has satisfied the positive region constraint. The ant[0] stops working and updates pheromone information.Stage 1.3 (pheromone updating).  After an artificial ant stops working, it will update the pheromone density of edges traveled by itself. In our algorithm, when one path is crossed by an ant, the pheromone diffusion will increase by one. Of course, the rule of pheromone updating may be designed in other methods. We adopt the way in this paper.In this example, the sequence attribute selection of ant[0] is Headache, Snivel, and Muscle pain. The pheromone of edge (0,3) and edge (3,1) increases by one.Each artificial ant runs in above three stages.Stage 2 (deletion stage).After all ants have stopped working, the algorithm will delete redundant attributes from the selected attribute sets of last few ants using positive region. If an attribute does not contribute to the positive region, we remove it.Stage 3 (filtration stage).After deletion, we choose the reduct obtained from last n ants with minimal test cost as the result.In this example, we adopt the centralized deletion strategy. Figure2 illustrates the distributed deletion strategy in an iteration. In the figure, the ant travels the paths 0, 1, 2, and 3. Suppose the attribute 1 is a redundant attribute, we remove it from the attribute subset and adjust the path. The adjusting process is shown in Figure 2(c).


## 3.1. The Problem Representation and Algorithm Framework

In order to apply ant colony optimization to tackle the minimal test cost reduction problem, we adopt the following model.Graph. The decision system is represented as a graph.Vertex. An attribute is represented as a vertex. Each feature vertex has information about test cost.Edge. There is an edge between any two vertexes. Each edge has information on pheromone density.Adjacent Matrix. The values of the matrix represent pheromone density of each edge.Because our objective is attribute reduction, not classification which produces rule sets, we do not adopt the tree structure such as the ant colony decision tree. In this section, we employ the first, simplified ant colony optimization—AS. We represent the general algorithm framework as follows.Stage 1 (addition stage).Compute the core of the dataset. K batches of ants are created with each batch containing k ants, therefore giving a total of K×k ants. Each ant takes the core as the starting position. From the initial positions, each ant traverses edges probabilistically until the stopping criterion is satisfied.Stage 2 (deletion stage).Delete redundant attributes and update the pheromone of the traveled path.Stage 3 (filtration stage).Gather the attribute subsets obtained by ants, and compute their test cost. Choose the attribute subset with minimal test cost, and output it as the result.Throughout the paper, the stopping criterion is the positive region condition. The selecting probability depends on the test cost of each adjacent attribute vertex and the pheromone of each adjacent edge. The probabilistic transition rule is(6)pijk=[τj]α[ηij]β∑l∈𝒩ik[τl]α[ηil]βj∈𝒩ik,
where k is the number of the ant, τj=1000/tc(aj) as the heuristic information, ηij means the pheromone density of the edge (i,j), α is the exponent of τj, β is the exponent of ηij, and 𝒩ik is the set of all unvisited adjacent vertexes of the attribute i.The difference between the centralized deletion strategy and the distributed strategy is the time to delete attributes. The deletion of redundant attributes follows the finish of the journey of all ants in the case of the centralized deletion strategy. When using the distributed deletion strategy, the algorithm deletes redundant attributes after each ant travels a path, and then it adjusts the path. In this situation, the pheromone of edges adjoining to redundant attribute vertexes will not be updated. The adjusting method is illustrated in Figure2. Monte Carlo methods are a broad class of computational algorithms that rely on repeated random sampling to obtain numerical results that is by running simulations many times over in order to calculate those same probabilities heuristically just like actually playing and recording your results in a real casino situation: hence the name. We employ Monte Carlo method to simulate the process that an artificial ant selects the next attribute vertex probabilistically.


## 3.2. The ACO with Distributed Deletion

We represent the substantial algorithm of the ant colony optimization with distributed deletion strategy as Algorithm1. The algorithm with centralized deletion strategy is similar.Algorithm 1:ACO with distributed deletion to MTR problem.
Input: S=(U,C,D,V,I,c)
Output: M, a reduct of S
Method: acomtr
(1)M=∅;//the best reduct obtained by ants
(2)for (i=0; i<K; i ++) do
(3)for (j=0; j<k; j++)  do
(4)    create ant[i×k+j];
(5)Bi×k+j=core(S);//takes the set of core attributes
(6)end for
//Each ant selects attributes
(7)for (each m∈[i×k,(i+1)×k] where ant[m] has not stopped) do
(8)   MonteCarlo(tc,ph);//select the next vertex
(9)if (POSBm(D)==POSC(D)) then
(10)    ant[m] stops; //Delete redundant vertexes.
(11)mtc=10000;//minimal test cost
(12)for  (j=|Bi|-1; j>=|Core(S)|; j- -) do
(13)if (POS(Bi-{aj})(D)==POSBi(D)) then
(14)Bi=Bi-{aj};
(15)end if
(16)end for
(17)    Adjust the traveled path
(18)    Update the pheromone of edges ant[m] traveled
(19)end if
(20)end for
(21)end for
//Choose the minimal test cost reduction from the lastn ants
(22)for (i=K×k-n; i<=K×k-1; i ++) do
(23)if  tc(Bi)≤mtc  then
(24)mtc=tc(Bi);
(25)M=Bi;
(26)end if
(27)end for
(28) returnM;We explain the algorithm in details as follows. Lines 1 through 8 correspond to Stage 1 in the algorithm framework. In lines 2 through 4,K batches of ants are created, with each batch containing k ants, therefore giving a total of K×k ants. Each ant takes the core, as indicated by line 5. Lines 9 through 21 represent Stage 2. Select vertexes until all ants in the batch meet the positive region condition. Whenever an ant stops, it removes the redundant attributes and adjusts the traveled path. In line 18, the ants release pheromone to the adjusted paths. After all ants finish their journey, select the best reduct as lines 22 through 27 have shown. Obviously, these lines relate to Stage 3. If we adopt the centralized deletion strategy, the algorithm deletes redundant attributes after all ants finish their journey.


## 3.3. A Running Example

The TCI-DS is given by Tables1 and 2. We illustrate the process of an ant in Figure 1, where attributes Headache, Muscle, Temperature, and  Snivel are represented by 0, 1, 2, and 3, respectively.(a) Start from the attribute Headache. (b) Select second attribute. (c) Select third attribute. (d) Stop.
(a)
(b)
(c)
(d)(a) The path which an ant traveled, (b) Delete the redundant attribute, (c) Reconstruct the path.
(a)
(b)
(c)After an ant travels a path, the redundant attributes are removed. After deletion, the algorithm reconstructs the path. This adjusting process is represented in Figure2.Each ant is placed on an attribute randomly. For simplicity, assumek=1; namely, there is one artificial ant in one batch; Consider Core(S)=∅; hence any ant initially takes no vertex.Stage 1 (addition stage)Stage 1.1 (artificial ants start). Assume ant[0] starts from the attribute[0].Stage 1.2 (artificial ants select attributes)Iteration 1.Ant[0] is at the attribute Headache. Ant[0] has three attributes to select, and they are Muscle pain, Temperature,and Snivel.According to (6), the denominator of selecting probability is (1000/25)2·12+(1000/40)2·12+(1000/20)2·12=1600+625+2500=4725.We compute the probability of three selections:(7)p01=1600/4725=0.339,p02=6254725=0.132,p03=25004725=0.529.Becausep03>p01>p02, the ant will select this as next feature with high probability. That is to say, the ant does not select the attribute necessarily. In this case, we assume the ant choose the feature Snivel. At the same time, the ant will add the selected attribute to the visited attribute set. We find the selected attribute set the ant has visited does not satisfy the positive region condition, so it will continue to select next attribute.Iteration 2.Ant[0] is at the attribute Snivel currently.The candidate features areMuscle pain and Temperature. According to the probability function, the denominator of selecting probability is (1000/25)2·12+(1000/40)2·12=1600+625=2225.Consider(8)p31=16002225=0.719,p32=6252225=0.281.More overp31>p32; then the artificial ant[0] almost selects the attribute Muscle pain.Now, the attribute set obtained by ant[0] has satisfied the positive region constraint. The ant[0] stops working and updates pheromone information.Stage 1.3 (pheromone updating).  After an artificial ant stops working, it will update the pheromone density of edges traveled by itself. In our algorithm, when one path is crossed by an ant, the pheromone diffusion will increase by one. Of course, the rule of pheromone updating may be designed in other methods. We adopt the way in this paper.In this example, the sequence attribute selection of ant[0] is Headache, Snivel, and Muscle pain. The pheromone of edge (0,3) and edge (3,1) increases by one.Each artificial ant runs in above three stages.Stage 2 (deletion stage).After all ants have stopped working, the algorithm will delete redundant attributes from the selected attribute sets of last few ants using positive region. If an attribute does not contribute to the positive region, we remove it.Stage 3 (filtration stage).After deletion, we choose the reduct obtained from last n ants with minimal test cost as the result.In this example, we adopt the centralized deletion strategy. Figure2 illustrates the distributed deletion strategy in an iteration. In the figure, the ant travels the paths 0, 1, 2, and 3. Suppose the attribute 1 is a redundant attribute, we remove it from the attribute subset and adjust the path. The adjusting process is shown in Figure 2(c).


## 4. Experiments

In this section, we try to answer the following questions by experimentation.(1)
How does the number of ants in the ant colony influence the result?(2)
How does the strategy of deleting redundant attributes influence the quality of the result?(3)
Does our ant colony optimization algorithm outperform the existing one?The UCI datasets we used are Zoo, Voting, Tic-tac-toe, and Mushroom. Since these datasets have no test cost settings, for statistical purposes, we apply three common distributions to generate random test cost. The three distributions are uniform distribution, normal distribution, and Pareto distribution. In this paper, the test cost is a random integer ranging from 1 to 100. The exponentβ in (6) is set to be 2, since under many circumstances this value is a good setup [28].We do not design the parameter learning mechanism. The employed competition approach has the selection mechanism of the parameter. Different applications use the same range, and users do not specify the value of the parameter. The strategy is more straightforward than many other parameter tuning strategies. However, we may design other strategies to save time. Finding optimal factor, minimal exceeding factor, and average exceeding factor are employed to measure the effectiveness of the algorithms. We run each algorithm on datasets with 1000 times. The results have statistical characteristics.


### 4.1. The Influence of Ant Counts on Experiment Results

In order to find the relationship between the number of ants and the quality of the result, we conduct this experiment. We run our algorithm with 100 ants, 150 ants, 200 ants … 400 ants. In this section, the number of experiments is set to be 100. The exponentα in (6) is set to be 2. Results are shown in Tables 3 and 4.Table 3
The finding optimal factor of ACO with different numbers of ants without competition method using centralized deletion.
Dataset
Distribution
Number of ants
100
150
200
250
300
350
400
Uniform
0.66
0.70
0.73
0.68
0.62
0.69
0.61
Iris
Normal
0.53
0.45
0.44
0.43
0.37
0.40
0.43
Pareto
0.95
0.93
0.94
0.96
0.87
0.93
0.94
Uniform
0.75
0.67
0.70
0.68
0.70
0.68
0.63
Zoo
Normal
0.66
0.48
0.38
0.29
0.30
0.33
0.34
Pareto
0.99
1.00
0.96
0.94
0.97
0.95
0.84
Uniform
0.90
0.94
0.84
0.80
0.86
0.81
0.83
Voting
Normal
1.00
0.94
0.91
0.95
0.95
0.90
0.93
Pareto
0.98
0.99
0.98
0.99
0.99
0.99
0.94
Uniform
0.96
0.91
0.89
0.78
0.68
0.69
0.68
Tic-tac-toe
Normal
1.00
0.96
0.96
0.96
0.89
0.86
0.83
Pareto
1.00
1.00
1.00
0.99
0.99
0.99
1.00
Uniform
0.75
0.72
0.69
0.64
0.59
0.60
0.49
Mushroom
Normal
0.61
0.55
0.52
0.51
0.39
0.41
0.31
Pareto
0.97
0.97
0.96
0.95
0.95
0.96
0.94Table 4
The finding optimal factor of ACO with different numbers of ants without competition method using distributed deletion.
Dataset
Distribution
Number of ants
100
150
200
250
300
350
400
Uniform
0.72
0.70
0.68
0.67
0.72
0.63
0.72
Iris
Normal
0.42
0.40
0.37
0.40
0.48
0.41
0.42
Pareto
0.90
0.95
0.94
0.93
0.92
0.91
0.97
Uniform
0.72
0.69
0.69
0.66
0.63
0.64
0.55
Zoo
Normal
0.59
0.50
0.39
0.33
0.33
0.36
0.27
Pareto
0.99
0.95
0.93
0.93
0.97
0.91
0.94
Uniform
0.93
0.84
0.87
0.86
0.81
0.79
0.83
Voting
Normal
1.00
1.00
0.97
0.98
0.96
0.96
0.97
Pareto
1.00
0.98
0.99
0.97
0.97
0.97
0.98
Uniform
0.90
0.78
0.74
0.77
0.71
0.74
0.66
Tic-tac-toe
Normal
0.99
0.99
0.92
0.90
0.90
0.91
0.77
Pareto
1.00
1.00
1.00
0.99
1.00
1.00
0.99
Uniform
0.80
0.71
0.70
0.62
0.53
0.59
0.58
Mushroom
Normal
0.62
0.54
0.53
0.34
0.33
0.27
0.29
Pareto
0.99
0.99
0.96
0.95
0.95
0.94
0.94


### 4.2. Comparison with Existing Heuristic Algorithms

According to the result of the above experiments, we find 100 is the optimal setting of the number of the ants. We conduct an empirical study to examine the effect of our algorithm. To improve the performance, we use the competition approach [1] to enhance our algorithm. The exponent α in (6) is set as integers ranging from 1 to 4. We illustrate the results among the information gain-based algorithm, GA-1, GA-2, ACO with centralized deletion and the ACO with distributed deletion, in Table 5. For clarity, the results on dataset Mushroom are shown in Figure 5.Table 5
Results ofλ-information gain algorithm and the test cost based ant colony optimization with centralized deletion and distributed deletion.
Dataset
Distribution
Information gain
GA-1
GA-2
ACO centralized
ACO distributed
FOF
MEF
AEF
FOF
MEF
AEF
FOF
MEF
AEF
FOF
MEF
AEF
FOF
MEF
AEF
Uniform
0.915
0.193
0.004
0.966
0.118
0.001
0.947
0.224
0.002
0.973
0.177
0.001
0.975
0.216
0.001
Zoo
Normal
0.833
0.028
0.001
0.970
0.109
0.001
0.924
0.101
0.002
0.849
0.040
0.002
0.844
0.048
0.002
Pareto
0.999
0.167
0.000
0.998
0.167
0.000
0.993
0.200
0.001
1.000
0.000
0.000
1.000
0.000
0.000
Uniform
0.998
0.053
0.000
1.000
0.000
0.000
0.999
0.009
0.000
0.997
0.065
0.000
1.000
0.000
0.000
Voting
Normal
1.000
0.000
0.000
1.000
0.000
0.000
1.000
0.000
0.000
1.000
0.000
0.000
1.000
0.000
0.000
Pareto
1.000
0.000
0.000
1.000
0.000
0.000
0.995
0.067
0.000
1.000
0.000
0.000
1.000
0.000
0.000
Uniform
0.877
0.055
0.002
1.000
0.000
0.000
1.000
0.000
0.000
1.000
0.000
0.000
1.000
0.000
0.000
Tic-tac-toe
Normal
0.408
0.018
0.003
1.000
0.000
0.000
1.000
0.000
0.000
1.000
0.000
0.000
1.000
0.000
0.000
Pareto
0.986
0.125
0.002
1.000
0.000
0.000
1.000
0.000
0.000
1.000
0.000
0.000
1.000
0.000
0.000
Uniform
0.762
0.523
0.023
0.448
1.000
0.067
0.766
0.543
0.019
0.948
0.085
0.002
0.946
0.144
0.002
Mushroom
Normal
0.176
0.306
0.010
0.321
0.401
0.063
0.551
0.340
0.022
0.958
0.020
0.000
0.966
0.020
0.000
Pareto
0.733
0.500
0.064
0.935
0.400
0.015
0.949
0.500
0.012
1.000
0.000
0.000
1.000
0.000
0.000Our algorithm is tested on the UCI datasets 1000 times, respectively. The new algorithm with different techniques is compared with the information gain-based algorithm [1] and the genetic algorithm [6]. Experimental results are listed in Figures 3, 4, and 5 and Tables 3, 4, and 5.ACO with centralized deletion on different distributions: (a) uniform; (b) normal; (c) pareto.
(a)
(b)
(c)ACO with distributed deletion on different distributions: (a) uniform; (b) normal; (c) Pareto.
(a)
(b)
(c)Result on Mushroom: (a) finding optimal factor; (b) maximal exceeding factor; (c) average exceeding factor.
(a)
(b)
(c)


### 4.3. Experimental Results

Now we can answer the questions proposed at the beginning of this section.(1)
Experiment results indicate that the effect of the algorithm becomes worse with the increment of the ant counts. We find that 100 is the optimal setting of the number of ants. So, we run our algorithm with 100 ants to compare with the existing heuristic algorithms.(2)
From Figures3 and 4, we infer that the distributed deletion outperforms the centralized one on the medium-sized dataset Mushroom without the competition approach. When we use the competition approach, the centralized deletion strategy and the distributed strategy produce the similar quality of results shown in Table 5 and Figure 5.(3)
The ant colony optimization provides better performance than the information gain one and the genetic one. Our algorithm overcomes the shortcoming of the genetic algorithm which performs badly on the larger dataset. On medium-sized dataset Mushroom, the ACO outperforms the genetic one significantly. For example, when the test cost distribution is normal, the FOF of the ant colony optimization using centralized deletion and distributed deletion is95.8%, 96.6%, respectively, outperforming the information gain one, GA-1, and GA-2 with 17.6%, 32.1%,  and 52.1%, respectively. Figure 5 gives a more intuitive understanding of results on the dataset Mushroom. One possible reason is that the ant colony optimization algorithm produces more diverse solutions than the existing one.


## 4.1. The Influence of Ant Counts on Experiment Results

In order to find the relationship between the number of ants and the quality of the result, we conduct this experiment. We run our algorithm with 100 ants, 150 ants, 200 ants … 400 ants. In this section, the number of experiments is set to be 100. The exponentα in (6) is set to be 2. Results are shown in Tables 3 and 4.Table 3
The finding optimal factor of ACO with different numbers of ants without competition method using centralized deletion.
Dataset
Distribution
Number of ants
100
150
200
250
300
350
400
Uniform
0.66
0.70
0.73
0.68
0.62
0.69
0.61
Iris
Normal
0.53
0.45
0.44
0.43
0.37
0.40
0.43
Pareto
0.95
0.93
0.94
0.96
0.87
0.93
0.94
Uniform
0.75
0.67
0.70
0.68
0.70
0.68
0.63
Zoo
Normal
0.66
0.48
0.38
0.29
0.30
0.33
0.34
Pareto
0.99
1.00
0.96
0.94
0.97
0.95
0.84
Uniform
0.90
0.94
0.84
0.80
0.86
0.81
0.83
Voting
Normal
1.00
0.94
0.91
0.95
0.95
0.90
0.93
Pareto
0.98
0.99
0.98
0.99
0.99
0.99
0.94
Uniform
0.96
0.91
0.89
0.78
0.68
0.69
0.68
Tic-tac-toe
Normal
1.00
0.96
0.96
0.96
0.89
0.86
0.83
Pareto
1.00
1.00
1.00
0.99
0.99
0.99
1.00
Uniform
0.75
0.72
0.69
0.64
0.59
0.60
0.49
Mushroom
Normal
0.61
0.55
0.52
0.51
0.39
0.41
0.31
Pareto
0.97
0.97
0.96
0.95
0.95
0.96
0.94Table 4
The finding optimal factor of ACO with different numbers of ants without competition method using distributed deletion.
Dataset
Distribution
Number of ants
100
150
200
250
300
350
400
Uniform
0.72
0.70
0.68
0.67
0.72
0.63
0.72
Iris
Normal
0.42
0.40
0.37
0.40
0.48
0.41
0.42
Pareto
0.90
0.95
0.94
0.93
0.92
0.91
0.97
Uniform
0.72
0.69
0.69
0.66
0.63
0.64
0.55
Zoo
Normal
0.59
0.50
0.39
0.33
0.33
0.36
0.27
Pareto
0.99
0.95
0.93
0.93
0.97
0.91
0.94
Uniform
0.93
0.84
0.87
0.86
0.81
0.79
0.83
Voting
Normal
1.00
1.00
0.97
0.98
0.96
0.96
0.97
Pareto
1.00
0.98
0.99
0.97
0.97
0.97
0.98
Uniform
0.90
0.78
0.74
0.77
0.71
0.74
0.66
Tic-tac-toe
Normal
0.99
0.99
0.92
0.90
0.90
0.91
0.77
Pareto
1.00
1.00
1.00
0.99
1.00
1.00
0.99
Uniform
0.80
0.71
0.70
0.62
0.53
0.59
0.58
Mushroom
Normal
0.62
0.54
0.53
0.34
0.33
0.27
0.29
Pareto
0.99
0.99
0.96
0.95
0.95
0.94
0.94


## 4.2. Comparison with Existing Heuristic Algorithms

According to the result of the above experiments, we find 100 is the optimal setting of the number of the ants. We conduct an empirical study to examine the effect of our algorithm. To improve the performance, we use the competition approach [1] to enhance our algorithm. The exponent α in (6) is set as integers ranging from 1 to 4. We illustrate the results among the information gain-based algorithm, GA-1, GA-2, ACO with centralized deletion and the ACO with distributed deletion, in Table 5. For clarity, the results on dataset Mushroom are shown in Figure 5.Table 5
Results ofλ-information gain algorithm and the test cost based ant colony optimization with centralized deletion and distributed deletion.
Dataset
Distribution
Information gain
GA-1
GA-2
ACO centralized
ACO distributed
FOF
MEF
AEF
FOF
MEF
AEF
FOF
MEF
AEF
FOF
MEF
AEF
FOF
MEF
AEF
Uniform
0.915
0.193
0.004
0.966
0.118
0.001
0.947
0.224
0.002
0.973
0.177
0.001
0.975
0.216
0.001
Zoo
Normal
0.833
0.028
0.001
0.970
0.109
0.001
0.924
0.101
0.002
0.849
0.040
0.002
0.844
0.048
0.002
Pareto
0.999
0.167
0.000
0.998
0.167
0.000
0.993
0.200
0.001
1.000
0.000
0.000
1.000
0.000
0.000
Uniform
0.998
0.053
0.000
1.000
0.000
0.000
0.999
0.009
0.000
0.997
0.065
0.000
1.000
0.000
0.000
Voting
Normal
1.000
0.000
0.000
1.000
0.000
0.000
1.000
0.000
0.000
1.000
0.000
0.000
1.000
0.000
0.000
Pareto
1.000
0.000
0.000
1.000
0.000
0.000
0.995
0.067
0.000
1.000
0.000
0.000
1.000
0.000
0.000
Uniform
0.877
0.055
0.002
1.000
0.000
0.000
1.000
0.000
0.000
1.000
0.000
0.000
1.000
0.000
0.000
Tic-tac-toe
Normal
0.408
0.018
0.003
1.000
0.000
0.000
1.000
0.000
0.000
1.000
0.000
0.000
1.000
0.000
0.000
Pareto
0.986
0.125
0.002
1.000
0.000
0.000
1.000
0.000
0.000
1.000
0.000
0.000
1.000
0.000
0.000
Uniform
0.762
0.523
0.023
0.448
1.000
0.067
0.766
0.543
0.019
0.948
0.085
0.002
0.946
0.144
0.002
Mushroom
Normal
0.176
0.306
0.010
0.321
0.401
0.063
0.551
0.340
0.022
0.958
0.020
0.000
0.966
0.020
0.000
Pareto
0.733
0.500
0.064
0.935
0.400
0.015
0.949
0.500
0.012
1.000
0.000
0.000
1.000
0.000
0.000Our algorithm is tested on the UCI datasets 1000 times, respectively. The new algorithm with different techniques is compared with the information gain-based algorithm [1] and the genetic algorithm [6]. Experimental results are listed in Figures 3, 4, and 5 and Tables 3, 4, and 5.ACO with centralized deletion on different distributions: (a) uniform; (b) normal; (c) pareto.
(a)
(b)
(c)ACO with distributed deletion on different distributions: (a) uniform; (b) normal; (c) Pareto.
(a)
(b)
(c)Result on Mushroom: (a) finding optimal factor; (b) maximal exceeding factor; (c) average exceeding factor.
(a)
(b)
(c)


## 4.3. Experimental Results

Now we can answer the questions proposed at the beginning of this section.(1)
Experiment results indicate that the effect of the algorithm becomes worse with the increment of the ant counts. We find that 100 is the optimal setting of the number of ants. So, we run our algorithm with 100 ants to compare with the existing heuristic algorithms.(2)
From Figures3 and 4, we infer that the distributed deletion outperforms the centralized one on the medium-sized dataset Mushroom without the competition approach. When we use the competition approach, the centralized deletion strategy and the distributed strategy produce the similar quality of results shown in Table 5 and Figure 5.(3)
The ant colony optimization provides better performance than the information gain one and the genetic one. Our algorithm overcomes the shortcoming of the genetic algorithm which performs badly on the larger dataset. On medium-sized dataset Mushroom, the ACO outperforms the genetic one significantly. For example, when the test cost distribution is normal, the FOF of the ant colony optimization using centralized deletion and distributed deletion is95.8%, 96.6%, respectively, outperforming the information gain one, GA-1, and GA-2 with 17.6%, 32.1%,  and 52.1%, respectively. Figure 5 gives a more intuitive understanding of results on the dataset Mushroom. One possible reason is that the ant colony optimization algorithm produces more diverse solutions than the existing one.


## 5. Conclusions

In this paper, we have pointed out the shortcoming of the existing heuristic algorithms including the information gain-based one and the genetic one. We have designed a new algorithm based on ant colony optimization to tackle the MTR problem. In deletion stage, our algorithm contains centralized deletion strategy and distributed deletion strategy. We have tested the new one with three representative test cost distributions on four UCI datasets. Experimental results show that the new algorithm outperforms the existing ones significantly, especially 1`on medium-sized dataset such as Mushroom. Moreover, when distribution is normal, the distributed deletion strategy obtains better results than the centralized one on Mushroom dataset.



---

*Source: 510167-2013-06-06.xml*