Tukivektorikone

Tukivektorikoneen yksinkertainen kaavio.

Tukivektorikone (engl. support-vector machine) on 1990-luvulla kehitetty lineaarinen luokitinmalli, joka soveltuu luokitteluun ja käyränsovitustehtävään. Tukivektorikone voidaan toteuttaa neuroverkolla.[1] Tukivektorikoneen yleistämiskyky on MLP-neuroverkkoon verrattuna parempi. Yleistämiskyky kuvaa luokittimen kykyä luokitella ennalta tuntemattomat näytteet oikein.

Perusteet

Päätöstaso (keskellä) ja marginaalitasot.

Tukivektorikoneen perusajatus on sovittaa kahden näytejoukon väliin sellainen taso, että sen kanssa yhdensuuntaisten marginaalitasojen välimatka on mahdollisimman suuri eikä yksikään näyte jää marginaalitasojen väliin. Marginaalitasojen välimatkaa rajoittavia näytevektoreita kutsutaan tukivektoreiksi. Luokittelun tulos riippuu ainoastaan näistä tukivektoreista.

Tukivektorikoneen tulee opetuksen jälkeen osata luokitella mielivaltainen näyte. Luokitin opetetaan opetusaineiston perusteella, ja siksi sen hyvyys osaltaan riippuu opetusaineiston hyvyydestä.

Kun opetusjoukot eivät ole separoituvia, käytetään tukivektorikonetoteutusta, jota kutsutaan joustavan marginaalin luokittimeksi. Menetelmä etsii opetusaineistosta ne näytevektorit, jotka määrittävät eri luokkien rajat.

Optimimarginaaliluokitin

Optimimarginaaliluokitin on tukivektorikoneen toteutus separoituville näytejoukoille. Optimimarginaaliluokitin määrää päätöstason w T x + b = 0 {\displaystyle \mathbf {w} ^{T}\mathbf {x} +b=0} kahden separoituvan näytejoukon välille. Lisäksi vaaditaan että kummankin luokan päätöstasoa lähimmän pisteen ja päätöstason välinen etäisyys on mahdollisimman suuri. Tämä tarkoittaa sitä että kummankin luokan lyhin etäisyys pintaan nähden on sama, koska kokonaisuuden kannalta lyhin etäisyys määrittyy lähempänä olevan luokan perusteella.

Luokat ovat lineaarisesti separoituvia silloin kun kahden eri luokan yksikäsitteinen erottaminen yhdellä tasolla on mahdollista.

Etsitään piirreavaruuden R N {\displaystyle R^{N}} tasoa ( w , b ) {\displaystyle (\mathbf {w} ,b)} siten, että näytevektorille x i {\displaystyle \mathbf {x} _{i}} pätee

{ w T x i + b > 0 , kun  x i  kuuluu luokkaan a w T x i + b < 0 , kun  x i  kuuluu luokkaan b {\displaystyle {\begin{cases}\mathbf {w} ^{T}\mathbf {x} _{i}+b>0\quad ,{\mbox{kun }}\mathbf {x} _{i}{\mbox{ kuuluu luokkaan a}}\\\mathbf {w} ^{T}\mathbf {x} _{i}+b<0\quad ,{\mbox{kun }}\mathbf {x} _{i}{\mbox{ kuuluu luokkaan b}}\end{cases}}} .

Sovitaan, että luokan a {\displaystyle a} pisteitä vastaa y {\displaystyle y} :n arvo + 1 {\displaystyle +1} ja luokan b {\displaystyle b} 1 {\displaystyle -1} . Muodostetaan yhtälö

y ( w T x + b ) > 0 {\displaystyle y(\mathbf {w} ^{T}\mathbf {x} +b)>0} ,

joka pätee kummankin luokan pisteille. Jos opetustieto koostuu l {\displaystyle l} eri pisteestä, jotka ovat x i {\displaystyle \mathbf {x} _{i}} , kun i = 1 , , l {\displaystyle i=1,\dots ,l} , ja y i {\displaystyle y_{i}} ovat niiden vastaavat luokat (joko + 1 {\displaystyle +1} tai 1 {\displaystyle -1} ), niin vaadimme että luokat erottava taso ( w , b ) {\displaystyle (\mathbf {w} ,b)} toteuttaa ehdon

y i ( w T x i + b ) > 0 i = 1 , , l {\displaystyle y_{i}(\mathbf {w} ^{T}\mathbf {x} _{i}+b)>0\quad \forall i=1,\dots ,l} .

Näytepisteen euklidinen etäisyys tasoon ( w , b ) {\displaystyle (\mathbf {w} ,b)} on

min i = 1 l | w T x i + b | w {\displaystyle \min _{i=1}^{l}{\frac {|\mathbf {w} ^{T}\mathbf {x} _{i}+b|}{\|\mathbf {w} \|}}} .

Päätöstason löytämiseksi on olemassa useita ratkaisutapoja. Kirjallisuudessa usein esiintyvä ratkaisutapa perustuu neliölliseen optimointiin. Toinen tapa on etsiä kummallekin luokalle konveksipeitteet, ja sen jälkeen etsiä konveksipeitteiden välinen lyhin etäisyys ja tätä vastaavat pisteet p a {\displaystyle p_{a}} ja p b {\displaystyle p_{b}} peitteiden pinnalla. Päätöstason normaalivektori w {\displaystyle \mathbf {w} } on pisteiden p a {\displaystyle p_{a}} ja p b {\displaystyle p_{b}} erotus vektorin suuntainen, ja b {\displaystyle b} on näiden pisteiden keskiarvo.

Ratkaiseminen neliöllisellä optimoinnilla

Optimimarginaaliluokitin määritetään ratkaisemalla seuraava neliöllinen optimointiongelma. Ratkaistaan päätöstaso minimoimalla lauseke 1 2 w 2 {\displaystyle {\frac {1}{2}}\|\mathbf {w} \|^{2}} tason normaalivektorin w {\displaystyle \mathbf {w} } suhteen ehdolla

y i ( w T x i + b ) 1 i = 1 , , l {\displaystyle y_{i}(\mathbf {w} ^{T}\mathbf {x} _{i}+b)\geq 1\quad \forall i=1,\dots ,l} .

Muuttuja b {\displaystyle b} ratkaistaan ehdosta.

Päätöstason yksiselitteiseen määrittämiseen tarvitsee tuntea kaksi tasoa lähinnä olevaa pistettä x n {\displaystyle \mathbf {x} _{n}} ja x m {\displaystyle \mathbf {x} _{m}} . Pisteen i {\displaystyle i} euklidinen etäisyys tasoon w {\displaystyle \mathbf {w} } on

δ i = | w T x i + b | w {\displaystyle \delta _{i}={\frac {|\mathbf {w} ^{T}\mathbf {x} _{i}+b|}{\|\mathbf {w} \|}}} ,

missä w {\displaystyle \|\mathbf {w} \|} on normalisointitermi.

Taso ( w , b ) {\displaystyle (\mathbf {w} ,b)} on yhdensuuntainen tason ( c w , c b ) {\displaystyle (c\mathbf {w} ,cb)} kanssa, kun c > 0 {\displaystyle c>0} , josta seuraa että c {\displaystyle c} voidaan valita siten, että | c w T x i + c b | = 1 {\displaystyle |c\mathbf {w} ^{T}\mathbf {x} _{i}+cb|=1} pisteille x n {\displaystyle \mathbf {x} _{n}} ja x m {\displaystyle \mathbf {x} _{m}} . Tällöin euklidisten etäisyyksien summaksi saadaan

δ n + δ m = | c w T x n + c b | c w + | c w T x m + c b | c w = 2 c w {\displaystyle \delta _{n}+\delta _{m}={\frac {|c\mathbf {w} ^{T}\mathbf {x} _{n}+cb|}{\|c\mathbf {w} \|}}+{\frac {|c\mathbf {w} ^{T}\mathbf {x} _{m}+cb|}{\|c\mathbf {w} \|}}={\frac {2}{\|c\mathbf {w} \|}}} .

Tämän seurauksena rajoittamalla reunaehdot annetun tehtävän mukaisesti ja minimoimalla w {\displaystyle \|\mathbf {w} \|} , marginaalitasoilla sijaitsevien pisteiden m ja n etäisyys maksimoituu.

Joustavan marginaalin luokitin

Joustavan marginaalin luokitin (C-SVM) optimimarginaaliluokittimen yleistys ei-separoituville näytejoukoille. Lähes kaikissa käytännön luokittelutehtävissä näytejoukot ovat ei-separoituvia. Ei-separoituvuus huomioidaan määrittämällä jokaiselle väärään luokkaan luokittuvalle näytteelle ns. slack-vakio γ i {\displaystyle \gamma _{i}} , mikä on näytteen etäisyys päätöspintaan. Vakio γ i = 0 {\displaystyle \gamma _{i}=0} , kun näyte luokittuu oikein.

Kahteen luokkaan luokitteleva joustavan marginaalin luokitin määritetään ratkaisemalla neliöllinen optimointiongelma

min 1 2 w 2 + C i = 1 l γ i {\displaystyle \min {\frac {1}{2}}\|\mathbf {w} \|^{2}+C\sum _{i=1}^{l}\gamma _{i}}

muuttujien w {\displaystyle \mathbf {w} } , b {\displaystyle b} ja γ i {\displaystyle \gamma _{i}} suhteen kun reunaehdot ovat

y i ( w T x i + b ) 1 γ i {\displaystyle y_{i}(\mathbf {w} ^{T}\mathbf {x} _{i}+b)\geq 1-\gamma _{i}}

ja

γ i 0 {\displaystyle \gamma _{i}\geq 0}

kun i = 1 , , l {\displaystyle i=1,\dots ,l} . Ongelma voidaan ratkaista ratkaisemalla Lagrangen menetelmällä konstruoitu duaaliongelma

min 1 2 i = 1 l j = 1 l α i α j K ( x i , x j ) C i = 1 l α i {\displaystyle \min {\frac {1}{2}}\sum _{i=1}^{l}\sum _{j=1}^{l}\alpha _{i}\alpha _{j}K(\mathbf {x_{i}} ,\mathbf {x_{j}} )-C\sum _{i=1}^{l}\mathbf {\alpha } _{i}}

missä C {\displaystyle C} on regularisointitermi ja reunaehdot ovat

α i 0 {\displaystyle \mathbf {\alpha } _{i}\geq 0}

ja

y T α = 0 {\displaystyle \mathbf {y} ^{T}\mathbf {\alpha } =0}

missä α i {\displaystyle \mathbf {\alpha } _{i}} ovat Lagrangen kertoimet ja y = [ y 1 , , y l ] T {\displaystyle \mathbf {y} =[y_{1},\dots ,y_{l}]^{T}} .

Epälineaarinen tukivektorikone

Coverin teoreeman mukaan kaksi ei-separoituvaa näytejoukkoa separoituvat suuremmalla todennäköisyydellä, kun ne kuvataan epälineaarisesti korkeampiulotteiseen avaruuteen. Tukivektorikoneen suorituskykyä voidaan parantaa kuvaamalla piirreavaruus kernelifunktiolla K ( x , y ) {\displaystyle K(\mathbf {x} ,\mathbf {y} )} , joka implisiittisesti määrittelee kuvauksen.

Mercerin teorian perusteella kernelifunktiolle K {\displaystyle K} on olemassa sisätuloavaruus jos matriisi M {\displaystyle M} on positiivisemidefiniitti ja K {\displaystyle K} on symmetrinen funktio, kun M i j = K ( x i , x j ) {\displaystyle M_{ij}=K(\mathbf {x} _{i},\mathbf {x} _{j})} kaikille i , j = 1 , , l {\displaystyle i,j=1,\dots ,l} . Luokittelu suoritetaan sisätuloavaruudessa lineaarisesti, mutta päätöspinta kuvautuu piirreavaruuteen epälineaariseksi, kuten ympyräksi tai ellipsiksi radiaalisella kernelifunktiolla. Mercerin teoria mahdollistaa epälineaarisen luokittimen määrittämisen intuitiivisesti.

Mielivaltainen piste x {\displaystyle \mathbf {x} } luokitellaan luokkaan + 1 {\displaystyle +1} , kun

d = i = 1 l α i K ( x i , x ) > 0 {\displaystyle d=\sum _{i=1}^{l}\alpha _{i}K(\mathbf {x} _{i},\mathbf {x} )>0}

ja muussa tapauksessa luokkaan 1 {\displaystyle -1} .

Havaitaan, että kun α i > 0 {\displaystyle \alpha _{i}>0} , niin x i {\displaystyle \mathbf {x} _{i}} on tukivektori. Näin ollen luokittelu riippuu ainoastaan tukivektoreista.

Kernelifunktiota

Yleisimpiä kernelifunktiota ovat

  • polynominen kernelifunktio K ( x i , x ) = ( x T x i + 1 ) p {\displaystyle K(\mathbf {x} _{i},\mathbf {x} )=(\mathbf {x} ^{T}\mathbf {x} _{i}+1)^{p}} , kun p {\displaystyle p} on positiivinen kokonaisluku
  • radiaalinen kernelifunktio K ( x i , x ) = exp ( 1 2 σ 2 x x i 2 ) {\displaystyle K(\mathbf {x} _{i},\mathbf {x} )=\exp(-{\frac {1}{2\sigma ^{2}}}\|\mathbf {x} -\mathbf {x} _{i}\|^{2})} , kun σ > 0 {\displaystyle \sigma >0}

Polynominen kernelifunktio erottelee piirteet lineaarisesti lähtöavaruudessa kun taas radiaalinen funktio erottelee piirteet lineaarisesti tietyssä korkeaulotteisessa avaruudessa. Lähtöavaruuden näkökulmasta radiaalinen kernelifunktio erottelee näytteet epälineaarisesti ja se voi näin ollen sen suorituskyky on joissain tilainteissa polynomista luokitinta parempi.

Tukivektorikonetoteutuksia

  • SVM light

Katso myös

  • Neuroverkot
  • Hahmontunnistus

Lähteet

  • Haykin, Simon: Neural Networks - A comprehensive foundation 2nd edition, Prentice-Hall, 1999

Viitteet

  1. C. Cortes, V. N. Vapnik: Support Vector Networks. Machine Learning, 1995, nro Volume 20, Number 3, s. 273–297.