AI Research · Symbolic Regression · Scientific Discovery25 Aug 2026 · Seoul
Comprehensive Survey · 2026/300+ Methods/Four Paradigms/Scientific Discovery

예측에서 방정식으로

A Comprehensive Survey on Symbolic Regression: State-of-the-Art Approaches, Key Applications, Benchmark Evaluations, and Future Research Directions

AuthorsV. Palakonda · S. Ghorbanpour · S. Yun · I.-M. Kim · J.-M. Kang · S. Moon
JournalArchives of Computational Methods in Engineering
Coverage1992–2025 · focus 2020–2025
DOI10.1007/s11831-026-10681-w
Central question

머신러닝이 정답에 가까운 숫자를 내놓는 데 성공했다면, 과학은 그 다음 질문을 던진다. “왜 그런가?” Symbolic Regression은 이 질문에 식으로 답하려 한다. 모델의 내부에 숨어 있는 관계를 사람이 읽을 수 있는 방정식으로 꺼내는 것이다.

Palakonda 등은 2026년 survey에서 300개가 넘는 방법을 deterministic, metaheuristic, neural-symbolic, hybrid의 네 패러다임으로 정리한다. 논문의 진짜 가치는 목록의 길이에 있지 않다. 서로 경쟁해 온 계산철학이 어디에서 만나기 시작했는지를 보여준다는 데 있다. physics-informed constraint, multi-objective optimization, semantic awareness, human collaboration이 패러다임을 가로지르는 공통축으로 떠오른다.

Reading position

이 글은 원 논문의 내용을 한국어로 기계 번역한 요약이 아니다. survey의 taxonomy, 평가표, benchmark, application map, limitations와 future directions를 빠짐없이 따라가되, “Symbolic Regression이 어떻게 자동 과학발견 시스템으로 변하고 있는가”라는 하나의 질문으로 재배열한 해설이다. 개별 방법의 성능 우열은 원 논문이 제시한 조건을 벗어나 일반화하지 않는다.

Part I · The Question

회귀가 아니라 구조를 찾는 문제

Symbolic Regression은 주어진 함수의 계수를 맞추는 문제가 아니라, 함수의 형태와 계수를 함께 찾아야 하는 조합최적화 문제다.

§1 · Definition

숫자를 맞히는 모델과 식을 발견하는 모델

전통적 회귀는 모델의 형태를 사람이 먼저 정한다. 선형회귀라면 \(y=ax+b\)라는 틀을 놓고 \(a,b\)를 찾는다. Symbolic Regression은 틀 자체를 탐색한다. 변수와 연산자를 조합해 함수구조를 만들고, 동시에 상수를 최적화한다.

\[f^*=\arg\min_{f\in\mathcal H}\;\mathcal L(f,\mathcal D)+\lambda\Omega(f)\]

여기서 \(\mathcal H\)는 symbolic expression의 hypothesis space, \(\mathcal L\)은 empirical loss, \(\Omega\)는 structural complexity, \(\lambda\)는 정확도와 간결성 사이의 균형을 뜻한다. 이 식 하나에 Symbolic Regression의 철학이 들어 있다. 잘 맞는 식만 찾는 것이 아니라, 이해할 수 있을 만큼 짧고 구조적으로 의미 있는 식을 찾는다.

원 논문의 Fig. 1은 이 차이를 시각적으로 분명히 한다. black-box model이 prediction을 주는 데 그친다면 Symbolic Regression은 structure discovery, transparency, explainability를 거쳐 explicit equation을 산출한다.

§2 · Why it is hard

식의 언어는 짧지만 탐색공간은 폭발한다

연산자 \(+,-,\times,\div,\sin,\exp\)와 변수 몇 개만 있어도 가능한 expression tree의 수는 깊이에 따라 조합적으로 증가한다. survey는 Symbolic Regression이 기본적인 산술 연산만으로도 NP-hard라는 이론적 결과를 짚는다. 결국 문제는 세 가지를 동시에 풀어야 하는 데 있다.

Structure

어떤 변수와 연산자를 어떤 순서로 연결할 것인가.

Constants

발견한 구조 안의 수치계수를 어떻게 정확히 최적화할 것인가.

Parsimony

오차를 낮추면서 불필요하게 긴 식과 과적합을 어떻게 막을 것인가.

그래서 Symbolic Regression은 단순 regression이라기보다 program search, numerical optimization, model selection, scientific prior가 만나는 문제다.

§3 · Survey scope

500편 이상을 찾고, 300편 이상을 남겼다

저자들은 Scopus, Web of Science, IEEE Xplore, Google Scholar에서 1992–2025년 문헌을 검색하되 2020년 이후 발전에 무게를 둔다. “symbolic regression”, “genetic programming symbolic”, “neural symbolic regression”, “sparse regression interpretable models”, “equation discovery from data” 등을 주요 검색어로 사용했다.

초기 후보는 500편이 넘었고, 방법론·benchmark·application·theory에 실질적인 기여를 한 300편 이상을 상세 검토 대상으로 남겼다. Symbolic Regression을 주변 구성요소로만 사용한 논문은 제외했다.

Contribution 1

300+ methods를 네 패러다임으로 통합하는 taxonomy.

Contribution 2

2020–2025 transformer, RL, LLM, RAG까지 현대적 범위를 포괄.

Contribution 3

50개가 넘는 hybrid method를 integration pattern 관점에서 분석.

Contribution 4

benchmark gap, failure mode, method-selection guideline을 함께 제시.

Part II · Four Paradigms

네 개의 계산철학이 같은 방정식을 찾는다

원 논문의 Fig. 3은 분야를 deterministic, metaheuristic, neural-symbolic, hybrid로 조직한다. 차이는 답의 형태가 아니라 답을 찾아가는 방식에 있다.

§4 · Comparative map

어떤 방법이 언제 유리한가

CriteriaDeterministicMetaheuristicNeural-SymbolicHybrid
Data sizeSmall–MediumAnyLarge preferredAny
Noise toleranceLowMediumMedium–HighHigh
DimensionalityLowLow–MediumHighHigh
InterpretabilityHighHighLow–MediumVaries
Computational costLowMediumHighMedium–High
HardwareCPU onlyCPUGPU preferredCPU/GPU
ScalabilityLimitedModerateHighModerate–High
Prior knowledgeRequiredOptionalNot requiredOptional
ExtrapolationHighMediumLowMedium–High
ImplementationLowMediumHighHigh

원 논문 Table 12의 저자 종합을 옮긴 것이다. 이는 특정 데이터셋에서 모든 구현에 보편적으로 성립하는 성능순위가 아니라 method-family 수준의 실무 가이드다.

§5 · Deterministic

탐색을 통제하고 보장을 얻는다

Deterministic 접근은 수학적 프로그래밍, exhaustive enumeration, sparse regression, structured representation을 이용해 search를 체계화한다. stochastic run마다 답이 흔들리지 않는 reproducibility와 이론적 convergence가 강점이다.

Exhaustive Search

ESR, expression DAG, probabilistic grammar로 가능한 식을 체계적으로 열거하고 중복을 줄인다.

MINLP

expression tree와 constants를 mixed-integer nonlinear programming으로 함께 표현하고 global optimization을 노린다.

Sparse Regression

SISSO 계열처럼 거대한 candidate feature space를 screening과 sparsification으로 압축한다.

Structured Representation

IT, TIR, block decomposition처럼 표현가능한 식의 형태를 제한해 search와 interpretability를 동시에 관리한다.

Bayesian / Physics-Informed

prior와 MCMC 또는 symmetry·separability·dimension 같은 domain structure를 이용한다.

Deterministic Optimization

DoME, DistilSR 등은 semantic representation과 constant optimization, parsimony를 체계화한다.

AI Feynman은 신경망으로 symmetry와 separability를 찾은 뒤 문제를 재귀적으로 분해한다. SISSO는 correlation screening과 sparse operator를 결합해 거대한 descriptor 후보군에서 저차원 표현을 찾는다. 공통점은 분명하다. 무작정 찾는 대신 문제구조를 이용해 search space를 줄인다.

§6 · Metaheuristic

정답을 보장하기보다 넓게 탐색한다

Koza의 Genetic Programming에서 출발한 metaheuristic 계열은 population을 진화시키며 expression tree를 탐색한다. predetermined structure가 필요 없고 다양한 domain에 적용하기 쉽지만, computational cost, run-to-run variance, constant optimization, bloat가 오래된 난제다.

FamilySurvey가 다루는 핵심 축대표적 아이디어
Core GPfitness·complexity, diversity, operators, specialized representationVC-dimension, correlation fitness, semantic operators, GP-GOMEA
Semantic GPexpression syntax가 아니라 output behavior를 search guide로 활용semantic crossover/backpropagation, geometric semantic GP, clustering
Multi-Objective GPaccuracy와 complexity·diversity·constraint를 동시 최적화NSGA-II, MOEA/D, Pareto dominance, shape constraints
Memetic GPglobal evolutionary search + local optimizationNelder–Mead, gradient refinement, Lamarckian/Baldwinian learning
Feature / Dimension고차원 변수 선택과 missing data 대응permutation importance, mRMR, NaN identity, Shapley-guided crossover
Transfer / Multitask관련 task 사이 subtree·feature·instance knowledge 재사용MFGP, TLGP, PST, MvSR, IRMTGP
GEPfixed-length chromosome → variable expression treeuniform design, self-learning functions, dimensional correction
Swarm / Otherbee, firefly, ant colony, immune, distribution, local searchABCP, Firefly Programming, ACO, GRASP, RILS-ROLS
Grammar-Basedformal grammar로 valid expression space를 정의SCFG-LGP, production ranking, dynamic pruning

survey의 metaheuristic 장이 길 수밖에 없는 이유는 GP가 단일 알고리즘이 아니라 하나의 생태계가 되었기 때문이다. selection, representation, crossover, mutation, local search, transfer, semantics, complexity control이 서로 독립적인 연구축을 이룬다.

Part III · The Neural-Symbolic Turn

식의 탐색을 생성 문제로 바꾸다

Transformer와 reinforcement learning은 Symbolic Regression을 “tree를 돌연변이시키는 문제”에서 “expression sequence를 생성하는 문제”로 다시 썼다.

§7 · Seq2Seq & Transformers

숫자 집합에서 expression token으로

Neural-symbolic 계열은 deep learning의 representation power를 사용하되 최종 결과를 symbolic expression으로 돌려준다. Seq2Seq-SR, NeSymReS, end-to-end transformer, SymFormer는 numerical observations를 expression skeleton과 constants로 번역하는 계열이다.

survey가 정리한 흐름은 세 단계로 읽을 수 있다. 먼저 structure와 constant를 분리해 찾았다. 다음에는 transformer가 둘을 함께 예측하기 시작했다. 최근에는 MMSR처럼 numeric data와 symbolic expression을 서로 다른 modality로 보고 contrastive alignment까지 도입한다.

검색 알고리즘이 식을 “조립”하던 시대에서, pretrained model이 식의 문법과 자주 나타나는 구조를 “기억하고 제안”하는 시대로 이동한 것이다.
§8 · Representative neural methods

생성, 동역학, 노이즈, 제약을 각각 다른 방식으로 다룬다

Method familyRepresentative methods in surveyKey idea
Seq2Seq coreSeq2Seq-SR, NeSymReS, End-to-End SR, SymFormer, MMSRskeleton/constant prediction, multimodal alignment
Dynamics extensionsNSODE, ODEFormer, PROSEtrajectory→ODE, operator learning + symbolic generation
Noise robustnessDN-CL, NRSRcontrastive clean/noisy views, noise-gated variable selection
RetrievalRAG-SRsemantic descent + retrieved symbolic tree library
RL generationDSR, uDSR, VSR-DPG, GraphDSR, CaMorisk-seeking policy gradient, modular substructure reuse
Alternative RLRL-GEP, GFN-SR, Deep-RL PDE, Sym-Qmulti-agent sampling, GFlowNet diversity, offline RL
Equation LearnerEQL÷, PEQL, SEQL/HEQLmathematical operators embedded as differentiable network units
Symbolic NNMetaSymNet, SyNISM, DDSR-NN, SymbolNetdynamic topology, sparse symbolic activations, differentiable distributions

survey는 PROSE가 25,600개 ODE system에서 99.9% validity를 보고했다고 정리한다. 이 숫자는 모든 symbolic task에서의 정확도를 뜻하지 않는다. validity는 syntactically/structurally valid output의 비율이라는 특정 문맥의 지표다.

§9 · Physics-Informed & Constrained

데이터에 맞는 식과 과학적으로 가능한 식은 다르다

PINN-SR은 scarce/noisy data에서 governing PDE를 찾기 위해 neural approximation과 sparse regression을 결합한다. Φ-SO는 dimensional analysis를 token generation에 넣어 dimensionally inconsistent expression을 애초에 만들지 않도록 한다. PINN-DPA는 PINN solution을 differentiable program과 pruning을 통해 interpretable expression으로 바꾼다.

\[\text{Numerical Fit}\quad\not\Rightarrow\quad\text{Scientific Validity}\]

survey 전반에서 physics-informed constraint가 반복해서 등장하는 이유다. dimensional homogeneity, conservation law, monotonicity, convexity, boundedness를 search에 넣으면 후보의 수를 줄이는 동시에 scientifically invalid한 식을 거를 수 있다. 물론 그 대가로 domain expert가 prior를 정확히 명시해야 한다.

Part IV · LLM, Retrieval & Interaction

언어모델은 방정식의 저자가 아니라 탐색의 조력자가 된다

2025년 이후 LLM은 functional form을 제안하고, concept library를 만들고, executable program을 생성한다. 동시에 hallucination이라는 새로운 실패모드도 가져온다.

§10 · LLM Integration

수학식도 언어처럼 생성할 수 있는가

survey는 MathLM, LA-DSO, SGA, ICSR, LASR, LLM-SR를 advanced neural paradigm으로 묶는다. MathLM은 Wikipedia에서 추출한 798,998개의 수학표현으로 language prior를 학습한다. 저자들은 이 대규모 corpus의 장점과 함께 notation inconsistency, incomplete expression, well-documented domain에 대한 selection bias도 경고한다.

ApproachLLM/LM의 역할Symbolic search와의 결합
MathLMmathematical expression priorDSR logits와 결합해 candidate generation을 유도
LA-DSOpretrained LM embedding유사 sequence가 가까워지는 latent search space
SGAdiscrete hypothesis generationphysics simulation이 continuous parameter를 최적화
ICSRmeta-prompted functional form proposalnonlinear least squares가 coefficients를 결정
LASRhypothesis + concept library evolutionzero-shot abstraction과 library refinement
LLM-SRequation as executable Python programmulti-island evolutionary strategy로 refinement

2025년 ICLR의 LLM-SR은 LLM이 만든 식을 실행가능한 Python program으로 표현하고 evolutionary island를 이용해 반복개선한다. 같은 해 RAG-SR은 pretraining 없이 retrieval-augmented semantic library를 이용한다. 이 두 흐름은 “큰 모델이 정답식을 외운다”보다 큰 모델이 search proposal을 만들고 symbolic machinery가 검증·수정한다는 방향이 더 지속가능함을 시사한다.

§11 · In-context reliability

예시 하나가 식의 방향을 바꿀 수 있다

survey는 future directions에서 in-context symbolic regression의 취약성을 따로 지적한다. demonstration selection이 생성식에 직접 영향을 주므로 잘못된 예시는 invalid expression이나 symbolic hallucination으로 이어질 수 있다. RDES, Meta-Sel 같은 demonstration-selection 방법이 이 문제를 줄이려는 시도다.

Interpretation

LLM이 Symbolic Regression에 들어오면서 search가 빨라질 수는 있다. 그러나 “그럴듯한 식을 쓰는 능력”과 “관측된 세계의 법칙을 찾아내는 능력”은 같은 것이 아니다. 후자는 dimensional check, constraint verification, held-out extrapolation, uncertainty, 실험 검증을 요구한다.

§12 · Human & Active Learning

전문가를 search loop 밖에 두지 않는다

Interactive-SR은 parameter adjustment, expression prioritization, diagnostics를 사람이 직접 조절하는 UI를 제공한다. RGB2-SR는 pairwise preference와 grammar-constrained suggestion을 search에 넣는다. QUOSR와 APPS는 반대로 시스템이 “어떤 데이터를 다음에 관측해야 가장 많이 배울 수 있는가”를 묻는다.

여기서 Symbolic Regression은 passive curve fitting을 벗어난다. 식 후보를 만들고, 불확실한 영역을 찾아, 필요한 데이터를 다시 요청하는 active scientific discovery loop가 된다.

Part V · Hybrid Convergence

좋은 알고리즘 하나보다 좋은 역할분담이 중요해진다

Hybrid 접근은 evolutionary global search, neural proposal, statistical fitting, formal constraint, numerical optimizer의 장점을 단계별로 배치한다.

§13 · Why hybrid

각 패러다임의 약점을 다른 패러다임으로 메운다

Evolutionary algorithm은 global exploration에 강하지만 constant refinement가 약하다. Neural model은 빠른 proposal을 만들지만 OOD와 formal guarantee가 약하다. Optimization solver는 정확하지만 search space가 커지면 힘들다. Domain constraint는 과학적 타당성을 높이지만 expert specification이 필요하다.

그래서 hybrid system은 역할을 나눈다. neural model이 promising region을 제안하고, evolutionary search가 구조를 변형하고, gradient/least-squares가 coefficient를 맞추고, physics/SMT가 식의 validity를 확인한다.

\[\text{Proposal}\rightarrow\text{Global Search}\rightarrow\text{Constant Refinement}\rightarrow\text{Constraint Verification}\]
§14 · Hybrid landscape

Survey가 포괄하는 주요 통합패턴

Integration patternExamples surveyedPurpose
Evolutionary + StatisticalSequential-SR, MRGP, GP-RVM, MLR-GPfeature engineering과 coefficient estimation 분업
Evolutionary + TransferGPTL, ITGP, multitree TL-GPdomain shift·incomplete data에서 knowledge transfer
Evolutionary + PhysicsPhysics-GP, SciMEDgoverning constraint와 data-driven discovery 결합
Neural + EvolutionDAE-GP, RNN-GP, SRNE, GENetslearned proposal/latent space + global evolutionary refinement
Transformer + GPPIGP, PGGP, SymbolicDPOpretrained initialization과 preference/evolutionary ranking
Multi-Objective / FormalCDSR, Shape-SR, TIRMOOaccuracy–complexity–constraint Pareto와 SMT counterexample
Planning + SearchTPSR, RSRM, SE-MCTSMCTS·RL로 expression construction의 lookahead 강화
Multi-ScaleNEEP-ADF, modular/hierarchical SRfine symbolic dependencies와 reusable modules 동시 학습
§15 · Practical selection

Hybrid가 항상 정답은 아니다

survey의 실무 가이드는 오히려 절제되어 있다. 변수가 5개 미만이고 noise가 적은 저차원 문제라면 single-paradigm method로 충분할 수 있다. dimensionality와 noise가 올라가거나 domain constraint를 강제해야 할 때 hybrid의 이점이 커진다.

  • High-dimensional: feature selection과 expression discovery를 같이 다루는 evolutionary-statistical integration이 유리할 수 있다.
  • Noisy / incomplete: transfer-learning, denoising hybrid가 robust한 선택이 된다.
  • Physics constrained: shape constraint나 SMT verification으로 domain consistency를 직접 확인한다.

반대로 neural-evolutionary hybrid는 search space 자체가 더 커질 수 있고, small data에서 ensemble overfitting이 생길 수 있으며, inter-paradigm hyperparameter가 새로운 민감도를 만든다. 혼합은 복잡성을 없애는 방법이 아니라 복잡성을 더 잘 배분하는 방법이다.

Part VI · Evaluation & Applications

방정식을 찾았다는 주장을 어떻게 믿을 것인가

좋은 Symbolic Regression 논문은 training fit만 보여주지 않는다. extrapolation, complexity, stochastic variance, runtime, reproducibility를 함께 보고해야 한다.

§16 · Software ecosystem

MATLAB toolbox에서 high-performance backend까지

FrameworkYearSurvey가 정리한 특징
Browser-based2007JavaScript interpreter, AJAX, distributed fitness evaluation
GPLab 3.02008bloat control, dynamic size limits, parsimony pressure
Pyevolve2009cross-platform, multiple chromosome representations
GPTIPS2010multigene SR, simplification, LaTeX conversion
DEAP2012modular architecture, distributed task management
GPOCL2013OpenCL acceleration, massive node evaluation throughput
HeuristicLab2014plugin architecture, visual algorithm design, interactive analysis
PonyGE22017grammatical evolution, genome/derivation-tree dual representation
Glyph2019client–server architecture, hardware integration
Parallel Computing2020MPI, GPU-CUDA, Apache Spark comparison
Operon2020compact linear tree encoding, automatic differentiation
Bingo2022component interchangeability, backend integration
PySR2023high-performance backend, multi-population evolution, runtime kernel fusion

도구의 역사는 연구문제의 변화를 그대로 보여준다. 초기에는 GP를 “돌릴 수 있는가”가 문제였다면, 이제는 distributed execution, automatic differentiation, backend modularity, reproducible workflow가 기본요건이 된다.

§17 · Benchmark genealogy

단순 polynomial에서 extrapolation과 noise까지

BenchmarkYear# VariablesComplexityRecommended use
Koza19921Lowbasic GP validation
Keijzer20031–2Low–Medinterval arithmetic / linear scaling
Vladislavleva20081–5Med–Highcomplexity-aware nonlinearity evaluation
Nguyen20101–2Low–Medsemantic GP evaluation
Korns20115Highcomplex ground-truth accuracy testing
White et al.2012benchmark protocol meta-analysis
AI Feynman20201–9Varied100 physics equations rediscovery
SRBench20211–100+Variedstandardized comparison + living leaderboard
iirsBenchmark2022VariedVariedinterpretability: stability / fidelity
SRSD2022Varied3 levelsphysics discovery: Easy / Medium / Hard
SRBench++2024VariedVariedrediscovery / extrapolation / noise robustness

원 논문 Table 14를 바탕으로 재구성했다. benchmark마다 목적이 다르므로 하나의 순위표처럼 사용하는 것은 부적절하다.

§18 · Minimum reporting protocol

“잘 맞았다”만으로는 과학적 발견이 아니다

Prediction. MSE/RMSE와 R² 같은 in-sample 지표뿐 아니라 held-out domain의 extrapolation error를 함께 보고한다.

Complexity. tree size 또는 expression length와 함께 structural, functional, statistical complexity 가운데 적어도 하나를 제시한다.

Stochastic variance. GP 등 stochastic method는 충분한 independent run과 run-to-run spread를 보고한다.

Benchmark fit. basic validation, complexity, dynamical system, extrapolation 등 연구질문에 맞는 benchmark를 선택한다.

Reproducibility. hardware, wall-clock runtime, hyperparameters, random seeds를 빠짐없이 공개한다.

§19 · Applications

식이 필요한 곳은 생각보다 넓다

원 논문의 Fig. 4는 Symbolic Regression 응용을 physical sciences, chemistry, energy, environmental science, healthcare & bioinformatics, information technology, economics & industry로 펼쳐 보인다. 각 분야는 같은 “interpretable equation”을 원하지만 이유는 서로 다르다.

Physical Sciences

constitutive law, material property, particle/nuclear/astrophysics, force law, symmetry와 governing equation 발견.

Chemistry

reaction kinetics, distillation monitoring, process dynamics, Michaelis–Menten, combustion, catalyst descriptor.

Engineering

structural response, fracture, manufacturing, fault detection, controller design, robotic locomotion.

Energy & Environment

wind/solar prediction, power systems, emissions, air quality, carbon modeling, climate/Earth PDE.

Healthcare & Bioinformatics

frailty, cancer biomarkers, EHR, pharmacometric covariates, protein dynamics, single-cell population equations.

Information Technology

neural-network approximation, LLM pruning metric, decision tree, RL reward, graph/network/channel models.

Economics & Industry

GDP/sentiment, commodity forecast, process optimization, drilling, building energy, petroleum correlations.

Scientific Discovery

domain knowledge를 넣을수록 단순 fit이 아니라 mechanism-oriented expression으로 이동한다.

survey의 cross-domain 분석은 중요한 차이를 짚는다. 물리과학은 exact symbolic recovery를 중시하고, engineering은 real-time robustness를, healthcare는 clinical interpretability와 regulatory acceptance를 더 중시한다. 따라서 method selection은 leaderboard보다 domain objective를 따라야 한다.

Part VII · Limits & Future

방정식이 읽힌다고 발견과정까지 투명한 것은 아니다

현재의 난제는 search complexity, OOD generalization, benchmark bias, uncertainty, deployment gap이다. 미래는 더 큰 search가 아니라 더 검증 가능한 discovery loop를 향한다.

§20 · Limitations

다섯 가지 벽

1 · Computation

operator와 depth가 늘수록 search space가 폭발한다. evolutionary는 많은 fitness evaluation을, neural은 pretraining/GPU를 요구한다.

2 · Generalization

noise와 spurious correlation에 과적합하고, bloat가 생기며, synthetic training distribution에서 OOD로 나가면 성능이 무너질 수 있다.

3 · Data & Domain

missing, heterogeneous, high-dimensional real data를 안정적으로 처리하기 어렵고, domain-agnostic fit이 scientifically invalid한 식을 만들 수 있다.

4 · Evaluation

low-dimensional synthetic benchmark 편중, inconsistent variance reporting, hyperparameter omission, standardized UQ/deployment protocol 부재가 남는다.

5 · Deployment

복잡한 evolutionary run은 hours-to-days가 걸릴 수 있고 stochastic spread가 크며, human-in-the-loop tools는 아직 prototype인 경우가 많다.

Mechanism gap

높은 predictive accuracy가 ground-truth structure recovery를 보장하지 않는다. 잘 맞는 식과 맞는 기전은 다르다.

Symbolic Regression의 역설은 여기에 있다. 결과식은 black box보다 읽기 쉽지만, 그 식이 왜 선택되었고 다른 후보가 왜 버려졌는지까지 자동으로 투명해지는 것은 아니다.
§21 · Timeline

저자들이 제시하는 연구개발 시간축

Near term · 1–2 years

advanced benchmark/evaluation protocol과 unified scalable architecture. 이미 존재하는 software·benchmark infrastructure를 정리하고 확장하는 단계다.

Medium term · 3–5 years

multi-objective complexity control, physics/domain-aware SR, hybrid neural-symbolic paradigm, robustness와 uncertainty quantification.

Long term · 5+ years

convergence/sample-complexity/generalization theory, formal analysis, interactive/explainable discovery, 새로운 과학 application과 autonomous scientific discovery.

§22 · Nine future directions

다음 세대 Symbolic Regression이 풀어야 할 아홉 문제

#Direction핵심 내용
1Unified & Scalable Architecturesdistributed training, hardware acceleration, streaming, standardized API/shared representation, software-engineering reproducibility, high-dimensional feature selection
2Multi-Objective Complexity Controlerror + description length + dimensional consistency + stability를 Pareto objective로 공동 최적화
3Physics / Domain Awarenessinvariance, conservation, dimensional analysis, monotonicity, convexity, boundedness, neural proposal + symbolic verification
4Advanced Benchmarkingrealistic noise, missingness, high dimension, complex dynamics, extrapolation, structural correctness, UQ
5Theoretical Foundationsconvergence, sample complexity, generalization bound, SMT, program synthesis, refinement-type verification, stopping criterion
6Hybrid Neural-Symbolictransformer/LM prior + evolutionary refinement, RAG semantic library, grammar-constrained decoding, multiscale modules
7Interactive & Explainable Systemsexpert constraints, operator guidance, candidate/search visualization, uncertainty display, natural-language interfaces
8Robustness & Trustworthinessstructural/coefficient uncertainty, Bayesian posterior, ensemble, noise filtering, distribution-shift test, formal verification, provenance
9Emerging Domainsmulti-output, time-varying systems, federated SR, PDEs, materials/drug/climate workflows, causal/interventional/counterfactual discovery
§23 · Convergence

분야가 수렴하는 세 개의 원칙

300개가 넘는 방법을 읽은 저자들이 마지막에 발견하는 공통점은 의외로 단순하다.

Physics-Informed

무엇이 가능한 식인지 domain law가 search를 제한한다.

Multi-Objective

accuracy 하나가 아니라 simplicity, stability, validity를 동시에 본다.

Semantic / Human

syntax만 탐색하지 않고 behavior, knowledge, expert feedback을 search signal로 쓴다.

Deterministic method의 guarantee, evolutionary method의 robust exploration, neural method의 speed와 prior, hybrid method의 역할분담이 서로 가까워진다. survey는 이를 단순한 알고리즘 유행이 아니라 robust, interpretable, scalable automated scientific discovery system으로의 수렴으로 읽는다.

§24 · What this survey really says

과학의 언어를 다시 기계가 쓰기 시작했다

오래전 과학자는 관측값을 보고 방정식을 썼다. 현대의 머신러닝은 그 사이에서 엄청난 예측력을 얻었지만, 종종 식을 잃었다. Symbolic Regression은 시간을 되돌리는 기술이 아니다. deep learning, optimization, evolution, formal methods를 이용해 예측력과 명시적 구조를 다시 같은 자리에 놓으려는 시도다.

Final synthesis

다음 세대 Symbolic Regression의 경쟁력은 “가장 복잡한 식을 가장 빨리 찾는 능력”보다 “어떤 식을 믿어도 되는지 설명하고, 어디에서 틀릴지 측정하고, 어떤 실험으로 반증할지 제안하는 능력”에서 결정될 가능성이 크다.

그 지점에서 Symbolic Regression은 interpretability tool을 넘어선다. 데이터에서 식을 만들고, 제약으로 검증하고, benchmark로 흔들어 보고, 전문가와 수정하고, 새 데이터를 요청하는 작은 과학방법론의 계산엔진이 된다. 이 survey가 그리는 미래는 자동화된 방정식 생성기가 아니라, 자동화와 해석가능성 사이의 긴장을 관리하는 discovery system이다.

Selected References & Resources

01
Palakonda et al. · Archives of Computational Methods in Engineering · 2026
이 글의 1차 자료. 300+ symbolic regression methods를 네 패러다임으로 통합하고 evaluation, applications, limitations, future directions를 분석한다.
02
Udrescu & Tegmark · Science Advances · 2020
symmetry와 separability를 이용한 physics-inspired deterministic decomposition의 대표 사례.
03
Petersen et al. · ICLR · 2021
risk-seeking policy gradient를 이용한 neural-guided expression generation의 대표 연구.
04
Tenachi et al. · The Astrophysical Journal · 2023
dimensional consistency를 generation constraint에 직접 넣는 physics-informed symbolic regression.
05
Shojaee et al. · ICLR · 2025
LLM-generated executable programs와 multi-island evolution을 결합한 LLM-SR.
06
Zhang et al. · ICLR · 2025
semantic library retrieval과 symbolic descent를 결합한 pretraining-free approach.
07
Han et al. · IEEE Computational Intelligence Magazine · 2025
pretrained transformer와 GP의 역할분담을 다루는 hybrid frontier.
08
Kumar et al. · Industrial & Engineering Chemistry Research · 2025
symbolic network와 parameter network를 결합한 chemical-process oriented neural-symbolic method.
09
Benchmark resource
standardized cross-method comparison, rediscovery, extrapolation, noise robustness를 위한 benchmark ecosystem.
10
Benchmark resource
physics-based equation discovery를 difficulty level에 따라 평가하는 benchmark.
11
Benchmark resource
Feynman Lectures 기반 100 equations를 이용한 physics-grounded rediscovery benchmark.
12
Authors' repository noted in survey
survey는 acceptance 후 complete BibTeX와 machine-readable comparison tables를 이 repository에 제공할 계획이라고 명시한다.