REVIEW 5 major objections 3 minor 30 references
Symbolically Regressing Fish Biomass Spectral Data: A Linear Genetic Programming Method with Tunable Primitives
T0 review · 5 major / 3 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read Adding tunable coefficients to genetic programming yields compact models that beat six regressors on fish biomass spectra.
desk verdict A genuine, useful extension of LGP with tunable primitives, but the significance testing on the fish data is anticonservative and needs a cluster-level correction before the empirical claims can be trusted. 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 load-bearing machinery is the tunable primitive: a function $T(X,W,\cdot)$ that takes immediate input values $X$, adjustable coefficients $W$, and, for terminals, two feature indices $\alpha,\beta$, and returns a vector meant to approximate the target $Y$. Three kinds are added to a linear genetic programming register machine: tunable terminals that fit linear regression or derived spectral statistics (average, standard deviation, fluctuation, slope, peak, valley, peak location) over a contiguous feature range $x_{\alpha:\beta}$; tunable functions that fit residual error of intermediate register values with linear, sine, exponential, or power forms; and a multivariate linear regression function appended at the end that maps several register values to the single output. Coefficients are tuned by least squares for the linear primitives and by normalized gradient descent for the trigonometric, exponential, and power primitives, so genetic search handles symbolic structure while local optimizers set the numbers.
What would settle it
Re-run the comparison with six-fold cross-validation in which each fish's three repeated spectra are explicitly assigned to the same fold, then compare mean test $R^2$ and the Friedman/Wilcoxon ranks. If LGP-TP no longer holds the best overall rank or the significant win counts reported in Fig. 2, the empirical claim fails. A second check is to run LGP-TP without the final multivariate linear regression function and see whether the mean-rank gain from LGP-TTF to LGP-TP disappears.
Extended reading notes
Core claim
The paper's central claim is that giving linear genetic programming three kinds of tunable primitives—tunable terminals that fit linear or shape-based features over a wavenumber range, tunable functions that fit residual error at intermediate program points, and a final multivariate linear regression over registers—lets evolution synthesize compact symbolic regressors that predict ten fish biomass components from Raman spectra better than PLSR, KNN, MLP, XGBoost, random forests, and FishCNN. On the paper's own terms, LGP-TP has the best overall mean rank (5.45) in ten runs of six-fold cross-validation, is statistically better than each compared method on at least four of the ten targets, and grows consistently in performance as the three primitive types are added one by one. The same setup transfers to three other spectral preprocessing treatments and, after removing spectral-specific terminals, ranks second in test $R^2$ among the 24 benchmark methods on the SRBench symbolic-regression benchmark while maintaining among the smallest model sizes and shortest training times.
Load-bearing premise
The reported test scores assume that the six cross-validation folds are split by fish rather than by individual measurement, so that spectra from the same fish never appear in both training and test sets; the paper never explicitly states this.
Editorial extensions
If this is right
- Fish production lines could estimate water, protein, fat, and seven other biomass components from Raman spectra with compact formulas instead of black-box models, across several spectral preprocessing protocols.
- The frequency heat maps identify the first three and the last peaks of the 1891.58–580.109 cm$^{-1}$ range as the informative regions, giving chemists specific wavenumbers to investigate.
- The ablation result implies that each of the three tunable-primitive types contributes to the gain, with the final multivariate linear regression function producing the largest jump in mean rank.
- The SRBench result suggests the tunable-primitive strategy generalizes beyond spectroscopy to tabular symbolic regression with competitive accuracy, smaller models, and faster training.
- Because the synthesized programs are compact and their coefficients are explicit, a domain scientist can read the model and reason about whether the highlighted features make chemical sense.
Reading between the lines
- Editorial inference: the performance advantage on small noisy datasets may come less from the evolutionary search itself than from the built-in least-squares coefficient fitting, which gives every candidate program a strong local head start; a testable extension would compare LGP-TP against a version with only random coefficients.
- Editorial inference: the paper never explicitly states that the six-fold split groups the three repeated measurements of each fish into the same fold, and the reported test $R^2$ would be optimistically biased if spectra from the same fish leaked across training and test folds; reporting the exact split procedure would settle this.
- Editorial inference: the same tunable-primitive recipe could be dropped into other high-dimensional, small-sample spectral domains such as soil, food, or pharmaceutical analysis, where interpretability is as valuable as accuracy.
- Editorial inference: the shared heat-map pattern across water, protein, and lipids yield suggests one multi-output model might exploit the same spectral peaks for all three targets; an explicit check is whether such a model retains the per-target accuracy of LGP-TP.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes LGP-TP, a linear genetic programming method with three kinds of tunable primitives (tunable terminals, tunable functions, and a multivariate linear regression function) for symbolic regression on fish biomass Raman spectra. The method is evaluated on ten biomass targets from 117 spectra of 39 fish using ten repeated six-fold cross-validation runs, and compared against six baselines, reporting test R2 values. The paper claims LGP-TP achieves the best overall prediction performance, with compact and interpretable models, and that the tunable primitives are responsible for significant improvements over basic LGP. Generality is further tested on three additional spectral data treatments and on the SRBench symbolic regression benchmark.
Significance. If the central claim holds, LGP-TP offers a practical, interpretable alternative for spectral regression with small, noisy datasets, and the tunable-primitive idea could transfer to other symbolic regression tasks. Strengths of the paper include the use of held-out test folds, a real-world spectral dataset, an ablation study, and an extension to a public benchmark with 122 problems. The authors provide a clear description of the algorithm and make the fish data available on reasonable request. However, the statistical inference supporting the superiority claims is weakened by pseudoreplication in the repeated cross-validation protocol, and several details of the method and comparison require clarification.
major comments (5)
- [Section 2.3] The Wilcoxon and Friedman tests in Section 2.3 treat the ten repeated six-fold cross-validation runs as ten independent observations. Because each run re-splits the same 39 fish (117 spectra), and the three measurements of each fish are kept together as implied by Section 2.1, the same fish appears in test folds multiple times across runs, so the R2 values entering the tests are correlated. This pseudoreplication makes the reported p-values (e.g., the ablation p=0.019 and the per-target p<0.0162 thresholds) anticonservative. The authors should re-analyze with a clustering-aware procedure (e.g., cluster bootstrap or permutation at the fish level), or justify why treating the ten runs as independent is valid.
- [Section 2.1] The cross-validation splitting is not fully specified: the text says each fold has approximately 20 instances (6 or 7 fish samples) but never states explicitly that the three repeated measurements of each fish are always placed in the same fold. If instances are split at the measurement level, spectra from the same fish would leak across training and test folds and inflate test R2. Please state explicitly that no fish appears in both training and test sets in any fold and describe how the folds are constructed.
- [Section 2.3] The description of the statistical adjustments is ambiguous: the sentence 'Friedman's test with a significance level of 0.05 and with the Bonferroni correction and the Wilcoxon rank-sum test with a significance level of 0.05' does not say which family of comparisons the Bonferroni correction applies to, and the threshold p<0.0162 used for per-target comparisons is not derived. Please specify the number of comparisons, the method of correction, and whether the pairwise Wilcoxon tests are adjusted across the ten targets, across the seven methods, or both.
- [Section 2.2, Table 1] The baseline methods KNN, MLP, XGB, and RF use default hyperparameters, whereas PLSR and FishCNN use recommended settings. Because the central claim is that LGP-TP achieves the best overall prediction performance, a comparison using untuned defaults for four of the six baselines may not reflect their true capabilities. The authors should either tune all baselines on the training folds (with nested cross-validation) or clearly state that the comparison is against default configurations and discuss the risk that tuned baselines would close the gap.
- [Section 4.2.2] The derivatives for the tunable functions contain an error: for ExpoRF, the derivative with respect to omega_1 includes the term (omega_1^2 + 1) x_j^{omega_1 - 1}, which is not the derivative of omega_0 + (omega_1^2 + 1) x_j + x_j with respect to omega_1; the correct derivative is 2 omega_1 x_j. Since the gradient descent update is used to tune these primitives, the method as described is not reproducible. Please correct the equation and, if the implementation differs, state the exact update used.
minor comments (3)
- [Section 2.1] The phrase 'wavenumber range of 1891.58-580.109 (cm−1)' lists the upper bound first, which is easy to misread; consider reordering to 580.109-1891.58 (cm−1). Also, Section 2.6 contains a typo '1064num' that should read '1064nm'.
- [Section 4.2.1] The definition of Fluctuate has a typo: it reads sum_{i=alpha+1}^{beta} |(x_j)_i - (x_j)_i|, likely intended as |(x_j)_i - (x_j)_{i-1}|, since otherwise the term is identically zero. Please correct the formula.
- [Section 2.7] The model size for LGP-TP in the SRBench comparison is computed by multiplying the number of effective instructions by 2.0, but it is not clear how this factor is chosen or how it aligns with the model-size metric used for other benchmark methods. Please clarify to make the Pareto plot in Figure 5b interpretable.
Circularity Check
No circularity: the central claim rests on out-of-sample cross-validation and an independent ablation, not on definitional equivalence or self-citation.
full rationale
The paper's central claim is empirical: LGP-TP predicts ten fish biomass targets and is compared with six baselines under six-fold cross-validation. All tunable-primitive parameters are fitted on training folds (Eqs. 1-3 and the least-squares solution for MVLR in Section 4.2) and the reported R2 values are computed on held-out test folds. No test target value is used to fit any coefficient, so the reported predictions do not reduce by construction to a fitted value. The ablation study compares basic LGP with crippled variants and LGP-TP under the same CV protocol; this is a genuine controlled comparison rather than a renamed fit. The paper does cite prior work by the same authors (Refs. 13, 16, 18), but those citations are used for contextual statements such as LGP having promising learning ability or for parameter settings; they are not load-bearing for the superiority claim, which is supported by the experiments. The reader-raised concern about fold grouping and the skeptic's pseudoreplication point are statistical validity or correctness risks, not circularity: even if the p-values are anticonservative, the test R2 values themselves are still out-of-sample estimates, and the claimed reduction of a prediction to its own input is absent. Therefore no circular step is identified.
Assumptions & free parameters
free parameters (6)
- LGP hyperparameters (fish data) =
popsize=250, n_gen=100, max_programsize=50, n_registers=30
- LGP hyperparameters (SRBench) =
popsize=500, n_gen=200, max_programsize=50, n_registers=8
- Gradient descent step size lambda =
0.1
- Tunable parameter bounds =
[-3, 3]
- Feature window cap =
50% of features for fish, 10% for SRBench
- Fitted primitive coefficients W =
varies per individual and run
assumptions (4)
- domain assumption The design matrices in the closed-form least-squares solutions (Eqs. 2-4) are full rank and numerically invertible after data augmentation.
- domain assumption The six-fold split keeps all three measurements of a fish in the same fold, so training and test sets do not share fish.
- domain assumption Data augmentation (spectral augmentation, mixup, Gaussian noise) preserves the true spectrum-to-biomass relationship.
- ad hoc to paper Standard LGP evolutionary search with the given primitive set and program-size limits can reach competitive models within 100 generations.
Cite this review
Pith. "Pith review of Symbolically Regressing Fish Biomass Spectral Data: A Linear Genetic Programming Method with Tunable Primitives." pith.science (2026). https://pith.science/paper/VD2FQVPU
@misc{pith2026250521901,
author = {Pith},
title = {Pith review of: Symbolically Regressing Fish Biomass Spectral Data: A Linear Genetic Programming Method with Tunable Primitives},
year = {2026},
howpublished = {\url{https://pith.science/paper/VD2FQVPU}},
note = {Machine review of arXiv:2505.21901}
}
read the original abstract
Machine learning techniques play an important role in analyzing spectral data. The spectral data of fish biomass is useful in fish production, as it carries many important chemistry properties of fish meat. However, it is challenging for existing machine learning techniques to comprehensively discover hidden patterns from fish biomass spectral data since the spectral data often have a lot of noises while the training data are quite limited. To better analyze fish biomass spectral data, this paper models it as a symbolic regression problem and solves it by a linear genetic programming method with newly proposed tunable primitives. In the symbolic regression problem, linear genetic programming automatically synthesizes regression models based on the given primitives and training data. The tunable primitives further improve the approximation ability of the regression models by tuning their inherent coefficients. Our empirical results over ten fish biomass targets show that the proposed method improves the overall performance of fish biomass composition prediction. The synthesized regression models are compact and have good interpretability, which allow us to highlight useful features over the spectrum. Our further investigation also verifies the good generality of the proposed method across various spectral data treatments and other symbolic regression problems.
Figures
Figures from the paper (3 more)
Reference graph
Works this paper leans on
-
[4]
Gade, A. et al. In-beam spectroscopy reveals competing nuclear shapes in the rare isotope 62Cr. Nat. Phys. (2024)
work page 2024
-
[11]
Ahmmed, F., Gordon, K. C., Killeen, D. P. & Fraser-Miller, S. J. Detection and Quantification of Adulteration in Krill Oil with Raman and Infrared Spectroscopic Methods. Molecules 28, 3695 (2023)
work page 2023
-
[1]
The State of World Fisheries and Aquaculture 2024 – Blue Transformation in action
FAO. The State of World Fisheries and Aquaculture 2024 – Blue Transformation in action. Tech. Rep., FAO, Rome, Italy (2024)
work page 2024
- [2]
-
[3]
Song, H. et al. Multimodal separation and cross fusion network based on Raman spectroscopy and FTIR spectroscopy for diagnosis of thyroid malignant tumor metastasis. Sci. Reports 14, 29125 (2024)
work page 2024
-
[5]
Usman, J. & Begum, S. SOC stocks prediction on the basis of spatial and temporal variation in soil properties by using partial least square regression. Sci. Reports 13, 7949 (2023)
work page 2023
-
[6]
Huang, L. et al. Rapid, label-free histopathological diagnosis of liver cancer based on Raman spectroscopy and deep learning. Nat. Commun. 14, 48 (2023)
work page 2023
-
[7]
Machine Learning for Raman Spectroscopy-based Cyber-Marine Fish Biochemical Composition Analysis
Zhou, Y .et al. Machine Learning for Raman Spectroscopy-based Cyber-Marine Fish Biochemical Composition Analysis. In Proceedings of International Conference on Neural Information Processing, 1–1 (2024)
work page 2024
Show all 30 references
-
[8]
H., Xue, B., Zhang, M
Wood, J., Nguyen, B. H., Xue, B., Zhang, M. & Killeen, D. Automated Fish Classification Using Unprocessed Fatty Acid Chromatographic Data: A Machine Learning Approach. In Proceedings of Advances in Artificial Intelligence, vol. 13728, 516–529 (Springer International Publishing...
2022
-
[9]
Herrero, A. M. Raman spectroscopy a promising technique for quality assessment of meat and fish: A review. Food Chem. 107, 1642–1651 (2008)
2008
-
[10]
Qu, C. et al. Raman spectroscopy for rapid fingerprint analysis of meat quality and security: Principles, progress and prospects. Food Res. Int. 161, 111805 (2022)
2022
-
[12]
& Banzhaf, W
Brameier, M. & Banzhaf, W. Linear Genetic Programming (Springer US, New York, NY , 2007)
2007
-
[13]
& Zhong, J
Huang, Z., Mei, Y . & Zhong, J. Semantic Linear Genetic Programming for Symbolic Regression.IEEE Transactions on Cybern. 54, 1321 – 1334 (2022). 13/15
2022
-
[14]
Linear Genetic Programming for Multi-class Classification Problems
Fogelberg, C. Linear Genetic Programming for Multi-class Classification Problems. Ph.D. thesis, Victoria University of Wellington (2005)
2005
-
[15]
Sotto, L. F. D. P., Rothlauf, F., de Melo, V . V . & Basgalupp, M. P. An Analysis of the Influence of Noneffective Instructions in Linear Genetic Programming. Evol. Comput. 30, 51–74 (2022)
2022
-
[16]
& Zhang, M
Huang, Z., Mei, Y ., Zhang, F. & Zhang, M. Toward Evolving Dispatching Rules With Flow Control Operations By Grammar-Guided Linear Genetic Programming. IEEE Transactions on Evol. Comput. 1–15 (2024). Early access
2024
-
[17]
& Ochoa, G
Banzhaf, W., Hu, T. & Ochoa, G. How the Combinatorics of Neutral Spaces Leads Genetic Programming to Discover Simple Solutions. In Genetic Programming Theory and Practice XX, 65–86 (Springer Nature Singapore, Singapore, 2024)
2024
-
[18]
& Zhang, M
Mei, Y ., Chen, Q., Lensen, A., Xue, B. & Zhang, M. Explainable Artificial Intelligence by Genetic Programming: A Survey. IEEE Transactions on Evol. Comput. 27, 621 – 641 (2022)
2022
-
[19]
& Zhang, M
Huang, Z., Mei, Y ., Zhang, F. & Zhang, M. Multitask Linear Genetic Programming with Shared Individuals and its Application to Dynamic Job Shop Scheduling. IEEE Transactions on Evol. Comput. 1–15 (2023). Early access
2023
-
[20]
J., Glahder, M
Bjerrum, E. J., Glahder, M. & Skov, T. Data Augmentation of Spectral Data for Convolutional Neural Network (CNN) Based Deep Chemometrics (2017). ArXiv:1710.01927 [cs]
2017 arXiv
-
[21]
Zhang, H., Cisse, M., Dauphin, Y . N. & Lopez-Paz, D. mixup: Beyond Empirical Risk Minimization (2018). ArXiv:1710.09412 [cs]
2018 arXiv
-
[22]
& Eriksson, L
Wold, S., Sjöström, M. & Eriksson, L. PLS-regression: a basic tool of chemometrics. Chemom. Intell. Lab. Syst. 58, 109–130 (2001)
2001
-
[23]
De Franca, F. O. et al. SRBench++: Principled Benchmarking of Symbolic Regression With Domain-Expert Interpretation. IEEE Transactions on Evol. Comput. 1–1 (2024). Early access
2024
-
[24]
ECJ then and now
Luke, S. ECJ then and now. In Proceedings of the Genetic and Evolutionary Computation Conference Companion, 1223 – 1230 (2017)
2017
-
[25]
& Zhang, H
Zhang, H., Zhou, A., Qian, H. & Zhang, H. PS-Tree: A piecewise symbolic regression tree. Swarm Evol. Comput. 71, 101061 (2022)
2022
-
[26]
Petersen, B. K. et al. Deep symbolic regression: Recovering mathematical expressions from data via risk-seeking policy gradients. In Proceedings of International Conference of Learning Representation (2021)
2021
-
[27]
Banzhaf, W., Nordin, P., Keller, R. E. & Francone, F. D. Genetic Programming: An Introduction On the Automatic Evolution of Computer Programs and Its Applications (Morgan Kaufmann, San Francisco, California, 1998)
1998
-
[28]
& Manzoni, L
Rovito, L., De Lorenzo, A. & Manzoni, L. Discovering Non-Linear Boolean Functions by Evolving Walsh Transforms with Genetic Programming. Algorithms 16, 499 (2023)
2023
-
[29]
Koza, J. R. Genetic Programming : On the Programming of Computers By Means of Natural Selection (Cambridge, MA, USA: MIT Press, 1992)
1992
-
[30]
& Punch, W
Topchy, A. & Punch, W. F. Faster Genetic Programming based on Local Gradient Search of Numeric Leaf Values.Proc. Genet. Evol. Comput. Conf. 155–162 (2001). Author contributions statement Z. H. proposed the algorithm in this manuscript, implemented the experiment, and wrote the...
2001
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.