While symbolic regression (SR) has been successfully used in science to discover new equations, its use in financial valuation is hindered by several limitations. Whereas the natural sciences provide objectively correct relationships, financial valuation constitutes a distinct class of symbolic discovery problems, as it admits multiple valid perspectives, operates under non-stationary market conditions, and involves noisy, continuous performance signals. In this work, we propose Multi-Agent Fundamental Analysis with Symbolic Adaptive learning (MUFASA), a hierarchical multi-agent framework for symbolic discovery in finance. MUFASA introduces (1) disentangled equation discovery via specialized agents representing distinct valuation perspectives, (2) a meta-coordinator that performs hierarchical-level reasoning over market context information, and (3) a memory mechanism that reasons over statistical performance summaries (e.g., accuracy, stability, and tail risk) to guide learning under noisy feedback. Experiments across datasets from multiple countries show that MUFASA achieves state-of-the-art performance on the valuation task compared to classical finance methods, financial large language models, and SR approaches, while simultaneously producing interpretable equations, which we share with the community. We also make publicly available the distilled learnings across evolution iterations and context-dependent strategy weights, which might offer useful insights for future research on financial fundamental analysis.
Figures & tables
Figure 1: Overview of the challenges in financial SR across different levels of the self-evolving loop.
Figure 2: Overview of the mufasa framework. mufasa uses disentangled multi-agent systems to encourage symbolic expression diversity. A meta-coordinator agent reasons over the non-stationary environment to integrate these expressions. Each agent maintains a memory to identify past statistical patterns.
Category
What is stored
Why
Accuracy
Mean MAPE, Median MAPE
Captures overall predictive accuracy while reducing sensitivity to outliers and skewed observations.
Stability
Standard Deviation, IQR, Coefficient of Variation
Measures whether performance is consistent or highly variable across samples or time periods.
Tail Risk
p90/p95 Error, Worst-Decile Mean
Exposes extreme failure modes that are often hidden by average performance metrics.
Complexity
Variable + Operator Count, Prediction Coverage
Reflects interpretability and numerical reliability of the expression in practical use.
Bias
Mean Signed Percentage Error (MSPE)
Detects systematic over- or under-prediction across observations rather than isolated errors.
Table 1: Statistical information collected by mufasa to enable learning under continuous feedback.
Model
S&P 500
Russell 2000
STOXX 600
FTSE APAC
CSI 300
All Data
Avg Δ
Machine Learning Based Methods
Random Forest
0.4812
0.6551
0.7490
0.8436
0.5605
0.6579
0.1498
HistGBM
0.5916
0.8121
0.9774
1.3771
0.6041
0.8725
0.3644
NeuralNet
0.6130
1.0302
0.8963
1.5292
0.5242
0.9186
0.4105
Financial LLMs
Fino1 ( Qian and others 2025 )
1.0230
2.1244
1.0958
1.2587
1.3388
1.3681
0.8600
Table 2: Performance comparison. The best baselines are underlined, and the best results are bolded.
Model Variant
S&P 500
Russell 2000
STOXX 600
FTSE APAC
CSI 300
S.P Agent
0.4971
0.8994
0.7608
0.6095
0.5928
M.P (Equal Weights)
0.4652
0.8768
0.5448
0.6442
0.5309
M.P (No Stats Memory)
0.4702
0.7525
0.5398
0.5852
0.4852
M.P (Best Single Agent)
0.4676
1.9173
0.5256
0.5833
0.4739
mufasa (Ours)
0.4465
0.5718
0.5064
0.5814
0.4343
Table 3: Ablation study on perspective variant. S.P and M.P denote the single-/multiple perspective, respectively. Best results are bolded, second best are underlined.
Table 9: Learnt equations for the CSI 300 dataset.
Figure 4: Learnt agent weighting across company sectors and market regimes for all datasets.
Model Variant
S&P 500
Russell 2000
STOXX 600
FTSE APAC
CSI 300
Binary Memory (Pass/Fail Signals Only)
Single-Perspective Agent
0.5699
0.7527
0.7220
0.6095
0.5735
Multi-Agent (Equal Weights)
0.5466
1.0459
0.6454
0.7088
0.5599
Multi-Agent (Best Single Agent)
0.5127
1.5160
0.6212
0.6083
0.5761
Multi-Agent (Meta-Coordinator)
0.4646
0.6445
0.5759
0.6086
0.5070
MAPE-Only Memory (No Additional Statistics)
Appendix
Table 10: Extended ablation study. We compare strategies under different memory settings.
Figure 5: Predicted vs. True log prices across Symbolic Regression models on the S&P 500 dataset.
Figure 6: Error distribution across sectors. The center lines show median errors, while the top and bottom edges of the box plots show the 25 th and 75 th percentile errors, respectively. Lower is better.
Model
25 th pct. ( ↓ )
Median ( ↓ )
75 th pct. ( ↓ )
PySR
0.4103
0.6249
0.8322
FunSearch
0.6136
0.7500
0.8283
LLM-SR
0.2663
0.4475
0.6086
mufasa
0.1725
0.3382
0.5556
Appendix
Figure 7: Aggregate APE distributions.
Figure 8: Mean Signed Percentage Error across valuation quartiles (Q1: Cheap to Q4: Expensive).
Department of Computer Science and Artificial Intelligence, Dongguk University-Seoul · Department of Industrial Engineering, Seoul National University · School of Finance, Soongsil University
Department of Information Science, University of North Texas, Denton, TX, United States · Department of Data Science, University of North Texas, Denton, TX, United States