Next Article in Journal
Target Image Mask Correction Based on Skeleton Divergence
Next Article in Special Issue
Planning the Schedule for the Disposal of the Spent Nuclear Fuel with Interactive Multiobjective Optimization
Previous Article in Journal
Solving Integer Linear Programs by Exploiting Variable-Constraint Interactions: A Survey
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

SVM-Based Multiple Instance Classification via DC Optimization

by
Annabella Astorino
1,*,†,
Antonio Fuduli
2,†,
Giovanni Giallombardo
3,† and
Giovanna Miglionico
3,†
1
Istituto di Calcolo e Reti ad Alte Prestazioni-C.N.R., 87036 Rende (CS), Italy
2
Dipartimento di Matematica e Informatica, Università della Calabria, 87036 Rende (CS), Italy
3
Dipartimento di Ingegneria Informatica, Modellistica, Elettronica e Sistemistica, Università della Calabria, 87036 Rende (CS), Italy
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 31 October 2019 / Revised: 19 November 2019 / Accepted: 20 November 2019 / Published: 23 November 2019

Abstract

:
A multiple instance learning problem consists of categorizing objects, each represented as a set (bag) of points. Unlike the supervised classification paradigm, where each point of the training set is labeled, the labels are only associated with bags, while the labels of the points inside the bags are unknown. We focus on the binary classification case, where the objective is to discriminate between positive and negative bags using a separating surface. Adopting a support vector machine setting at the training level, the problem of minimizing the classification-error function can be formulated as a nonconvex nonsmooth unconstrained program. We propose a difference-of-convex (DC) decomposition of the nonconvex function, which we face using an appropriate nonsmooth DC algorithm. Some of the numerical results on benchmark data sets are reported.

1. Introduction

Multiple instance learning (MIL) is a recent machine learning paradigm [1,2,3], which consists of classifying sets of points. Each set is called bag, while the points inside the bags are called instances. The main characteristic of an MIL problem is that in the learning phase the instance labels are hidden and only the labels of the bags are known.
An MIL seminal paper is [4], where a drug-design problem has been faced. Such a problem consists of determining whether a drug molecule (bag) is active or non-active. A molecule provides the desired drug effect (positive label) if, and only if, at least one of its conformations (instances) binds to the target site. The crucial question is that it is not known a priori which conformation makes the molecule active.
Some MIL applications are image classification [5,6,7,8], drug discovery [9,10], classification of text documents [11], bankruptcy prediction [12], and speaker identification [13].
For this kind of problems, there are various solutions in the literature that fall into three different classes: instance-space approaches, bag-space approaches, and embedding-space approaches. In instance-space approaches, classification is performed at the instance level, finding a separation surface directly in the instance space, without looking at the global structure of the bags; the label of each bag is determined as an aggregation of the labels of its corresponding instances. Vice-versa, in bag-space approaches (for example, see [14,15,16]), the separation is performed at a global level, considering the bag as a whole entity. A compromise between these two kinds of approaches is constituted by embedding-space techniques, where each bag is represented by one feature vector and the classification is consequently performed in the instance space. An example of an embedding-space approach is presented in [17].
The method we propose uses the instance-space approach and provides a separation hyperplane for the binary case, where the objective is to discriminate between positive and negative bags. We start from the standard MIL assumption stating that a bag is positive if, and only if, at least one of its instances is positive and it is negative whenever all its instances are negative.
Some examples of linear instance-space MIL classifiers can be found in [18,19,20,21,22]. In particular, in [18], two different models have been proposed. The first one is a mixed-integer nonlinear optimization problem solved using a heuristic technique based on the block coordinate descent method [23] and faced in [19] using a Lagrangian relaxation technique. The second model, which will be the objective of our analysis in the next section, is a nonsmooth nonconvex optimization problem, solved in [21] using the bundle type method described in [24]. In [20], a semi-proximal support vector machine (SVM) approach is used, coming from a compromise between the classical SVM [25] and the proximal approach proposed in [26] for supervised classification. Finally, an optimization problem with bilinear constraints is analyzed in [22], where each positive bag is expressed as a convex combination of its instances and a local solution is obtained by solving successive linear programs.
Recently, nonlinear instance-space MIL classifiers have also been proposed in the literature, such as in [27] and in [28], where a spherical separation approach is adopted: in particular, in the former a variable neighborhood search method [29] is used, while in the latter a DC (difference of convex) model is solved using an appropriate DC algorithm [30]. In passing, we stress that many DC models have been introduced in machine learning, in the supervised [31,32,33,34,35], semisupervised [36,37] and unsupervised cases [38,39,40].
In this work, we propose a DC optimization model providing a linear classifier for binary MIL problems. The solution method we adopt is the proximal bundle method introduced in [30] for the minimization of nonsmooth DC functions. The paper is organized as follows. In the next two sections, we describe, respectively, the DC optimization model and the corresponding nonsmooth solution algorithm. Finally, in Section 4, we report the results of our numerical experimentation performed on some data sets drawn from the literature.

2. A DC Decomposition of the SVM-Based MIL

We tackle a binary MIL problem whose objective is to discriminate between m positive bags and k negative ones using a hyperplane
H ( w , b ) { x n | w T x + b = 0 } ,
where w n and b . Indicating by J i + , i = 1 , , m , the index set of the instances belonging to the ith positive bag and by J i , i = 1 , , k , the index set of the instances belonging to the ith negative bag, we recall that, on the basis of the standard MIL assumption, a bag is positive if, and only if, at least one of its instances is positive and it is negative vice-versa. As a consequence, while a positive bag is allowed to, possibly, straddle the hyperplane, the negative bags should lie completely on the negative side.
More formally, indicating by x j n the jth instance of a positive or negative bag, the hyperplane H ( w , b ) performs a correct separation if, and only if, the following conditions hold:
w T x j + b 1 , for   at   least   an   index   j J i + and   for   all   i = 1 , , m w T x j + b 1 , for   all   j J i and   for   all   i = 1 , , k .
As a consequence (see Figure 1 and Figure 2), a positive bag J i + , i = 1 , , m , is misclassified if
max j J i + ( w T x j + b ) < 1
and a negative one J i , i = 1 , , k , is misclassified if
max j J i ( w T x j + b ) > 1 .
Then, we come out with the following error function, already introduced in [18]:
f ( w , b ) 1 2 w 2 + C i = 1 m max { 0 , 1 max j J i + ( w T x j + b ) } + i = 1 k max { 0 , 1 + max j J i ( w T x j + b ) } ,
where C > 0 represents the trade-off between two objectives: the maximization of the separation margin, characterizing the classical SVM [25] approach, and the minimization of the classification error.
To minimize function f, we propose a DC decomposition based on the following formula:
max { 0 , 1 h ( y ) } = max { 1 , h ( y ) } h ( y ) ,
where h is a convex function. By applying Equation (2) to our case, we can write f in the form:
f ( w , b ) = f 1 ( w , b ) f 2 ( w , b ) ,
where
f 1 ( w , b ) 1 2 w 2 + C i = 1 k max { 0 , 1 + max j J i ( w T x j + b ) } + C i = 1 m max { 1 , max j J i + ( w T x j + b ) }
and
f 2 ( w , b ) C i = 1 m max j J i + ( w T x j + b )
are convex functions. Hence, we come up with the following nonconvex nonsmooth optimization problem, DC-MIL,
min w , b [ f 1 ( w , b ) f 2 ( w , b ) ] .

3. Solving DC-MIL Using a Nonsmooth DC Algorithm

We start by recalling some preliminary property of the DC optimization problem, by adopting the same notation as above. Given the DC optimization problem
min y [ f 1 ( y ) f 2 ( y ) ]
where both f 1 and f 2 are convex nonsmooth functions, we say that a point y * is a local minimizer if f 1 ( y * ) f 2 ( y * ) is finite and there exists a neighborhood N of y * such that
f 1 ( y * ) f 2 ( y * ) f 1 ( y ) f 2 ( y ) , y N .
Considering that, in general, the Clarke subdifferential calculus cannot be used to compute subgradients of the DC function since
c l f ( y ) f 1 ( y ) f 2 ( y ) ,
where c l f ( · ) denotes Clarke’s subdifferential, different stationary points can be defined for nonsmooth DC functions. A point y * is called inf-stationary for problem Equation (4) if
f 2 ( y * ) f 1 ( y * ) .
Furthermore, a point y * is called Clarke stationary for problem Equation (4) if
0 c l f ( y * ) ,
while, it is called a critical point of f if
f 2 ( y * ) f 1 ( y * ) .
Denoting the set of inf-stationary points by S i n f , the set of Clarke stationary points by  S c l , and the set of critical points of the function f by S c r , the following inclusions hold
S i n f S c l S c r
as shown in (Proposition 3, [30]).
Nonsmooth DC functions have attracted the interest of several researchers, both, from the theoretical and from the algorithmic viewpoint. Focusing in particular on the algorithmic side, the most relevant contribution has probably been provided by the methods based on the linearization of function f 2 (see, [41] and references therein), where the problem is tackled via successive convexifications of function f. In the last years, nonsmooth-tailored DC programming has experienced a lot of attention as it has a lot of practical applications (see [28,42]). In fact, several nonsmooth DC algorithms have been developed ([30,43,44,45,46,47]).
Here, we adopt the algorithm DCPCA, a bundle-type method introduced in [30] to solve problem Equation (4), which is based on a model function built by combining two convex piecewise approximations, each related to one component function. More in details, a simplified version of DCPCA works as follows:
  • It iteratively builds two separate piecewise-affine approximations of the component functions, grouping the corresponding information in two separate bundles.
  • It combines the two convex piecewise-affine approximations and generates a DC piecewise-affine model.
  • The DC (hence, nonconvex) model is locally approximated using an auxiliary quadratic program, whose solution is used to certify approximate criticality, or to generate a descent search-direction to be explored via a backtracking line-search approach.
  • Whenever no descent is achieved along the search direction, the bundle of the first function is enriched, thus, obtaining a better model function with this being the fundamental feature of any cutting plane algorithm.
In fact, the DCPCA is based on constructing a model function as the pointwise maximum of several concave piecewise-affine pieces. To construct this model, starting from some cutting-plane ideas, the information coming from the two component functions are kept separate in two bundles. We denote the stability center by z (i.e., an estimate of the minimizer), and by I and L, the index sets of the points generated by the algorithm where the information of function f 1 and f 2 have been evaluated, respectively. Therefore, we denote the two bundles of information as
B 1 = { ( g i ( 1 ) , α i ( 1 ) ) : i I }
and
B 2 = { ( g l ( 2 ) , α l ( 2 ) ) : l L }
where, for every i I , g i ( 1 ) f 1 ( y i ) with
α i ( 1 ) = f 1 ( z ) f 1 ( y i ) + g i ( 1 ) T ( z y i ) ,
and, for every l L , g l ( 2 ) f 2 ( y l ) with
α l ( 2 ) = f 2 ( z ) f 2 ( y l ) + g l ( 2 ) T ( z y l ) .
We remark that both component functions, along with their subgradients, could be evaluated at some iterate-point, and, indeed, we assume that ( g ( 1 ) ( z ) , 0 ) B 1 and ( g ( 2 ) ( z ) , 0 ) B 2 , where  g ( 1 ) ( z ) f 1 ( z ) and g ( 2 ) ( z ) f 2 ( z ) .
To approximate the difference function
f 1 ( z + d ) f 2 ( z + d ) f 1 ( z ) f 2 ( z )
at a given iteration k the following nonconvex model function Γ k ( d ) is introduced
Γ k ( d ) max i I min l L g i ( 1 ) g l ( 2 ) T d α i ( 1 ) + α l ( 2 ) ,
which is defined as the maximum of finitely many concave piecewise-affine functions. The model- function Γ k is used to state a sufficient descent condition of the type
f 1 ( z + d ) f 2 ( z + d ) f 1 ( z ) f 2 ( z ) m Γ k ( d )
where m ( 0 , 1 ) . The interesting property of such a model-function is that whenever the sufficient descent is not achieved at points that are close to the stability center, say z + d ¯ , then an improved cutting-plane model can be obtained by only updating the bundle of f 1 with the appropriate information related to the point z + d ¯ . On the other hand, it looks obviously difficult to adopt the minimization of the model-function Γ k as a building block of any algorithm, given its nonconvexity. In fact, DCPCA does not require the direct minimization of Γ k ( d ) , but the search direction can be obtained by solving the following auxiliary quadratic problem:
min d R n , v R v + 1 2 d 2 v ( g i ( 1 ) g l ¯ ( 2 ) ) T d α i ( 1 ) i I Q P ( I )
where l ¯ L ( 0 ) { l L : α l ( 2 ) = 0 } . We observe that L ( 0 ) as B 2 is assumed to contain the information about the current stability center. More precisely, DCPCA works by forcing L ( 0 ) to be a singleton, hence by letting g l ¯ ( 2 ) = g ( 2 ) ( z ) . Denoting the unique optimal solution of Equation (QP(I)) by ( d ¯ , v ¯ ) , a standard duality argument ensures that
d ¯ = i I λ ¯ i ( g i ( 1 ) g l ¯ ( 2 ) )
v ¯ = i I λ ¯ i ( g i ( 1 ) g l ¯ ( 2 ) ) 2 i I λ ¯ i α i ( 1 )
where λ ¯ i 0 , i I , are the optimal variables of the dual of Q P ( I ) , with  i I λ ¯ i = 1 .
Given that any starting point z = y 0 , DCPCA returns an approximate critical point z * , see (Theorem 1, [30]). The following parameters are adopted: the optimality parameter θ > 0 , the subgradient threshold η > 0 , the linearization-error threshold ε > 0 , the approximate line-search parameter m ( 0 , 1 ) , and the step-size reduction parameter σ ( 0 , 1 ) . In Algorithm 1, we report an algorithmic scheme of the main iteration, namely of the set of steps where the stability center is unchanged. An exit from the main iteration is obtained as soon as a stopping criterion is satisfied or whenever the stability center is updated. To make the presentation clearer, without impairing convergence properties, we skip the description of some rather technical steps, which are strictly related to the management of bundle B 2 . Details can be found in [30].
Algorithm 1 DCPCA Main Iteration
1: Solve QP(I) and obtain ( d ¯ , v ¯ ) ▹ Find the search-direction and the predicted-reduction
2: if | v ¯ | θ then▹ Stopping test
3:     set z * = z and exit▹ Return the approximate critical point z *
4: end if
5: Set t = 1 ▹ Start the line-search
6: if f ( z + t d ¯ ) f ( z ) m t v ¯ then▹ Descent test
7:     set z : = z + t d ¯ ▹ Make a serious step
8:     calculate g + ( 1 ) f 1 ( z ) and g + ( 2 ) f 2 ( z ) ▹   |       
9:     update α i ( 1 ) for all i I and α l ( 2 ) for all l L ▹   |       
10:    set B 1 = B 1 \ { ( g i ( 1 ) , α i ( 1 ) ) : α i ( 1 ) > ε , i I } { ( g + ( 1 ) , 0 ) } ▹   |       
11:    set B 2 = B 2 { ( g + ( 2 ) , 0 ) } ▹   |       
12:    update appropriately I and L, and go to 1▹   |       
13: else if t d ¯ > η then▹ Closeness test
14:    set t = σ t and go to 6▹ Reduce the step-size and iterate the line-search
15: end if
16: Calculate g + ( 1 ) f 1 ( z + t d ¯ ) ▹ Make a null step
17: calculate α + ( 1 ) = f 1 ( z ) f 1 ( z + t d ¯ ) + t g + ( 1 ) d ¯ ▹   |    
18: set B 1 = B 1 { ( g + ( 1 ) , α + ( 1 ) ) } , update appropriately I, and go to 1▹   |    
We remark that the stopping condition v ¯ θ , checked at Step 2 of the DCPCA, is an approximate θ -criticality condition for z * . Indeed, taking into account Equation (12), the stopping condition ensures that
i I λ i * g i ( 1 ) g l ¯ ( 2 ) θ and i I λ i * α i ( 1 ) θ ,
which in turn implies that g * ( 1 ) θ f 1 ( z * ) and g * ( 2 ) f 2 ( z * ) such that
g * ( 1 ) g * ( 2 ) 2 θ ,
namely, that
dist θ f 1 ( z * ) , f 2 ( z * ) θ ,
an approximate θ -criticality condition for z * , see Equation (9).

4. Results

We tested the performance of the algorithm DCPCA applied to the DC-MIL formulation (3) by adopting two sets of medium- and large-size problems extracted from [18]. The relevant characteristics of each problem are reported in Table 1 and Table 2, where we list the problem dimension n, the number of instances, and the number of bags.
The two-level cross-validation protocol has been used to tune C and to train the classifier. Before proceeding with the training phase, the model-selection phase is aimed at finding a promising value of parameter C in the set { 2 7 , 2 6 , , 1 , , 2 6 , 2 7 } , using a lower-level cross-validation protocol on each training set. The selected C value, for each training set, is the one returning the highest average test-correctness in the model-selection phase.
Choosing a good starting point is a critical phase to ensure good performance for a local optimization algorithm like DCPCA. For each training set, denoting the barycenter of all the instances belonging to positive bags by w ¯ + and the barycenter of all the instances belonging to negative bags by w ¯ , we have selected the starting point ( w 0 , b 0 ) by setting
w 0 = w ¯ + w ¯
and choosing b 0 such that the corresponding hyperplane correctly classifies all the positive bags.
We adopted the Java implementation of algorithm DCPCA by running the computational experiments on a 3.50 GHz Intel Core i7 computer. We limited the computational budget for every execution of DCPCA to 500 and 200 evaluations of the objective function for medium-size and large-size problems, respectively, and we restricted the size of the bundle to 100 elements adopting a restart strategy, as soon as, the bundle size exceeds the threshold and a new stability center is obtained. The QP solver of IBM ILOG CPLEX 12.8 has been used to solve quadratic subprograms. The following set of parameters, according to the notation introduced in [45], has been selected: the optimality parameter θ = 0.7 , the subgradient threshold η = 0.7 , the approximate linesearch parameter m = 0.01 , the step-size reduction parameter σ = 0.01 , and the linearization-error threshold ϵ = 0.95 .
We compare our DC-MIL approach against the algorithms mi-SVM [18], MI-SVM [18], MICA [22], MIL-RL [19], and for medium-size problems also against the MIC B u n d l e [21] and DC-SMIL [28]. All such methods have been briefly surveyed in the introduction section.
To analyze the reliability of our approach, in Table 3 and Table 4, we report the numerical results in terms of the percentage test-correctness averaged over 10 folds, with the best performance being underlined. We remark that some data are not reported in Table 5 as the corresponding results are obtained by adopting only nonlinear kernels in [18,22]. Moreover, to provide some insight into the efficiency of DC-MIL, we report in Table 5 and Table 6, the average train-correctness (train, %), the average cpu time (cpu, sec), the average number of function evaluations (nF), and the average number of subgradient evaluations of the two functions (nG1 and nG2). The reliability results show a good and balanced performance of the DC-MIL approach equipped with DCPCA, both, for the medium-size problems, where in one case DC-MIL slightly outperforms the other approaches, and for the large-size problems. Moreover, we observe that our approach looks strongly efficient as it manages to achieve high train-correctness in reasonably small execution times even for large-size problems.

5. Conclusions

We have considered a multiple instance learning problem consisting of classifying sets instead of single points. The resulting binary classification problem, addressed by a support vector machine approach, is formulated as an unconstrained nonsmooth optimization problem for which an original DC decomposition is presented. The problem is solved by a proximal bundle-type method, specialized for nonsmooth DC optimization, which is tested on some benchmark datasets against a set of state-of-the-art approaches. The numerical results in terms of reliability show, on one hand, that there are no outperforming methods on all the test problems, on the other hand, that our method achieves comparable performance with other approaches. Moreover, the encouraging results obtained in terms of efficiency show that there is room for improvement by further investigating the parameter settings in relation to specific test problems.

Author Contributions

Methodology, A.A., A.F., G.G., G.M.; software, A.A., A.F., G.G., G.M.; writing–review & editing, A.A., A.F., G.G., G.M.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
MILMultiple instance learning
SVMSupport vector machine
DCDifference of convex

References

  1. Amores, J. Multiple instance classification: Review, taxonomy and comparative study. Artif. Intell. 2013, 201, 81–105. [Google Scholar] [CrossRef]
  2. Carbonneau, M.; Cheplygina, V.; Granger, E.; Gagnon, G. Multiple instance learning: A survey of problem characteristics and applications. Pattern Recognit. 2018, 77, 329–353. [Google Scholar] [CrossRef]
  3. Herrera, F.; Ventura, S.; Bello, R.; Cornelis, C.; Zafra, A.; Sanchez-Tarrago, D.; Vluymans, S. Multiple Instance Learning. Foundations and Algorithms; Springer: Berlin/Heidelberg, Germany, 2016; pp. 1–233. [Google Scholar]
  4. Dietterich, T.G.; Lathrop, R.H.; Lozano-Pérez, T. Solving the multiple instance problem with axis-parallel rectangles. Artif. Intell. 1997, 89, 31–71. [Google Scholar] [CrossRef]
  5. Astorino, A.; Fuduli, A.; Gaudioso, M.; Vocaturo, E. Multiple Instance Learning Algorithm for Medical Image Classification. CEUR Workshop Proceedings. 2019, Volume 2400. Available online: http://ceur-ws.org/Vol-2400/paper-46.pdf (accessed on 25 September 2019).
  6. Astorino, A.; Fuduli, A.; Veltri, P.; Vocaturo, E. Melanoma detection by means of multiple instance learning. Interdiscip. Sci. Comput. Life Sci. 2019. [Google Scholar] [CrossRef] [PubMed]
  7. Astorino, A.; Gaudioso, M.; Fuduli, A.; Vocaturo, E. A multiple instance learning algorithm for color images classification. In ACM International Conference Proceeding Series; ACM: New York, NY, USA, 2018; pp. 262–266. [Google Scholar]
  8. Quellec, G.; Cazuguel, G.; Cochener, B.; Lamard, M. Multiple-Instance Learning for Medical Image and Video Analysis. IEEE Rev. Biomed. Eng. 2017, 10, 213–234. [Google Scholar] [CrossRef]
  9. Fu, G.; Nan, X.; Liu, H.; Patel, R.Y.; Daga, P.R.; Chen, Y.; Wilkins, D.E.; Doerksen, R.J. Implementation of multiple-instance learning in drug activity prediction. BMC Bioinform. 2012, 13. [Google Scholar] [CrossRef]
  10. Zhao, Z.; Fu, G.; Liu, S.; Elokely, K.M.; Doerksen, R.J.; Chen, Y.; Wilkins, D.E. Drug activity prediction using multiple-instance learning via joint instance and feature selection. BMC BioInform. 2013, 14. [Google Scholar] [CrossRef]
  11. Liu, B.; Xiao, Y.; Hao, Z. A selective multiple instance transfer learning method for text categorization problems. Knowl.-Based Syst. 2018, 141, 178–187. [Google Scholar] [CrossRef]
  12. Kotsiantis, S.; Kanellopoulos, D. Multi-instance learning for bankruptcy prediction. In Proceedings of the 2008 Third International Conference on Convergence and Hybrid Information Technology, Busan, Korea, 11–13 November 2008; Volume 1, pp. 1007–1012. [Google Scholar]
  13. Briggs, F.; Lakshminarayanan, B.; Neal, L.; Fern, X.Z.; Raich, R.; Hadley, S.J.K.; Hadley, A.S.; Betts, M.G. Acoustic classification of multiple simultaneous bird species: A multi-instance multi-label approach. J. Acoust. Soc. Am. 2012, 131, 4640–4650. [Google Scholar] [CrossRef]
  14. Gärtner, T.; Flach, P.A.; Kowalczyk, A.; Smola, A.J. Multi-instance kernels. In Proceedings of the 19th International Conference on Machine Learning, Sydney, Australia, 8–12 July 2002; pp. 179–186. [Google Scholar]
  15. Wang, J.; Zucker, J.D. Solving the multiple-instance problem: A lazy learning approach. In Proceedings of the Seventeenth International Conference on Machine Learning, Stanford, CA, USA, 29 June–2 July 2000; Morgan Kaufmann: San Francisco, CA, USA, 2000; pp. 1119–1126. [Google Scholar]
  16. Wen, C.; Zhou, M.; Li, Z. Multiple instance learning via bag space construction and ELM. In Proceedings of the International Society for Optical Engineering, Shanghai, China, 15–17 August 2018; Volume 10836. [Google Scholar]
  17. Wei, X.; Wu, J.; Zhou, Z. Scalable Algorithms for Multi-Instance Learning. IEEE Trans. Neural Netw. Learn. Syst. 2017, 28, 975–987. [Google Scholar] [CrossRef]
  18. Andrews, S.; Tsochantaridis, I.; Hofmann, T. Support vector machines for multiple-instance learning. In Advances in Neural Information Processing Systems; Becker, S., Thrun, S., Obermayer, K., Eds.; MIT Press: Cambridge, UK, 2003; pp. 561–568. [Google Scholar]
  19. Astorino, A.; Fuduli, A.; Gaudioso, M. A Lagrangian relaxation approach for binary multiple instance classification. IEEE Trans. Neural Netw. Learn. Syst. 2019, 30, 2662–2671. [Google Scholar] [CrossRef] [PubMed]
  20. Avolio, M.; Fuduli, A. A semi-proximal support vector machine approach for binary multiple instance learning. 2019; submitted. [Google Scholar]
  21. Bergeron, C.; Moore, G.; Zaretzki, J.; Breneman, C.; Bennett, K. Fast bundle algorithm for multiple instance learning. IEEE Trans. Pattern Anal. Mach. Intell. 2012, 34, 1068–1079. [Google Scholar] [CrossRef] [PubMed]
  22. Mangasarian, O.; Wild, E. Multiple instance classification via successive linear programming. J. Optim. Theory Appl. 2008, 137, 555–568. [Google Scholar] [CrossRef]
  23. Tseng, P. Convergence of a block coordinate descent method for nondifferentiable minimization. J. Optim. Theory Appl. 2001, 109, 475–494. [Google Scholar] [CrossRef]
  24. Fuduli, A.; Gaudioso, M.; Giallombardo, G. Minimizing nonconvex nonsmooth functions via cutting planes and proximity control. SIAM J. Optim. 2004, 14, 743–756. [Google Scholar] [CrossRef]
  25. Vapnik, V. The Nature of the Statistical Learning Theory; Springer: New York, NY, USA, 1995. [Google Scholar]
  26. Fung, G.; Mangasarian, O. Proximal support vector machine classifiers. In Proceedings of the Seventh ACM Sigkdd International Conference on Knowledge Discovery and Data Mining, San Francisco, CA, USA, 26–29 August 2001; Provost, F., Srikant, R., Eds.; ACM: New York, NY, USA, 2001; pp. 77–86. [Google Scholar]
  27. Plastria, F.; Carrizosa, E.; Gordillo, J. Multi-instance classification through spherical separation and VNS. Comput. Oper. Res. 2014, 52, 326–333. [Google Scholar] [CrossRef]
  28. Gaudioso, M.; Giallombardo, G.; Miglionico, G.; Vocaturo, E. Classification in the multiple instance learning framework via spherical separation. Soft Comput. 2019. [Google Scholar] [CrossRef]
  29. Hansen, P.; Mladenović, N.; Moreno Pérez, J.A. Variable neighbourhood search: Methods and applications. 4OR 2008, 6, 319–360. [Google Scholar] [CrossRef]
  30. Gaudioso, M.; Giallombardo, G.; Miglionico, G.; Bagirov, A.M. Minimizing nonsmooth DC functions via successive DC piecewise-affine approximations. J. Glob. Optim. 2018, 71, 37–55. [Google Scholar] [CrossRef]
  31. Astorino, A.; Fuduli, A.; Gaudioso, M. DC models for spherical separation. J. Glob. Optim. 2010, 48, 657–669. [Google Scholar] [CrossRef]
  32. Astorino, A.; Fuduli, A.; Gaudioso, M. Margin maximization in spherical separation. Comput. Optim. Appl. 2012, 53, 301–322. [Google Scholar] [CrossRef]
  33. Astorino, A.; Gaudioso, M.; Seeger, A. Conic separation of finite sets. I. The homogeneous case. J. Convex Anal. 2014, 21, 1–28. [Google Scholar]
  34. Astorino, A.; Gaudioso, M.; Seeger, A. Conic separation of finite sets. II. The non-homogeneous case. J. Convex Anal. 2014, 21, 819–831. [Google Scholar]
  35. Le Thi, H.A.; Le, H.M.; Pham Dinh, T.; Van Huynh, N. Binary classification via spherical separator by DC programming and DCA. J. Glob. Optim. 2013, 56, 1393–1407. [Google Scholar] [CrossRef]
  36. Astorino, A.; Fuduli, A. Semisupervised spherical separation. Appl. Math. Model. 2015, 39, 6351–6358. [Google Scholar] [CrossRef]
  37. Wang, J.; Shen, X.; Pan, W. On efficient large margin semisupervised learning: Method and theory. J. Mach. Learn. Res. 2009, 10, 719–742. [Google Scholar]
  38. Bagirov, A.M.; Taheri, S.; Ugon, J. Nonsmooth DC programming approach to the minimum sum-of-squares clustering problems. Pattern Recognit. 2016, 53, 12–24. [Google Scholar] [CrossRef]
  39. Karmitsa, N.; Bagirov, A.M.; Taheri, S. New diagonal bundle method for clustering problems in large data sets. Eur. J. Oper. Res. 2017, 263, 367–379. [Google Scholar] [CrossRef]
  40. Khalaf, W.; Astorino, A.; D’Alessandro, P.; Gaudioso, M. A DC optimization-based clustering technique for edge detection. Optim. Lett. 2017, 11, 627–640. [Google Scholar] [CrossRef]
  41. Le Thi, H.; Pham Dinh, T. The DC (difference of convex functions) programming and DCA revisited with DC models of real world nonconvex optimization problems. J. Glob. Optim. 2005, 133, 23–46. [Google Scholar]
  42. Astorino, A.; Miglionico, G. Optimizing sensor cover energy via DC programming. Optim. Lett. 2016, 10, 355–368. [Google Scholar] [CrossRef]
  43. De Oliveira, W. Proximal bundle methods for nonsmooth DC programming. J. Glob. Optim. 2019, 75, 523–563. [Google Scholar] [CrossRef]
  44. De Oliveira, W.; Tcheou, M.P. An inertial algorithm for DC programming. Set-Valued Var. Anal. 2019, 27, 895–919. [Google Scholar] [CrossRef]
  45. Gaudioso, M.; Giallombardo, G.; Miglionico, G. Minimizing piecewise-concave functions over polytopes. Math. Oper. Res. 2018, 43, 580–597. [Google Scholar] [CrossRef]
  46. Joki, K.; Bagirov, A.M.; Karmitsa, N.; Mäkelä, M.M. A proximal bundle method for nonsmooth DC optimization utilizing nonconvex cutting planes. J. Glob. Optim. 2017, 68, 501–535. [Google Scholar] [CrossRef]
  47. Joki, K.; Bagirov, A.M.; Karmitsa, N.; Mäkelä, M.M.; Taheri, S. Double bundle method for finding Clarke stationary points in nonsmooth DC programming. Siam J. Optim. 2018, 28, 1892–1919. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Positive bag J i + .
Figure 1. Positive bag J i + .
Algorithms 12 00249 g001
Figure 2. Negative bag J i .
Figure 2. Negative bag J i .
Algorithms 12 00249 g002
Table 1. Medium-size test problems.
Table 1. Medium-size test problems.
Data SetDimensionInstancesBags
Elephant2301320200
Fox2301320200
Tiger2301220200
Musk-116647692
Musk-21666598102
Table 2. Large-size test problems.
Table 2. Large-size test problems.
Data SetDimensionInstancesBags
TST166683224400
TST268423344400
TST365683246400
TST466263391400
TST770373367400
TST969823300400
TST1070733453400
Table 3. Average test-correctness (%) for medium-size problems.
Table 3. Average test-correctness (%) for medium-size problems.
Data SetDC-MILMIL-RLDC-SMILmi-SVMMI-SVMMICAMIC Bundle
Elephant84.083.084.582.281.480.580.5
Fox57.054.556.058.257.858.758.3
Tiger84.575.081.078.484.082.679.1
Musk-174.580.076.7---75.6
Musk-274.073.079.0---76.8
Underlined means the best performance being.
Table 4. Average test-correctness (%) for large-size problems.
Table 4. Average test-correctness (%) for large-size problems.
Data SetDC-MILMIL-RLmi-SVMMI-SVMMICA
TST0194.395.593.693.994.5
TST0280.085.578.284.585.0
TST0386.586.887.082.286.0
TST0486.079.882.882.487.7
TST0779.883.581.378.078.9
TST0968.368.867.560.261.4
TST1078.077.579.679.582.3
Underlined means the best performance being.
Table 5. DC-MIL average efficiency. Medium-size test problems.
Table 5. DC-MIL average efficiency. Medium-size test problems.
Data SetTrainCpunFnG1nG2
Elephant91.03.14500243208
Fox79.93.055008180
Tiger95.52.83500237197
Musk-196.91.29500197177
Musk-293.56.52500174167
Table 6. DC-MIL average efficiency. Large-size test problems.
Table 6. DC-MIL average efficiency. Large-size test problems.
Data SetTrainCpunFnG1nG2
TST01100.070.222009391
TST0294.269.872008382
TST0399.664.772008281
TST0493.567.582008483
TST0799.274.112008584
TST0994.467.992008281
TST1091.972.242008180

Share and Cite

MDPI and ACS Style

Astorino, A.; Fuduli, A.; Giallombardo, G.; Miglionico, G. SVM-Based Multiple Instance Classification via DC Optimization. Algorithms 2019, 12, 249. https://0-doi-org.brum.beds.ac.uk/10.3390/a12120249

AMA Style

Astorino A, Fuduli A, Giallombardo G, Miglionico G. SVM-Based Multiple Instance Classification via DC Optimization. Algorithms. 2019; 12(12):249. https://0-doi-org.brum.beds.ac.uk/10.3390/a12120249

Chicago/Turabian Style

Astorino, Annabella, Antonio Fuduli, Giovanni Giallombardo, and Giovanna Miglionico. 2019. "SVM-Based Multiple Instance Classification via DC Optimization" Algorithms 12, no. 12: 249. https://0-doi-org.brum.beds.ac.uk/10.3390/a12120249

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop