26subset of the data and trains <strong>on</strong> that. The rest of the training data is tested <strong>on</strong> the resultingclassier, and a list of the errors is c<strong>on</strong>structed, sorted by how far <strong>on</strong> the wr<strong>on</strong>g side of themargin they lie (i.e. how egregiously the KKT c<strong>on</strong>diti<strong>on</strong>s are violated). The next chunk isc<strong>on</strong>structed from the rst N of these, combined with the N S support vectors already found,where N + N S is decided heuristically (a chunk size that is allowed to grow tooquickly ortoo slowly will result in slow overall c<strong>on</strong>vergence). Note that vectors can be dropped fromachunk, and that support vectors in <strong>on</strong>e chunk may not appear in the nal soluti<strong>on</strong>. Thisprocess is c<strong>on</strong>tinued until all data points are found to satisfy the KKT c<strong>on</strong>diti<strong>on</strong>s.The above method requires that the number of support vectors N S be small enough so thata Hessian of size N S by N S will t in memory. An alternative decompositi<strong>on</strong> algorithm hasbeen proposed which overcomes this limitati<strong>on</strong> (Osuna, Freund and Girosi, 1997b). Again,in this algorithm, <strong>on</strong>ly a small porti<strong>on</strong> of the training data is trained <strong>on</strong> at a given time, andfurthermore, <strong>on</strong>ly a subset of the support vectors need be in the \working set" (i.e. that setof points whose 's are allowed to vary). This method has been shown to be able to easilyhandle a problem with 110,000 training points and 100,000 support vectors. However, itmust be noted that the speed of this approach relies <strong>on</strong> many of the support vectors havingcorresp<strong>on</strong>ding Lagrange multipliers i at the upper bound, i = C.These training algorithms maytakeadvantage of parallel processing in several ways. First,all elements of the Hessian itself can be computed simultaneously. Sec<strong>on</strong>d, each elementoften requires the computati<strong>on</strong> of dot products of training data, which could also be parallelized.Third, the computati<strong>on</strong> of the objective functi<strong>on</strong>, or gradient, which isaspeedbottleneck, can be parallelized (it requires a matrix multiplicati<strong>on</strong>). Finally, <strong>on</strong>e can envisi<strong>on</strong>parallelizing at a higher level, <strong>for</strong> example by training <strong>on</strong> dierent chunks simultaneously.Schemes such as these, combined with the decompositi<strong>on</strong> algorithm of (Osuna, Freund andGirosi, 1997b), will be needed to makevery large problems (i.e. >> 100,000 support vectors,with many not at bound), tractable.5.1.2. Testing In test phase, <strong>on</strong>e must simply evaluate Eq. (61), which will requireO(MN S ) operati<strong>on</strong>s, where M is the number of operati<strong>on</strong>s required to evaluate the kernel.For dot product and RBF kernels, M is O(d L ), the dimensi<strong>on</strong> of the data vectors. Again,both the evaluati<strong>on</strong> of the kernel and of the sum are highly parallelizable procedures.In the absence of parallel hardware, <strong>on</strong>e can still speed up test phase by a large factor, asdescribed in Secti<strong>on</strong> 9.6. The VC Dimensi<strong>on</strong> of <strong>Support</strong> <strong>Vector</strong> <strong>Machines</strong>We now show that the VC dimensi<strong>on</strong> of SVMs can be very large (even innite). We willthen explore several arguments as to why, in spite of this, SVMs usually exhibit goodgeneralizati<strong>on</strong> per<strong>for</strong>mance. However it should be emphasized that these are essentiallyplausibility arguments. Currently there exists no theory which guarantees that a givenfamily of SVMs will have high accuracy <strong>on</strong> a given problem.We will call any kernel that satises Mercer's c<strong>on</strong>diti<strong>on</strong> a positive kernel, and the corresp<strong>on</strong>dingspace H the embedding space. We will also call anyembedding space with minimaldimensi<strong>on</strong> <strong>for</strong> a given kernel a \minimal embedding space". We have the followingTheorem 3 Let K be a positive kernel which corresp<strong>on</strong>ds to a minimal embedding spaceH. Then the VC dimensi<strong>on</strong> of the corresp<strong>on</strong>ding support vector machine (where the errorpenalty C in Eq. (44) is allowed to take all values) is dim(H)+1.
27Proof: If the minimal embedding space has dimensi<strong>on</strong> d H ,thend H points in the image ofL under the mapping can be found whose positi<strong>on</strong> vectors in H are linearly independent.From Theorem 1, these vectors can be shattered by hyperplanes in H. Thus by eitherrestricting ourselves to SVMs <strong>for</strong> the separable case (Secti<strong>on</strong> 3.1), or <strong>for</strong> which the errorpenalty C is allowed to take all values (so that, if the points are linearly separable, a C canbe found such that the soluti<strong>on</strong> does indeed separate them), the family of support vectormachines with kernel K can also shatter these points, and hence has VC dimensi<strong>on</strong> d H +1.Let's look at two examples.6.1. The VC Dimensi<strong>on</strong> <strong>for</strong> Polynomial KernelsC<strong>on</strong>sider an SVM with homogeneous polynomial kernel, acting <strong>on</strong> data in R dL :K(x 1 x 2 )=(x 1 x 2 ) p x 1 x 2 2 R dL (78)As in the case when d L = 2 and the kernel is quadratic (Secti<strong>on</strong> 4), <strong>on</strong>e can explicitlyc<strong>on</strong>struct the map . Letting z i = x 1i x 2i ,sothatK(x 1 x 2 )=(z 1 + + z dL ) p ,we see thateach dimensi<strong>on</strong> of H corresp<strong>on</strong>ds to a term with given powers of the z i in the expansi<strong>on</strong> ofK. In fact if we choose to label the comp<strong>on</strong>ents of (x) inthismanner,we can explicitlywrite the mapping <strong>for</strong> any p and d L : r1r 2r dL(x) =This leads tos p!r 1 !r 2 ! r dL !x r11 xr22 xr d Ld Ld L Xi=1r i = p r i 0 (79)Theorem 4 If the space in which the data live has dimensi<strong>on</strong> d L (i.e. L = R dL ), thedimensi<strong>on</strong> of the minimal embedding space, <strong>for</strong> ; homogeneous polynomial kernels of degree p(K(x 1 x 2 )=(x 1 x 2 ) p x 1 x 2 2 R dL d), is L+p;1 p .; dL+p;1p(The proof is in the Appendix). Thus the VC dimensi<strong>on</strong> of SVMs with these kernels is+1. As noted above, this gets very large very quickly.6.2. The VC Dimensi<strong>on</strong> <strong>for</strong> Radial Basis Functi<strong>on</strong> KernelsTheorem 5 C<strong>on</strong>sider the class of Mercer kernels <strong>for</strong> which K(x 1 x 2 ) ! 0 as kx 1 ; x 2 k!1, and <strong>for</strong> which K(x x) is O(1), and assume that the data can be chosen arbitrarily fromR d . Then the family of classiers c<strong>on</strong>sisting of support vector machines using these kernels,and <strong>for</strong> which the error penalty is allowed to take all values, has innite VC dimensi<strong>on</strong>.Proof: The kernel matrix, K ij K(x i x j ), is a Gram matrix (a matrix of dot products:see (Horn, 1985)) in H. Clearly we can choose training data such that all o-diag<strong>on</strong>alelements K i6=j can be made arbitrarily small, and by assumpti<strong>on</strong> all diag<strong>on</strong>al elements K i=jare of O(1). The matrix K is then of full rank hence the set of vectors, whose dot productsin H <strong>for</strong>m K, are linearly independent (Horn, 1985) hence, by Theorem 1, the points can beshattered by hyperplanes in H, and hence also by support vector machines with sucientlylarge error penalty. Since this is true <strong>for</strong> any nite number of points, the VC dimensi<strong>on</strong> ofthese classiers is innite.