
Në këtë artikull, ne do të shqyrtojmë teoritë rreth transformimit funksioni i regresionit linear në funksionin e transformimit logit invers (ndryshe i njohur si funksioni i përgjigjes logjistike). Më pas, duke shfrytëzuar arsenalin metodës së maksimumit të mundësisë, sipas modelit të regresionit logjistik, do të nxjerrim funksionin e humbjes Humbja Logjistike, ose, me fjalë të tjera, do të përcaktojmë funksionin që, në modelin e regresionit logjistik, përshtat parametrat e vektorit të peshave
.
Plani i artikullit:
- Do të përkujtojmë mbi lidhjen e drejtpërdrejtë midis dy variablave
- Do të identifikojmë nevojën për transformimin funksioni i regresionit linear
në funksionin e përgjigjes logjistike 
- Do të kryejmë transformime dhe do të nxjerrim funksionin e përgjigjes logjistike
- Do të përpiqemi të kuptojmë se përse metoda e katrorëve të vogla është e papërshtatshme për përshtatjen e parametrave
funksione Humbja Logjistike - Përdorim metoda e maksimumit të mundësisë për të identifikuar funksionet e përshtatjes së parametrave
:5.1. Rasti 1: funksioni Humbja Logjistike për objektet me etiketat e klasave 0 dhe 1:

5.2. Rasti 2: funksioni Humbja Logjistike për objektet me etiketat e klasave -1 dhe +1:

Artikulli Ă«shtĂ« i mbushur me shembuj tĂ« thjeshtĂ«, ku tĂ« gjitha llogaritjet mund tĂ« kryhen lehtĂ«sisht gojarisht ose nĂ« letĂ«r, ndonjĂ«herĂ« mund tĂ« nevojitet njĂ« kalkulator. Prandaj, pĂ«rgatituni đ
Ky artikull është kryesisht i orientuar ndaj shkencëtarëve të të dhënave me njohuri fillestare në parimet e mësimit të makinerive.
NĂ« artikull do tĂ« jepen gjithashtu kodet pĂ«r vizatimin e grafikĂ«ve dhe llogaritjet. TĂ« gjitha kodet janĂ« tĂ« shkruara nĂ« gjuhĂ«n python 2.7. Le ta sqaroj pĂ«r "novitetin" e versionit tĂ« pĂ«rdorur â kjo Ă«shtĂ« njĂ« nga kĂ«rkesat pĂ«r tĂ« kaluar kursin e njohur nga Yandex nĂ« njĂ« nga platformat mĂ« tĂ« njohura tĂ« arsimit online Coursera, dhe, siç mund tĂ« supozohet, materiali Ă«shtĂ« pĂ«rgatitur nĂ« frymĂ«n e kĂ«tij kursi.
01. Lidhja e drejtpërdrejtë
Ă«shtĂ« plotĂ«sisht e arsyeshme tĂ« bĂ«jmĂ« pyetjen â çfarĂ« ka tĂ« bĂ«jĂ« lidhja e drejtpĂ«rdrejtĂ« me regresionin logjistik?
Të gjitha janë të thjeshta! Regresioni logjistik është një nga modelet që i përket klasifikatorit linear. Thjesht, detyra e klasifikatorit linear është parashikimi i vlerave të synuara
nga variablat (regresorët)
. Në këtë rast, supozohet se lidhja midis karakteristikave
dhe vlerave të synuara
Ă«shtĂ« lineare. Prandaj quhet klasifikatori i tillĂ« â linear. NĂ«se duam tĂ« pĂ«rgjithĂ«sojmĂ«, nĂ« themel tĂ« modelit tĂ« regresionit logjistik Ă«shtĂ« supozimi pĂ«r prani tĂ« njĂ« lidhjeje lineare midis karakteristikave
dhe vlerave të synuara
. Ja, ky është lidhja.
NĂ« studio Ă«shtĂ« shembulli i parĂ«, dhe ai, saktĂ«sisht, pĂ«r lidhjen e drejtpĂ«rdrejtĂ« tĂ« variablave tĂ« hulumtuar. NĂ« procesin e pĂ«rgatitjes sĂ« artikullit, kam hasur nĂ« njĂ« shembull qĂ« Ă«shtĂ« bĂ«rĂ« tashmĂ« i njohur pĂ«r shumĂ« â lidhja midis forcĂ«s sĂ« rrymĂ«s dhe tensionit ("Analiza e Regresionit tĂ« Aplikuar", N.Dreiper, G.Smith). KĂ«tĂ« do ta shqyrtojmĂ« kĂ«tu.
Sipas me ligjin e Ohmit:
, ku
â forca e rrymĂ«s,
â tensioni,
â rezistenca.
Po të mos e dinim ligjin e Ohmit, mund të gjenim lidhjen empirike, duke ndryshuar
dhe duke matur
, duke mbajtur në të njëjtën kohë
konstant. Atëherë do të shihnim se grafiku i varësisë
nga
tregon një linjë më të vogël ose më shumë, e cila kalon përmes origjinës së koordinatave. Thame "më shumë ose më pak", sepse, megjithëse lidhja është në të vërtetë e saktë, masat tona mund të kenë gabime të vogla, dhe për këtë arsye pikat në grafik, mund të mos bien saktësisht mbi vijë, por do të shpërndahen rreth saj në mënyrë të rastësishme.
Grafiku 1 "Lidhja
nga
»

Kodi për vizatimin e grafikut
import matplotlib.pyplot as plt
%matplotlib inline
import numpy as np
import random
R = 13.75
x_line = np.arange(0,220,1)
y_line = []
for i in x_line:
y_line.append(i/R)
y_dot = []
for i in y_line:
y_dot.append(i+random.uniform(-0.9,0.9))
fig, axes = plt.subplots(figsize = (14,6), dpi = 80)
plt.plot(x_line,y_line,color = 'purple',lw = 3, label = 'I = U/R')
plt.scatter(x_line,y_dot,color = 'red', label = 'Rezultatet aktuale')
plt.xlabel('I', size = 16)
plt.ylabel('U', size = 16)
plt.legend(prop = {'size': 14})
plt.show()02. Nevoja për transformime të ekuacionit të regresionit linear
Le të shqyrtojmë një shembull tjetër. Le të parafytyrojmë se ne punojmë në një bankë dhe para nesh qëndron detyra për të përcaktuar probabilitetin e shpërndarjes së një kredi nga kredimarrësi në varësi të disa faktorëve. Për të thjeshtuar detyrën, le të shqyrtojmë vetëm dy faktorë: paga mujore e kredimarrësit dhe shuma e pagesës mujore për shlyerjen e kredisë.
Detyra është shumë hipotetike, por me këtë shembull do të mund të kuptojmë pse për zgjidhjen e saj nuk është e mjaftueshme aplikimi funksioni i regresionit linear, si dhe do të zbulojmë cilat transformime kërkohen për funksionin.
Kthehemi në shembullin. E qartë se sa më e lartë të jetë paga, aq më shumë mund të drejtojë borxhi çdo muaj për shlyerjen e kredisë. Megjithatë, për një farë diapazoni pagash, kjo varësi do të jetë mjaft lineare. Për shembull, le të marrim një diapazon pagash nga 60.000R deri në 200.000R dhe të supozojmë se në këtë diapazon pagash, varësia e pagesës mujore ndaj pagës është lineare. Le të supozojmë se në këtë diapazon pagash është konstatuar se raporti i pagës ndaj pagesës nuk mund të jetë më i vogël se 3 dhe gjithashtu borxhli duhet të ketë ende 5.000R brenda. Vetëm në këtë rast, ne do të mendojmë se borxhli do ta kthejë kredinë në bankë. Atëherë, ekuacioni i regresionit linear do të marrë këtë formë:

ku
,
,
,
â paga
-t e borxhlisë,
â pagesa e kredisĂ«
-t e borxhlisë.
Duke futur në ekuacionin pagën dhe pagesën e kredisë me parametra të fiksuar
mund të merret një vendim për dhënien ose refuzimin e kredisë.
Duke u shkuar përpara, vëmë re se, me parametrat e dhëna
funkcioni i regresionit linear, i përdorur në funkcioni i reagimit logjik do të jep rezultate të mëdha, të cilat do të vështirësojnë kryerjen e llogaritjeve për të përcaktuar probabilitetet e shlyerjes së kredisë. Prandaj, propozohet që të zvogëlohen koeficientët tanë, le të themi, me 25.000 herë. Nga ky transformim në koeficientë, vendimi për dhënien e kredisë nuk do të ndryshojë. Ta mbajmë këtë moment në mend për të ardhmen, dhe tani, për ta bërë edhe më të qartë, le të shqyrtojmë situatën me tre borxhlinj potencialë.
Tabela 1 "Borxhlinj potencialë"

Kodi për formimin e tabelës
import pandas as pd
r = 25000.0
w_0 = -5000.0/r
w_1 = 1.0/r
w_2 = -3.0/r
data = {'Borxhli':np.array(['Vasya', 'Fedya', 'Lesha']),
'Paga':np.array([120000,180000,210000]),
'Pagesa':np.array([3000,50000,70000])}
df = pd.DataFrame(data)
df['f(w,x)'] = w_0 + df['Paga']*w_1 + df['Pagesa']*w_2
decision = []
for i in df['f(w,x)']:
if i > 0:
dec = 'Miratuar'
decision.append(dec)
else:
dec = 'Refuzim'
decision.append(dec)
df['Vendimi'] = decision
df[['Borxhli', 'Paga', 'Pagesa', 'f(w,x)', 'Vendimi']]Sipas të dhënave të tabelës, Vasya me një pagë prej 120.000R dëshiron të marrë një kredi të tillë, në mënyrë që ta paguajë atë me 3.000R çdo muaj. Ne kemi përcaktuar se për të miratuar kredinë, niveli i pagës së Vasya duhet të tejkalojë tre herë shumën e pagesës, dhe që të mbetet ende 5.000R. Ky kërkesë Vasya e përmbush:
. Mbetej edhe 106.000R. Pavarësisht se në llogaritjen tonë
ne zvogëluam koeficientët
me 25.000 herĂ«, rezultati ishte po ai â kredia mund tĂ« miratohet. Fedya gjithashtu do tĂ« marrĂ« kredit, por Lesha, pavarĂ«sisht se merr mĂ« shumĂ« se tĂ« tjerĂ«t, do tĂ« duhet tĂ« kufizojĂ« apetitat e tij.
Të vizatojmë një grafik për këtë rast.
Grafiku 2 "Klasifikimi i borxhlinjve"

Kodi për vizatimin e grafikut
salary = np.arange(60000,240000,20000)
payment = (-w_0-w_1*salary)/w_2
fig, axes = plt.subplots(figsize = (14,6), dpi = 80)
plt.plot(salary, payment, color = 'grey', lw = 2, label = '$f(w,x_i)=w_0 + w_1x_{i1} + w_2x_{i2}$')
plt.plot(df[df['Vendimi'] == 'Miratuar']['Paga'], df[df['Vendimi'] == 'Miratuar']['Pagesa'],
'o', color ='green', markersize = 12, label = 'Vendimi - Kredia e miratuar')
plt.plot(df[df['Vendimi'] == 'Refuzim']['Paga'], df[df['Vendimi'] == 'Refuzim']['Pagesa'],
's', color = 'red', markersize = 12, label = 'Vendimi - Refuzimi i kredisë')
plt.xlabel('Paga', size = 16)
plt.ylabel('Pagesa', size = 16)
plt.legend(prop = {'size': 14})
plt.show()Prandaj, linea jonë, e ndërtuar sipas funksionit
, ndan "borxhlinjtĂ« e kĂ«qij" nga "borxhlinjtĂ« e mirĂ«". Ata borxhlinj qĂ« dĂ«shirat e tyre nuk pĂ«rputhen me mundĂ«sitĂ« janĂ« mbi linjĂ« (Lesha), ata qĂ« janĂ« nĂ« gjendje, sipas parametra tonĂ« modeli, tĂ« kthejnĂ« kredisĂ«, janĂ« nĂ«n linjĂ« (Vasya dhe Fedya). NĂ«se mund ta themi ndryshe â linea jonĂ« ndan borxhlinjtĂ« nĂ« dy klasa. Tâi quajmĂ« kĂ«shtu: nĂ« klasĂ«n
do të përfshijmë borxhlinjtë që do të mund të kthejnë kredisë, ndërsa në klasën
ose
do të përfshijmë borxhlinjtë që ndoshta nuk do të mund ta kthejnë kredinë.
Të përmbledhim konkluzionet nga kjo shembull i thjeshtë. Le të marrim pikën
dhe, duke futur koordinatat e pikës në ekuacionin përkatës të linjës
, të shqyrtojmë tri variante:
- Nëse pika ndodhet nën linjë, dhe ne e klasifikojmë atë në klasën
, atëherë vlera e funksionit
do të jetë pozitive nga
deri në
. Kështu ne mund të mendojmë se probabiliteti i shlyerjes së kredisë, është brenda
. Sa më e lartë të jetë vlera e funksionit, aq më e lartë është probabiliteti. - Nëse pika ndodhet mbi linjë dhe ne e klasifikojmë atë në klasën
ose
, atëherë vlera e funksionit do të jetë negative nga
deri në
. Atëherë do të mendojmë se probabiliteti i shlyerjes së borxhit është brenda
dhe, sa më i lartë në modulin e vlerës së funksionit, aq më e lartë është siguria jonë. - Pika ndodhet në linjë, në kufirin midis dy klasave. Në këtë rast vlera e funksionit
do të jetë e barabartë me
dhe probabiliteti i shlyerjes së kredisë është i barabartë me
.
Tani, imagjinoni se ne nuk kemi dy faktorë, por disa, borxhlinj nuk janë tre, por mijëra. Atëherë në vend të linjës ne do të kemi pllakë m-dimensional edhe koeficientët
Ne do të marrim të dhëna rastësore, por do t'i nxjerrim sipas rregullave, bazuar në të dhënat e grumbulluara për huamarrësit që kanë kthyer ose nuk kanë kthyer kreditë. Në të vërtetë, vërejeni, tani jemi duke përzgjedhur huamarrësit me koeficientë të njohur.
Sidoqoftë, objekti i modelit të regresionit logjistik është pikërisht të përcaktojë parametrat
, për të cilat vlera e funksionit të humbjeve Humbja Logjistike do të përpiqet të minimizohet. Por se si llogaritet vektori
, do ta mësojmë në seksionin e pestë të artikullit. Ndërkohë, le të kthehemi te banka jonë dhe tre klientët e saj.
Falë funksionit
Ne e dimĂ« se kujt mund t'i jepet kredi dhe kujt duhet t'i refuzohet. Por me kĂ«tĂ« informacion, nuk mund tĂ« shkojmĂ« te drejtori, sepse ata donin tĂ« merrnin probabilitetin e kthimit tĂ« kreditĂ«ve nga çdo huamarrĂ«s. ĂfarĂ« tĂ« bĂ«jmĂ«? PĂ«rgjigjja Ă«shtĂ« e thjeshtĂ« â na nevojitet tĂ« transformojmĂ« funksionin
, vlerat e të cilit bien në gamën
në një funksion, vlerat e të cilit do të bien në gamën
. Dhe një funksion i tillë ekziston, e quajnë funksionin e përgjigjes logjistike ose transformimin invers-logit. Prezantim:

Le të shohim hap pas hapi si formohet funksioni i përgjigjes logjistike. Vërejmë se do të ecim në anën e kundërt, dmth. do të supozojmë se na është dhënë vlera e probabilitetit, e cila bie brenda kufijve nga
deri në
dhe më pas do ta "zgjerojmë" këtë vlerë në të gjithë intervalin e numrave nga
deri në
.
03. Formojmë funksionin e përgjigjes logjistike
Hapi 1. Të përkthejmë vlerat e probabilitetit në gamën 
Për kohën e transformimit të funksionit
në funksionin e përgjigjes logjistike
do ta lëmë pa përmendur analistin tonë të kredive, dhe përkundrazi do të shqyrtojmë bastet në zyrat e basteve. Jo, natyrisht, ne nuk do të vendosim bast, ajo që na intereson është vetëm kuptimi i shprehjes, për shembull, probabiliteti 4 me 1. Shanset, të njohura për të gjithë lojtarët e basteve, përfaqësojnë raportin e "sukseseve" ndaj "dështimeve". Nga pikëpamja e probabiliteteve, shanset janë probabiliteti i shfaqjes së ngjarjes, i ndarë me probabilitetin që ngjarja nuk do të ndodhë. Le të shkruajmë formulën e probabilitetit të shfaqjes së ngjarjes
:

, ku
â probabiliteti i shfaqjes sĂ« ngjarjes,
â probabiliteti qĂ« ngjarja NUK do tĂ« ndodhĂ«
Për shembull, nëse probabiliteti që një kalë i ri, i fortë dhe energjik i quajtur "Erë" të kalojë kalë të vjetër dhe të dobët me emrin "Matilda" është
, atëherë shanset për suksesin e "Erës" do të jenë
për
dhe anasjelltas, duke ditur shanset, nuk do të na jetë e vështirë të llogarisim probabilitetin
:

Në këtë mënyrë, ne kemi mësuar të "përkthejmë" probabilitetin në shanse, të cilat pranojnë vlera nga
deri në
. Të bëjmë një hap tjetër dhe të mësojmë "të përkthejmë" probabilitetin në të gjithë vijën numerike nga
deri në
.
Hapi 2. Të përkthejmë vlerat e probabilitetit në gamën 
Ky hap Ă«shtĂ« shumĂ« i thjeshtĂ« â do ta logarizojmĂ« shanset me bazĂ«n e numrit tĂ« Euler-it
dhe do të marrim:

Tani e dimë se nëse
, atëherë llogaritja e vlerës
do të jetë shumë e thjeshtë dhe, më tepër, duhet të jetë pozitive:
. Kështu është.
Për kuriozitet, le të kontrollojmë se nëse
, atëherë presim të shohim një vlerë negative
. Kontrollojmë:
. E saktë.
Tani e dimë si të përkthejmë vlerën e probabilitetit nga
deri në
në të gjithë vijën numerike nga
deri në
. Në hapin tjetër do të bëjmë gjithçka në të kundërt.
Ndërkohë, vërejmë se sipas rregullave të logaritmimit, duke ditur vlerën e funksionit
, mund të llogarisim shanset:

Ky metodë për përcaktimin e shanseve do të na ndihmojë në hapin tjetër.
Hapi 3. Do të nxjerrim formulën për përcaktimin e 
Pra, kemi mësuar të gjejmë vlerat e funksionit duke pasur
. Por, nĂ« tĂ« vĂ«rtetĂ«, ajo qĂ« na nevojitet Ă«shtĂ« e kundĂ«rta â duke ditur vlerĂ«n
për të gjetur
. Për këtë, do t'i referohemi konceptit të funksionit invers të shanseve, sipas të cilit:
Në këtë artikull nuk do të nxjerrim formulën e mësipërme, por do ta testojmë me numra nga shembulli më parë. E dimë se kur shanset janë 4 me 1 (

), probabiliteti i shfaqjes së ngjarjes është 0.8 (
). Do të bëjmë një zëvendësim:
. Kjo përputhet me llogaritjet tona të bëra më parë. Le të vazhdojmë.
Në hapin e kaluar nxorrëm se
, prandaj mund të bëjmë një zëvendësim në funksionin invers të shanseve. Do të marrim:
Të ndajmë si numëruesin dhe emëruesin me

, atëherë:
Për çdo rast, për t'u siguruar se nuk kemi bërë ndonjë gabim, do të bëjmë një kontroll të vogël. Në hapin 2, ne përcaktuam se

. Atëherë, duke zëvendësuar vlerën
në funksionin e përgjigjes logjistike, presim të marrim
. Zëvendësojmë dhe marrim:
Urime, lexues i nderuar, sapo kemi nxjerrë dhe testuar funksionin e përgjigjes logjistike. Le të shikojmë grafikun e funksionit.
Grafiku 3 "Funksioni i përgjigjes logjistike" 
Urime ju, lexues i nderuar, sapo kemi nxjerrë dhe testuar funksionin e përgjigjes logistike. Le të shohim grafikën e funksionit.
Grafiku 3 «Funksioni i përgjigjes logistike»

Kodi për vizatimin e grafikut
import math
def logit (f):
return 1/(1+math.exp(-f))
f = np.arange(-7,7,0.05)
p = []
for i in f:
p.append(logit(i))
fig, axes = plt.subplots(figsize = (14,6), dpi = 80)
plt.plot(f, p, color = 'grey', label = '$ 1 / (1+e^{-w^Tx_i})$')
plt.xlabel('$f(w,x_i) = w^Tx_i$', size = 16)
plt.ylabel('$p_{i+}$', size = 16)
plt.legend(prop = {'size': 14})
plt.show()Në literaturë mund të takoni gjithashtu emrin e kësaj funksioni si funksioni sigmoid. Nga grafiku, është e qartë se ndryshimi kryesor i probabilitetit të përkatësisë së objektit në klasë ndodh në një gamë relativisht të vogël
, diku nga
deri në
.
Propozoj qĂ« t'i kthehemi analistit tonĂ« tĂ« kreditit dhe t'i ndihmojmĂ« atij nĂ« llogaritjen e probabilitetit tĂ« shlyerjes sĂ« kredive, pĂ«rndryshe ai rrezikon tĂ« mbetet pa bonus đ
Tabela 2 "Huamarrësit potencialë"

Kodi për formimin e tabelës
proba = []
for i in df['f(w,x)']:
proba.append(round(logit(i),2))
df['Probability'] = proba
df[['Huamarrësi', 'Paga', 'Pagesa', 'f(w,x)', 'Vendimi', 'Probabiliteti']]Pra, probabiliteti i kthimit të kredive është përcaktuar. Në përgjithësi, kjo duket e arsyeshme.
Në të vërtetë, probabiliteti që Vasja, me një pagë prej 120.000R, të mund të kthejë çdo muaj 3.000R në bankë është afër 100%. Për më tepër, duhet të kuptojmë se banka mund t'i japë një kredi edhe Lleshës, në rast se politika e bankës parashikon, për shembull, të kreditohen klientët me një probabilitet kthimi të kredive më shumë se, le të themi, 0.3. Thjesht në këtë rast, banka do të formojë një rezervë më të madhe për humbjet e mundshme.
Gjithashtu, duhet të theksohet se raporti midis pagës dhe pagesës, jo më pak se 3 dhe me një rezerva prej 5.000R, ishte marrë nga ajri. Prandaj, ne nuk mund ta kishim përdorur në formën e saj fillestare vektorët e peshave
. Ne kishim nevojë që t'i zvogëlojmë ndjeshëm koeficientët dhe në këtë rast ndamë çdo koeficient në 25.000, domethënë, në thelb e përshtatëm rezultatin. Por kjo u bë qëllimisht për të thjeshtuar kuptimin e materialit në fazën fillestare. Në jetën reale, na nevojitet të gjejmë koeficientët, jo t'i shpikim ose t'i përshtasim. Në seksionet e ardhshme të artikullit, do të nxjerrim ekuacionet me ndihmën e të cilave përcaktohen parametrat
.
04. Metoda e katërzuar e minimit të kuadrateve për përcaktimin e vektorit të peshave
në funksionin e përgjigjes logjistike
Metoda e përcaktimit të vektorit të peshave është tashmë e njohur për ne
, si metoda e katërzuar e minimit (MKK) dhe, përse të mos e përdorim atë në problemet e klasifikimit binar? Në të vërtetë, asgjë nuk e ndalon përdorimin MNMK, por kjo metodë në problemet e klasifikimit jep rezultate më pak të sakta se Humbja Logjistike. Kjo ka një bazë teorike. Le të shohim fillimisht një shembull të thjeshtë.
Supozoni se modelet tona (që përdorin MSE dhe Humbja Logjistike) kanë filluar të përshtatin vektorin e peshave
dhe ne e ndaluam llogaritjen në një hap të caktuar. Ka pak rëndësi nëse është në mes, në fund ose në fillim, ajo që është e rëndësishme është se ne tashmë kemi disa vlera të vektorit të peshave dhe le të supozojmë se në këtë hap, vektori i peshave
për të dyja modelet nuk ka ndryshime. Pastaj të marrim peshat e marra dhe t'i vendosim ato në funksionin e përgjigjes logjistike (
) për ndonjë objekt që i përket klasës
. Le tĂ« shqyrtojmĂ« dy raste, kur sipas vektorit tĂ« pĂ«rshtatur tĂ« peshave, modeli ynĂ« bĂ«n njĂ« gabim tĂ« madh dhe pĂ«rkundrazi â modeli Ă«shtĂ« shumĂ« i sigurt se objekti i pĂ«rket klasĂ«s
. Le të shohim se çfarë dënimesh do të "jepen" në përdorimin e MNMK dhe Humbja Logjistike.
Kodi për llogaritjen e dënimeve në varësi të funksionit të humbjes së përdorur
# ĐșлаŃŃ ĐŸĐ±ŃĐ”ĐșŃа
y = 1
# ĐČĐ”ŃĐŸŃŃĐœĐŸŃŃŃ ĐŸŃĐœĐ”ŃĐ”ĐœĐžŃ ĐŸĐ±ŃĐ”ĐșŃа Đș ĐșлаŃŃŃ ĐČ ŃĐŸĐŸŃĐČĐ”ŃŃŃĐČОО Ń ĐżĐ°ŃĐ°ĐŒĐ”ŃŃĐ°ĐŒĐž w
proba_1 = 0.01
MSE_1 = (y - proba_1)**2
print 'ĐšŃŃĐ°Ń MSE ĐżŃĐž ĐłŃŃĐ±ĐŸĐč ĐŸŃОбĐșĐ” =', MSE_1
# ĐœĐ°ĐżĐžŃĐ”ĐŒ ŃŃĐœĐșŃĐžŃ ĐŽĐ»Ń ĐČŃŃĐžŃĐ»Đ”ĐœĐžŃ f(w,x) ĐżŃĐž ОзĐČĐ”ŃŃĐœĐŸĐč ĐČĐ”ŃĐŸŃŃĐœĐŸŃŃĐž ĐŸŃĐœĐ”ŃĐ”ĐœĐžŃ ĐŸĐ±ŃĐ”ĐșŃа Đș ĐșлаŃŃŃ +1 (f(w,x)=ln(odds+))
def f_w_x(proba):
return math.log(proba/(1-proba))
LogLoss_1 = math.log(1+math.exp(-y*f_w_x(proba_1)))
print 'ĐšŃŃĐ°Ń Log Loss ĐżŃĐž ĐłŃŃĐ±ĐŸĐč ĐŸŃОбĐșĐ” =', LogLoss_1
proba_2 = 0.99
MSE_2 = (y - proba_2)**2
LogLoss_2 = math.log(1+math.exp(-y*f_w_x(proba_2)))
print '**************************************************************'
print 'ĐšŃŃĐ°Ń MSE ĐżŃĐž ŃОлŃĐœĐŸĐč ŃĐČĐ”ŃĐ”ĐœĐœĐŸŃŃĐž =', MSE_2
print 'ĐšŃŃĐ°Ń Log Loss ĐżŃĐž ŃОлŃĐœĐŸĐč ŃĐČĐ”ŃĐ”ĐœĐœĐŸŃŃĐž =', LogLoss_2Rasti i gabimit tĂ« rĂ«ndĂ« â modeli e kategorizon objektin nĂ« klasĂ«n
me probabilitet 0.01
Dënimi duke përdorur MNMK do të jetë:

Dënimi duke përdorur Humbja Logjistike do të jetë:

Rasti me siguri tĂ« fortĂ« â modeli e kategorizon objektin nĂ« klasĂ«n
me probabilitet 0.99
Dënimi duke përdorur MNMK do të jetë:

Dënimi duke përdorur Humbja Logjistike do të jetë:

Ky shembull ilustron mirë se në rastin e një gabimi të rëndë, funksioni i humbjes Log Loss dënon modelin shumë më fort se MSE. Le të kuptojmë tani se cilat janë parakushtet teorike për përdorimin e funksionit të humbjes Log Loss në problemet e klasifikimit.
05. Metoda e maksimalizimit të ndershmërisë dhe regresioni logjistik
Siç kishim premtuar nĂ« fillim, artikulli Ă«shtĂ« plot me shembuj tĂ« thjeshtĂ«. NĂ« studion tonĂ« kemi njĂ« shembuj tĂ« rinj dhe mysafirĂ«t e vjetĂ«r â huamarrĂ«sit e bankĂ«s: Vasja, Fedja dhe Llesha.
Për çdo rast, para se të zhvillojmë shembullin, le të kujtojmë se në jetë kemi të bëjmë me një mostër trajnimi nga mijëra ose miliona objekte me dhjetëra ose qindra karakteristika. Megjithatë, këtu shifrat janë marrë kështu që ato lehtë përputhen në mendjen e një të saposhtun në fushën e të dhënave.
Po kthehemi në shembullin. Le të supozojmë se drejtori i bankës vendosi të jep kredi të gjithëve që kanë nevojë, pavarësisht se algoritmi sugjeronte të mos ia jepte atij Leshës. Pas një kohe, mësojmë se cili nga tre heronjtë e pagoi kredinë dhe cili jo. Siç u prit, Vasia dhe Fedja e paguan kredinë, ndërsa Lesha jo. Tani le të imagjinojmë që ky rezultat do të shërbejë si një grup i ri edukativ për ne dhe, në të njëjtën kohë, sikur janë zhdukur të gjitha të dhënat mbi faktorët që ndikojnë në probabilitetin e shlyerjes së kredisë (paga e huamarrësit, shuma e pagesës mujore). Atëherë intuitivisht mund të mendojmë se çdo huamarrës i tretë nuk e kthen kredinë, ose në fjalë të tjera probabiliteti i kthimit të kredisë nga huamarrësi i ardhshëm
. Ky supozim intuitiv ka një konfirmim teorik dhe është i bazuar në metodën e maksimalit të mundësisë, shpesh në literaturë quhet principi i maksimalit të mundësisë.
Në fillim, le të njihemi me terminologjinë.
MundĂ«sia e grupit â Ă«shtĂ« probabiliteti i marrjes sĂ« kĂ«tij grupi specifik, i marrjes sĂ« kĂ«tyre vĂ«zhgimeve / rezultateve, domethĂ«nĂ« produkti i probabiliteteve tĂ« marrjes sĂ« secilit nga rezultatet e grupit (pĂ«r shembull, paguar ose jo paguar kredia nga Vasia, Fedja dhe Lesha nĂ« tĂ« njĂ«jtĂ«n kohĂ«).
Funksioni i mundësisë lidhet me mundësinë e grupit dhe vlerat e parametrave të shpërndarjes.
Në rastin tonë, grupi edukativ përbën një skemë të përgjithësuar Bernuli, në të cilën gjendja e rastësishme merr vetëm dy vlera:
ose
. Ndërkaq, mundësia e grupit mund të shkruhet si një funksion mundësie në varësi të parametrin
si nga:


Shkruan regjistrimin e mësipërm mund të interpretohet kështu. Probabiliteti i bashkët që Vasia dhe Fedja të paguajnë kredinë është i barabartë me
, probabiliteti që Lesha Tà MOS e paguajë kredinë është
(pasi ndodhi pikërisht MOS pagimi i kredisë), kështu që probabiliteti bashkëkombëtar i të tre ngjarjeve është
.
Metoda e maksimalit tĂ« mundĂ«sisĂ« â Ă«shtĂ« njĂ« metodĂ« e vlerĂ«simit tĂ« parametrave tĂ« panjohur duke maksimizuar funksionin e mundĂ«sisĂ«. NĂ« rastin tonĂ«, kĂ«rkohet tĂ« gjejmĂ« njĂ« vlerĂ«
, për të cilin
arrin maksimumin.
Prej nga vjen ideja â tĂ« kĂ«rkojmĂ« vlerĂ«n e parametrin tĂ« panjohur, pĂ«r tĂ« cilin funksioni i mundĂ«sisĂ« arrin maksimumin? Origjinat e ideve burojnĂ« nga pĂ«rfaqĂ«simi se grupi Ă«shtĂ« burimi i vetĂ«m i njohurive tona pĂ«r popullatĂ«n e pĂ«rgjithshme. Gjithçka qĂ« dimĂ« pĂ«r popullatĂ«n e pĂ«rgjithshme paraqitet nĂ« grup. Prandaj, gjithçka qĂ« mund tĂ« themi Ă«shtĂ« se grupi Ă«shtĂ« pasqyrimi mĂ« i saktĂ« i popullatĂ«s sĂ« pĂ«rgjithshme qĂ« kemi. Prandaj, na nevojitet njĂ« parametr i tillĂ«, pĂ«r tĂ« cilin grupi ekzistues bĂ«het mĂ« e mundshme.
Në mënyrë të qartë, kemi të bëjmë me një problem optimizimi, në të cilin duhet të gjejmë pikën ekstreme të funksionit. Për të gjetur pikën ekstreme, nevojitet të shqyrtojmë kushtin e rendit të parë, do të thotë të barazojmë derivatin e funksionit me zero dhe të zgjidhim ekuacionin në lidhje me parametrin e kërkuar. Megjithatë, kërkimi i derivatit të produktit të shumë faktorëve mund të jetë një punë e gjatë, për të shmangur këtë, ekziston një teknikë speciale - kalimi në logaritem. funksionin e mundësisë. Pse është i mundshëm një kalim i tillë? Le të vëmë re se ne jemi në kërkim jo të vetë ekstremit të funksionit
, por të pikës ekstreme, do të thotë atë vlerë të parametrin të panjohur
, për të cilin
arrin maksimumin. Kur kalojmë në logaritem, pika ekstreme nuk ndryshon (ndonëse vetë ekstremi do të ndryshojë), pasi logaritemi është një funksion monoton.
Tani, në përputhje me më sipër, le të vazhdojmë të zhvillojmë shembullin tonë me kreditë për Vasin, Fedin dhe Leshën. Në fillim, le të kalojmë në logaritemin e funksionit të mundësisë:

Tani mund të lehtësojmë diferencimin e shprehjes në
:

Dhe përfundimisht, le të shqyrtojmë kushtin e rendit të parë - do ta barazojmë derivatin e funksionit me zero:

Si rezultat, vlerësimi ynë intuitiv për probabilitetin e shlyerjes së kredisë
ishte teorikisht i justifikuar.
Shumë mirë, por çfarë duhet të bëjmë tani me këtë informacion? Nëse do të mendojmë se çdo huamarrës i tretë nuk do ta kthejë paratë në bankë, banka përfundimisht do të falimentojë. E vërteta është, megjithatë, se kur vlerësojmë probabilitetin e shlyerjes së kredisë si
ne nuk e kemi marrë parasysh faktorët që ndikojnë në kthimin e kredisë: paga e huamarrësit dhe shuma e pagesës mujore. Le të kujtojmë se më parë llogaritëm probabilitetin e shlyerjes së kredisë nga çdo klient duke marrë parasysh këta faktorë.
.
Le të përcaktojmë plausibilitetin e mostrave:
Kodi për llogaritjet e plausibilitetit të mostrave
from functools import reduce
def likelihood(y,p):
line_true_proba = []
for i in range(len(y)):
ltp_i = p[i]**y[i]*(1-p[i])**(1-y[i])
line_true_proba.append(ltp_i)
likelihood = []
return reduce(lambda a, b: a*b, line_true_proba)
y = [1.0,1.0,0.0]
p_log_response = df['Probability']
const = 2.0/3.0
p_const = [const, const, const]
print 'Plausibiliteti i mostrës me vlerën e konstantës p=2/3:', round(likelihood(y,p_const),3)
print '****************************************************************************************************'
print 'Plausibiliteti i mostrës me vlerën e llogaritur p:', round(likelihood(y,p_log_response),3)Plausibiliteti i mostrës me vlerën e konstantës
:

Plausibiliteti i mostrës në llogaritjen e probabilitetit për shlyerjen e kredisë duke marrë parasysh faktorët
:


Plausibiliteti i mostrĂ«s me probabilitetin e llogaritur nĂ« varĂ«si tĂ« faktorĂ«ve doli tĂ« ishte mĂ« i lartĂ« se plausibiliteti me vlerĂ«n e konstantĂ«s. ĂfarĂ« do tĂ« thotĂ« kjo? Kjo do tĂ« thotĂ« se njohuritĂ« mbi faktorĂ«t lejuan pĂ«rshtatjen mĂ« tĂ« saktĂ« tĂ« probabilitetit tĂ« shlyerjes sĂ« kredisĂ« pĂ«r secilin klient. Prandaj, kur jepet njĂ« kredi, Ă«shtĂ« mĂ« e saktĂ« tĂ« pĂ«rdoret modeli i vlerĂ«simit tĂ« probabilitetit tĂ« shlyerjes sĂ« borxhit, i sugjeruar nĂ« fund tĂ« seksionit tĂ« tretĂ« tĂ« artikullit.
Por atëherë, nëse na nevojitet maksimalizimi i funksionit të plausibilitetit të mostrës, atëherë pse të mos përdorim ndonjë algoritëm që do të jepte probabilitetet për Vasilin, Fedin dhe Leshin, për shembull, si 0.99, 0.99 dhe 0.01 përkatësisht. Ndoshta një algoritëm i tillë do të performonte mirë në mostrat e trajnimit, pasi do ta afronte vlerën e plausibilitetit të mostrës me
, por, sĂ« pari, njĂ« algoritĂ«m i tillĂ« do tĂ« ketĂ« me siguri vĂ«shtirĂ«si nĂ« kapacitetin e pĂ«rgjithshĂ«m, dhe sĂ« dyti, ky algoritĂ«m me siguri nuk do tĂ« jetĂ« linear. Dhe nĂ«se metodat pĂ«r luftimin e mbipopullimit (po aq sa dhe kapaciteti i dobĂ«t i pĂ«rgjithshĂ«m) nuk janĂ« padyshim nĂ« planin e kĂ«tij artikulli, atĂ«herĂ« pĂ«r pikĂ«n e dytĂ« le tĂ« kalojmĂ« mĂ« nĂ« hollĂ«si. PĂ«r kĂ«tĂ«, mjafton tĂ« pĂ«rgjigjemi nĂ« njĂ« pyetje tĂ« thjeshtĂ«. A mund tĂ« ketĂ« tĂ« njĂ«jtĂ«n probabilitet pĂ«r shlyerjen e kredisĂ« Vasil dhe Fed? Nga pikĂ«pamja e logjikĂ«s sĂ« shĂ«ndoshĂ«, sigurisht qĂ« jo, nuk mundet. Pra, pĂ«r shlyerjen e kredisĂ« Vasil do tĂ« jepte 2.5% tĂ« pagĂ«s sĂ« tij çdo muaj, ndĂ«rsa Fedi â pothuajse 27.8%. Gjithashtu, nĂ« grafikun 2 "Klasifikimi i klientĂ«ve" ne shohim se Vasil Ă«shtĂ« shumĂ« mĂ« larg nga line, qĂ« ndan klasat, se sa Fedi. Dhe pĂ«r mĂ« tepĂ«r, ne e dimĂ« se funksioni
për Vasilin dhe Fedin merr vlera të ndryshme: 4.24 për Vasilin dhe 1.0 për Fedin. Tani, po të thoshim se Fedi për shembull, fitonte shumë më tepër ose kërkonte një kredi më të vogël, atëherë probabilitetet e shlyerjes së kredisë për Vasilin dhe Fedin do të ishin të ngjashme. Me fjalë të tjera, nuk e mashtrohet varësia lineare. Dhe nëse ne do të kishim llogaritur me të vërtetë koeficientët
, e jo t'i kishim marrë nga tavani, atëherë do të mund të deklaronim se vlerat tona
ofrojnë më së miri vlerësimin e probabilitetit të shlyerjes së kredisë nga çdo huamarrës, por duke marrë parasysh se ne e kemi vendosur se përcaktimi i koeficientëve
Ă«shtĂ« kryer sipas rregullave tĂ« gjitha, kĂ«shtu qĂ« edhe ne do tĂ« konsiderojmĂ« â koeficientĂ«t tanĂ« lejojnĂ« tĂ« ofrojnĂ« vlerĂ«simin mĂ« tĂ« mirĂ« tĂ« probabilitetit đ
Megjithatë, ne u shpëtuam. Në këtë seksion ne duhet të kuptojmë se si përcaktohet vektori i peshave
, i cili është i nevojshëm për vlerësimin e probabilitetit të kthimit të kredisë nga çdo huamarrës.
Përmbledhim shkurtimisht se me çfarë arsenali do të dalim në kërkim të koeficientëve
:
1. Ne supozojmĂ« se varĂ«sia midis variablit tĂ« synuar (vlefta e parashikuar) dhe faktorĂ«ve qĂ« ndikojnĂ« nĂ« rezultat â Ă«shtĂ« lineare. PĂ«r kĂ«tĂ« arsye aplikohet funkcioni i regresionit linear nĂ« formĂ«
, linja e së cilës ndan objektet (klientët) në klasa
dhe
ose
(klientë, të cilët janë në gjendje të shlyejnë kredinë dhe ata që s'kanë mundësi). Në rastin tonë, ekuacioni ka formën
.
2. Ne përdorim funksionin e kthimit të logit-transformimit në formë
për të vlerësuar probabilitetin e përkatësisë së objektit në klasë
.
3. Ne e konsiderojmë mostrën tonë të trajnimit si një realizim të përgjeneralizuar të skemës së Bernulit, domethënë për çdo objekt gjenerohet një variabël e rastit, e cila me probabilitetin
(me probabilitetin e saj për çdo objekt) merr vlerën 1 dhe me probabilitetin
â 0.
4. Ne e dimë se na nevojitet maksimalizimi funksionit të plausibilitetit të mostrës funksionit të plausibilitetit për të bërë mostrat ekzistuese sa më plausible. Me fjalë të tjera, ne duhet të përshtatim parametra të tillë, ku mostrat do të jenë më të plausibla. Në rastin tonë, parametri që përshtatet është probabiliteti i shlyerjes së kredisë
, i cili nga ana e tij varet nga koeficientët e panjohur
. Pra, na nevojitet të gjejmë një vektor peshash
, për të cilin plausibiliteti i mostrës do të jetë maksimal.
5. Ne e dimë se për maksimalizimin funksionit të plausibilitetit të mostrës mund të përdoret metoda e maksimumit të mundësisë. Dhe ne dimë të gjitha hilet për të punuar me këtë metodë.
Kjo duket si njĂ« lĂ«vizje shumĂ«dimensionale đ
Tani le të kujtojmë se në fillim të artikullit doja të nxirrnim dy lloje të funksioneve të humbjes Humbja Logjistike në përputhje me mënyrën se si shënohen klasat e objekteve. Kështu që është e zakonshme që në detyrat e klasifikimit me dy klasa, klasat shënohen si
dhe
ose
. Në përputhje me shenjat, rezultati do të jetë funksioni përkatës i humbjes.
Rasti 1. Klasifikimi i objekteve në
dhe 
Më parë, kur e përcaktonim kodin e mundësisë së mostrës, ku probabiliteti i shlyerjes së borxhit nga huamarrësi është llogaritur në bazë të faktorëve dhe koeficientëve të caktuar
, ne aplikuam formulën:

Në të vërtetë
â ky Ă«shtĂ« vlera e funksionit tĂ« pĂ«rgjigjes logjistike
me një vektor peshe të caktuar 
Atëherë asgjë nuk na pengon ta shënojmë funksionin e mundësisë së mostrës kështu:

Ndonjëherë ndodh që ndonjë analist i ri ka vështirësi të kuptojë menjëherë si funksionon ky funksion. Le të shqyrtojmë 4 shembuj të shkurtër që do t'i ndriçojnë të gjitha:
1. Nëse
(dmth. në përputhje me mostrat e trajnimit, objekti i përket klasës +1), dhe algoritmi ynë
përcakton probabilitetin e përfshirjes së objektit në klasën
të barabartë me 0.9, atëherë ky copë i mundësisë së mostrës do të llogaritet kështu:

2. Nëse
, dhe
, kështu që llogaritja do të jetë kështu:

3. Nëse
, dhe
, kështu që llogaritja do të jetë kështu:

4. Nëse
, dhe
, kështu që llogaritja do të jetë kështu:

E qartĂ« se funksioni i mundĂ«sisĂ« do tĂ« maksimizohet nĂ« rastet 1 dhe 3 ose nĂ« pĂ«rgjithĂ«si â kur vlerat e probabilitetit pĂ«rshkrimit tĂ« objektit nĂ« klasĂ« janĂ« tĂ« sakta
.
Duke pasur parasysh se në përcaktimin e probabilitetit të përfshirjes së objektit në klasë
na mungojnë vetëm koeficientët
, kështu që ne do t'i kërkojmë ata. Siç u tha më parë, kjo është një detyrë optimizimi, ku fillimisht na nevojitet të gjejmë derivatën e funksionit të mundësisë në vektorin e peshave
. Megjithatë, më parë ka kuptim ta thjeshtojmë detyrën: do të kërkojmë derivatën nga logaritmi funksionin e mundësisë.

Pse pas logaritmimit, në funksionin e gabimit logjistik, ne ndërruam shenjën nga
në
. ĂshtĂ« e thjeshtĂ«, pasi nĂ« detyrat e vlerĂ«simit tĂ« cilĂ«sisĂ« sĂ« modelit Ă«shtĂ« zakon qĂ« tĂ« minimizosh vlerĂ«n e funksionit, kĂ«shtu qĂ« ne e shumezuam anĂ«n e djathtĂ« tĂ« shprehjes me
dhe për rrjedhojë, tani minimizojmë funksionin.
Pra, tani, para jush Ă«shtĂ« nxjerrĂ« me vuajtje funksioni i humbjes â Humbja Logjistike pĂ«r mostrĂ«n e trajnuar me dy klasa:
dhe
.
Tani, për të gjetur koeficientët, mjafton të gjejmë derivatën funksionin e gabimit logjistik dhe më pas, duke përdorur metoda numerike optimizimi, si rënia e gradienteve ose rënia e gradienteve stohastike, të gjejmë koeficientët më optimal.
. Por, duke marrë parasysh, tashmë jo një volum të vogël të artikullit, propozohet të bëni diferencimin vetë ose, ndoshta, kjo do të jetë tema e artikullit të ardhshëm me më shumë matematikë pa shembuj kaq të detajuar.
Rasti 2. Klasifikimi i objekteve në
dhe 
Qasja këtu do të jetë e ngjashme me ato të klasave
dhe
, por vetĂ« rruga pĂ«r tĂ« nxjerrĂ« funksionin e humbjes Humbja Logjistike, do tĂ« jetĂ« mĂ« e ndĂ«rlikuar. TĂ« fillojmĂ«. Do tĂ« pĂ«rdorim operatorin "nĂ«seâŠ, atĂ«herĂ«âŠ". DomethĂ«nĂ«, nĂ«se
-i objekt i përket klasës
, atëherë për llogaritjen e mundësisë së mostrës përdorim probabilitetin
, nëse objekti i përket klasës
, atëherë në mundësi zëvendësojmë
. Kështu duket funksioni i mundësisë:

Për ta thjeshtuar, le të shohim si funksionon. Le të shqyrtojmë 4 raste:
1. Nëse
dhe
, atëherë në mundësinë e mostrës "do të shkojë" 
2. Nëse
dhe
, atëherë në mundësinë e mostrës "do të shkojë" 
3. Nëse
dhe
, atëherë në mundësinë e mostrës "do të shkojë" 
4. Nëse
dhe
, atëherë në mundësinë e mostrës "do të shkojë" 
E qartë se në rastin 1 dhe 3, kur probabilitetet ishin përcaktuar saktë nga algoritmi, funksioni i mundësisë do të maksimizohet, pra, pikërisht këtë donim të arrijmë. Megjithatë, kjo qasje është mjaft voluminoze dhe më tej do të shqyrtojmë një shkrim më kompakt. Por për fillim, do ta logaritmojmë funksionin e mundësisë me ndërrimin e shenjës, pasi tani do të minimizojmë atë.

Le të zëvendësojmë me
shprehjen
:

Le të thjeshtojmë pjesën e djathtë nën logaritm, duke përdorur metoda të thjeshta aritmetike dhe të marrim:

Dhe tani ka ardhur koha tĂ« heqim operatorin "nĂ«seâŠ, atĂ«herĂ«âŠ". VĂ«rejmĂ« se kur objekti
i përket klasës
, atëherë në shprehjen nën logaritm, në emëruesin,
ngrihet në fuqi
, nëse objekti i përket klasës
, kështu që $e$ ngrihet në fuqi
. Si rezultat, shprehja e fuqisĂ« mund tĂ« thjeshtohet â tĂ« dy rastet tĂ« bashkohen nĂ« njĂ«:
. Atëherë funksioni i gabimit logjistik do të marrë formën:

Sipas rregullave të logaritmimit, ne e kthejmë fraksionin përmbys dhe nxjerrim shenjën "
" (minus) përpara logaritmit, duke marrë:

Para jush është funksioni i humbjes logistic Loss, i cili përdoret në mostrat e trajnuara me objekte që i përkasin klasave:
dhe
.
ĂfarĂ« do, nĂ« kĂ«tĂ« moment po e lĂ« dhe mbyllim artikullin.
Puna e mĂ«parshme e autorit â "Kthimi i ekuacionit tĂ« regresionit linear nĂ« formĂ« matricore"
Materiale ndihmëse
1. Literatura
1) Analiza regresionare e aplikimeve / N. Draper, G. Smith â botimi i dytĂ« â M.: Financa dhe Statistika, 1986 (pĂ«rkthim nga anglishtja)
2) Teoria e probabilitetit dhe statistika matematike / V.E. GmurmĂĄn â botimi i nĂ«ntĂ« â M.: Shkolla e LartĂ«, 2003
3) Teoria e probabilitetit / N.I. Ăernova â Novosibirsk: Universitete i Shtetit tĂ« Novosibirskut, 2007
4) Biznes-analitika: nga tĂ« dhĂ«nat nĂ« njohuri / Paklin N. B., Oreshkov V. I. â botimi i dytĂ« â ShĂ«n Petersburg: Piter, 2013
5) Data Science Shkenca e tĂ« dhĂ«nave nga e para / Joel Gras â ShĂ«n Petersburg: BHK Petersburg, 2017
6) Statistika praktike pĂ«r specialistĂ«t e Data Science / P.Bruce, E.Bruce â ShĂ«n Petersburg: BHK Petersburg, 2018
2. Lectures, courses (video)
1)
2)
3)
4)
5)
3. Internet sources
1)
2)
3)
4)
5)
6)
7)
8)
Burimi: habr.com

në funksionin e përgjigjes logjistike 
funksione Humbja Logjistike
:

, atëherë vlera e funksionit
do të jetë pozitive nga
deri në
. Kështu ne mund të mendojmë se probabiliteti i shlyerjes së kredisë, është brenda
. Sa më e lartë të jetë vlera e funksionit, aq më e lartë është probabiliteti.
ose
, atëherë vlera e funksionit do të jetë negative nga
deri në
. Atëherë do të mendojmë se probabiliteti i shlyerjes së borxhit është brenda
dhe, sa më i lartë në modulin e vlerës së funksionit, aq më e lartë është siguria jonë.
do të jetë e barabartë me
dhe probabiliteti i shlyerjes së kredisë është i barabartë me
.


në funksionin e përgjigjes logjistike
dhe 
dhe 