REVIEW 4 major objections 5 minor 46 references
Forecasting VIX using interpretable Kolmogorov-Arnold networks
T0 review · 4 major / 5 minor · reviewed 2026-08-09 · deepseek-v4-flash
Pith's one-line read The paper claims that a minimal Kolmogorov-Arnold network turns VIX forecasting into a closed-form linear equation that matches much larger neural networks and exposes mean reversion and the leverage effect.
desk verdict A competent empirical KAN-for-VIX study with believable forecasting results, but the mean-reversion "insight" is a relabeling of the fitted line, not a finding. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The object that carries the argument is the Kolmogorov-Arnold network, a neural network whose trainable parameters are univariate spline activation functions on the edges rather than linear weights, with the Kolmogorov-Arnold representation theorem as its theoretical basis. Here it takes a two-layer form with two hidden nodes, trained on squared error, then pruned and symbolified—the learned splines are replaced by the closest elementary functions, here $y=x$ and $y=0$, followed by affine fine-tuning. The mechanism that matters is that the trained splines are nearly linear, so the whole network collapses into a closed-form linear equation; pruning also drops irrelevant inputs such as $V_{t-21}$ and the quarterly average $V_q$, turning the model into an interpretable variable-selection device.
What would settle it
Take the Dataset 3 closed-form expression for any period, compute the residual $\hat{\varepsilon}_t = \Delta V_t - \hat{\kappa}(\hat{\theta} - V_{t-1})$ from the fitted coefficients, then regress $\hat{\varepsilon}_t$ on $V_{t-1}$ or bin $V_{t-1}$ and test the mean within each bin; if the conditional mean of $\hat{\varepsilon}_t$ is not zero across levels of $V_{t-1}$, the mean-reversion interpretation in equation (3.4) does not follow from the fitted model.
Extended reading notes
Core claim
The paper's central claim is that a deliberately minimal KAN is a competitive, interpretable VIX forecaster. Using only lagged VIX values (five consecutive days, four non-uniform lags, or yesterday's VIX plus weekly/monthly/quarterly averages), a two-layer KAN with two hidden nodes is trained out-of-sample; the B-spline activations consistently train to near-linear shapes, and after pruning and symbolification the network reduces to explicit formulas such as $\hat{V}_t = 0.8584 V_{t-1} + 0.0206 V_{t-2} + 0.0632 V_{t-3} - 0.0374 V_{t-4} + 0.0801 V_{t-5} + 0.2793$. These formulas match the forecasting accuracy of MLPs with 20,901 parameters and LSTMs with roughly 63,500 to 63,800 parameters while using only 60 to 72 parameters, and symbolification improves $R^2$ in eight of nine setups. The paper interprets the Dataset 3 formulas as a typical mean-reverting process and, after augmenting the model with the lagged S&P 500 excess return, as evidence of the leverage effect.
Load-bearing premise
The load-bearing premise is that the leftover term $\varepsilon_t$ in equation (3.4) can be treated as zero-mean noise even though it is allowed to depend on $V_{t-1}$; if that premise fails, the mean-reversion reading of the fitted equation is just a relabeling.
Editorial extensions
If this is right
- VIX forecasting can be delivered as a transparent equation rather than a black-box prediction, with explicit coefficients for each lag and average.
- A very small KAN (60 to 72 parameters) is enough to reach the accuracy of MLPs and LSTMs that are hundreds to thousands of times larger, so interpretability does not have to come at the cost of accuracy.
- Symbolification slightly improves out-of-sample $R^2$ in eight of the nine data/period combinations, suggesting that parsimonious linear structure helps in noisy financial time series.
- The lagged S&P 500 excess return systematically enters with a negative coefficient except in one anomalous setup, giving a learned closed-form expression of the leverage effect.
- For one-day-ahead forecasts, very long-horizon averages are pruned, so the useful memory lies in the most recent VIX level and short/medium-term averages rather than monthly or quarterly averages.
Reading between the lines
- Because every trained spline is near-linear, the KAN is effectively discovering a linear model in this application; the practical gain is data-driven lag selection and a closed-form output rather than a nonlinear functional form.
- A direct test of whether the residual in equation (3.4) has zero mean conditional on yesterday's VIX would turn the mean-reversion reading from a rearrangement of the fitted line into an empirical finding.
- The same pruning-and-symbolification recipe could be applied to other persistent volatility series, where a natural check is whether the linear activations and the specific pruned horizons recur or whether nonlinear forms appear.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper applies two-layer Kolmogorov-Arnold Networks (KANs) to one-day-ahead VIX forecasting, using three lag/aggregate datasets and three train–validation–test splits. After pruning and symbolification, the KAN collapses to closed-form linear forecasting equations. The paper reports that these KANs achieve accuracy comparable to or slightly better than MLP, LSTM, ARMA/ARIMA, and HAR baselines while using far fewer parameters, and it claims that the closed forms reveal two interpretable facts about the VIX: mean reversion and the leverage effect.
Significance. If the empirical claims hold, the paper is a useful contribution to interpretable financial machine learning: it demonstrates that a very small symbolic KAN can match the out-of-sample accuracy of much larger neural networks on a standard volatility-forecasting task, and the closed-form output is genuinely transparent. The paper's strengths include a systematic comparison across nine configurations, an honest report that symbolification usually slightly improves R-squared, and a candid acknowledgment that some symbolification replacements for the leverage effect are not strictly accurate. However, the mean-reversion interpretation in Section 3.2.2 is not established by the derivation given, the forecasting comparisons lack any uncertainty quantification, and a few coefficient inconsistencies need correction. These issues affect the paper's advertised interpretability contribution rather than the raw forecasting tables.
major comments (4)
- [3.2.2, Eqs. (3.3)–(3.4), Table 2] The mean-reversion interpretation is not supported. First, Eq. (3.3) uses coefficients (0.8290, 0.1472, 0.4866) that do not match the D3P1 closed form in Table 2 (0.8297, 0.1477, 0.4756), so the calculation is based on an expression that is not the fitted model. Second, defining epsilon_t as whatever remains after extracting kappa(theta - V_{t-1}) makes Eq. (3.4) a relabeling: the epsilon_t in this case is -0.0238V_{t-1}+0.4866, whose conditional expectation depends on V_{t-1}, so it is not a zero-mean innovation in any standard sense. The claim that the KAN-derived equation 'can be interpreted as representing a typical mean-reverting process' therefore does not follow. Since the abstract and conclusion advertise mean reversion as a key interpretable insight, this point should be either substantiated (e.g., by testing E[epsilon_t | V_{t-1}] = 0 and showing the decomposition is not purely arbitrary) or removed from the claims.
- [3.2.5, Tables 5–7] All forecasting comparisons are point estimates without error bars, confidence intervals, or significance tests. The differences between KAN and the HAR or forward-filling baselines are very small (for example, Period 3 MSE is 1.9575 for KAN-D3 versus 1.9579 for HAR(4), and Period 1 MSE is 4.3957 for KAN-D3 versus 4.4337 for HAR(3)), so the statement that KAN 'tends to excel in MSE' is not statistically supported. The paper should report standard errors, Diebold-Mariano tests, or at least multiple random seeds with mean +/- std. In addition, the '# params' column should state exactly which parameters are counted; if it counts only the affine parameters after symbolification, the comparison with MLP/LSTM parameter counts is not a like-for-like measure of model complexity during training.
- [3.2.1, ARMA/ARIMA order selection] The order selection for ARMA and ARIMA is described only as 'the optimal values of p and q are selected based on AIC' and 'd is determined using the KPSS test,' without specifying whether this is done on the training split of each period or on the full sample. If the same ARMA(1,1) and ARIMA(1,1,1) specifications are used for all three periods after a single in-sample fit, those baselines are not strictly out-of-sample, and the comparison in Tables 5–7 is unfair to the neural networks, which are refitted per period. Please clarify the selection procedure and, if necessary, redo the baseline estimation within each training window.
- [3.2.4, Table 4] The quantitative leverage-effect coefficients in Table 4 are obtained by manually replacing the learned activation function for R^e_{t-1} with y = -x, which the authors acknowledge is not strictly accurate. The reported coefficients therefore describe a constrained auxiliary model, not the fitted KAN. The visual downward slope in Figure 2 is suggestive, but the paper should either quantify the approximation error of the y = -x replacement or present the leverage-effect conclusion as qualitative only. The current wording ('the coefficients for R^e_{t-1} range from -0.02 to -0.07... this result suggests a nonnegligible negative impact') overstates the evidential weight.
minor comments (5)
- [Abstract and Section 3.2.2] There are several typos and spacing errors: 'V olatility' in the abstract, 'Moverover' and 'the the value' in Section 3.2.2, and the repeated 'a ffine' should all be corrected to 'affine'.
- [2.2] The layer definition uses q for the input index and p for the output index, which is the reverse of the convention used in the original KAN paper [26]; consider aligning the notation to avoid confusion.
- [3.2.2] The text says 'the coefficient of Vt-1 is estimated at 0.8290' while Table 2 gives 0.8297 for D3P1; this is part of the inconsistency noted in Major Comment 1 and should be corrected independently.
- [3.2.3] It should be stated explicitly whether the Mincer-Zarnowitz and Durbin-Watson tests are applied to out-of-sample forecast errors; currently this is only implicit.
- [Data availability] The data availability statement is vague ('available from DataGuide... restrictions apply'). Please add a code availability statement or a detailed description of the KAN implementation (library version, hyperparameters, pruning thresholds) to enable reproducibility.
Circularity Check
Mean-reversion insight is a relabeling: Eq. (3.4) is made true by defining ε_t as the V_{t-1}-dependent residual, so this advertised interpretable result is not empirically demonstrated; the forecasting comparisons themselves remain independent.
-
self definitional
[Section 3.2.2, Eqs. (3.3)–(3.4)]
"If we let (−0.0238Vt−1 + 0.4866) in Period 1 be ϵt, then the equation (3.3) becomes identical to (3.4). Consequently, the closed-form expressions derived by KAN with Dataset 3 can be interpreted as representing a typical mean-reverting process reverting to its medium- or long-term average."
Eq. (3.4) is imposed rather than tested: the paper defines ϵt as the leftover term, explicitly allowing ϵt to depend on Vt−1. Since the residual −0.0238Vt−1 + 0.4866 is a function of Vt−1, E[ϵt|Vt−1] ≠ 0, so the conditional expectation of ΔVt is not κ(θ − Vt−1). The actual Vt−1 coefficient in (3.3) is −0.1472 − 0.0238 = −0.171, not κ = 0.1472, and the 'mean level' θ would have to be the time-varying weekly average Vw rather than the constant θ in (3.4). Any persistent linear autoregression of this form can be rewritten as (3.4) by absorbing Vt−1 dependence into ϵt, so the claimed empirical demonstration of mean reversion is a relabeling, not a derivation.
full rationale
The forecasting-performance core of the paper is a self-contained benchmark comparison: KAN is compared against forward filling, ARMA, ARIMA, HAR, MLP, and LSTM on out-of-sample data in Tables 5–7, and the parameter-count comparison is direct. Those claims are not circular. However, the paper's other advertised headline insight, mean reversion, is circular: Eq. (3.4) is satisfied by construction because the paper permits ϵt to depend on Vt−1 and defines it as the residual term in Eq. (3.3). The 'typical mean-reverting process' conclusion is therefore a relabeling of the fitted linear equation rather than an empirical finding. Because this is a central interpretability contribution, though not the forecasting comparison itself, the circularity is partial rather than total.
Assumptions & free parameters
free parameters (14)
- D1P1 regression coefficients =
0.8482*Vt-1 + 0.0536*Vt-2 + ... + 0.3779
- D1P2 regression coefficients =
0.8785*Vt-1 + ... + 0.2345
- D1P3 regression coefficients =
0.8584*Vt-1 + ... + 0.2793
- D2P1 regression coefficients =
0.8661*Vt-1 + 0.0634*Vt-5 + 0.0655*Vt-10 + 0.0712
- D2P2 regression coefficients =
0.8800*Vt-1 + 0.0733*Vt-5 + 0.0359*Vt-10 + 0.1950
- D2P3 regression coefficients =
0.8840*Vt-1 + 0.0742*Vt-5 + 0.0344*Vt-10 + 0.1099
- D3P1 regression coefficients =
0.8297*Vt-1 + 0.1477*Vw + 0.4756
- D3P2 regression coefficients =
0.8189*Vt-1 + 0.1660*Vw + 0.2749
- D3P3 regression coefficients =
0.8255*Vt-1 + 0.1292*Vw + 0.0377*Vm + 0.1025
- Leverage model affine coefficients (Table 4) =
e.g., 0.9974*Vhat - 0.0662*R_{t-1} - 0.0393 for D1P1
- B-spline grid size =
3
- B-spline order =
3
- Regularization coefficient lambda =
0
- Learning rates and early stopping =
0.04 then 0.0004
assumptions (6)
- standard math Kolmogorov-Arnold representation theorem
- standard math B-spline basis functions can approximate the learned activations
- domain assumption The VIX is a function of the chosen explanatory variables plus noise, and out-of-sample evaluation reflects predictive ability
- ad hoc to paper Symbolification with only y=x and y=0 is sufficient to represent the true relationship
- domain assumption The residual from the VIX-only model separates additively into a function of excess returns plus irreducible noise
- ad hoc to paper epsilon_t in Eq. (3.4) can have zero mean while depending on V_{t-1}
Cite this review
Pith. "Pith review of Forecasting VIX using interpretable Kolmogorov-Arnold networks." pith.science (2026). https://pith.science/paper/CHSJT4V5
@misc{pith2026250200980,
author = {Pith},
title = {Pith review of: Forecasting VIX using interpretable Kolmogorov-Arnold networks},
year = {2026},
howpublished = {\url{https://pith.science/paper/CHSJT4V5}},
note = {Machine review of arXiv:2502.00980}
}
read the original abstract
This paper presents the use of Kolmogorov-Arnold Networks (KANs) for forecasting the CBOE Volatility Index (VIX). Unlike traditional MLP-based neural networks that are often criticized for their black-box nature, KAN offers an interpretable approach via learnable spline-based activation functions and symbolification. Based on a parsimonious architecture with symbolic functions, KAN expresses a forecast of the VIX as a closed-form in terms of explanatory variables, and provide interpretable insights into key characteristics of the VIX, including mean reversion and the leverage effect. Through in-depth empirical analysis across multiple datasets and periods, we show that KANs achieve competitive forecasting performance while requiring significantly fewer parameters compared to MLP-based neural network models. Our findings demonstrate the capacity and potential of KAN as an interpretable financial time-series forecasting method.
Figures
Reference graph
Works this paper leans on
-
[1]
LeCun, Y
Y . LeCun, Y . Bengio, G. Hinton, Deep learning, Nature 521 (2015) 436–444
2015
-
[2]
G. Cybenko, Approximation by superpositions of a sigmoidal function, Mathematics of Control, Signals and Systems 2 (1989) 303–314
work page 1989
- [3]
-
[4]
R. Cont, Empirical properties of asset returns: stylized facts and statistical issues, Quantitative Finance 1 (2001) 223
work page 2001
-
[5]
R. T. Baillie, Long memory processes and fractional integration in econometrics, Journal of Econometrics 73 (1996) 5–59
work page 1996
-
[6]
S. J. Taylor, Modelling financial time series, World Scientific, 2008
work page 2008
-
[7]
S. Satchell, J. Knight, Forecasting volatility in the financial markets, Elsevier, 2011
work page 2011
-
[8]
G. E. Box, G. M. Jenkins, G. C. Reinsel, G. M. Ljung, Time series analysis: forecasting and control, John Wiley & Sons, 2015
2015
Show all 46 references
-
[9]
K. R. French, G. W. Schwert, R. F. Stambaugh, Expected stock returns and volatility, Journal of financial Economics 19 (1987) 3–29
1987
-
[10]
Fernandes, M
M. Fernandes, M. C. Medeiros, M. Scharth, Modeling and predicting the CBOE market volatility index, Journal of Banking & Finance 40 (2014) 1–10
2014
-
[11]
Corsi, A simple approximate long-memory model of realized volatility, Journal of Financial Econometrics 7 (2009) 174–196
F. Corsi, A simple approximate long-memory model of realized volatility, Journal of Financial Econometrics 7 (2009) 174–196
2009
-
[12]
U. A. M ¨uller, M. M. Dacorogna, R. D. Dav´e, O. V . Pictet, R. B. Olsen, J. R. Ward, Fractals and intrinsic time-a challenge to econometricians, SSRN, 2008
2008
-
[13]
T. G. Andersen, T. Bollerslev, F. X. Diebold, Roughing it up: Including jump components in the measurement, modeling, and forecasting of return volatility, The Review of Economics and Statistics 89 (2007) 701–720
2007
-
[14]
A. J. Patton, K. Sheppard, Optimal combinations of realised volatility estimators, International Journal of Forecasting 25 (2009) 218–238
2009
-
[15]
Izzeldin, M
M. Izzeldin, M. K. Hassan, V . Pappas, M. Tsionas, Forecasting realised volatility using ARFIMA and HAR models, Quantitative Finance 19 (2019) 1627–1638
2019
-
[16]
Bollerslev, Generalized autoregressive conditional heteroskedasticity, Journal of Econometrics 31 (1986) 307–327
T. Bollerslev, Generalized autoregressive conditional heteroskedasticity, Journal of Econometrics 31 (1986) 307–327
1986
-
[17]
W. Ge, P. Lalbakhsh, L. Isai, A. Lenskiy, H. Suominen, Neural network–based financial volatility forecasting: A systematic review, ACM Computing Surveys (CSUR) 55 (2022) 1–30
2022
-
[18]
Mostafa, T
F. Mostafa, T. Dillon, E. Chang, Computational intelligence approach to capturing the implied volatility, in: Artificial Intelligence in Theory and Practice IV: 4th IFIP TC 12 International Conference on Artificial Intelligence, IFIP AI 2015, Held as Part of WCC 2015, Daejeon,...
2015
-
[19]
Kumar, S
H. Kumar, S. B. Patil, Estimation & forecasting of volatility using arima, arfima and neural network based techniques, in: 2015 IEEE International Advance Computing Conference (IACC), IEEE, 2015, pp. 992–997
2015
-
[20]
S. Pyo, J. Lee, Exploiting the low-risk anomaly using machine learning to enhance the black–litterman framework: Evidence from south korea, Pacific-Basin Finance Journal 51 (2018) 1–12. 17
2018
-
[21]
Petneh ´azi, J
G. Petneh ´azi, J. G ´all, Exploring the predictability of range-based volatility estimators using recurrent neural networks, Intelligent Systems in Accounting, Finance and Management 26 (2019) 109–116
2019
-
[22]
H. Y . Kim, C. H. Won, Forecasting the volatility of stock price index: A hybrid model integrating lstm with multiple garch-type models, Expert Systems with Applications 103 (2018) 25–37
2018
-
[23]
Liu, Novel volatility forecasting using deep learning–long short term memory recurrent neural networks, Expert Systems with Applications 132 (2019) 99–109
Y . Liu, Novel volatility forecasting using deep learning–long short term memory recurrent neural networks, Expert Systems with Applications 132 (2019) 99–109
2019
-
[24]
Lahmiri, A variational mode decompoisition approach for analysis and forecasting of economic and financial time series, Expert Systems with Applications 55 (2016) 268–273
S. Lahmiri, A variational mode decompoisition approach for analysis and forecasting of economic and financial time series, Expert Systems with Applications 55 (2016) 268–273
2016
-
[25]
Litimi, A
H. Litimi, A. BenSaida, L. Belkacem, O. Abdallah, Chaotic behavior in financial market volatility, Journal of Risk, Forthcoming (2018)
2018
-
[26]
Z. Liu, Y . Wang, S. Vaidya, F. Ruehle, J. Halverson, M. Solja ˇci´c, T. Y . Hou, M. Tegmark, KAN: Kolmogorov-Arnold networks, arXiv preprint arXiv:2404.19756 (2024)
2024 arXiv
-
[27]
A. N. Kolmogorov, On the representation of continuous functions of many variables by superposition of continuous functions of one variable and addition, in: Doklady Akademii Nauk, volume 114, Russian Academy of Sciences, 1957, pp. 953–956
1957
-
[28]
A. N. Kolmogorov, On the representation of continuous functions of several variables by superpositions of continuous functions of a smaller number of variables, American Mathematical Society, 1961
1961
-
[29]
Y . Peng, M. He, F. Hu, Z. Mao, X. Huang, J. Ding, Predictive modeling of flexible ehd pumps using kolmogorov-arnold networks, arXiv preprint arXiv:2405.07488 (2024)
2024 arXiv
-
[30]
C. J. Vaca-Rubio, L. Blanco, R. Pereira, M. Caus, Kolmogorov-arnold networks (kans) for time series analysis, arXiv preprint arXiv:2405.08790 (2024)
2024
-
[31]
Alves, F
D. Alves, F. Mendonc ¸a, S. S. Mostafa, F. Morgado-Dias, On the use of kolmogorov–arnold networks for adapting wind numerical weather forecasts with explainability and interpretability: application to madeira international airport, Environmental Research Communications 6 (2024) 105008
2024
-
[32]
K. Xu, L. Chen, S. Wang, Kolmogorov-arnold networks for time series: Bridging predictive power and interpretability, arXiv preprint arXiv:2406.02496 (2024)
2024 arXiv
-
[33]
Inzirillo, R
H. Inzirillo, R. Genet, Sigkan: Signature-weighted kolmogorov-arnold networks for time series, arXiv preprint arXiv:2406.17890 (2024)
2024 arXiv
-
[34]
Genet, H
R. Genet, H. Inzirillo, Tkan: Temporal kolmogorov-arnold networks, arXiv preprint arXiv:2405.07344 (2024)
2024 arXiv
-
[35]
Genet, H
R. Genet, H. Inzirillo, A temporal kolmogorov-arnold transformer for time series forecasting, arXiv preprint arXiv:2406.02486 (2024)
2024 arXiv
-
[36]
R. E. Whaley, Understanding the vix, Journal of Portfolio Management 35 (2009) 98–105
2009
-
[37]
Z. Ding, C. W. Granger, R. F. Engle, A long memory property of stock market returns and a new model, Journal of Empirical Finance 1 (1993) 83–106
1993
-
[38]
Kwiatkowski, P
D. Kwiatkowski, P. C. Phillips, P. Schmidt, Y . Shin, Testing the null hypothesis of stationarity against the alternative of a unit root: How sure are we that economic time series have a unit root?, Journal of Econometrics 54 (1992) 159–178
1992
-
[39]
Goodfellow, Deep learning, 2016
I. Goodfellow, Deep learning, 2016
2016
-
[40]
F. X. Diebold, Elements of forecasting, Citeseer, 1998
1998
-
[41]
J. D. Hamilton, Time series analysis, Princeton university press, 2020
2020
-
[42]
S. L. Heston, A closed-form solution for options with stochastic volatility with applications to bond and currency options, The Review of Financial Studies 6 (1993) 327–343
1993
-
[43]
L. O. Scott, Option pricing when the variance changes randomly: Theory, estimation, and an application, Journal of Financial and Quantitative Analysis 22 (1987) 419–438
1987
-
[44]
J. A. Mincer, V . Zarnowitz, The evaluation of economic forecasts, in: Economic forecasts and expectations: Analysis of forecasting behavior and performance, NBER, 1969, pp. 3–46
1969
-
[45]
Durbin, G
J. Durbin, G. S. Watson, Testing for serial correlation in least squares regression. I, in: Breakthroughs in Statistics: Methodology and Distribution, Springer, 1992, pp. 237–259
1992
-
[46]
Durbin, G
J. Durbin, G. S. Watson, Testing for serial correlation in least squares regression. II, in: Breakthroughs in Statistics: Methodology and Distribution, Springer, 1992, pp. 260–266. 18
1992
Reviewed August 9, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.