\[
% MathJax has no bm package; redefine \bm in terms of \boldsymbol, which it supports natively
\newcommand{\bm}[1]{\boldsymbol{#1}}
% --- Operators -----------------------------------------------------------
% Expectation operator
\DeclareMathOperator{\E}{E}
% Variance operator
\DeclareMathOperator{\Var}{Var}
% Covariance operator
\DeclareMathOperator{\Cov}{Cov}
% Correlation operator
\DeclareMathOperator{\Corr}{Corr}
% Rank operator (Spearman rank correlation)
\DeclareMathOperator{\Rank}{Rank}
% Skewness operator
\DeclareMathOperator{\Skewness}{Skewness}
% Kurtosis operator
\DeclareMathOperator{\Kurtosis}{Kurtosis}
% Difference/differential operator (upright d per ISO 80000-2)
\newcommand{\Diff}{\mathrm{d}}
% --- Risk Measures -------------------------------------------------------
% Value at Risk
\DeclareMathOperator{\VaR}{VaR}
% Expected Shortfall
\DeclareMathOperator{\ES}{ES}
% Marginal Value at Risk
\DeclareMathOperator{\MVaR}{MVaR}
% Component Value at Risk
\DeclareMathOperator{\CompVaR}{CVaR}
% Incremental Value at Risk
\DeclareMathOperator{\IVaR}{IVaR}
% Component Expected Shortfall
\DeclareMathOperator{\CompES}{CES}
% Weighted sensitivity for risk factor in FRTB SBM bucket
\DeclareMathOperator{\WS}{WS}
% Hedge-benefit ratio for FRTB DRC bucket
\DeclareMathOperator{\HBR}{HBR}
% Stress scenario risk measure for FRTB NMRF capital add-on
\DeclareMathOperator{\SES}{SES}
% Jump-to-default exposure for obligor (FRTB DRC)
\DeclareMathOperator{\JTD}{JTD}
% Loss given default for obligor (FRTB DRC)
\DeclareMathOperator{\LGD}{LGD}
% Default risk weight for obligor (FRTB DRC)
\DeclareMathOperator{\RW}{RW}
% ES ratio function
\newcommand{\ESratio}{\lambda}
% Arbitrary risk measure function (coherence axioms)
\newcommand{\RiskMeasure}{\varphi}
% Asset position (coherence axioms)
\newcommand{\Asset}{A}
% First named asset in examples
\newcommand{\AssetA}{A}
% Second named asset in examples
\newcommand{\AssetB}{B}
% Third named asset in examples
\newcommand{\AssetC}{C}
% Constant (risk measure axioms)
\newcommand{\Constant}{c}
% Bucket-level aggregate sensitivity in FRTB SBM
\newcommand{\BucketAgg}{S}
% Cross-bucket correlation in FRTB SBM
\newcommand{\CrossBucketCorr}{\gamma}
% Notional amount for an instrument (FRTB DRC and RRAO)
\newcommand{\Notional}{\text{Notional}}
% --- Distributions -------------------------------------------------------
% Binomial distribution
\DeclareMathOperator{\Binomial}{Binomial}
% Uniform distribution
\DeclareMathOperator{\Uniform}{Uniform}
% Normal distribution
\newcommand{\NormalDist}{\mathcal{N}}
% Student-t CDF
\newcommand{\StudentCDF}{t}
% Student-t PDF
\newcommand{\StudentPDF}{f}
% --- Returns -------------------------------------------------------------
% Price
\newcommand{\Price}{P}
% Simple (arithmetic) return
\newcommand{\SimpleReturns}{R}
% Log (compound) return
\newcommand{\CompoundReturns}{y}
% Matrix of compound returns
\newcommand{\ReturnMatrix}{\bm{y}}
% Dividend payment
\newcommand{\dividend}{d}
% --- Volatility ----------------------------------------------------------
% Realised variance
\DeclareMathOperator{\RealisedVar}{RV}
% Realised volatility (square root of realised variance)
\DeclareMathOperator{\RealisedVol}{RVol}
% Volatility
\newcommand{\Vol}{\sigma}
% GARCH intercept parameter
\newcommand{\GARCHconst}{\omega}
% ARCH coefficient (weight on past squared returns)
\newcommand{\ARCHcoeff}{\alpha}
% GARCH coefficient (weight on past variance)
\newcommand{\GARCHcoeff}{\beta}
% EWMA decay factor (typically 0.94)
\newcommand{\EWMAdecay}{\lambda}
% Standardized residual / error term
\newcommand{\StdNormal}{\epsilon}
% Leverage parameter in apARCH
\newcommand{\APARCHleverage}{\zeta}
% Power parameter in apARCH
\newcommand{\APARCHpower}{\delta}
% Degrees of freedom (Student-t)
\newcommand{\DOF}{\nu}
% GJR-GARCH leverage parameter
\newcommand{\GJRleverage}{\gamma}
% GARCH-X external regressor coefficient
\newcommand{\RegressorCoeff}{\xi}
% Mean (first moment)
\newcommand{\Mean}{\mu}
% Lag order in volatility models
\newcommand{\Lag}{L}
% Dummy variable (indicator in regressions)
\newcommand{\DummyVar}{D}
% --- Portfolio -----------------------------------------------------------
% Portfolio weight (scalar)
\newcommand{\weight}{w}
% Portfolio weight vector
\newcommand{\weights}{\bm{w}}
% Covariance matrix
\newcommand{\CovMatrix}{\bm{\Sigma}}
% Correlation matrix
\newcommand{\CorrMatrix}{\bm{C}}
% Correlation coefficient
\newcommand{\correlation}{\rho}
% Number of assets in portfolio
\newcommand{\NumberAssets}{K}
% Portfolio value
\newcommand{\PortfolioValue}{\vartheta}
% Asset beta (CAPM)
\newcommand{\AssetBeta}{\beta}
% Ledoit-Wolf shrinkage intensity
\newcommand{\ShrinkageIntensity}{\delta}
% Sample covariance matrix
\newcommand{\SampleCov}{\bm{S}}
% --- Time Series ---------------------------------------------------------
% Sample size
\newcommand{\SampleSize}{T}
% Generic count (number of lags, payments, etc.)
\newcommand{\Count}{N}
% Estimation window length
\newcommand{\EstWindow}{W_E}
% Testing window length
\newcommand{\TestWindow}{W_T}
% Stress window length
\newcommand{\StressWindow}{W_S}
% AR coefficient
\newcommand{\ARcoeff}{\phi}
% MA coefficient
\newcommand{\MAcoeff}{\psi}
% --- Probability ---------------------------------------------------------
% Probability level (e.g. 0.01 for 1% VaR)
\newcommand{\probability}{p}
% Cumulative sorted scenario weight
\newcommand{\CumScenarioWeight}{\Omega}
% p-quantile of the P and L distribution
\newcommand{\Quantile}{q}
% Quantile function
\newcommand{\QuantileFunction}{Q}
% Probability density function
\newcommand{\PDF}{f}
% Cumulative distribution function
\newcommand{\CDF}{F}
% Standard normal CDF
\newcommand{\NormalCDF}{\Phi}
% Inverse standard normal (quantile function)
\newcommand{\NormalQuantile}{\Phi^{-1}}
% Standard normal density function
\newcommand{\NormalPDF}{\phi}
% Likelihood function
\newcommand{\lik}{\mathcal{L}}
% Log-likelihood
\newcommand{\LogLikelihood}{\ell}
% Parameter vector (MLE estimation)
\newcommand{\ParamSet}{\theta}
% Parameter space
\newcommand{\ParamSpace}{\Theta}
% Fisher information matrix
\newcommand{\FisherInfo}{\mathcal{I}}
% Statistical power (Type II error rate)
\newcommand{\StatPower}{\beta}
% Significance level (CI coverage complement)
\newcommand{\SignifLevel}{\gamma}
% --- Options -------------------------------------------------------------
% Call option label
\newcommand{\CallOption}{\text{call}}
% Put option label
\newcommand{\PutOption}{\text{put}}
% Strike price
\newcommand{\Strike}{X}
% Risk-free interest rate
\newcommand{\RiskFree}{r_f}
% Option maturity (time to expiration)
\newcommand{\OptionMaturity}{\tau}
% Option delta
\newcommand{\OptionDelta}{\Delta}
% Option gamma
\newcommand{\OptionGamma}{\Gamma}
% Option vega
\newcommand{\Vega}{\mathcal{V}}
% Holding period
\newcommand{\HoldingPeriod}{H}
% Simulation profit/loss (signed; negative is a loss)
\newcommand{\ProfitLoss}{\Pi}
% Futures price
\newcommand{\Futures}{F}
% Hedge ratio
\newcommand{\HedgeRatio}{h}
% Units of basic asset held
\newcommand{\StockHolding}{x^b}
% Units of options held
\newcommand{\OptionHolding}{x^o}
% --- Interest Rates ------------------------------------------------------
% Bond convexity
\newcommand{\Convexity}{C}
% Interest rate / yield
\newcommand{\Yields}{r}
% Key rate change, basis points (scalar tenor component); used as \KeyRateChange_k or \KeyRateChange_{t,k}. EXCEPTION to the no-macro-to-macro-indirection rule, decided 2026-08-09: composed from \Yields on purpose, so a future change to the rate glyph propagates here without a follow-up edit. \Yields is defined earlier in generated output (same interest_rates category, TOML order), so expansion order is safe; verify this still holds if either entry is reordered.
\newcommand{\KeyRateChange}{\Delta \Yields^{\mathrm{bp}}}
% Key rate changes, basis points (vector); used bare or as \KeyRateChanges_t. EXCEPTION to the no-macro-to-macro-indirection rule, decided 2026-08-09: composed from \Yields on purpose, so a future change to the rate glyph propagates here without a follow-up edit. \bm{\Yields} depends on \bm/\boldsymbol correctly expanding a macro argument rather than a bare letter --- confirm this renders correctly the first time either consumer compiles it.
\newcommand{\KeyRateChanges}{\Delta \bm{\Yields}^{\mathrm{bp}}}
% Dollar value of a basis point
\newcommand{\DV}{\text{DV01}}
% Mean reversion speed (interest rate models)
\newcommand{\MeanReversion}{\varkappa}
% Modified duration
\newcommand{\ModDur}{D}
% Macaulay duration
\newcommand{\MacDur}{D_{\text{Mac}}}
% Effective duration
\newcommand{\EffDur}{D_{\text{eff}}}
% DV01 vector
\newcommand{\DVvec}{\bm{d}}
% Long-run mean rate (Vasicek/CIR/Hull-White)
\newcommand{\LongRunRate}{\theta}
% Hull-White time-varying drift
\newcommand{\HullWhiteDrift}{\varphi}
% Cash flow (coupon/principal payment)
\newcommand{\CashFlow}{c}
% --- Extreme Value -------------------------------------------------------
% Tail index (EVT, Pareto-type tails)
\newcommand{\TailIndex}{\iota}
% Shape parameter (xi = 1/iota)
\newcommand{\ShapeParam}{\varsigma}
% Threshold value (EVT)
\newcommand{\Threshold}{u}
% Extremal index (Leadbetter 1983)
\newcommand{\ExtremalIndex}{\theta}
% Multivariate extremal coefficient, theta in [1,d]. theta=1 complete tail dependence; theta=d asymptotic independence. Glyph collides with ExtremalIndex and ParamSet; disambiguated at source by macro name.
\newcommand{\ExtremalCoef}{\theta}
% GEV distribution function
\newcommand{\GEV}{\mathcal{H}}
% GPD distribution function
\newcommand{\GPD}{\mathcal{G}}
% GPD scale parameter
\newcommand{\GPDscale}{\beta}
% GPD scale parameter at a threshold u, distinct from the generic GPD scale
\newcommand{\GPDscaleAtThreshold}{\beta_u}
% Number of threshold exceedances (EVT)
\newcommand{\TailCount}{C}
% Sample maximum (block maxima, EVT)
\newcommand{\SampleMax}{M}
% Kurtosis value (realized, e.g. sample kurtosis)
\newcommand{\KurtosisVal}{\kappa}
% Pareto scaling constant (EVT)
\newcommand{\ParetoConst}{\mathcal{A}}
% Generic constant (EVT asymptotics)
\newcommand{\RemainderConst}{\mathcal{C}}
% Little-o asymptotic notation
\newcommand{\LittleO}{o}
% Moment order (EVT)
\newcommand{\MomentOrder}{\mathscr{m}}
% Doubled sample size (EVT block maxima)
\newcommand{\DoubleSample}{D}
% GEV normalizing location constant
\newcommand{\GEVloc}{a}
% GEV normalizing scale constant
\newcommand{\GEVscale}{b}
% Negative return (loss, EVT context)
\newcommand{\NegativeReturn}{L}
% --- Copulas -------------------------------------------------------------
% Copula function
\newcommand{\Copula}{\mathcal{C}}
% Copula density
\newcommand{\CopulaDensity}{c}
% Generator function (Archimedean copulas)
\newcommand{\generator}{\varphi}
% Copula dependence parameter
\newcommand{\CopulaParam}{\theta}
% Copula marginal parameter set
\newcommand{\MarginalParamSet}{\eta}
% Lower tail dependence coefficient
\newcommand{\LowerTailDep}{\lambda_L}
% Upper tail dependence coefficient
\newcommand{\UpperTailDep}{\lambda_U}
% Marginal CDF (second variable)
\newcommand{\MarginalCDF}{G}
% Joint density function
\newcommand{\JointDensity}{h}
% Marginal density function
\newcommand{\MarginalDensity}{g}
% Joint distribution function
\newcommand{\JointCDF}{H}
% Gaussian copula correlation parameter
\newcommand{\CopulaCorr}{\rho}
% --- Multivariate --------------------------------------------------------
% Diagonal matrix of conditional volatilities
\newcommand{\DiagVolD}{\bm{D}}
% Idiosyncratic variance matrix (diagonal, factor models)
\newcommand{\IdioVarMatrix}{\bm{\Psi}}
% Covariance matrix entry (two indices)
\newcommand{\CovElement}{\sigma}
% Auxiliary matrix in DCC dynamics
\newcommand{\DCCauxQ}{\bm{Q}}
% BEKK constant matrix
\newcommand{\BEKKconst}{\bm{\Omega}}
% BEKK ARCH parameter matrix
\newcommand{\BEKKarch}{\bm{A}}
% BEKK GARCH parameter matrix
\newcommand{\BEKKgarch}{\bm{B}}
% DCC rescaling diagonal matrix
\newcommand{\DCCrescaleZ}{\bm{Z}}
% DCC auxiliary matrix element
\newcommand{\DCCelement}{q}
% Factor loadings matrix (loadings/eigenvectors of the covariance matrix)
\newcommand{\FactorLoadings}{\bm{\Lambda}}
% Factor subscript label
\newcommand{\factor}{\text{factor}}
% PCA eigenvalue
\newcommand{\Eigenvalue}{\lambda}
% Number of factors in factor model
\newcommand{\NumberFactors}{m}
% Factor return vector
\newcommand{\FactorReturn}{\bm{g}}
% PCA factor score (scalar; not an eigenvector)
\newcommand{\PCAFactor}{F}
% BEKK cross-covariance coefficient
\newcommand{\BEKKcross}{\delta}
% Factor loading element
\newcommand{\FactorLoading}{\Lambda}
% DCC correlation persistence parameter
\newcommand{\DCCxi}{\xi}
% DCC news coefficient (weight on recent shocks)
\newcommand{\DCCzeta}{\zeta}
% --- Simulation ----------------------------------------------------------
% Number of simulation paths
\newcommand{\NumberSims}{B}
% Cholesky factor of covariance matrix
\newcommand{\Cholesky}{\bm{L}}
% Uniform random number
\newcommand{\UniformDraw}{\tilde{u}}
% --- Backtesting ---------------------------------------------------------
% Violation ratio (observed/expected violations)
\DeclareMathOperator{\ViolRatio}{VR}
% Stressed Value at Risk
\DeclareMathOperator{\StressedVaR}{SVaR}
% QLIKE loss function for variance forecasts
\DeclareMathOperator{\QLIKE}{QLIKE}
% Mean squared error (forecast loss function)
\DeclareMathOperator{\MSE}{MSE}
% Mean absolute error (forecast loss function)
\DeclareMathOperator{\MAE}{MAE}
% Markov transition probability
\newcommand{\TransProb}{\pi}
% Exception indicator: 1 if y_t < -VaR_t
\newcommand{\ExceptionInd}{\eta}
% Exception count
\newcommand{\ExceptionCount}{\upsilon}
% Probability Integral Transform
\newcommand{\ProbIntegralTransform}{\hat{u}}
% Quantile score (pinball loss function)
\newcommand{\QuantileScore}{\mathcal{S}}
% Kolmogorov-Smirnov test statistic
\newcommand{\KSstat}{D}
% Test statistic (generic base letter)
\newcommand{\TestStat}{J}
% --- General -------------------------------------------------------------
% Regression intercept
\newcommand{\RegressionIntercept}{a}
% Regression slope
\newcommand{\RegressionSlope}{\beta}
% Risk factor
\newcommand{\RiskFactor}{x}
% Forecast horizon
\newcommand{\ForecastHorizon}{h}
% Block length (bootstrap)
\newcommand{\BlockLength}{\ell}
% Capital tau variant
\newcommand{\DeliveryTime}{\Upsilon}
% Today's calendar time in years (option pricing input; distinct from trading-date index t)
\newcommand{\CalendarTime}{t^*}
% Probability measure
\newcommand{\ProbMeasure}{\mathbb{P}}
% Risk-neutral measure
\newcommand{\RiskNeutral}{\mathbb{Q}}
% Pricing function (Black-Scholes, bond pricing, etc.)
\newcommand{\PricingFn}{V}
% Indicator function
\newcommand{\Indicator}{\bm{1}}
% Vector of ones
\newcommand{\OnesVector}{\bm{1}}
% Risk factor sensitivity (first-order P&L)
\newcommand{\Sensitivity}{\delta}
% Second-order sensitivity (convexity/gamma)
\newcommand{\ConvexitySens}{\gamma}
% Output floor percentage (Basel)
\newcommand{\OutputFloor}{\alpha}
% --- Subscript Labels ----------------------------------------------------
% Subscript label: annualized
\newcommand{\Annual}{a}
% Subscript label: implied
\newcommand{\Implied}{I}
% Subscript label: portfolio
\newcommand{\Portfolio}{\pi}
% --- Text Abbreviations --------------------------------------------
% S&P 500 index
\newcommand{\SP}{\text{S\&P-500}}
% Student-t distribution (text)
\newcommand{\St}{\text{Student-t}}
% Geopolitical Risk index (Caldara and Iacoviello)
\newcommand{\GPR}{\text{GPR}}
\]
6Getting started
There is no risk forecasting without code. Every method in these notes is implemented three times over, in R, Python and Julia. Work in whichever you know.
6.1 Installation
6.1.1 R and RStudio
To work with R, you need two components:
R, the language and statistical computing engine.
RStudio, an integrated development environment (IDE) that makes working with R easier.
If prompted, install R first — you will be directed to the correct version for your operating system.
Note: If you are using a Windows ARM-based machine, you may need to search separately for the compatible R version, as it is not always the default option.
Once both are installed, launch RStudio. You should see a console window where you can start entering R code immediately.
6.1.2 Positron
Positron is Posit’s newer IDE. R and Python are both first-class runtimes in it, each with its own console, variables pane and plots pane, so it suits work that moves between the two. It is free for personal and academic use.
Install it as you would any other application — on Windows it also needs the Microsoft Visual C++ Redistributable, which most machines already have;
Open the folder holding your work with File > Open Folder. Positron works on a folder rather than a file, and your data files are read relative to it;
Click Start Session at the top right and choose the interpreter you want to use.
The last step matters more than it looks. A machine usually carries several Pythons — one belonging to the operating system, others installed by other software, another from a distribution such as Anaconda. Positron finds all of them and starts the newest unless told otherwise, which need not be the one holding your packages. The interpreter in use is named at the top right whenever a session is running, and Positron remembers the choice for that folder.
Quarto comes with Positron, so rendering to HTML or Word needs nothing further. PDF is the exception, because Quarto produces it through LaTeX and supplies none. If you have no LaTeX distribution already, run quarto install tinytex once in the terminal.
6.1.3 Python
To work with Python, we recommend uv, a fast package and environment manager that can also install Python itself:
Download and install the current stable release for your operating system;
Alternatively, use juliaup for managing Julia versions.
Once installed, launch Julia to access the REPL (interactive prompt) where you can start entering code immediately.
6.1.5 Packages used in this book
Installing R, Python and Julia gives you the languages, not the packages this book’s code imports. The commands below install those packages at whatever version is current when you run them, not the exact versions used to write the book, so results may differ slightly from what is shown here. We do not distribute a locked, version-pinned environment for any of the three languages.
A few chapters that fetch data directly from the EODHD API use additional author-maintained packages not listed here.
Package installation is the last of the setup. The rest of this chapter is about the environment you work in. Chapter 29 collects the learning material, documentation and community links for all three languages.
6.2 User interfaces
The interface matters. Some tools help beginners learn quickly. Others are better for speed, automation or mixed-language work.
6.2.1 RStudio
For most beginners, RStudio is the practical choice. It keeps the editor, console, help and plots in one place.
6.2.2 Posit Cloud
Posit Cloud provides RStudio through a web browser, without any local installation. It runs the same RStudio interface on remote servers.
Posit Cloud helps when you cannot install software on your computer, when you move between machines or when a class or team needs identical R environments. It sets up R, packages and file storage automatically. Projects are easy to share, and you can reach your work from any device with an internet connection.
The free tier covers learning and small projects but has usage limits and needs an internet connection throughout. For extensive analysis or sensitive data, a local RStudio installation is more appropriate.
6.2.3 Quarto
Quarto combines code, output and prose in one file. That is why we use it for these notes. Chapter 17 has more information.
6.2.4 Jupyter
Jupyter notebooks mix code and prose in a cell-based document that runs one cell at a time and exports to HTML or PDF. Each language needs a support package, known as a kernel (IRkernel for R, ipykernel for Python, IJulia for Julia).
Jupyter is convenient for sharing analysis with colleagues who use different languages, since the same interface serves R, Python and Julia, and it suits teaching and presenting analysis step by step. It is less comfortable for long scripts, refactoring and any task that depends on seeing the whole program at once.
6.2.5 VSCode
Visual Studio Code or VSCode is a general-purpose editor with solid support for R and Quarto. It has become a bit bloated and unwieldy, however, and we have moved to Zed for writing these notes.
VSCode supports AI-assisted coding through extensions such as GitHub Copilot. Cursor is a separate editor built on VSCode with AI features throughout.
6.2.6 Zed
Zed is a newer, faster editor than VSCode, with a useful extension for R and AI assistance built in. It does not yet support Quarto documents, so we edit .qmd files in Zed as plain markdown and render them from the command line.
Zed favours performance and simplicity over breadth. Its R support covers syntax highlighting and basic language features, but falls short of RStudio or VSCode, and its extension ecosystem is smaller. We use Zed. Beginners should start with RStudio.
6.2.7 Basic editors
Each language comes with a basic interface:
R: The base R editor gives a simple text-editing window alongside the R console, with basic syntax highlighting.
Python: Python’s IDLE gives a simple interactive shell and editor.
Julia: The Julia REPL (Read-Eval-Print Loop) is an interactive command-line interface.
None of them offers code completion, a plot pane or any notion of a project. Use them when the machine is too old or too locked down to install anything else. Otherwise use RStudio or VSCode.
6.2.8 Command line
All three languages can be accessed from the command line, which is useful for remote servers, high-performance computing clusters or SSH sessions:
R: Type R in a terminal to start the R interpreter.
Python: Type python or python3 to start the Python interpreter.
Julia: Type julia to start the Julia REPL.
The command line draws no graphics, so it stays usable over an SSH connection that would make an IDE unbearable, and it is how most production work runs. The cost is that errors arrive as text, with no plot pane to look at and no variable explorer to poke around in, which makes debugging harder.
6.2.9 Non-interactive scripts
For automated workflows and batch processing, all three languages can run scripts non-interactively:
R: Rscript filename.R
Python: python filename.py
Julia: julia filename.jl
This is how a risk report gets rebuilt at six every morning with nobody present. A script run this way has to be self-contained. It loads its own data, sets its own parameters and writes its output to a file, because there is no one to answer a prompt or notice a warning scrolling past. That is also why error handling matters more here than in interactive work — an unhandled failure at six in the morning is discovered when someone asks for the report.
6.3 Which language
The next three chapters cover the same ground — assignment, vectors, matrices, dictionaries, data frames, packages and source files — once in R, once in Python and once in Julia. Read the one you work in and skip the other two. They are alternatives, not a sequence, and nothing later in the book assumes you have read more than one of them.