Máy vectơ hỗ trợ

Máy vectơ hỗ trợ (SVM - viết tắt tên tiếng Anh support vector machine) là một khái niệm trong thống kê và khoa học máy tính cho một tập hợp các phương pháp học có giám sát liên quan đến nhau để phân loại và phân tích hồi quy. SVM dạng chuẩn nhận dữ liệu vào và phân loại chúng vào hai lớp khác nhau. Do đó SVM là một thuật toán phân loại nhị phân. Với một bộ các ví dụ luyện tập thuộc hai thể loại cho trước, thuật toán luyện tập SVM xây dựng một mô hình SVM để phân loại các ví dụ khác vào hai thể loại đó. Một mô hình SVM là một cách biểu diễn các điểm trong không gian và lựa chọn ranh giới giữa hai thể loại sao cho khoảng cách từ các ví dụ luyện tập tới ranh giới là xa nhất có thể. Các ví dụ mới cũng được biểu diễn trong cùng một không gian và được thuật toán dự đoán thuộc một trong hai thể loại tùy vào ví dụ đó nằm ở phía nào của ranh giới.

== Tổng quan về máy vectơ hỗ trợ ==
Một máy vectơ hỗ trợ xây dựng một siêu phẳng hoặc một tập hợp các siêu phẳng trong một không gian nhiều chiều hoặc vô hạn chiều, có thể được sử dụng cho phân loại, hồi quy, hoặc các nhiệm vụ khác. Một cách trực giác, để phân loại tốt nhất thì các siêu phẳng nằm ở càng xa các điểm dữ liệu của tất cả các lớp (gọi là hàm Biên) càng tốt, vì nói chung Biên càng lớn thì sai số tổng quát hóa của thuật toán phân loại càng bé.

Trong nhiều trường hợp, không thể phân chia các lớp dữ liệu một cách tuyến tính trong một không gian ban đầu được dùng để mô tả một vấn đề. Vì vậy, nhiều khi cần phải ánh xạ các điểm dữ liệu trong không gian ban đầu vào một không gian mới nhiều chiều hơn, để việc phân tách chúng trở nên dễ dàng hơn trong không gian mới. Để việc tính toán được hiệu quả, ánh xạ sử dụng trong thuật toán SVM chỉ đòi hỏi tích vô hướng của các vectơ dữ liệu trong không gian mới có thể được tính dễ dàng từ các tọa độ trong không gian cũ. Tích vô hướng này được xác định bằng một hàm hạt nhân "K"("x","y") phù hợp. Một siêu phẳng trong không gian mới được định nghĩa là tập hợp các điểm có tích vô hướng với một vectơ cố định trong không gian đó là một hằng số. Vectơ xác định một siêu phẳng sử dụng trong SVM là một tổ hợp tuyến tính của các vectơ dữ liệu luyện tập trong không gian mới với các hệ số "α". Với siêu phẳng lựa chọn như trên, các điểm "x" trong không gian đặc trưng được ánh xạ vào một siêu mặt phẳng là các điểm thỏa mãn:

Ghi chú rằng nếu "K"("x","y") nhận giá trị ngày càng nhỏ khi "y" xa dần khỏi "x" thì mỗi số hạng của tổng trên được dùng để đo độ tương tự giữa "x" với điểm "x" tương ứng trong dữ liệu luyện tập. Như vậy, tác dụng của tổng trên chính là so sánh khoảng cách giữa điểm cần dự đoán với các điểm dữ liệu đã biết. Lưu ý là tập hợp các điểm "x" được ánh xạ vào một siêu phẳng có thể có độ phức tạp tùy ý trong không gian ban đầu, nên có thể phân tách các tập hợp thậm chí không lồi trong không gian ban đầu.

== Lịch sử ==
Thuật toán SVM ban đầu được tìm ra bởi Vladimir N. Vapnik và dạng chuẩn hiện nay sử dụng Biên mềm được tìm ra bởi Vapnik và Corinna Cortes năm 1995.

== Đặt vấn đề ==
Phân loại thống kê là một nhiệm vụ phổ biến trong học máy. Trong mô hình học có giám sát, thuật toán được cho trước một số điểm dữ liệu cùng với nhãn của chúng thuộc một trong hai lớp cho trước. Mục tiêu của thuật toán là xác định xem một điểm dữ liệu "mới" sẽ được thuộc về lớp nào. Mỗi điểm dữ liệu được biểu diễn dưới dạng một vector p-chiều, và ta muốn biết liệu có thể chia tách hai lớp dữ liệu bằng một siêu phẳng "p" − 1 chiều. Đây gọi là phân loại tuyến tính. Có nhiều siêu phẳng có thể phân loại được dữ liệu. Một lựa chọn hợp lý trong chúng là siêu phẳng có Biên lớn nhất giữa hai lớp.

== SVM tuyến tính ==
Ta có một tập huấn luyện formula_1 gồm "n" điểm có dạng

với "y" mang giá trị 1 hoặc −1, xác định lớp của điểm formula_3. Mỗi formula_4 là một vectơ thực "p"-chiều. Ta cần tìm siêu phẳng có Biên lớn nhất chia tách các điểm có formula_5 và các điểm có formula_6. Mỗi siêu phẳng đều có thể được viết dưới dạng một tập hợp các điểm formula_7 thỏa mãn

với formula_9 ký hiệu cho tích vô hướng và formula_10 là một vectơ pháp tuyến của siêu phẳng. Tham số formula_11 xác định khoảng cách giữa gốc tọa độ và siêu phẳng theo hướng vectơ pháp tuyến formula_10.

Chúng ta cần chọn formula_10 và formula_14 để cực đại hóa Biên, hay khoảng cách giữa hai siêu mặt song song ở xa nhau nhất có thể trong khi vẫn phân chia được dữ liệu. Các siêu mặt ấy được xác định bằng

và

Để ý rằng nếu dữ liệu huấn luyện có thể được chia tách một cách tuyến tính, thì ta có thể chọn hai siêu phẳng của Biên sao cho không có điểm nào ở giữa chúng và sau đó tăng khoảng cách giữa chúng đến tối đa có thể. Bằng hình học, ta tìm được khoảng cách giữa hai siêu phẳng là formula_17. Vì vậy ta muốn cực tiểu hóa giá trị formula_18. Để đảm bảo không có điểm dữ liệu nào trong Biên, ta thêm vào các điều kiện sau:

Với mỗi formula_19 ta có

hoặc

Có thể viết gọn lại như sau với mọi formula_22:

Tóm lại, ta có bài toán tối ưu hóa sau:

Cực tiểu hóa (theo formula_24)

với điều kiện (với mọi formula_26)

=== Dạng ban đầu ===
Bài toán tối ưu ở mục trên tương đối khó giải vì hàm mục tiêu phụ thuộc vào ||w||, là một hàm có khai căn. Tuy nhiên có thể thay ||w|| bằng hàm mục tiêu formula_28 (hệ số 1/2 để tiện cho các biến đổi toán học sau này) mà không làm thay đổi lời giải (lời giải của bài toán mới và bài toán ban đầu có cùng w và b). Đây là một bài toán quy hoạch toàn phương. Cụ thể hơn:

Cực tiểu hóa (theo formula_24)

với điều kiện (với mọi formula_26)

Bằng cách thêm các nhân tử Lagrange formula_33, bài toán trên trở thành

nghĩa là ta cần tìm một điểm yên ngựa. Khi đó, tất cả các điểm không nằm trên Biên, nghĩa là formula_35 đều không ảnh hưởng đến giá trị hàm mục tiêu vì ta có thể chọn formula_36 bằng không.

Có thể giải bài toán này bằng các kĩ thuật thông thường cho quy hoạch toàn phương. Theo điều kiện Karush–Kuhn–Tucker, lời giải có thể được viết dưới dạng tổ hợp tuyến tính của các vectơ luyện tập

Chỉ có một vài formula_36 nhận giá trị lớn hơn 0. Các điểm formula_39 tương ứng là các "vectơ hỗ trợ" nằm trên Biên và thỏa mãn formula_40. Từ điều kiện này, ta nhận thấy
từ đó ta suy ra được giá trị formula_14. Trên thực tế, một cách thức tốt hơn để tính formula_14 là tính giá trị trung bình từ tất cả formula_44 vectơ hỗ trợ:

=== Dạng đối ngẫu ===
Nếu viết điều kiện phân loại dưới dạng đối ngẫu không điều kiện thì sẽ dễ dàng nhận thấy siêu phẳng với Biên lớn nhất, và do đó nhiệm vụ phân loại, chỉ phụ thuộc vào các điểm luyện tập nằm trên Biên, còn gọi là các "vectơ hỗ trợ".

Vì formula_46 và formula_47, ta nhận thấy bài toán đối ngẫu của SVM là chính là bài toán tối ưu hóa sau:

Cực đại hóa (theo formula_36)

với điều kiện (với mọi formula_26)

và điều kiện sau ứng với việc cực tiểu hóa theo formula_52

Ở đây hàm hạt nhân được định nghĩa là formula_54.

Sau khi giải xong, có thể tính formula_55 từ các giá trị formula_56 tìm được như sau:

== Biên mềm ==
Năm 1995, Corinna Cortes và Vladimir N. Vapnik đề xuất một ý tưởng mới cho phép thuật toán gán nhãn sai cho một số ví dụ luyện tập. Nếu không tồn tại siêu phẳng nào phân tách được hai lớp dữ liệu, thì thuật toán "Biên mềm" sẽ chọn một siêu phẳng phân tách các ví dụ luyện tập tốt nhất có thể, và đồng thời cực đại hóa khoảng cách giữa siêu phẳng với các ví dụ được gán đúng nhãn. Phương pháp này sử dụng các biến bù formula_58, dùng để đo độ sai lệch của ví dụ formula_59

Hàm mục tiêu có thêm một số hạng mới để phạt thuật toán khi formula_58 khác không, và bài toán tối ưu hóa trở thành việc trao đổi giữa Biên lớn và mức phạt nhỏ. Nếu hàm phạt là tuyến tính thì bài toán trở thành:

với điều kiện (với mọi formula_63)

Có thể giải bài toán trên bằng nhân tử Lagrange tương tự như trường hợp cơ bản ở trên. Bài toán cần giải trở thành:

với formula_66.

=== Dạng đối ngẫu ===
Cực đại hóa (theo formula_36)

với điều kiện (với mọi formula_26)
và

Ưu điểm của việc dùng hàm phạt tuyến tính là các biến bù biến mất khỏi bài toán đối ngẫu, và hằng số "C" chỉ xuất hiện dưới dạng một chặn trên cho các nhân tử Lagrange. Cách đặt vấn đề trên đã mang lại nhiều thành quả trong thực tiễn, và Cortes và Vapnik đã nhận được giải Paris Kanellakis của ACM năm 2008 cho đóng góp này. Các hàm phạt phi tuyến cũng được sử dụng, đặc biệt là để giảm ảnh hưởng của các trường hợp ngoại lệ, tuy nhiên nếu không lựa chọn hàm phạt cẩn thận thì bài toán trở thành không lồi, và việc tìm lời giải tối ưu toàn cục thường là rất khó.

== Xem thêm ==
- In situ adaptive tabulation
- Máy hạt nhân
- Predictive analytics
- Relevance vector machine, Một mô hình máy hạt nhân thưa xác suất có dạnghàm số giống như SVM.
- Tối ưu hóa nhỏ nhất tuần tự

== Liên kết ngoài ==
- * www.kernel-machines.org "(thông tin tổng quan và danh sách các bài báo nghiên cứu)"
- www.support-vector-machines.org "(Bài báo nghiên cứu, đánh giá, phần mềm, liên kết có liên quan đến máy vectơ hỗ trợ)"
- videolectures.net "(Video bài giảng về SVM)"
- Phim ngắn: Minh họa SVM sử dụng hàm hạt nhân đa thức.
- Một hướng dẫn sử dụng SVM cho người mới học bởi Tristan Fletcher .
- www.shogun-toolbox.org "(Shogun (hộp công cụ) gồm khoảng 20 thư viện lập trình SVM)"
- libsvm libsvm là một thư viện lập trình SVM
- liblinear liblinear là một thư viện lập trình gồm nhiều thuật toán phân loại tuyến tính, trong đó có SVM
- flssvm flssvm là một thư viện lập trình svm bình phương nhỏ nhất viết bằng fortran
- Shark Shark là một thư viện học máy viết bằng C++ có chứa nhiều loại SVM
- dlib dlib là một thư viện C++ cho máy hạt nhân và SVM
- SVM light là một bộ phần mềm cho học máy và phân loại bằng SVM.

==Tham khảo==
- Sergios Theodoridis and Konstantinos Koutroumbas "Pattern Recognition", 4th Edition, Academic Press, 2009, ISBN 978-1-59749-272-0
- Nello Cristianini and John Shawe-Taylor. "An Introduction to Support Vector Machines and other kernel-based learning methods". Cambridge University Press, 2000. ISBN 0-521-78019-5 "( SVM Book)"
- Huang T.-M., Kecman V., Kopriva I. (2006), Kernel Based Algorithms for Mining Huge Data Sets, Supervised, Semi-supervised, and Unsupervised Learning, Springer-Verlag, Berlin, Heidelberg, 260 pp. 96 illus., Hardcover, ISBN 3-540-31681-7
- Vojislav Kecman: "Learning and Soft Computing — Support Vector Machines, Neural Networks, Fuzzy Logic Systems", The MIT Press, Cambridge, MA, 2001.
- Bernhard Schölkopf and A. J. Smola: "Learning with Kernels". MIT Press, Cambridge, MA, 2002. "(Partly available on line: .)" ISBN 0-262-19475-9
- Bernhard Schölkopf, Christopher J.C. Burges, and Alexander J. Smola (editors). "Advances in Kernel Methods: Support Vector Learning". MIT Press, Cambridge, MA, 1999. ISBN 0-262-19416-3.
- John Shawe-Taylor and Nello Cristianini. "Kernel Methods for Pattern Analysis". Cambridge University Press, 2004. ISBN 0-521-81397-2 "( Kernel Methods Book)"
- Ingo Steinwart and Andreas Christmann. "Support Vector Machines". Springer-Verlag, New York, 2008. ISBN 978-0-387-77241-7 "( SVM Book)"
- P.J. Tan and D.L. Dowe (2004), MML Inference of Oblique Decision Trees, Lecture Notes in Artificial Intelligence (LNAI) 3339, Springer-Verlag, pp1082-1088. (This paper uses minimum message length (MML) and actually incorporates probabilistic support vector machines in the leaves of decision trees.)
- Vladimir Vapnik. "The Nature of Statistical Learning Theory". Springer-Verlag, 1995. ISBN 0-387-98780-0
- Vladimir Vapnik, S.Kotz "Estimation of Dependences Based on Empirical Data" Springer, 2006. ISBN 0-387-30865-2, 510 pages [this is a reprint of Vapnik's early book describing philosophy behind SVM approach. The 2006 Appendix describes recent development].
- Dmitriy Fradkin and Ilya Muchnik "Support Vector Machines for Classification" in J. Abello and G. Carmode (Eds) "Discrete Methods in Epidemiology", DIMACS Series in Discrete Mathematics and Theoretical Computer Science, volume 70, pp. 13–20, 2006. . Succinctly describes theoretical ideas behind SVM.
- Kristin P. Bennett and Colin Campbell, "Support Vector Machines: Hype or Hallelujah?", SIGKDD Explorations, 2,2, 2000, 1–13. . Excellent introduction to SVMs with helpful figures.
- Ovidiu Ivanciuc, "Applications of Support Vector Machines in Chemistry", In: "Reviews in Computational Chemistry", Volume 23, 2007, pp. 291–400. Reprint available:
- Catanzaro, Sundaram, Keutzer, "Fast Support Vector Machine Training and Classification on Graphics Processors", In: "International Conference on Machine Learning", 2008