Rizikos analizės metodai ir rodikliai

Mokslo vadovėlis: Aukštojo mokslo studijoms skirta medžiaga

Anotacija#

Šis skyrius pereina nuo rizikos koncepcijų (3 skyrius) prie konkrečių rizikos analizės metodų ir jų taikymo verslo planavime. Pristatoma keturių pagrindinių metodų sistema: jautrumo analizė atsako į klausimą „kas svarbu?“, scenarijų analizė — „kas galėtų nutikti?“, Monte Carlo simuliacija — „kokia tikimybė?“, o rizikos rodikliai — „kiek rizikuojame?“. Šie metodai nėra konkurentai — jie sudaro laipsniškai gilėjančią analitinę seką, kurioje kiekvienas lygmuo papildo ankstesnį ir atsako į vis sudėtingesnį klausimą.

Skyriaus centrinė tezė — rizikos analizė yra ne vienas universalus įrankis, o struktūruota sistema, kurios kiekvienas elementas turi aiškias galimybes ir ribas. Jautrumo analizė identifikuoja kritinius kintamuosius, tačiau ignoruoja jų tarpusavio sąveiką. Scenarijų analizė apčiuopia struktūrinius lūžius, tačiau priklauso nuo subjektyvių tikimybių. Monte Carlo simuliacija kiekybiškai integruoja neapibrėžtumą, tačiau propaguoja prielaidų klaidas su statistiniu tikslumu. Rizikos rodikliai (VaR, CVaR, sverto koeficientai) suteikia standartizuotas metrikas, tačiau gali sukurti klaidingo saugumo iliuziją. Efektyvi rizikos analizė reikalauja derinti visus keturis metodus, sąmoningai suvokiant kiekvieno ribas ir taikant proporcingumo principą — analizės sudėtingumas turi atitikti sprendimo mastą.

Teorinėje dalyje atsekama metodų raida nuo Hertzo (1964) pionierinio Monte Carlo taikymo investicijoms, per Kahno ir Wacko scenarijų revoliuciją iki Artzner ir kolegų (1999) koherentiškų rizikos matų teorijos, pakeitusios bankų reguliavimo paradigmą. Analitinėje dalyje pateikiamos trylika formulių (4.1–4.13) su skaitiniais pavyzdžiais, apimančios elastingumo koeficientą, scenarijų tikėtinę vertę ir dispersiją, Monte Carlo pasikliautinąjį intervalą, lūžio tašką, operacinio ir finansinio sverto laipsnius, parametrinį VaR, CVaR, sprendimų medžio tikėtinę piniginę vertę, Omega koeficientą, didžiausią nuosmukį (MDD) ir Sobol’ globalios jautrumo analizės indeksą. Skyrius tiesiogiai remiasi 3 skyriaus rizikos koncepcijomis ir baziniais matais (standartinis nuokrypis, beta, Sortino, variacijos koeficientas), juos gilindamas ir operacionalizuodamas. Taip pat rengiamas pagrindas 5 skyriui (diversifikacijos sprendimai), 11 skyriui (WACC jautrumo analizė) ir 14 skyriui (realieji opcionai ir sprendimų medžiai).

4.1 Teorinė dalis#

4.1.1 Istorinis ir konceptualus kontekstas#

Nuo intuicijos prie struktūruoto neapibrėžtumo valdymo#

Verslo sprendimuose rizika visada buvo vertinama — tačiau ilgą laiką tai buvo veikiau intuicijos, patirties ir nuojautos sritis nei struktūruoto analitinio proceso objektas. Dar 3 skyriuje aptartas Knighto (1921) skyrimas tarp rizikos (kur tikimybės žinomos arba apskaičiuojamos) ir tikro neapibrėžtumo (kur nei tikimybės, nei galimi scenarijai nėra aiškūs) ilgą laiką buvo filosofinė distinkcija, neturinti praktinių įrankių. Šio skyriaus metodai — jautrumo analizė, scenarijai, Monte Carlo, rizikos rodikliai — yra būtent tie įrankiai, kurie per septynių dešimtmečių evoliuciją leido Knighto skirtį paversti operacionalia analitine sistema.

Pirmieji žingsniai: inžinerija ir karas (1940–1960)#

Monte Carlo metodas — šiandien vienas galingiausių rizikos analizės įrankių — gimė ne finansuose, o fizikoje. 1946–1947 m. matematikas Stanisławas Ulamas, sveikdamas po ligos ir žaisdamas pasjansą, suvokė, kad sudėtingų sistemų elgseną galima tirti ne analitiškai sprendžiant lygtis, o generuojant daugybę atsitiktinių bandymų ir stebint rezultatų pasiskirstymą (Eckhardt, 1987). Kartu su Johnu von Neumannu jis formalizavo šią idėją ir pritaikė ją neutronų sklaidos skaičiavimams Manhattan projekte — JAV branduolinio ginklo kūrimo programoje. Pavadinimas „Monte Carlo“ buvo pasirinktas kaip nuoroda į garsųjį Monako kazino — žaidimų ir atsitiktinumo simbolį (Metropolis & Ulam, 1949).

Lygiagrečiai, karo ir pokario laikotarpiu, RAND korporacijoje Hermanas Kahnas pradėjo taikyti scenarijų metodą — struktūruotą mąstymą apie alternatyvias ateitis. Kahnas naudojo scenarijus ne prognozavimui, o branduolinio karo strateginių galimybių sisteminiam tyrinėjimui — jo knyga On Thermonuclear War (1960) sukėlė viešą polemiką, tačiau įtvirtino scenarijų planavimą kaip rimtą analitinę discipliną. Esminis Kahno įnašas buvo epistemologinis: scenarijai nėra prognozės — jie yra struktūruoti hipotetiniai pasakojimai, skirti ne atspėti ateitį, o išplėsti sprendimų priėmėjų mentalinį horizontą (Kahn & Wiener, 1967).

Finansų revoliucija: Hertz ir investicijų rizikos analizė (1964)#

Lūžis finansų srityje įvyko 1964 m., kai Davidas B. Hertzas Harvard Business Review žurnale paskelbė straipsnį „Risk Analysis in Capital Investment“ — pirmąjį sistemingą Monte Carlo metodo pritaikymą investiciniams sprendimams. Hertzas argumentavo, kad tradicinis investicijų vertinimas, pagrįstas vienu „labiausiai tikėtinu“ scenarijumi, yra fundamentaliai klaidingas: jis ignoruoja rezultatų sklaidą ir sukuria netikro tikrumo iliuziją. Vietoj to Hertzas pasiūlė kiekvienam DCF modelio kintamajam priskirti tikimybinį skirstinį ir generuoti tūkstančius simuliacijos iteracijų, gaunant ne vieną NPV skaičių, o visą NPV pasiskirstymą (Hertz, 1964).

Hertzo straipsnis buvo revoliucinis ne tik metodologiškai, bet ir epistemologiškai. Jis pirmasis aiškiai suformulavo mintį, kuri šiandien atrodo akivaizdi: vieno skaičiaus vertinimas yra ne tikslumas, o tikslumas be informacijos apie neapibrėžtumą. Ši idėja tiesiogiai siejasi su 2 skyriuje aptartu klaidingo tikslumo fenomenu — sudėtingas modelis, pateikiantis vieną konkretų rezultatą, gali būti mažiau informatyvus nei paprastesnis modelis, parodantis galimų rezultatų intervalą.

Tačiau Hertzo idėjos dešimtmečius liko daugiau teorinės nei praktinės — ne dėl konceptualių trūkumų, o dėl technologinių ribų. 1964 m. Monte Carlo simuliacija reikalavo mainframe kompiuterio; praktinis jos taikymas verslo planavime tapo įmanomas tik aštuntajame–devintajame dešimtmetyje, kai atsirado personaliniai kompiuteriai ir skaičiuoklių priedai, tokie kaip @RISK (1987).

Scenarijų revoliucija: Shell ir naftos krizė (1970-aisiais)#

Scenarijų metodas iš karinės strategijos į verslo pasaulį persikėlė per vieną iškalbingiausių atvejų valdymo istorijoje. Pierre’as Wackas, Royal Dutch Shell planavimo skyriaus analitikas, septintojo dešimtmečio pabaigoje pradėjo taikyti struktūruotus scenarijus naftos rinkos analizei. Tuo metu dominavo viena prielaida: naftos kaina liks stabili ir prognozuojama, nes OPEC šalys neturės motyvacijos drastiškai keisti tiekimo. Wackas sukūrė alternatyvų scenarijų, kuriame OPEC panaudoja naftą kaip geopolitinį ginklą — ir kai 1973 m. naftos embargas tapo realybe, Shell buvo vienintelė didelė naftos bendrovė, turėjusi parengtą veiksmų planą (Wack, 1985).

Wacko indėlis buvo ne tik praktinis, bet ir konceptualus. Jis suformulavo „mentalinių modelių laužymo“ (angl. breaking mental models) principą: scenarijų tikroji vertė yra ne prognozuoti ateitį, o priversti vadovus kvestionuoti savo giliai įsišaknijusias prielaidas apie tai, kaip veikia pasaulis. Tai tiesiogiai siejasi su 3 skyriuje aptartais kognityviniais šališkumais — ypač inkaravimo efektu ir pertekliniu pasitikėjimu.

Peteris Schwartzas, buvęs Shell strategas, vėliau demokratizavo scenarijų metodą, įkurdamas Global Business Network (GBN) ir parašydamas knygą The Art of the Long View (1991), kuri tapo praktinio scenarijų planavimo vadovu. Schoemaker (1995) suteikė metodui akademinį pagrindą, aprašydamas struktūruotą scenarijų kūrimo procesą.

Jautrumo analizė: nuo inžinerijos iki finansų#

Jautrumo analizė — paprasčiausias, bet neretai praktiškiausias rizikos analizės įrankis — neturi vieno „tėvo“ ar aiškaus gimimo momento. Ji natūraliai kilo iš inžinerinės ir mokslinės praktikos: bet kuris modelis reikalauja suprasti, kurie įvesties parametrai labiausiai veikia rezultatą. Finansuose jautrumo analizė tapo neatsiejama modeliavimo dalimi kartu su skaičiuoklių revoliucija (žr. 1 ir 2 skyrius): kai VisiCalc ir Excel leido interaktyviai keisti prielaidas, „kas, jeigu?“ (angl. what-if) klausimai tapo kasdiene analitiko praktika.

Akademinį pagrindą jautrumo analizei suteikė Hillier (1963), kuris pirmasis pasiūlė investicijų vertinime naudoti NPV tikimybinį skirstinį, ir Saltelli su kolegomis (2008), sukūrusiais globalios jautrumo analizės (angl. Global Sensitivity Analysis, GSA) metodologiją, kuri įvertina ne tik pavienių kintamųjų, bet ir jų sąveikų poveikį modelio rezultatui.

Rizikos matų revoliucija: nuo VaR iki koherentiškų matų (1994–1999)#

Paskutinis didelis metodologinis lūžis įvyko rizikos rodiklių srityje. 1994 m. JP Morgan viešai paskelbė RiskMetrics — standartizuotą VaR (angl. Value at Risk) skaičiavimo metodologiją, kuri, kaip aptarta 3 skyriuje, pirmą kartą pavertė sudėtingą rizikos modelį viešai prieinamu standartu. VaR greitai tapo bankų reguliavimo pagrindu (Basel II), tačiau 2008 m. finansų krizė atskleidė fundamentalų jo trūkumą: VaR parodo, kokio didžiausio nuostolio galima tikėtis „normaliomis“ sąlygomis, bet nieko nepasako apie tai, kas nutinka, kai tos sąlygos pažeidžiamos.

Matematinį atsakymą į šią kritiką pateikė Artzner, Delbaen, Eber ir Heath (1999) savo įtakingame straipsnyje „Coherent Measures of Risk“, kuriame suformulavo keturias aksiomas, kurias turi tenkinti „geras“ rizikos matas: monotoninumo, subadityvumo, teigiamo homogeniškumo ir transliacinio invariantiškumo. VaR pažeidžia subadityvumo aksiomą — tai reiškia, kad dviejų portfelių sujungimas teoriškai gali padidinti VaR, nors diversifikacija turėtų riziką mažinti. CVaR (angl. Conditional Value at Risk), dar vadinamas tikėtinu trūkumu (angl. Expected Shortfall, ES), šias aksiomas tenkina ir todėl Bazilio komitetas nuo Basel III / FRTB perėjo nuo VaR prie ES kaip pagrindinio reguliacinio rizikos mato.

Sintezė: keturių metodų sistema#

4.1 pav. Rizikos analizės metodų evoliucija ir tarpusavio ryšiai

graph LR
    subgraph "1940–1960"
        A["Ulam & von Neumann<br/><i>Monte Carlo metodas</i><br/>1946–1949"]
        B["Kahn (RAND)<br/><i>Scenarijų planavimas</i><br/>1960-aisiais"]
    end
    subgraph "1960–1990"
        C["Hertz<br/><i>MC investicijoms</i><br/>1964"]
        D["Wack (Shell)<br/><i>Verslo scenarijai</i><br/>1970-aisiais"]
        E["Hillier / Saltelli<br/><i>Jautrumo analizė</i><br/>1963–2008"]
        F["@RISK, Crystal Ball<br/><i>Skaičiuoklių priedai</i><br/>1987"]
    end
    subgraph "1990–dabar"
        G["JP Morgan<br/><i>RiskMetrics / VaR</i><br/>1994"]
        H["Artzner et al.<br/><i>CVaR / ES</i><br/>1999"]
        I["Integruota<br/>rizikos analizė"]
    end

    A --> C
    B --> D
    C --> F
    D --> I
    E --> I
    F --> I
    G --> H
    H --> I

    style I fill:#4a90d9,color:#fff

Ši schema parodo, kad šiandieninė integruota rizikos analizė nėra vienas metodas, o keturių nepriklausomai kilusių tradicijų konvergencija: fizikos (Monte Carlo), karinės strategijos (scenarijai), inžinerijos (jautrumas) ir finansų reguliavimo (rizikos matai). Kiekvienas metodas atsako į skirtingą klausimą, ir jų galia atsiskleidžia tik naudojant juos kartu — tai detaliau aptariama 4.1.2 poskyryje.

4.1.2 Pagrindinės teorijos ir modeliai#

Jautrumo analizė: kas labiausiai veikia rezultatą?#

Jautrumo analizė (angl. sensitivity analysis) yra paprasčiausias ir dažniausiai taikomas rizikos analizės metodas. Jos esmė — sistemingai keisti vieną ar kelis modelio įvesties kintamuosius ir stebėti, kaip tai paveikia rezultatą (paprastai — NPV, pelną ar pinigų srautą). Tai leidžia identifikuoti kritinius kintamuosius — tuos, kurių nedideli pokyčiai sukelia neproporcingai didelius rezultato svyravimus.

Vienkryptė jautrumo analizė (angl. one-way sensitivity analysis) keičia vieną kintamąjį vienu metu, laikydama visus kitus pastoviais (ceteris paribus). Pavyzdžiui, DCF modelyje galima atskirai tirti, kaip NPV kinta priklausomai nuo WACC (nuo 8 % iki 14 %), pajamų augimo tempo (nuo 2 % iki 8 %) ar EBIT maržos (nuo 10 % iki 20 %). Rezultatas paprastai vaizduojamas tornado diagrama — horizontaliu stulpeliu grafiku, kuriame kintamieji surikiuojami pagal jų poveikio dydį, nuo didžiausio iki mažiausio. Tornado diagrama yra itin efektyvus komunikacijos įrankis: ji akimirksniu parodo vadovams, „kurie trys veiksniai gali mus „nužudyti“ greičiausiai” — ir todėl plačiai naudojama investicinių komitetų prezentacijose ir due diligence ataskaitose.

Kitas vizualizacijos būdas — vorų diagrama (angl. spider diagram), kurioje kiekvieno kintamojo poveikis vaizduojamas kaip linija, einanti nuo bazinio scenarijaus. Statesnis linijos nuolydis reiškia didesnį jautrumą. Voratinklio diagrama ypač naudinga, kai norima lyginti ne tik poveikio dydį, bet ir jo kryptį (teigiamą ar neigiamą) bei tiesinumą — ar poveikis yra proporcingas, ar paspartėjantis.

Kiekybiškai jautrumas matuojamas elastingumo koeficientu — santykiu tarp procentinio rezultato pokyčio ir procentinio kintamojo pokyčio (formulė (4.1), žr. 4.2 poskyrį). Elastingumas leidžia palyginti skirtingų kintamųjų poveikį vienodais matais, nepriklausomai nuo jų absoliutinių reikšmių.

Daugiakryptė jautrumo analizė (angl. multi-way sensitivity analysis) keičia kelis kintamuosius vienu metu. Dažniausias variantas — dviejų kintamųjų lentelė (angl. two-way data table), pavyzdžiui, WACC × augimo tempas, generuojanti NPV matricą. Tai itin naudinga DCF vertinime, nes parodo, kaip du labiausiai neapibrėžti parametrai — diskonto norma ir ilgalaikis augimas — kartu veikia vertę. Tačiau jau trijų kintamųjų vienu metu varijavimas ($5^3 = 125$ kombinacijos) tampa sunkiai vizualizuojamas, o šešių ($5^6 = 15\ 625$) — nebeįmanomas interpretuoti be automatizuotų įrankių.

Šią vienkryptės ir dvikryptės analizės ribą sprendžia globali jautrumo analizė (angl. Global Sensitivity Analysis, GSA), kurią formalizavo Saltelli ir kolegos (2008). Skirtingai nuo lokalios analizės, kuri keičia kintamuosius po vieną aplinkui bazinį tašką, GSA vertina kintamųjų poveikį visame jų galimų reikšmių diapazone, įskaitant tarpusavio sąveikas (angl. interactions). Tai ypač svarbu, kai modelyje yra netiesiškumų ar koreliacijų, dėl kurių vienkryptė analizė gali klaidinti.

Sobol’ indeksai — populiariausias GSA matas — remiasi dispersijos dekompozicija. Pirmos eilės Sobol’ indeksas $S_i$ parodo, kokią bendros rezultato dispersijos dalį paaiškina vienas kintamasis $X_i$:

kur $\mathbb{E}[V \mid X_i]$ yra sąlyginė tikėtina modelio rezultato vertė, žinant $X_i$ reikšmę. Jei $S_i = 0{,}45$, tai reiškia, kad 45 % viso modelio neapibrėžtumo kyla iš šio vieno kintamojo. Totalinis Sobol’ indeksas $S_{T_i}$ papildomai apima visas kintamojo sąveikas su kitais kintamaisiais — skirtumas $S_{T_i} - S_i$ parodo sąveikų svarbą. Kai $\sum S_i \approx 1$, modelis yra beveik adityvus ir sąveikos nereikšmingos; kai $\sum S_i \ll 1$, sąveikos dominuoja ir vienkryptė analizė fundamentaliai nepakankama.

Naujesni tyrimai siūlo Shapley efektus (angl. Shapley effects) kaip alternatyvą Sobol’ indeksams, ypač kai kintamieji yra tarpusavyje koreliuoti — tokiu atveju klasikinė dispersijos dekompozicija tampa dviprasmiška, o Shapley efektai, paremti kooperatyvinės žaidimų teorijos logika, leidžia gauti stabilesnę ir aiškiau interpretuojamą įtakos metriką (Palar et al., 2023).

Praktikoje GSA reikalauja tūkstančių ar dešimčių tūkstančių modelio paleidimų ir todėl dažniausiai derinama su Monte Carlo simuliacija (žr. toliau). Šiuolaikiniai įrankiai (@RISK, Crystal Ball, Python SALib biblioteka) automatizuoja šį procesą, tačiau analitikui svarbu suprasti, kad GSA rezultatai priklauso nuo pasirinktų kintamųjų diapazonų ir skirstinių — t. y. nuo to, kiek plačiai „leidžiame“ kintamiesiems svyruoti.

Jautrumo analizės vertė verslo planavime yra dvejopa. Analitiškai ji identifikuoja, kurioms prielaidoms verta skirti daugiausia dėmesio ir išteklių — jei 70 % NPV neapibrėžtumo kyla iš pajamų augimo tempo, būtent šiai prielaidai turi būti skirta giliausia rinkos analizė. Komunikaciškai tornado diagrama yra vienas efektyviausių būdų perduoti rizikos struktūrą ne finansų specialistams — ji leidžia valdybos nariams ar investuotojams akimirksniu suvokti, kas modelyje yra svarbiausia.

Scenarijų analizė: kas galėtų nutikti?#

Scenarijų analizė (angl. scenario analysis) yra kokybiškai kitokio pobūdžio rizikos analizės metodas nei jautrumo analizė. Jei jautrumo analizė klausia „kiek pasikeistų rezultatas, jei šis skaičius būtų kitoks?“, tai scenarijų analizė klausia „kaip atrodytų visas verslo kontekstas, jei pasaulis pasuktų kita linkme?“. Tai ne parametrų varijavimas, o alternatyvių ateičių konstravimas — struktūruotas mąstymas apie kokybiškai skirtingas galimybes.

Epistemologinė distinkcija. Scenarijai nėra prognozės. Prognozė siekia vieno „teisingo“ atsakymo — labiausiai tikėtino rezultato. Scenarijai sąmoningai atsisako šio tikslo ir pripažįsta, kad ateitis yra fundamentaliai daugiavariantė, ypač Knighto neapibrėžtumo zonoje (žr. 3 skyrių), kur nei tikimybės, nei galimi scenarijai nėra stabilūs. Kaip suformulavo Wackas (1985): scenarijų tikslas — ne atspėti ateitį, o priversti sprendimų priėmėjus paruošti savo protus kelioms kokybiškai skirtingoms ateičių versijoms.

Klasikinis trijų scenarijų modelis — optimistinis, bazinis, pesimistinis — yra plačiausiai paplitęs verslo planuose, tačiau turi rimtų metodologinių trūkumų. Pirma, jis sukuria inkaravimo efektą: bazinis scenarijus psichologiškai tampa „tikruoju“, o kiti du — tik formaliais priedais. Antra, trys scenarijai dažnai konstruojami paprasčiausiai keičiant parametrų dydį (pajamos +10 % / 0 % / –10 %), iš esmės atliekant jautrumo analizę su kitais pavadinimais, o ne kuriant kokybiškai skirtingus pasaulio vaizdus.

Šiuolaikinis struktūruotas požiūris naudoja 2×2 neapibrėžtumo matricą: identifikuojamos dvi svarbiausios ir mažiausiai prognozuojamos jėgos (pvz., reguliacinė aplinka ir technologinė disrupcija), sukuriamos keturios jų kombinacijos, ir kiekvienai kuriamas nuoseklus, vidujai logiškas pasakojimas — scenarijus. Keturi scenarijai nėra „blogas–vidutinis–geras“ spektras; jie yra keturi kokybiškai skirtingi pasauliai, iš kurių kiekvienas reikalauja kitokios strategijos.

Scenarijų kiekybinė integracija. Nors scenarijai yra pirmiausia naratyvai, jie gali būti kiekybiškai integruojami priskiriant kiekvienam scenarijui subjektyvią tikimybę ir skaičiuojant tikėtinę vertę $E(V) = \sum p_i V_i$ bei dispersiją tarp scenarijų (formulės (4.2) ir (4.3), žr. 4.2 poskyrį). Tačiau būtina suvokti esminę epistemologinę problemą: šios tikimybės yra subjektyvios ekspertinio vertinimo išraiškos, o ne objektyvios dažnių statistikos — todėl jas reikia traktuoti kaip struktūruotus spėjimus, o ne kaip tikslius duomenis.

Scenarijai verslo planavime. Praktikoje scenarijai ypač vertingi ilgalaikiam strateginiam planavimui, kur kiekybiniai modeliai yra silpniausi. Lietuvos ir Baltijos kontekste scenarijų analizė yra itin aktuali dėl mažos atviros ekonomikos pažeidžiamumo išoriniams šokams: energijos kainų svyravimai (po Ignalinos AE uždarymo, SGD terminalo, NordBalt ir LitPol Link atsiradimo), ES reguliaciniai pokyčiai (Žaliasis kursas, CBAM, CSRD), geopolitinė rizika (Rusijos–Ukrainos karas), eksporto rinkų dinamika. Šiose srityse vieno „bazinio“ scenarijaus naudojimas yra ne supaprastinimas, o struktūrinis aplaidumas.

Monte Carlo simuliacija: kokia tikimybė?#

Monte Carlo simuliacija (angl. Monte Carlo simulation) yra galingiausias kiekybinis rizikos analizės metodas, leidžiantis iš deterministinio modelio sukurti tikimybinį — t. y. gauti ne vieną rezultatą, o visą galimų rezultatų pasiskirstymą. Jei jautrumo analizė keičia kintamuosius po vieną, o scenarijų analizė konstruoja kelis kokybiškai skirtingus pasaulius, tai Monte Carlo vienu metu keičia visus kintamuosius pagal jiems priskirtus tikimybinius skirstinius ir kartoja šį procesą tūkstančius ar dešimtis tūkstančių kartų.

Metodologijos esmė. Monte Carlo simuliacijos procesas apima keturis žingsnius:

1. Modelio identifikavimas. Nustatomas deterministinis modelis (pvz., DCF) su aiškiai apibrėžtais įvesties kintamaisiais ir rezultato formule.
2. Skirstinių priskyrimas. Kiekvienam neapibrėžtam kintamajam priskiriamas tikimybinis skirstinys, atspindintis jo neapibrėžtumo mastą ir formą:
   • Normalusis — kai kintamasis svyruoja simetriškai apie vidurkį (pvz., makroekonominiai rodikliai).
   • Trikampis — kai žinomi minimumo, maksimumo ir labiausiai tikėtino scenarijaus taškai, bet duomenų mažai (pvz., ekspertiniai vertinimai).
   • Lognormalusis — kai kintamasis negali būti neigiamas ir turi dešinę asimetriją (pvz., akcijų kainos, žaliavų kainos).
   • Beta / PERT — kai kintamasis svyruoja apibrėžtame intervale su žinoma moda (pvz., projekto trukmė). PERT skirstinys (angl. Program Evaluation and Review Technique) yra ypač patogus verslo planavime, nes reikalauja tik trijų ekspertinių įverčių: minimumo $a$, modos $c$ (labiausiai tikėtina reikšmė) ir maksimumo $b$. Iš jų apskaičiuojamas vidurkis $\mu = \frac{a + 4c + b}{6}$ ir standartinis nuokrypis $\sigma \approx \frac{b - a}{6}$. Pasvėrimas koeficientu 4 suteikia modai didesnį svorį nei kraštinėms reikšmėms, todėl PERT skirstinys yra mažiau jautrus ekstremumams nei trikampis skirstinys ir geriau atspindi ekspertų lūkesčius. Praktikoje PERT plačiai naudojamas projektų trukmės, kaštų ir pardavimų apimčių modeliavime, kai istorinių duomenų trūksta, bet ekspertai gali pateikti tris orientyrus (Vose, 2008).
3. Simuliacija. Atsitiktinių skaičių generatorius (paprastai Mersenne Twister algoritmas) kiekvienoje iteracijoje parenka kintamųjų reikšmes pagal priskirtus skirstinius, ir modelis apskaičiuoja rezultatą. Po $N$ iteracijų (paprastai 10 000–100 000) gaunamas rezultato pasiskirstymas.
4. Analizė. Rezultatų histograma leidžia apskaičiuoti tikėtinę vertę, standartinį nuokrypį, pasikliautinuosius intervalus ir tikimybes (pvz., $P(NPV < 0)$ — tikimybę, kad projektas bus nuostolingas).

Konvergencija ir iteracijų skaičius. Monte Carlo standartinė klaida mažėja proporcingai $1/\sqrt{N}$ — tai reiškia, kad norint sumažinti klaidą 10 kartų, reikia 100 kartų daugiau iteracijų. Praktikoje 10 000 iteracijų paprastai pakanka vidurkio ir standartinio nuokrypio stabilumui, tačiau uodegos tikimybėms (pvz., $P(NPV < 0)$, kai ta tikimybė maža) gali prireikti 50 000–100 000 iteracijų. Pasikliautinojo intervalo formulė pateikiama kaip (4.4) formulė 4.2 poskyryje.

Koreliacijų modeliavimas. Viena dažniausių Monte Carlo klaidų praktikoje — kintamųjų nepriklausomumo prielaida. Realybėje daugelis kintamųjų yra koreliuoti: pajamų augimas ir EBIT marža dažnai turi neigiamą koreliaciją (greitas augimas mažina maržą dėl investicijų poreikio), o palūkanų normos ir nekilnojamojo turto kainos — taip pat neigiamą. Ignoruojant koreliacijas, Monte Carlo simuliacija gali generuoti nerealias kintamųjų kombinacijas ir iškraipyti rezultatų pasiskirstymą. Techninis sprendimas — Cholesky dekompozicija (angl. Cholesky decomposition), matematinis metodas, kuris transformuoja nepriklausomus atsitiktinius skaičius į koreliuotus pagal nurodytą koreliacijos matricą (Vose, 2008). Alternatyva — Iman-Conover rango koreliacijų metodas, naudojamas @RISK ir panašiuose įrankiuose.

Latin Hypercube Sampling (LHS). Paprastas atsitiktinis atrankos metodas gali „praleisti“ tam tikras skirstinio sritis, ypač uodegas. LHS — patobulintas metodas, kuris padalina kiekvieno kintamojo skirstinį į $N$ lygių intervalų ir užtikrina, kad kiekvienoje iteracijoje būtų paimta po vieną reikšmę iš kiekvieno intervalo (McKay, Beckman, & Conover, 1979). Tai žymiai pagerina konvergenciją ir yra standartinis pasirinkimas šiuolaikiniuose simuliacijos įrankiuose.

Monte Carlo ir scenarijai: papildantys, ne konkuruojantys metodai. Monte Carlo ir scenarijų analizė atlieka fundamentaliai skirtingas funkcijas. Scenarijai kuria diskrečias, kokybiškai skirtingas ateitis su naratyvine logika — jie geriausiai tinka struktūriniams lūžiams, kur praeities statistika nebetinka. Monte Carlo kuria tolydų tikimybinį pasiskirstymą pagal parametrinius skirstinius — jis geriausiai tinka parametriniam neapibrėžtumui kiekybiškai vertinti. Efektyvi praktika — derinti abu: naudoti scenarijus struktūriniams lūžiams identifikuoti, o kiekvieno scenarijaus viduje taikyti Monte Carlo parametriniam neapibrėžtumui kiekybiškai vertinti.

Rizikos rodikliai: kiek rizikuojame?#

3 skyriaus 3.2 poskyryje pristatyti baziniai rizikos matai — standartinis nuokrypis ($\sigma$, formulė (3.1)), beta koeficientas ($\beta$, formulė (3.2)), Sortino koeficientas (formulė (3.5)), variacijos koeficientas (CV, formulė (3.6)) ir šalies rizikos premija (formulė (3.7)). Šie matai vertina bendrą rizikos lygį, tačiau neatsako į specifiškesnius klausimus: koks didžiausias galimas nuostolis? Kokia rizika uodegose? Kaip svertas didina jautrumą? Šiame poskyryje pristatomi specializuoti rodikliai, papildantys 3 skyriaus bazinį rinkinį.

VaR (rizikos vertė). VaR (angl. Value at Risk) atsako į konkretų klausimą: „kokio didžiausio nuostolio galima tikėtis per nurodytą laikotarpį su nurodyta tikimybe?“. Pavyzdžiui, „dienos VaR 99 % pasikliovimo lygiu yra €150 000“ reiškia: 99 dienų iš 100 nuostolis neviršys €150 000. VaR skaičiuojamas trimis metodais:

• Parametrinis (variacijos-kovariacijos): daroma prielaida, kad grąžos pasiskirsčiusios pagal normalųjį skirstinį, ir VaR apskaičiuojamas tiesiai iš vidurkio ir standartinio nuokrypio (formulė (4.8), žr. 4.2 poskyrį).
• Istorinė simuliacija: naudojami faktiniai istoriniai duomenys, nerenkant jokios skirstinio prielaidos — VaR yra tiesiog atitinkamas procentilis istorinių nuostolių eilutėje.
• Monte Carlo VaR: generuojamas nuostolių pasiskirstymas per Monte Carlo simuliaciją ir iš jo imamas atitinkamas procentilis.

VaR tapo bankų reguliavimo pagrindu per Basel II (žr. 3 skyrių), tačiau turi fundamentalų trūkumą: jis yra nekoherentiškas matas — pažeidžia subadityvumo aksiomą (Artzner et al., 1999). Tai reiškia, kad dviejų portfelių VaR suma gali būti mažesnė nei sujungto portfelio VaR, sukuriant paradoksą, kur diversifikacija teoriškai „padidina“ riziką. Be to, VaR nieko nepasako apie nuostolių dydį, kai VaR riba viršijama — jis parodo ribą, bet ne tai, kas už jos.

CVaR / Expected Shortfall. CVaR (angl. Conditional Value at Risk), dar vadinamas tikėtinu trūkumu (angl. Expected Shortfall, ES), atsako būtent į šį klausimą: „jei nuostolis viršija VaR ribą, koks yra vidutinis nuostolio dydis?“ Matematiškai tai yra sąlyginė tikėtina vertė $E[L \mid L > VaR_\alpha]$ (formulė (4.9)). CVaR tenkina visas keturias Artzner ir kolegų koherentiškumo aksiomas, todėl Bazilio komitetas nuo FRTB (angl. Fundamental Review of the Trading Book) perėjo prie CVaR kaip pagrindinio reguliacinio rizikos mato.

Lūžio taškas (break-even). Paprasčiausias, bet vienas naudingiausių rizikos rodiklių verslo planuose. Lūžio taškas parodo, kiek vienetų ar kokia apimtis reikalinga, kad pajamos padengtų visas sąnaudas (formulė (4.5)). Jo atvirkštinis rodiklis — saugumo riba (angl. margin of safety) — parodo, kiek dabartiniai pardavimai viršija lūžio tašką, t. y. kokį pardavimų sumažėjimą verslas gali atlaikyti nepradėdamas nuostoliauti.

Operacinio ir finansinio sverto laipsniai. Svertas (angl. leverage) yra rizikos „stiprintuvas“ — jis padidina ir pelną augimo metu, ir nuostolius nuosmukio metu. Operacinio sverto laipsnis DOL (angl. Degree of Operating Leverage) parodo, kiek procentų pakinta EBIT, kai pardavimai pakinta 1 % (formulė (4.6)). Finansinio sverto laipsnis DFL (angl. Degree of Financial Leverage) parodo analogišką grynojo pelno jautrumą EBIT pokyčiams (formulė (4.7)). Bendras sverto laipsnis DTL = DOL × DFL parodo kumuliacinį efektą: jei DOL = 3 ir DFL = 2, tai 10 % pardavimų sumažėjimas sukelia 60 % grynojo pelno sumažėjimą. Tai ypač svarbu vertinant kapitalo imlius ar stipriai finansuotus skolintomis lėšomis verslus.

Stresiniai testai. Stresinis testas (angl. stress test) yra scenarijų analizės forma, kurioje analizuojamos ekstremalios, bet įmanomos situacijos. Istoriniai stresiniai testai naudoja praeities krizių parametrus (pvz., 2008 m. finansų krizės metu akcijų rinkos smuko ~50 %, kredito maržos išsiplėtė 300–500 bazinių punktų). Hipotetiniai stresiniai testai kuria scenarijus, kurių praeityje nebuvo (pvz., kibernetinė ataka, paralyžiuojanti visą tiekimo grandinę). Atvirkštinis stresinis testas (angl. reverse stress test) veikia priešinga kryptimi: vietoj klausimo „kas nutiktų, jei…?“ jis klausia „kokios sąlygos priverstų mus tapti nemokiais?“ — ir nuo šio nesėkmės taško analizuoja atgal, kokių įvykių ar jų kombinacijų reikėtų, kad toks rezultatas taptų realybe. Metodologiškai tai yra trijų žingsnių procesas: (1) apibrėžiamas kritinis slenkstis (pvz., nuosavas kapitalas < 0, likvidumo koeficientas < 1,0, skolos aptarnavimo rodiklis < 1,2); (2) identifikuojamos kintamųjų kombinacijos, galinčios privesti prie šio slenksčio; (3) vertinama, ar tokios kombinacijos yra realios, remiantis istoriniais precedentais ir ekspertiniais vertinimais. Atvirkštinis stresinis testas yra itin vertingas, nes priverčia organizaciją tiesiogiai susidurti su savo pažeidžiamumais, o ne juos ignoruoti — jis demaskuoja „akląsias zonas“, kurių tradicinis stresinis testas gali neaptikti, nes pastarasis testuoja tik iš anksto numatomus scenarijus. EBA (Europos bankininkystės institucija) ir PRA (Jungtinės Karalystės reguliuotojas) reikalauja atvirkštinio stresinio testavimo iš finansų institucijų; verslo planavime ši logika ne mažiau vertinga, nors rečiau taikoma.

Omega koeficientas ir Maximum Drawdown. Omega koeficientas (Keating & Shadwick, 2002) yra alternatyvus rizikos-grąžos matas, kuris, skirtingai nuo Sharpe’o ar Sortino koeficientų, apima visą grąžų pasiskirstymą — ne tik pirmuosius du momentus (vidurkį ir dispersiją), bet ir asimetriją bei ekscesą. Maximum Drawdown (MDD) — didžiausias nuosmukis nuo piko iki dugno — yra intuityvus matas, ypač svarbus investuotojams: –30 % nuosmukis reikalauja +42,8 % atsigavimo vien tam, kad grįžtų į pradinį lygį. Calmar koeficientas (metinė grąža / |MDD|) leidžia palyginti investicijas pagal rizikos-grąžos santykį atsižvelgiant į didžiausią nuosmukį.

Integruota rizikos analizė: nuo metodų prie sistemos#

Aukščiau aprašyti keturi metodai nėra konkurentai — jie sudaro laipsniškai gilėjančią sistemą, kurioje kiekvienas lygmuo papildo ankstesnį:

4.2 pav. Integruotos rizikos analizės darbo srautas

graph TD
    A["1. Jautrumo analizė<br/><i>Tornado diagrama</i><br/>Klausia: kas svarbu?"] --> B["2. Scenarijų analizė<br/><i>2×2 matrica / 3 scenarijai</i><br/>Klausia: kas galėtų nutikti?"]
    B --> C["3. Monte Carlo<br/><i>Tikimybinis skirstinys</i><br/>Klausia: kokia tikimybė?"]
    C --> D["4. Rizikos rodikliai<br/><i>VaR, CVaR, DOL/DFL</i><br/>Klausia: kiek rizikuojame?"]
    D --> E["5. Sprendimas<br/><i>Priimti / Koreguoti / Atmesti</i>"]
    E -.->|Iteracija| A

    style A fill:#50b5a9,color:#fff
    style B fill:#f5a623,color:#fff
    style C fill:#4a90d9,color:#fff
    style D fill:#e8705a,color:#fff
    style E fill:#7b8a8e,color:#fff

Rizikos registras (angl. risk register) yra centrinis integruotos rizikos analizės dokumentas. Jis apima kiekvienos identifikuotos rizikos aprašymą, kategoriją (strateginė, operacinė, finansinė, reguliacinė), tikimybės ir poveikio vertinimą, rizikos prioriteto skaičių $RPN = P \times I$ (kur $P$ — tikimybė balais, $I$ — poveikis balais), rizikos valdymo strategiją ir atsakingą asmenį. Rizikos registras nėra statinis dokumentas — jis reguliariai atnaujinamas per visą projekto ar verslo plano gyvavimo ciklą.

Rizikos žemėlapis (angl. risk heat map) yra vizuali rizikos registro matrica, kurioje horizontalioji ašis vaizduoja tikimybę, o vertikalioji — poveikį. Rizikos spalvinamos pagal zonas: žalia (mažas RPN — priimtina rizika), geltona (vidutinis — stebėti ir valdyti), oranžinė (reikšminga — būtinas valdymo planas), raudona (kritinė — automatinis eskalacijos signalas). Rizikos žemėlapis yra efektyvus komunikacijos įrankis valdybai ir investuotojams, tačiau turi būti suvokiamas kaip diskusijos pradžios taškas, o ne galutinis atsakymas — tikimybių ir poveikio vertinimai yra subjektyvūs ir gali būti veikiami tų pačių kognityvinių šališkumų, kurie aptarti 3 skyriuje.

Sprendimų medžiai (angl. decision trees) yra metodas, leidžiantis modeliuoti nuoseklius sprendimus su neapibrėžtumu kiekviename etape. Sprendimų medis susideda iš sprendimų mazgų (kur priimamas sprendimas), tikimybės mazgų (kur atsitiktinis įvykis nulemia rezultatą) ir galinių mazgų (rezultatai). Kiekvieno kelio tikėtina piniginė vertė EMV (angl. Expected Monetary Value) apskaičiuojama kaip svertinis tikimybių ir rezultatų vidurkis (formulė (4.10), žr. 4.2 poskyrį), ir sprendimas priimamas atgaline indukcija — nuo galinių mazgų „atgal“ iki pirmojo sprendimo.

Sprendimų medžiai ypač svarbūs, kai investicijos yra stadijinės (angl. staged investment) — pavyzdžiui, pirmoji stadija yra pilotinis projektas, o antroji — pilno masto plėtra. Tai tiesiogiai siejasi su realių opcionų logika (plačiau — 14 skyriuje): stadijinis investavimas sukuria atsitraukimo opcioną (galimybę sustabdyti projektą, jei pilotas nepasiteisina) ir plėtros opcioną (galimybę padidinti investiciją, jei pilotas sėkmingas). Standartinis statinis NPV šių opcionų vertės neapima — todėl sprendimų medžiai ir realių opcionų mąstymas gali fundamentaliai pakeisti investicinio sprendimo logiką.

Šią laipsniškai gilėjančią seką galima vizualizuoti kaip vientisą analitinę grandinę:

4.2 pav. Integruota rizikos analizės seka

graph LR
    A["1. Jautrumo analizė<br/>Kas svarbu?<br/>(tornado diagrama)"] --> B["2. Scenarijų analizė<br/>Kas galėtų nutikti?<br/>(2×2 matrica)"]
    B --> C["3. Monte Carlo<br/>Kokia tikimybė?<br/>(tikimybinis skirstinys)"]
    C --> D["4. Rizikos rodikliai<br/>Kiek rizikuojame?<br/>(VaR, CVaR, DOL/DFL)"]
    D --> E["5. Sprendimas<br/>Ką darome?<br/>(priimti / atidėti / stadijinis)"]

    style A fill:#50b5a9,color:#fff
    style B fill:#f5a623,color:#fff
    style C fill:#4a90d9,color:#fff
    style D fill:#e8705a,color:#fff
    style E fill:#7b8a8e,color:#fff

Kiekvienas lygmuo prideda informacijos, bet ir sudėtingumo — todėl ne kiekvienam sprendimui reikalingi visi penki žingsniai.

Proporcingumo principas. Integruota rizikos analizė yra galinga, tačiau gali tapti „analizės paralyžiaus“ (angl. analysis paralysis) šaltiniu. Chapman ir Ward (2003) perspėja: rizikos analizės sudėtingumas turi būti proporcingas sprendimo mastui ir organizacijos gebėjimams. Mažai įmonei, vertinančiai €50 000 investiciją, pakanka jautrumo analizės ir kelių scenarijų. Dideliam infrastruktūros projektui (€50 mln.) reikalinga pilna Monte Carlo simuliacija su stresiniais testais ir rizikos registru. Taikyti pilną arsenalą kiekvienam sprendimui yra ne kruopštumas, o išteklių švaistymas. Pareto principas galioja ir čia: paprastai 80 % rizikos vertės kyla iš 20 % rizikos veiksnių — tornado diagrama juos identifikuoja, ir būtent jiems turi būti skiriamas didžiausias analitinis dėmesys.

4.1.3 Kritinė analizė ir diskusija#

Klaidingo tikslumo problema#

Kiekvienas šiame skyriuje pristatytas metodas turi bendrą pagundą — sukurti tikslumo iliuziją. Jautrumo analizė pateikia konkrečius skaičius lentelėje, ir skaitytojas linkęs juos traktuoti kaip „atsakymus“, nors tai tik deterministiniai taškai be tikimybių. Monte Carlo simuliacija generuoja pasiskirstymus su statistinėmis charakteristikomis (vidurkis, standartinis nuokrypis, procentiliai), kurie atrodo moksliškai griežti — tačiau jie yra tik tiek patikimi, kiek patikimos jiems priskirtos prielaidos. Kaip aptarta 2 skyriuje, sudėtingesnis modelis nebūtinai yra tikslesnis — jis tiesiog sukuria stipresnę „moksliškumo“ regimybę.

Ši problema ypač aštri Monte Carlo kontekste. Simuliacija gali generuoti NPV pasiskirstymą su trimis skaitmenimis po kablelio, tačiau jei pajamų augimo tempas priskirtas kaip „normalusis skirstinys su vidurkiu 5 % ir standartiniu nuokrypiu 3 %“ remiantis vieno analitiko spėjimu, o ne empiriniais duomenimis, tai visas pasiskirstymas yra ne kas kita, kaip statistiškai preciziškas to spėjimo atvaizdas. „Garbage in — garbage out“ principas Monte Carlo kontekste transformuojasi į subtilesnę formą: „garbage in — beautifully distributed garbage out“.

VaR: nekoherentiškumo paradoksas ir reguliacinė pamoka#

VaR kritika iliustruoja platesnę problemą — kaip plačiai priimtas standartas gali sukurti sisteminę riziką. Artzner ir kolegų (1999) parodyta subadityvumo aksiomų pažaida reiškia ne tik matematinę eleganciją — ji turi konkrečias praktines pasekmes. Jei bankas valdo du portfelius, kurių kiekvieno VaR yra priimtinas, jų sujungimas gali formaliai „padidinti“ VaR — tai reiškia, kad diversifikacija, kuri turėtų riziką mažinti, reguliaciniame rėmelyje gali ją „padidinti“. Šis paradoksas buvo viena priežasčių, dėl kurių Bazilio komitetas perėjo nuo VaR prie CVaR/ES.

Tačiau gilesnė pamoka yra ne techninė, o epistemologinė. VaR, net jei jis būtų matematiškai tobulas, turi struktūrinę ribą: jis nieko nepasako apie tai, kas vyksta už nurodytos ribos. 99 % VaR parodo „blogiausią normalią dieną“ — bet 2008 m. krizė parodė, kad būtent tas 1 % atvejų, apie kuriuos VaR tyli, nulemia ar institucija išgyvena, ar žlunga. Tai Talebo (2007) „juodosios gulbės“ kritikos esmė: modeliai, sukalibruoti pagal „gero oro“ duomenis, sistemingai nuvertina ekstremalių įvykių riziką.

Storų uodegų problema ir normalumo prielaida#

Normaliojo skirstinio prielaida — nuo kurios priklauso parametrinis VaR, daugelis Monte Carlo konfigūracijų ir standartinė portfelio teorija — yra viena labiausiai kritikuojamų prielaidų šiuolaikiniuose finansuose. Kaip aptarta 3 skyriuje, Mandelbrotas dar 1963 m. parodė, kad finansinių grąžų skirstiniai turi „storesnes uodegas“ (angl. fat tails) nei numato normalusis skirstinys — tai reiškia, kad ekstremalūs įvykiai (krizės, krachai) nutinka žymiai dažniau, nei prognozuoja standartiniai modeliai.

Monte Carlo kontekste ši problema pasireiškia dvejopai. Pirma, jei kintamiesiems priskiriami normalieji skirstiniai, simuliacija sistemingai nuvertins ekstremalių rezultatų tikimybę. Antra, krizinėmis sąlygomis koreliacijos tarp aktyvų ir kintamųjų staigiai padidėja (vadinamasis correlation breakdown efektas) — modelis, sukalibruotas pagal „ramaus laiko“ koreliacijas, gali generuoti per daug optimistišką įvairovę kriziniais periodais. Tai nereiškia, kad Monte Carlo yra nenaudingas — tai reiškia, kad jo rezultatus reikia interpretuoti sąmoningai, suprantant, kokių prielaidų ribose jie galioja.

Scenarijų subjektyvumo problema#

Scenarijų analizė pretenduoja valdyti struktūrinį neapibrėžtumą — tačiau pati yra subjektyvumo produktas. Kelios konkrečios problemos:

Tikimybių priskyrimo subjektyvumas. Kai scenarijams priskiriamos tikimybės (pvz., optimistinis — 25 %, bazinis — 50 %, pesimistinis — 25 %), šios tikimybės yra ekspertiniai spėjimai, o ne statistiniai dažniai. Skirtingi ekspertai tiems patiems scenarijams priskirs skirtingas tikimybes, ir nėra objektyvaus būdo nustatyti, kuris teisesnis. Formulė $E(V) = \sum p_i V_i$ sukuria matematinio tikslumo regimybę, tačiau yra tik tiek tikra, kiek tikros yra įvesties tikimybės.

Siauro rėmo pavojus (angl. framing bias). Sukurti trys ar keturi scenarijai neišvengiamai atspindi kūrėjų mentalinį horizontą. Scenarijai, kurie atrodo „per ekstremalūs“ arba „per keisti“, dažnai atmetami kaip nerealūs — tačiau būtent jie gali apimti tikrąsias „juodąsias gulbes“. 2020 m. COVID-19 pandemija buvo tokio pobūdžio įvykis: daugelis įmonių turėjo scenarijus ekonominiam nuosmukiui, bet labai mažai kas turėjo scenarijų globaliam tiekimo grandinių paralyžiui dėl epidemiologinio šoko.

Inkaravimas baziniame scenarijuje. Net kai sukuriami keli scenarijai, psichologiškai bazinis scenarijus tampa „tikruoju“ — sprendimų priėmėjai linkę vertinti optimistinį ir pesimistinį scenarijus kaip „mažai tikėtinus priedus“, o ne kaip lygiavertes galimybes. Tai tiesiogiai siejasi su 3 skyriuje aptartu inkaravimo šališkumu (angl. anchoring bias) ir yra viena priežasčių, kodėl 2×2 matricos požiūris (keturi kokybiškai skirtingi scenarijai be „bazinio“) yra metodologiškai pranašesnis.

Elgsenos šališkumai rizikos analizės procese#

3 skyriuje pristatyti septyni kognityviniai šališkumai pasireiškia ne tik rizikos suvokime, bet ir pačiame rizikos analizės procese:

• Perteklinis pasitikėjimas (angl. overconfidence) skatina analitikus priskirti Monte Carlo kintamiesiems per siaurus skirstinius — „aš tikrai žinau, kad augimas bus tarp 3 % ir 7 %“ — nors empirinė patirtis rodo, kad ekspertai sistemingai per siaurai vertina galimų rezultatų intervalą.
• Prieinamumo euristika (angl. availability heuristic) lemia, kad rizikos registre dominuoja neseniai patirtos arba žiniasklaidoje matomos rizikos, o „nematomos“ — bet potencialiai reikšmingesnės — rizikos praleidžiamos.
• Grupinis mąstymas (angl. groupthink) rizikos komitetuose gali slopinti disidentinių scenarijų kūrimą — jei dominuojanti nuomonė yra optimistinė, pesimistinis scenarijus gali būti „sušvelnintas“ iki socialiai priimtino lygio.
• Nuostolių eskalavimas (angl. escalation of commitment) reiškia, kad jau priimtam projektui rizikos analizė gali būti taikoma ne kaip objektyvus vertinimas, o kaip patvirtinimo paieška (angl. confirmation bias) — ieškant argumentų tęsti, o ne sustoti.

Šių šališkumų valdymo strategijos — pre-mortem analizė, red team, reference class forecasting — aptartos 3 skyriuje. Šio skyriaus kontekste svarbu pabrėžti, kad struktūruoti metodai (jautrumo analizė, scenarijai, Monte Carlo) yra patys savyje šališkumų mažinimo instrumentai — jie priverčia eksplicitiškai apibrėžti prielaidas, jas testuoti ir kvestionuoti. Tačiau jie nėra imunitetai nuo šališkumų — prielaidos, kurios įvedamos į modelį, vis tiek yra žmogaus sprendimai.

ESG ir dvigubos materialumo problema#

ESG (angl. Environmental, Social, and Governance — aplinkosauginiai, socialiniai ir valdysenos) veiksniai sukuria specifinę rizikos analizės problemą — dvigubą materialumą (angl. double materiality). Tai reiškia, kad ESG rizika turi dvi kryptis: rizika verslui (fizinė klimato rizika, reguliaciniai pokyčiai, tiekimo grandinės sutrikimai) ir rizika iš verslo (anglies emisijos, darbo sąlygos, valdymo skaidrumas). ES reguliacinis karkasas (ypač CSRD — angl. Corporate Sustainability Reporting Directive) jau reikalauja vertinti abi kryptis.

Praktikoje ESG integruojamas į rizikos analizę keliais būdais: kaip papildomas kintamasis jautrumo analizėje (pvz., anglies mokesčio poveikis EBIT maržai), kaip scenarijų ašis (pvz., „2°C scenarijus“ vs „business as usual“), ir kaip atskira kategorija rizikos registre. TCFD (angl. Task Force on Climate-related Financial Disclosures, 2017) rekomendacijos tapo de facto standartu klimato scenarijų analizei — jos reikalauja, kad įmonės analizuotų savo verslą bent dviejuose klimato scenarijuose (≤2°C ir >2°C globalaus atšilimo).

Dvigubo skaičiavimo pavojus. Integruojant ESG į DCF vertinimą, iškyla esminis metodologinis spąstas: ta pati ESG rizika gali būti „suskaičiuota du kartus“ — ir koreguojant pinigų srautus (pvz., sumažinant prognozes dėl galimo anglies mokesčio), ir padidinant diskonto normą (pvz., pridedant ESG rizikos premiją prie WACC). Cornell ir Damodaran (2020) perspėja, kad tai yra viena dažniausių vertinimo klaidų ESG kontekste: jei klimato rizika jau atspindėta pesimistiškesnėse pinigų srautų prognozėse, papildoma diskonto normos korekcija yra perteklinė, ir atvirkščiai. Analitikui būtina pasirinkti vieną integravimo kanalą — arba koreguoti pinigų srautus (skaidresnis ir labiau rekomenduojamas metodas), arba koreguoti diskonto normą — bet ne abu vienu metu. Šis principas taikytinas ne tik ESG, bet ir bet kuriai rizikos kategorijai, integruojamai į DCF modelį (pvz., šalies rizika, reguliacinė rizika).

Modelio rizika: struktūrinė klaida pavojingesnė nei parametrinė#

Galutinė ir galbūt svarbiausia kritinė pastaba: visi šiame skyriuje aptarti metodai vertina parametrinį neapibrėžtumą — t. y. klausimą „kas nutiktų, jei šie skaičiai būtų kitokie?“. Tačiau jie neapsaugo nuo modelio rizikos (angl. model risk) — galimybės, kad pats modelis yra struktūriškai neteisingas. Jei DCF modelis neapima esminio pajamų šaltinio, jokia jautrumo analizė to neaptiks — ji testuos tik tuos kintamuosius, kurie modelyje yra. Jei scenarijuose nenumatyta kokybiškai nauja rizikos kategorija (pvz., pandemija), Monte Carlo simuliacija neparodys jos tikimybės — ji simuliuos tik tas rizikas, kurios buvo eksplicitiškai apibrėžtos.

Tai Knighto neapibrėžtumo (žr. 3 skyrių) praktinė manifestacija: kiekybiniai metodai gali valdyti riziką (kur tikimybės ir scenarijai yra bent iš dalies žinomi), tačiau tikro neapibrėžtumo zonoje jie yra bejėgiai. Būtent todėl integruota rizikos analizė turi būti papildyta kokybinėmis strategijomis — nuolatiniu aplinkos stebėjimu, organizaciniu mokymusi, strateginiu lankstumu ir realių opcionų mąstymu (plačiau — 14 skyriuje). Kaip taikliai suformulavo Danielssonas (2011): „rizikos modeliai geriausiai veikia tada, kai jų labiausiai nereikia — ramiomis dienomis“.

4.1.4 Technologijų perspektyva#

Rizikos analizės metodų praktinis taikymas visada buvo neatsiejamai susietas su skaičiavimo technologijomis. Kiekviena technologinė banga ne tik suteikė naujų įrankių, bet ir pakeitė tai, kokie metodai tapo praktiškai prieinami ir kas juos gali taikyti.

Kompiuterizacijos era (1970–1995): nuo mainframe prie skaičiuoklės#

Hertzo (1964) Monte Carlo idėjos pirmaisiais dešimtmečiais liko daugiausia mainframe kompiuterių srityje — simuliacijos reikalavo programavimo įgūdžių ir brangios aparatūros. Lūžis įvyko devintajame dešimtmetyje, kai @RISK (Palisade Corporation, 1987) ir „Crystal Ball“ (Decisioneering, 1988) pavertė Monte Carlo simuliaciją skaičiuoklių priedu — analitikas galėjo vykdyti tūkstančius iteracijų tiesiai Microsoft Excel aplinkoje, be jokio programavimo. Tai demokratizavo simuliaciją: ji tapo prieinama ne tik investiciniams bankams, bet ir korporacijų finansų skyriams bei konsultacinėms įmonėms.

Jautrumo analizė irgi transformavosi kartu su skaičiuoklėmis. Excel Data Table funkcija (vienkryptė ir dvikryptė) leido automatizuoti „kas, jeigu?“ skaičiavimus — analitikas galėjo per minutes sugeneruoti tornado diagramai reikalingus duomenis, kai anksčiau tam reikėdavo valandų rankinio darbo. Sprendimų medžiams atsirado specializuota programinė įranga — DPL (Decision Programming Language) ir TreePlan, integruojamas į Excel.

Interneto ir duomenų era (1995–2020): nuo įrankio prie ekosistemos#

Šiuo laikotarpiu rizikos analizės programinė įranga perėjo konsolidacijos etapą: Oracle 2007 m. įsigijo Crystal Ball, integruodama jį į savo verslo analitikos platformą; Palisade (dabar Lumivero) tęsė @RISK plėtrą; Vose Software 2008 m. išleido ModelRisk kaip alternatyvą. Šie įrankiai tapo pramonės standartu — naftos ir dujų sektoriuje, farmacijoje ir infrastruktūros projektuose Monte Carlo simuliacija su @RISK ar Crystal Ball virto privalomu investicijų pagrindimo komponentu. Įvairios rinkos tyrimų įmonės (Grand View Research, Allied Market Research) vertina, kad globalinė rizikos analitikos rinka 2020-aisiais siekė ~$10–11 mlrd. ir sparčiai auga, atspindėdama augantį institucinį rizikos valdymo poreikį.

Šią erą taip pat transformavo dvi kitos jėgos. Pirma, atviro kodo ekosistema: Python (su NumPy, SciPy, pandas bibliotekomis), R ir kitos aplinkos leido vykdyti Monte Carlo simuliacijas su šimtais tūkstančių iteracijų, modeliuoti sudėtingas koreliacijos struktūras ir taikyti pažangius metodus (Latin Hypercube Sampling, Copula modelius) be komercinių licencijų kainos. Tai ypač svarbu akademiniam ir mažų įmonių sektoriui. Python SALib biblioteka padarė Sobol’ ir kitų GSA metodų skaičiavimą prieinamą vos keliomis kodo eilutėmis.

Antra, verslo analitikos (BI) platformos — Tableau, Power BI ir analogiški įrankiai — perkėlė rizikos vizualizaciją iš statinių lentelių į interaktyvius valdymo skydus (angl. dashboards). Rizikos žemėlapiai, tornado diagramos ir Monte Carlo histogramos tapo realiu laiku atnaujinamos vizualizacijos, prieinamos ne tik analitikams, bet ir vadovams. Bankų sektoriuje EBA (Europos bankininkystės institucija) ir ECB automatizuoti stresų testai — reguliariai vykdomi visoje Europos bankų sistemoje — pareikalavo sudėtingos IT infrastruktūros ir standartizuotų duomenų formatų, paverčiant rizikos analizę institucine funkcija, o ne pavienio analitiko užduotimi.

DI ir DI agentų era (2020→): rizikos analizė kaip dialogas#

DI agentai (angl. AI agents) — autonomiškos programinės sistemos, gebančios vykdyti daugiažingsnes užduočių sekas su minimaliu žmogaus įsikišimu.

Dirbtinis intelektas keičia rizikos analizę keliais lygmenimis. Rizikos identifikavimas — NLP (angl. Natural Language Processing — natūralios kalbos apdorojimo) modeliai geba automatiškai skenuoti naujienų srautus, reguliacinius dokumentus ir tiekimo grandinių duomenis, identifikuodami naujas rizikos kategorijas anksčiau nei tradiciniai metodai. Tai ypač vertinga ESG rizikos stebėsenai, kur reguliacinė aplinka kinta greitai.

Monte Carlo automatizavimas — generatyviniai DI modeliai jau geba kurti Monte Carlo simuliacijos konfigūracijas iš tekstinių prielaidų aprašymų, parinkti tikėtinus skirstinius ir interpretuoti rezultatus. Tai mažina barjerą — analitikas, nemokantis programuoti ar naudoti @RISK, gali gauti tikimybinį vertinimą per pokalbį su DI asistentu.

Rizikos registrų automatizavimas — DI gali analizuoti praeitinius projektus, identifikuoti pasikartojančias rizikas ir siūlyti standartizuotus rizikos registrus naujiems projektams, sumažindamas „tuščio lapo“ problemą.

Pažangūs DI metodai — mašininio mokymosi modeliai (gradientinio stiprinimo medžiai, neuroniniai tinklai) naudojami kredito rizikos vertinime, nes gali aptikti netiesinius ryšius tarp šimtų kintamųjų, kurių tradiciniai logistinės regresijos modeliai neaprėpia (Lessmann et al., 2015). Tačiau šie modeliai kelia „juodosios dėžės“ problemą — reguliuotojai (ypač EBA, ECB) reikalauja modelio paaiškinamumo (angl. explainability), todėl SHAP (angl. Shapley Additive Explanations — metodas, paaiškinantis kiekvieno kintamojo indėlį į modelio prognozę; Lundberg & Lee, 2017) ir LIME metodai tampa standartiniu DI modelių validavimo komponentu.

Tačiau DI era sustiprina jau aptartą „automation bias“ riziką: jei Monte Carlo simuliacijos konfigūraciją sugeneravo DI, o ne žmogus-ekspertas, gali kilti pernelyg didelis pasitikėjimas rezultatais, nesuprantant, kokios prielaidos buvo padarytos. Lietuvos ir Baltijos kontekste, kur fintech sektorius sparčiai auga (žr. 3 skyriaus technologijų perspektyvą), ši įtampa tarp prieinamumo ir kritiško supratimo yra ypač aktuali — Lietuvos banko priežiūros praktikoje rizikos valdymo reikalavimai fintech subjektams nuolat griežtėja, o Nasdaq Vilnius duomenys leidžia vykdyti Monte Carlo simuliacijas su realiais Baltijos rinkos parametrais. Technologijų ir rizikos analizės santykio gilesnė analizė pateikiama 19 skyriuje.

4.2 Analitiniai modeliai ir formulės#

Šiame poskyryje pateikiamos trylika pagrindinių rizikos analizės formulių su kintamųjų paaiškinimais ir skaitiniais pavyzdžiais. Formulės tęsia 3 skyriaus numeraciją (3.1–3.7) ir apima jautrumo analizės, scenarijų, Monte Carlo, verslo rizikos rodiklių ir integruotos analizės metodus.

Elastingumo koeficientas — jautrumo matas#

Elastingumo koeficientas (angl. elasticity) matuoja santykinį modelio rezultato jautrumą santykiniam kintamojo pokyčiui:

kur:

• $\varepsilon_{Y,X}$ — rezultato $Y$ elastingumas kintamojo $X$ atžvilgiu,
• $Y_0$ — bazinė rezultato reikšmė,
• $X_0$ — bazinė kintamojo reikšmė,
• $\Delta Y$, $\Delta X$ — atitinkami pokyčiai.

Skaitinis pavyzdys. Įmonės DCF modelyje bazinis NPV = €500 000 esant WACC = 10 %. Padidinus WACC iki 12 % (pokytis +20 %), NPV sumažėja iki €350 000 (pokytis –30 %).

Elastingumas –1,5 reiškia: 1 % WACC padidėjimas sumažina NPV maždaug 1,5 %. Neigiamas ženklas rodo atvirkštinį ryšį. Jei analogiškai apskaičiuotas pajamų augimo elastingumas yra +0,8, tai WACC yra beveik dvigubai svarbesnis nei pajamų augimas — ir tornado diagramoje jis bus viršuje.

Tikėtina vertė su scenarijų tikimybėmis#

Kai sukurti keli scenarijai su priskirtomis tikimybėmis, tikėtina vertė skaičiuojama kaip svertinis vidurkis:

kur:

• $E(V)$ — tikėtina vertė,
• $p_i$ — scenarijaus $i$ tikimybė ($\sum p_i = 1$),
• $V_i$ — rezultato vertė scenarijuje $i$,
• $n$ — scenarijų skaičius.

Skaitinis pavyzdys. Lietuvos eksporto įmonė vertina investicinį projektą trimis scenarijais:

Tikėtinas NPV yra teigiamas (€337 500), tačiau vien šis skaičius nepakanka sprendimui — būtina įvertinti ir rezultatų sklaidą.

Dispersija tarp scenarijų#

Scenarijų dispersija parodo, kiek rezultatai gali varijuoti apie tikėtinę vertę:

Skaitinis pavyzdys (tęsinys). Naudojant aukščiau pateiktus duomenis ($E(V) = 337\ 500$):

Standartinis nuokrypis €337 200 yra beveik lygus tikėtinai vertei €337 500. Variacijos koeficientas (žr. 3 skyriaus formulę (3.6)): $CV = 337\ 200 / 337\ 500 \approx 1{,}00$ — tai rodo labai aukštą santykinę riziką. Pesimistiniame scenarijuje projektas yra nuostolingas — sprendimų priėmėjas turi nuspręsti, ar 25 % nuostolingumo tikimybė yra priimtina.

Monte Carlo pasikliautinasis intervalas#

Po $N$ Monte Carlo iteracijų rezultato vidurkio tikslumas vertinamas standartiniu nuokrypiu:

kur:

• $\bar{Y}$ — simuliacijos rezultatų vidurkis,
• $s_Y$ — simuliacijos rezultatų standartinis nuokrypis,
• $N$ — iteracijų skaičius,
• $1{,}96$ — standartinio normaliojo skirstinio $z$-reikšmė 95 % pasikliovimo lygiui.

Skaitinis pavyzdys. Po 10 000 Monte Carlo iteracijų DCF modelis generuoja: $\bar{Y} = 420\ 000$ €, $s_Y = 310\ 000$ €.

Vidurkio 95 % pasikliautinasis intervalas: [€413 924; €426 076]. Tai reiškia, kad pats vidurkis yra nustatytas gana tiksliai (±€6 076). Tačiau atskirų rezultatų sklaida yra žymiai didesnė — 95 % rezultatų patenka į intervalą $\bar{Y} \pm 1{,}96 \times s_Y$, t. y. maždaug [–€187 600; €1 027 600]. Šis intervalas parodo tikrąjį neapibrėžtumo mastą.

Lūžio taškas — bazinis rizikos rodiklis#

Lūžio taškas (angl. break-even point) parodo minimalią pardavimų apimtį, reikalingą sąnaudoms padengti:

kur:

• $Q_{BE}$ — lūžio taško pardavimų kiekis (vienetais),
• $FC$ — pastoviosios sąnaudos (angl. fixed costs),
• $P$ — pardavimo kaina už vienetą,
• $VC$ — kintamosios sąnaudos už vienetą (angl. variable costs).

Skaitinis pavyzdys. Lietuvos gamybinė įmonė: pardavimo kaina €50/vnt., kintamosios sąnaudos €30/vnt., pastoviosios sąnaudos €120 000/metus.

Jei dabartiniai pardavimai yra 8 000 vnt., saugumo riba = $(8\ 000 - 6\ 000) / 8\ 000 = 25\%$ — tai reiškia, kad pardavimai gali sumažėti 25 % prieš įmonei pradedant nuostoliauti. Lūžio taškas yra intuityvus ir galingas komunikacijos įrankis — jis verčia klausimą „ar šis verslas pelningas?“ į „kiek vnt. / kokią apimtį turime parduoti?“ — tai klausimas, į kurį operatyviniai vadovai gali atsakyti geriau nei finansų analitikai.

Operacinio sverto laipsnis (DOL)#

Sverto matai — DOL ir DFL — buvo formaliai apibrėžti Van Horne’o (1968) vadovėlyje Financial Management and Policy, o Mandelker ir Rhee (1984) empiriškai parodė jų ryšį su sistemine rizika. Šie rodikliai tapo standartine verslo rizikos vertinimo priemone korporacijų finansuose.

DOL (angl. Degree of Operating Leverage) parodo, kiek procentų pakinta EBIT, kai pardavimai pakinta 1 %:

kur:

• $Q$ — pardavimų kiekis,
• $P - VC$ — ribinis pelnas (angl. contribution margin) už vienetą,
• $FC$ — pastoviosios sąnaudos.

Skaitinis pavyzdys (tęsinys). Esant 8 000 vnt. pardavimams:

DOL = 4,0 reiškia: 10 % pardavimų sumažėjimas sukels 40 % EBIT sumažėjimą. Kuo arčiau lūžio taško veikia įmonė, tuo didesnis DOL — ir tuo didesnė operacinė rizika. Tai ypač aktualu kapitalo imlioms pramonės šakoms (gamyba, energetika, logistika), kur pastoviosios sąnaudos sudaro didelę dalį.

Finansinio sverto laipsnis (DFL)#

DFL (angl. Degree of Financial Leverage) parodo, kaip grynasis pelnas reaguoja į EBIT pokyčius:

kur:

• $EBIT$ — pelnas prieš palūkanas ir mokesčius,
• $I$ — palūkanų sąnaudos.

Skaitinis pavyzdys. Įmonės EBIT = €40 000 (iš ankstesnio pavyzdžio), metinės palūkanos €15 000.

DFL = 1,6 reiškia: 10 % EBIT sumažėjimas sukels 16 % grynojo pelno sumažėjimą. Bendras sverto laipsnis DTL = DOL × DFL = 4,0 × 1,6 = 6,4. Tai reiškia, kad 10 % pardavimų sumažėjimas galiausiai sukels 64 % grynojo pelno sumažėjimą. Ši „sverto kaskada“ paaiškina, kodėl įmonės su aukštomis pastoviomis sąnaudomis ir didelėmis skolomis yra itin pažeidžiamos ekonominiam nuosmukiui — ir kodėl svertas yra vienas svarbiausių veiksnių rizikos analizėje (plačiau apie kapitalo struktūrą — 11 skyriuje).

Parametrinis VaR#

Parametrinis VaR (angl. Value at Risk) su normaliojo skirstinio prielaida:

kur:

• $VaR_{\alpha}$ — didžiausias tikėtinas nuostolis pasikliovimo lygiu $\alpha$,
• $\mu$ — tikėtina grąža (arba nulis, jei skaičiuojama per trumpą laikotarpį),
• $z_{\alpha}$ — standartinio normaliojo skirstinio kvantilis (pvz., $z_{0{,}99} = 2{,}326$; $z_{0{,}95} = 1{,}645$),
• $\sigma$ — grąžos standartinis nuokrypis.

Skaitinis pavyzdys. Investicinio portfelio vertė €1 000 000, dienos grąžos standartinis nuokrypis $\sigma = 1{,}5\%$, $\mu \approx 0$ (per vieną dieną tikėtina grąža artima nuliui, todėl formulė supaprastėja iki $VaR_{\alpha} \approx z_{\alpha} \cdot \sigma \cdot V$).

Dienos VaR 99 % lygiu = €34 890. Tai reiškia: 99 dienų iš 100 dienos nuostolis neviršys €34 890. Tačiau tą vieną dieną iš 100 nuostolis gali būti bet kokio dydžio — VaR apie tai nieko nepasako. Būtent todėl CVaR yra informatyvesnis matas.

CVaR / Expected Shortfall#

CVaR (angl. Conditional Value at Risk) arba tikėtinas trūkumas (angl. Expected Shortfall, ES):

Normalusis artinys:

kur $\phi(z_{\alpha})$ — standartinės normaliosios tankio funkcijos reikšmė taške $z_{\alpha}$.

Skaitinis pavyzdys (tęsinys). Esant $z_{0{,}99} = 2{,}326$, $\phi(2{,}326) \approx 0{,}0264$:

CVaR = €39 600 — tai reiškia, kad vidutinis nuostolis tomis dienomis, kai VaR riba viršijama (t. y. tą 1 dieną iš 100), yra €39 600. CVaR yra didesnis nei VaR (€39 600 > €34 890), nes jis apima visą uodegos sritį, ne tik ribinį tašką.

Tikėtina piniginė vertė sprendimų medyje (EMV)#

EMV (angl. Expected Monetary Value) leidžia palyginti sprendimo alternatyvas atsižvelgiant į neapibrėžtumą:

kur:

• $p_k$ — $k$-ojo rezultato tikimybė,
• $C_k$ — $k$-ojo rezultato piniginė vertė,
• $I_0$ — pradiniai investicijos kaštai,
• $m$ — galimų rezultatų skaičius.

Skaitinis pavyzdys. Įmonė svarsto du variantus: (A) Pilna investicija €500 000 iš karto, arba (B) Pilotinis projektas €150 000, po kurio — sprendimas dėl pilnos plėtros (papildomi €400 000).

Variantas A (pilna investicija):

• Sėkmė (p = 0,6): grąža €900 000
• Nesėkmė (p = 0,4): grąža €200 000

Variantas B (stadijinis):

• Pilotas atskleidžia informaciją; jei pilotas sėkmingas (p = 0,6): investuojama pilnai, sėkmės tikimybė padidėja iki 0,85
  • $EMV_{plėtra} = 0{,}85 \times 900\ 000 + 0{,}15 \times 200\ 000 - 400\ 000 = 765\ 000 + 30\ 000 - 400\ 000 = 395\ 000$ €
• Jei pilotas nesėkmingas (p = 0,4): sustabdoma, nuostolis tik €150 000

Nors variantas A turi didesnį EMV (€120 000 vs €87 000), variantas B turi esmingai mažesnę riziką: blogiausiu atveju prarandama tik €150 000, o ne €300 000. Tai realių opcionų logika praktikoje — stadijinė investicija turi vertę net kai jos EMV yra mažesnis, nes ji suteikia informacijos vertę ir atsitraukimo galimybę (plačiau — 14 skyriuje).

Omega koeficientas#

Omega koeficientas (Keating & Shadwick, 2002) yra alternatyvus rizikos-grąžos matas, apimantis visą grąžų pasiskirstymą — ne tik vidurkį ir dispersiją, bet ir asimetriją bei ekscesą:

kur:

• $F(r)$ — grąžų kumuliacinė pasiskirstymo funkcija,
• $\theta$ — pasirinktas slenkstinis grąžų lygis (angl. threshold).

Skaitiklis reiškia „laimėjimų plotą“ (tikėtiną perteklių virš $\theta$), vardiklis — „praradimų plotą“ (tikėtiną trūkumą žemiau $\theta$). Kai $\Omega > 1$, investicija „laimi“ daugiau nei „pralaimi“ prie pasirinkto slenksčio. Omega pranašumas prieš Sharpe’o koeficientą — jis nereikalauja normalumo prielaidos ir natūraliai apima visus pasiskirstymo momentus.

Skaitinis pavyzdys. Turime du investicinius fondus, kurių mėnesinės grąžos per 12 mėnesių (%) yra:

• Fondas A: 2, 4, –1, 3, –2, 5, 1, –3, 6, 2, –1, 4 (vidurkis ≈ 1,67 %)
• Fondas B: 1, 1, 1, 2, 0, 1, 2, 1, 0, 1, 1, 1 (vidurkis ≈ 1,00 %)

Pasirinkus slenkstį $\theta = 0\%$ (t. y. „ar fondas uždirba?“):

• Fondas A: 9 mėnesiai virš 0 %, 3 mėnesiai žemiau. Laimėjimų suma = 27 %, praradimų suma = 6 %. $\Omega_A(0) \approx 27/6 = 4{,}5$.
• Fondas B: 10 mėnesių virš 0 %, 2 mėnesiai lygiai 0 %. Laimėjimų suma = 12 %, praradimų suma ≈ 0 %. $\Omega_B(0) \to \infty$ (praktiškai labai didelis).

Nors fondas A turi aukštesnį vidurkį, fondas B turi aukštesnę Omega — jis beveik niekada nepraranda. Tai iliustruoja, kodėl Omega koeficientas gali atskleisti rizikos-grąžos aspektus, kurių vidurkis ir standartinis nuokrypis neparodo — ypač kai grąžų pasiskirstymas yra asimetriškas.

Maximum Drawdown (MDD) ir Calmar koeficientas#

Maximum Drawdown — didžiausias nuosmukis nuo piko iki dugno — yra intuityvus ir psichologiškai reikšmingas rizikos matas:

kur $V_{pikas,t}$ — didžiausia pasiekta portfelio vertė iki momento $t$, $V_{dugnas,t}$ — mažiausia vertė po piko.

Praktinė interpretacija: –30 % nuosmukis reikalauja +42,9 % atsigavimo vien tam, kad grįžtų į pradinį lygį ($1/0{,}7 - 1 \approx 0{,}429$). Calmar koeficientas (angl. Calmar ratio) normalizuoja grąžą pagal MDD:

Calmar leidžia palyginti investicijas ar verslo projektus pagal rizikos-grąžos santykį atsižvelgiant į didžiausią nuosmukį, o ne tik į dispersiją.

Sobol’ pirmos eilės indeksas — globalios jautrumo analizės matas#

Formulė (4.1) matuoja vietinį jautrumą (vienas kintamasis, vienas pokytis). Sobol’ indeksas (Saltelli et al., 2008) yra globalus jautrumo matas, vertinantis kintamojo poveikį visame jo galimų reikšmių diapazone:

kur:

• $S_i$ — pirmos eilės Sobol’ indeksas kintamajam $X_i$,
• $\mathbb{E}[V \mid X_i]$ — sąlyginė tikėtina modelio rezultato vertė, žinant $X_i$ reikšmę,
• $\text{Var}(V)$ — bendra modelio rezultato dispersija.

$S_i = 0{,}45$ reiškia: 45 % viso modelio neapibrėžtumo kyla iš kintamojo $X_i$. Totalinis indeksas $S_{T_i}$ papildomai apima visas kintamojo sąveikas su kitais kintamaisiais. Skirtumas $S_{T_i} - S_i$ parodo sąveikų svarbą.

Skaitinis pavyzdys. BaltTech DCF modelyje (žr. 4.4) GSA su 10 000 iteracijų duoda: pardavimų apimties $S_1 = 0{,}42$, kainos $S_2 = 0{,}28$, žaliavų $S_3 = 0{,}15$, WACC $S_4 = 0{,}08$, darbo sąnaudų $S_5 = 0{,}04$. Suma $\sum S_i = 0{,}97 \approx 1$ — modelis beveik adityvus, sąveikos nereikšmingos. Tai patvirtina tornado analizės išvadą: pardavimai ir kaina yra du svarbiausi veiksniai. Tačiau jei $\sum S_i$ būtų, pvz., 0,65 — tai reikštų, kad 35 % neapibrėžtumo kyla iš kintamųjų sąveikų, kurių vienkryptė jautrumo analizė visiškai nematė.

Formulių suvestinė#

4.3 Apibendrinimas#

Šis skyrius pristatė keturis pagrindinius rizikos analizės metodus ir jų integraciją į vientisą sistemą. Pagrindinės išvados:

1. Jautrumo analizė identifikuoja kritinius kintamuosius ir yra efektyviausias komunikacijos įrankis (tornado diagrama), tačiau ignoruoja kintamųjų sąveikas ir neparodo tikimybių. Saltelli ir kolegų (2008) globali jautrumo analizė sprendžia sąveikų problemą, bet reikalauja Monte Carlo infrastruktūros.

2. Scenarijų analizė yra vienintelis metodas, gebantis apčiuopti struktūrinius lūžius — kokybiškai naujus pasaulio vaizdus, kurių kiekybiniai modeliai neapima. Tačiau ji priklauso nuo kūrėjų mentalinio horizonto ir yra pažeidžiama inkaravimo bei siauro rėmo šališkumų. Šiuolaikinis 2×2 matricinis požiūris pranašesnis už klasikinį optimistinį/bazinį/pesimistinį modelį.

3. Monte Carlo simuliacija transformuoja deterministinį modelį į tikimybinį, leisdama kiekybiškai vertinti neapibrėžtumą per pasiskirstymo funkciją. Jos galia ir jos silpnybė yra ta pati: ji propaguoja prielaidas su statistiniu tikslumu — jei prielaidos geros, rezultatai informatyvūs; jei prielaidos blogos, rezultatai klaidingi, bet atrodo moksliški.

4. Rizikos rodikliai (VaR, CVaR, Omega, MDD, DOL/DFL, lūžio taškas, stresiniai testai) suteikia standartizuotas metrikas sprendimų priėmimui. VaR nekoherentiškumo problema ir Bazilio perėjimas prie CVaR iliustruoja platesnę pamoką: kiekybinis matas yra tiek naudingas, kiek gerai suprantamos jo ribos. Omega koeficientas ir Maximum Drawdown papildo tradicinius dispersijos matus, atsižvelgdami į visą pasiskirstymo formą ir didžiausius nuostolius.

5. Integruota sistema sujungia visus metodus į laipsniškai gilėjančią seką: tornado (kas svarbu?) → scenarijai (kas galėtų nutikti?) → Monte Carlo (kokia tikimybė?) → rodikliai (kiek rizikuojame?) → sprendimas. Tačiau ši sistema turi būti taikoma proporcingai — „analizės paralyžius“ yra reali rizika.

6. Modelio rizika — galimybė, kad pats modelis yra struktūriškai neteisingas — yra svarbiausia ir sunkiausiai valdoma rizikos dimensija. Jokie kiekybiniai metodai negali apsaugoti nuo nežinomų nežinomųjų — todėl kiekybinė analizė turi būti papildyta kokybinėmis strategijomis (aplinkos stebėjimas, organizacinis mokymasis, strateginis lankstumas, realių opcionų mąstymas).

Šio skyriaus metodai rengia pagrindą tolesnėms knygos temoms: 5 skyriuje diversifikacijos sprendimai remiasi portfelio rizikos matais; 9 skyriuje finansinių ataskaitų analizė naudoja sverto rodiklius; 11 skyriuje WACC jautrumo analizė parodo, kaip diskonto normos neapibrėžtumas veikia visą vertinimą; 14 skyriuje realieji opcionai ir sprendimų medžiai taikomi startuolių ir M&A vertinimui.

4.4 Atvejo analizė ir praktiniai pavyzdžiai#

Atvejis 1: UAB „BaltTech“ — integruota gamybos investicijos rizikos analizė#

Tai hipotetinis mokomasis pavyzdys, paremtas tipine Lietuvos baldų gamybos sektoriaus situacija. Skaičiai iliustraciniai, atspindintys vidutinės eksportuojančios gamybos įmonės proporcijas.

Situacija#

UAB „BaltTech“ yra Lietuvos baldų gamybos įmonė, eksportuojanti ~60 % produkcijos į Vokietiją ir Skandinaviją. Įmonė svarsto €1,5 mln. investiciją į automatizuotą gamybos liniją, kuri turėtų sumažinti darbo sąnaudas 30 % ir padidinti gamybos pajėgumus 40 %. Projekto trukmė — 7 metai. Bazinis WACC = 11 %, bazinis NPV = €420 000.

Valdyba prašo finansų direktoriaus atlikti pilną rizikos analizę prieš priimant sprendimą.

1 žingsnis: Jautrumo analizė (tornado diagrama)#

Finansų direktorius identifikuoja penkis kritinius kintamuosius ir atlieka vienkryptę jautrumo analizę, kiekvieną kintamąjį keisdamas ±20 % nuo bazinės reikšmės:

Tornado diagrama atskleidžia: pardavimų apimtis (€1 020 000 intervalas) ir produkcijos kaina (€750 000) yra du kritiniai veiksniai. WACC ir darbo sąnaudos turi mažesnį poveikį. Elastingumo koeficientai: pardavimų apimtis $\varepsilon \approx 1{,}21$; kaina $\varepsilon \approx 0{,}89$; žaliavos $\varepsilon \approx -0{,}71$.

Išvada: didžiausias dėmesys turi būti skiriamas pardavimų prognozei ir kainų strategijai.

2 žingsnis: Scenarijų analizė#

Remiantis tornado rezultatais, kuriami trys scenarijai, keičiant visus kintamuosius koordinuotai (ne tik po vieną):

$CV = 358\ 000 / 402\ 500 \approx 0{,}89$ — aukšta santykinė rizika. Pesimistiniame scenarijuje projektas nuostolingas.

3 žingsnis: Monte Carlo simuliacija#

Penkiems kintamiesiems priskiriami tikimybiniai skirstiniai:

• Pardavimai: trikampis (8 000; 10 000; 12 500) — ekspertinis vertinimas
• Kaina: normalusis ($\mu = 150$, $\sigma = 12$) — istoriniai duomenys
• Žaliavos: lognormalusis ($\mu = 60$, $\sigma = 10$) — žaliavų rinkos asimetrija
• WACC: normalusis ($\mu = 11\%$, $\sigma = 1{,}5\%$)
• Darbo sąnaudos: trikampis (20; 25; 32) — su automatizacijos neapibrėžtumu

Koreliacija: pardavimai ↔ kaina: $\rho = +0{,}3$ (geresnė paklausa → galimybė didinti kainą); pardavimai ↔ žaliavos: $\rho = +0{,}2$ (bendras ekonomikos ciklas).

Po 10 000 iteracijų rezultatai:

• Vidutinis NPV: €405 000
• Standartinis nuokrypis: €310 000
• $P(NPV < 0) = 18{,}5\%$
• 5-asis procentilis: –€180 000
• 95-asis procentilis: €920 000

Interpretacija: beveik vienas iš penkių scenarijų yra nuostolingas. Blogiausiu tikėtinu atveju (5-asis procentilis) nuostolis siekia –€180 000.

4 žingsnis: Sprendimas — stadijinė investicija#

Atsižvelgdama į 18,5 % nuostolingumo tikimybę, valdyba nusprendžia taikyti stadijinį investavimą: pirmiausia investuoti €500 000 į pilotinę liniją (6 mėnesių bandomasis periodas), ir tik pasitvirtinus rezultatams — investuoti likusius €1 000 000. Tai sukuria:

• Atsitraukimo opcioną: jei pilotas nepasiteisina, nuostolis ribojamas €500 000 (o ne €1 500 000)
• Informacijos vertę: pilotas atskleidžia realius pardavimų ir kaštų parametrus, mažindamas Monte Carlo prielaidų neapibrėžtumą

Sprendimų medžio analizė (pagal formulę (4.10)) patvirtina, kad stadijinis požiūris, nors ir turintis šiek tiek mažesnį EMV, turi fundamentaliai geresnį rizikos profilį.

Atvejis 2: Lietuvos energetikos sektoriaus scenarijų 2×2 matrica#

Mišrus mokomasis pavyzdys: scenarijų kontekstas paremtas realiais Lietuvos energetikos sektoriaus įvykiais (Ignalinos uždarymas, SGD terminalas, NordBalt/LitPol jungtys, 2022 m. energetikos krizė), tačiau konkreti įmonė ir investicijos duomenys yra hipotetiniai.

Situacija#

Atsinaujinančios energetikos bendrovė planuoja €5 mln. investiciją į saulės elektrinių parką Lietuvoje. Du esminiai neapibrėžtumo šaltiniai: (1) elektros kainų režimas (aukštas vs žemas, priklausantis nuo geopolitikos ir energijos rinkos) ir (2) reguliacinė aplinka (palanki vs ribojanti, priklausanti nuo ES politikos ir nacionalinių sprendimų).

2×2 matrica#

Keturi scenarijai nėra „geras–blogas“ spektras — jie reikalauja skirtingų strategijų:

• „Žaliasis bumas“ → agresyvi plėtra, ilgalaikės PPA sutartys
• „Tyli transformacija“ → fokusas į reguliacinių dotacijų administravimą
• „Laisvoji banga“ → orientacija į didmeninę rinką, mažesnė priklausomybė nuo valstybės
• „Ramybė prieš audrą“ → minimalus investavimas, laukimo opcionas

Lietuvos kontekstas. Šie scenarijai nėra abstraktūs — jie atspindi realias Lietuvos energetikos sektoriaus trajektorijas. Po Ignalinos AE uždarymo (2009), SGD terminalo paleidimo (2014), NordBalt ir LitPol Link jungčių atsiradimo — Lietuvos energetikos rinka fundamentaliai pasikeitė. 2022 m. Rusijos invazija į Ukrainą sukėlė energijos kainų šuolį, kuris vienų investuotojų projektus pavertė itin pelningais, o kitų — nuostolingais. Scenarijų analizė būtų leidusi pasiruošti abiem trajektorijoms, užuot statant viską ant vieno „bazinio“ scenarijaus.

4.5 Excel modelio konstravimas#

Tikslas#

Sukurti integruotą rizikos analizės modelį, apimantį jautrumo analizę (tornado diagrama), scenarijų analizę ir supaprastintą Monte Carlo simuliaciją. Modelis taikomas investicinio projekto NPV vertinimui — analogiškai „BaltTech“ atvejui (žr. 4.4).

Failo struktūra#

Žingsnis po žingsnio instrukcija#

1 žingsnis: Lapas „Duomenys“#

Sukuriama bazinių kintamųjų lentelė:

2 žingsnis: Lapas „DCF“#

Sukuriamas 7 metų pinigų srautų modelis:

• Eilutė „Pajamos“: =Duomenys!B4*Duomenys!B5
• Eilutė „Kintamosios sąnaudos“: =Duomenys!B4*Duomenys!B6
• Eilutė „Ribinis pelnas“: =Pajamos-Kintamosios
• Eilutė „EBIT“: =Ribinis pelnas-Duomenys!B7
• Eilutė „Grynasis pelnas“: =EBIT*(1-Duomenys!B9)
• Eilutė „FCF“: =Grynasis pelnas + Nusidėvėjimas (nusidėvėjimas = Investicija/7)
• NPV langelis: =NPV(Duomenys!B8, FCF_range) - Duomenys!B2

3 žingsnis: Lapas „Jautrumas“#

Vienkryptė analizė (tornado duomenys):

Kiekvienoje eilutėje pakeičiama viena Duomenų reikšmė, perskaičiuojamas NPV ir grąžinama bazinė reikšmė. Tornado diagramai: stulpelinė diagrama su kintamaisiais Y ašyje, NPV intervalais X ašyje.

Dvikryptė analizė:

• Naudoti Excel duomenų lentelę (angl. Data Table; meniu: What-If Analysis → Data Table): WACC eilutėje (8 %–14 %), augimo tempas stulpelyje (0 %–6 %) → NPV matrica.

4 žingsnis: Lapas „Scenarijai“#

• E(NPV): =SUMPRODUCT(tikimybės, NPV_eilutė)
• σ²: =SUMPRODUCT(tikimybės, (NPV_eilutė - E_NPV)^2)
• σ: =SQRT(σ²)
• CV: =σ / E_NPV

5 žingsnis: Lapas „MonteCarlo“#

Supaprastinta MC su Excel RAND():

Rezultatų analizė:

• =AVERAGE(E:E) — vidutinis NPV
• =STDEV(E:E) — standartinis nuokrypis
• =COUNTIF(E:E,"<0")/1000 — $P(NPV < 0)$
• =PERCENTILE(E:E, 0.05) — 5-asis procentilis
• Histograma: Įterpimas → Diagrama → Histograma (angl. Insert → Chart → Histogram) iš E stulpelio

Pastaba: Excel RAND() generuoja tik nepriklausomus atsitiktinius skaičius. Koreliacijoms modeliuoti reikia @RISK, Python ar specialių VBA makrokomandų. Šiame modelyje koreliacijos nemodelizuojamos — tai supaprastinimas, kurį studentas turėtų suvokti ir kritiškai vertinti.

6 žingsnis: Lapas „Rodikliai“#

Formulių sąrašas#

Failai:

• Excel/CH04/Ch04_Formula_Guide.md — žingsnis-po-žingsnio instrukcija (vienas tiesos šaltinis)
• Excel/CH04/Ch04_Model_Master.xlsm — pilnas modelis su formulėmis ir diagramomis (dėstytojui)
• Excel/CH04/Ch04_Model_Student.xlsx — šablonas su tuščiais formulių langeliais (studentui pildyti pagal Formula Guide)

4.6 Klausimai ir praktinės užduotys#

Diskusiniai klausimai#

1. Kodėl trijų scenarijų modelis (optimistinis/bazinis/pesimistinis) gali būti pavojingesnis nei vieno scenarijaus naudojimas? Aptarkite inkaravimo efektą ir siauro rėmo problemą.

2. Monte Carlo simuliacija generuoja NPV pasiskirstymą su tiksliais statistiniais parametrais. Kodėl šis „tikslumas“ gali būti klaidinantis? Paaiškinkite „garbage in — beautifully distributed garbage out“ problemą.

3. VaR buvo bankų reguliavimo pagrindas (Basel II), tačiau Bazilio komitetas perėjo prie CVaR (Basel III / FRTB). Kodėl? Paaiškinkite subadityvumo aksiomą ir jos praktinę reikšmę.

4. Kaip „analizės paralyžius“ pasireiškia praktikoje? Aptarkite proporcingumo principą: kokio lygio rizikos analizė tinkama (a) €50 000 investicijai, (b) €5 mln. investicijai, (c) €50 mln. infrastruktūros projektui?

5. Danielssonas (2011) teigia: „rizikos modeliai geriausiai veikia tada, kai jų labiausiai nereikia.“ Ką tai reiškia ir kokios praktinės pasekmės rizikos analizės procesui?

Skaičiavimo užduotys#

1. Jautrumo analizė. Investicinio projekto bazinis NPV = €200 000. Atlikus vienkryptę jautrumo analizę, nustatyta: (a) padidinus WACC nuo 10 % iki 12 %, NPV sumažėja iki €120 000; (b) padidinus pajamas 10 %, NPV padidėja iki €280 000. Apskaičiuokite elastingumo koeficientus ir nustatykite, kuris kintamasis svarbesnis.

2. Scenarijų analizė. Trys scenarijai: optimistinis (NPV = €500 000, p = 0,20), bazinis (NPV = €180 000, p = 0,55), pesimistinis (NPV = –€80 000, p = 0,25). Apskaičiuokite: (a) tikėtinę vertę $E(V)$; (b) standartinį nuokrypį $\sigma$; (c) variacijos koeficientą $CV$; (d) tikimybę, kad projektas bus nuostolingas.

3. Svertas. Įmonės duomenys: pardavimai 5 000 vnt., kaina €80/vnt., kintamosios sąnaudos €50/vnt., pastoviosios sąnaudos €90 000, palūkanos €20 000. Apskaičiuokite: (a) lūžio tašką; (b) saugos maržą; (c) DOL; (d) DFL; (e) DTL. Kaip 15 % pardavimų sumažėjimas paveiktų grynąjį pelną?

4. VaR ir CVaR. Portfelio vertė €2 000 000, dienos grąžos standartinis nuokrypis 2 %, tikėtina dienos grąža ≈ 0 %. Apskaičiuokite: (a) dienos VaR 95 % lygiu ($z_{0{,}95} = 1{,}645$); (b) dienos VaR 99 % lygiu ($z_{0{,}99} = 2{,}326$); (c) CVaR 99 % lygiu ($\phi(2{,}326) \approx 0{,}0264$). Paaiškinkite, kodėl CVaR > VaR.

Savarankiško darbo užduotys#

1. Excel modelis. Sukurkite investicinio projekto rizikos analizės modelį pagal 4.5 poskyrio instrukciją. Naudodami savo sukurtą modelį, atlikite: (a) tornado analizę 5 kintamiesiems; (b) trijų scenarijų analizę su tikėtina verte; (c) 1 000 iteracijų Monte Carlo simuliaciją. Pateikite rezultatų interpretaciją ir rekomendaciją — investuoti ar ne.

2. Scenarijų 2×2 matrica. Pasirinkite Lietuvos verslo sektorių (pvz., fintech, gamyba, logistika, nekilnojamasis turtas). Identifikuokite du svarbiausius neapibrėžtumo šaltinius. Sukurkite 2×2 scenarijų matricą su keturiais kokybiškai skirtingais scenarijais. Kiekvienam scenarijui aprašykite: (a) naratyvą (kas nutiko?); (b) pagrindinius finansinius parametrus; (c) rekomenduojamą strategiją. Paaiškinkite, kuo šis požiūris pranašesnis už klasikinį optimistinį/bazinį/pesimistinį modelį.

Literatūros sąrašas#

Artzner, P., Delbaen, F., Eber, J.-M., & Heath, D. (1999). Coherent measures of risk. Mathematical Finance, 9(3), 203–228.

Brealey, R. A., Myers, S. C., & Allen, F. (2020). Principles of Corporate Finance (13th ed.). McGraw-Hill.

Chapman, C., & Ward, S. (2003). Project Risk Management: Processes, Techniques and Insights (2nd ed.). Wiley.

Chermack, T. J., & Lynham, S. A. (2002). Definitions and outcome variables of scenario planning. Human Resource Development Review, 1(3), 366–383. https://doi.org/10.1177/1534484302013006

Chermack, T. J., & Nimon, K. (2008). The effects of scenario planning on participant decision-making style. Human Resource Development Quarterly, 19(4), 351–372. https://doi.org/10.1002/hrdq.1245

Cornell, B., & Damodaran, A. (2020). Valuing ESG: Doing good or sounding good? The Journal of Impact and ESG Investing, 1(1), 76–93. https://doi.org/10.3905/jesg.2020.1.1.076

COSO. (2017). Enterprise Risk Management — Integrating with Strategy and Performance. Committee of Sponsoring Organizations of the Treadway Commission.

Damodaran, A. (2012). Investment Valuation: Tools and Techniques for Determining the Value of Any Asset (3rd ed.). Wiley.

Danielsson, J. (2011). Financial Risk Forecasting. Wiley.

Eckhardt, R. (1987). Stan Ulam, John von Neumann, and the Monte Carlo method. Los Alamos Science, 15, 131–137.

Haeffner, M., Leone, D., Coons, L., & Chermack, T. (2012). The effects of scenario planning on perceptions of learning organization characteristics. Human Resource Development Quarterly, 23(4), 519–542. https://doi.org/10.1002/hrdq.21147

Hertz, D. B. (1964). Risk analysis in capital investment. Harvard Business Review, 42(1), 95–106.

Hillier, F. S. (1963). The derivation of probabilistic information for the evaluation of risky investments. Management Science, 9(3), 443–457.

Iman, R. L., & Conover, W. J. (1982). A distribution-free approach to inducing rank correlation among input variables. Communications in Statistics — Simulation and Computation, 11(3), 311–334. https://doi.org/10.1080/03610918208812265

Jorion, P. (2007). Value at Risk: The New Benchmark for Managing Financial Risk (3rd ed.). McGraw-Hill.

Kahn, H., & Wiener, A. J. (1967). The Year 2000: A Framework for Speculation on the Next Thirty-Three Years. Macmillan.

Keating, C., & Shadwick, W. F. (2002). A universal performance measure. Journal of Performance Measurement, 6(3), 59–84.

Knight, F. H. (1921). Risk, Uncertainty and Profit. Houghton Mifflin.

Lessmann, S., Baesens, B., Seow, H.-V., & Thomas, L. C. (2015). Benchmarking state-of-the-art classification algorithms for credit scoring: An update of research. European Journal of Operational Research, 247(1), 124–136. https://doi.org/10.1016/j.ejor.2015.05.030

Lundberg, S. M., & Lee, S.-I. (2017). A unified approach to interpreting model predictions. Advances in Neural Information Processing Systems 30 (pp. 4765–4774). Neural Information Processing Systems Foundation. https://proceedings.neurips.cc/paper/7062-a-unified-approach-to-interpreting-model-predictions.pdf

Mandelker, G. N., & Rhee, S. G. (1984). The impact of the degrees of operating and financial leverage on systematic risk of common stock. Journal of Financial and Quantitative Analysis, 19(1), 45–57. https://doi.org/10.2307/2331000

McKay, M. D., Beckman, R. J., & Conover, W. J. (1979). A comparison of three methods for selecting values of input variables in the analysis of output from a computer code. Technometrics, 21(2), 239–245. https://doi.org/10.1080/00401706.1979.10489755

Metropolis, N., & Ulam, S. (1949). The Monte Carlo method. Journal of the American Statistical Association, 44(247), 335–341.

Palar, P. S., Zuhal, L. R., & Shimoyama, K. (2023). Enhancing the explainability of regression-based polynomial chaos expansion by Shapley additive explanations. Reliability Engineering & System Safety, 232, 109045.

Phadnis, S., Caplice, C., Sheffi, Y., & Singh, M. (2015). Effect of scenario planning on field experts’ judgment of long-range investment decisions. Strategic Management Journal, 36(9), 1401–1411. https://doi.org/10.1002/smj.2293

Rubinstein, R. Y., & Kroese, D. P. (2017). Simulation and the Monte Carlo Method (3rd ed.). Wiley.

Saltelli, A., Ratto, M., Andres, T., Campolongo, F., Cariboni, J., Gatelli, D., Saisana, M., & Tarantola, S. (2008). Global Sensitivity Analysis: The Primer. Wiley.

Schoemaker, P. J. H. (1995). Scenario planning: A tool for strategic thinking. Sloan Management Review, 36(2), 25–40.

Schwartz, P. (1991). The Art of the Long View: Planning for the Future in an Uncertain World. Currency Doubleday.

Taleb, N. N. (2007). The Black Swan: The Impact of the Highly Improbable. Random House.

TCFD. (2017). Recommendations of the Task Force on Climate-related Financial Disclosures. Financial Stability Board.

Van Horne, J. C. (1968). Financial Management and Policy. Prentice-Hall.

van der Heijden, K. (2005). Scenarios: The Art of Strategic Conversation (2nd ed.). Wiley.

Vose, D. (2008). Risk Analysis: A Quantitative Guide (3rd ed.). Wiley.

Wack, P. (1985). Scenarios: Uncharted waters ahead. Harvard Business Review, 63(5), 73–89.