Po co zawracać sobie głowę podwójnym problemem przy montażu SVM?


50

Biorąc pod uwagę punkty danych i etykiety , podstawowym problemem z twardym marginesem SVM jestx1,,xnRdy1,,yn{1,1}

minimizew,w012wTw
s.t.i:yi(wTxi+w0)1

który jest programem kwadratowym ze zmiennymi , które należy zoptymalizować dla ograniczeń i . Podwójnyd+1i

maximizeαi=1nαi12i=1nj=1nyiyjαiαjxiTxj
s.t.i:αi0i=1nyiαi=0
to program kwadratowy z zmiennymi do optymalizacji i nierównościami i ograniczeniami równości.n+1nn

Kiedy wdrażam SVM z twardym marginesem, dlaczego miałbym rozwiązać podwójny problem zamiast pierwotnego? Pierwotny problem wydaje mi się bardziej „intuicyjny” i nie muszę się martwić luką w dualności, stanem Kuhna-Tuckera itp.

Rozsądne byłoby rozwiązanie podwójnego problemu, jeśli , ale podejrzewam, że istnieją lepsze powody. Czy tak jest w przypadku?dn


26
Krótka odpowiedź to jądra. Długa odpowiedź to keeerneeels (-;

Najważniejszą kwestią podwójnego problemu jest wprowadzenie sztuczki jądra, która ma na celu zmapowanie oryginalnych danych w przestrzeń o większym wymiarze.
BigeyeDestroyer

Odpowiedzi:


40

Na podstawie notatek z wykładu, o których mowa w odpowiedzi @ user765195 (dzięki!), Najbardziej widocznymi przyczynami wydają się:

Rozwiązując pierwotny problem, uzyskujemy optymalne , ale nic nie wiemy o . Aby sklasyfikować punkt zapytania , musimy jawnie obliczyć iloczyn skalarny , co może być kosztowne, jeśli jest duże.wαixwTxd

Rozwiązując podwójny problem, otrzymujemy (gdzie dla wszystkich oprócz kilku punktów - wektorów wsparcia). Aby sklasyfikować punkt zapytania , obliczamyα i = 0 xαiαi=0x

wTx+w0=(i=1nαiyixi)Tx+w0=i=1nαiyixi,x+w0

Ten termin jest bardzo skutecznie obliczany, jeśli istnieje tylko kilka wektorów pomocniczych. Ponadto, ponieważ teraz mamy produkt skalarny obejmujący wyłącznie wektory danych , możemy zastosować sztuczkę jądra .


6
Poczekaj poczekaj. Załóżmy, że masz dwa wektory pomocnicze x1 i x2. Nie możesz mieć mniej niż dwóch, prawda? Czy mówisz, że obliczenia <x1, x> i <x2, x> są szybsze niż <w, x>?
Leo

1
@Leo: Pamiętaj, że używam <x1, x>i wTx. Pierwszy z nich jest używany jako symbol oceny K (x1, x) jądra, która rzutuje x1 i x na przestrzeń o bardzo dużych wymiarach i domyślnie oblicza iloczyn skalarny rzutowanych wartości. To ostatnie jest to normalny iloczyn skalarny, to wi xmuszą być rzutowany bezpośrednio, a następnie oblicza się iloczyn skalarny wyraźnie. W zależności od wyboru jądra, jedno jawne obliczenie może wymagać znacznie więcej obliczeń niż wiele ocen jądra.
blubb

1
Jak rozumiem pierwotny problem, to mnożniki Lagrange'a, więc dlaczego nie możemy rozwiązać pierwotnego problemu, aby znaleźć ? To znaczy, że prawdopodobnie nie trzeba uciekać się do dualnego, aby dowiedzieć się „s, prawda? α αααα
awokado

2
„Ponadto, ponieważ obecnie mamy produkt skalarny obejmujący wyłącznie wektory danych, możemy zastosować sztuczkę jądra”. - Dotyczy to również pierwotnego sformułowania.
Firebug

2
Jeśli ludzie chcą więcej szczegółów na temat komentarza z @Firebug ... sprawdź równania 10-12 lib.kobe-u.ac.jp/repository/90001050.pdf (która jest nieograniczoną wersją pierwotną).
MrDrFenner


3

Oto jeden z powodów, dla których podwójne sformułowanie jest atrakcyjne z punktu widzenia optymalizacji numerycznej. Szczegóły można znaleźć w następującym artykule :

Hsieh, C.-J., Chang, K.-W., Lin, C.-J., Keerthi, SS i Sundararajan, S., „Metoda zejścia z podwójną współrzędną dla SVM o dużej skali” 25. międzynarodowa konferencja nt. Uczenia maszynowego, Helsinki, 2008.

Podwójne sformułowanie obejmuje jedno ograniczenie równości afinicznej i n ograniczenia.

1. Ograniczenie równości afinicznej można „wyeliminować” z podwójnego sformułowania.

Można to zrobić po prostu patrząc na swoje dane w R ^ (d + 1) poprzez osadzenie R ^ d w R ^ (d + 1), rezygnując z dodania pojedynczej współrzędnej „1” do każdego punktu danych, tj. R ^ d ----> R ^ (d + 1): (a1, ..., reklama) | ---> (a1, ..., reklama, 1).

Wykonanie tego dla wszystkich punktów w zestawie treningowym przekształca problem liniowej separowalności w R ^ (d + 1) i eliminuje stały składnik w0 z twojego klasyfikatora, co z kolei eliminuje ograniczenie równości afinicznej z podwójnego.

2. W punkcie 1, dual można łatwo rzutować jako wypukły kwadratowy problem optymalizacji, którego ograniczenia są tylko ograniczeniami związanymi.

3. Podwójny problem można teraz skutecznie rozwiązać, tj. Za pomocą algorytmu podwójnego współrzędnego opadania, który daje optymalne dla epsilonu rozwiązanie w O (log (1 / epsilon)).

Odbywa się to poprzez zauważenie, że naprawienie wszystkich alf z wyjątkiem jednego daje rozwiązanie w formie zamkniętej. Następnie możesz kolejno przełączać wszystkie alfy (np. Wybierając jedną losowo, naprawiając wszystkie pozostałe alfy, obliczając rozwiązanie formy zamkniętej). Można pokazać, że w ten sposób uzyskasz prawie optymalne rozwiązanie „dość szybko” (patrz Twierdzenie 1 we wspomnianym artykule).

Istnieje wiele innych powodów, dla których podwójny problem jest atrakcyjny z punktu widzenia optymalizacji, niektóre z nich wykorzystują fakt, że ma tylko jedno ograniczenie równości afinicznej (wszystkie pozostałe ograniczenia są ograniczeniami związanymi), podczas gdy inne wykorzystują obserwację, że w rozwiązaniu podwójnego problemu „często większość alf” ma wartość zero (niezerowe alfy odpowiadające wektorom wspierającym).

Dobry przegląd zagadnień optymalizacji numerycznej dla maszyn SVM można znaleźć w prezentacji Stephena Wrighta w Computational Learning Workshop (2009).

PS: Jestem tu nowy. Przepraszamy za niedostateczne wykorzystanie notacji matematycznej na tej stronie.


1
Informacje na temat korzystania z składu matematycznego znajdują się tutaj: math.meta.stackexchange.com/questions/5020/…
Przywróć Monikę

-5

Moim zdaniem w notatkach z wykładu Andrew ng wyraźnie zaznaczono, że pierwotny problem 1 / || w || jest problemem niewypukłym. Podwójny jest problemem wypukłym i zawsze łatwo jest znaleźć optymalną funkcję wypukłą.


1
Pierwotna SVM, jak podano powyżej, jest wypukła.
Dougal
Korzystając z naszej strony potwierdzasz, że przeczytałeś(-aś) i rozumiesz nasze zasady używania plików cookie i zasady ochrony prywatności.
Licensed under cc by-sa 3.0 with attribution required.