ALL Metrics
-
Views
Get PDF
Get XML
Cite
Export
Track
Research Article

Employing Fixed-Point Theory for Fuzzy Regression Analysis: Methodology and Empirical Application

[version 1; peer review: awaiting peer review]
PUBLISHED 23 Dec 2025
Author details Author details
OPEN PEER REVIEW
REVIEWER STATUS AWAITING PEER REVIEW

This article is included in the Fallujah Multidisciplinary Science and Innovation gateway.

Abstract

Background

Traditional fuzzy regression approaches, such as Tanaka’s fuzzy minimum method and fuzzy least squares, often lack theoretical guarantees of existence, uniqueness, and numerical stability. These limitations are of paramount importance in engineering applications involving uncertainty, such as predicting the compressive strength of concrete. This study addresses these issues by presenting a mathematically rigorous fuzzy regression model, based on fixed-point theory, formulated within the full metric space of trapezoidal fuzzy numbers using the metric scale d ∞ .

Methodology

We define a shrinkage coefficient on the trapezoidal fuzzy coefficient vector space and prove, using Banach’s fixed-point theory, the existence and uniqueness of the regression solution. An iterative algorithm is constructed to estimate the coefficients, using alpha computation and Lipschitz continuity to ensure convergence. The University of California, Irvine concrete compressive strength dataset was amplified using ASTM and ACI-based uncertainty coefficients, and the proposed fixed-point model was evaluated against the Tanaka method and least-squares fuzzy regression. Performance was assessed by mean squared error (MSE), coefficient ambiguity, convergence behavior, and toughness under ±5% noise.

Results

The proposed method demonstrated consistent geometric convergence with an average of 12.3 iterations and zero divergence across all experiments. It also reduced the overall mean squared error by an average of 12.5%, and by up to 25.1% at best, compared to the comparator methods. Coefficient ambiguity—measured by ambiguity width—was reduced by 18.3% compared to the Tanaka method and by 13.3% compared to ambiguity squares. Under noise perturbation, the model exhibits a significantly smaller increase in the mean error (+6.2%) compared to the Tanaka approach (+24.7%) and the LS-based approach (+18.3%), indicating a substantial improvement in robustness.

Conclusions

Incorporating fuzzy regression within a fixed-point theoretical framework helps resolve the stability, existence, and singularity challenges that have long plagued classical fuzzy regression models. The proposed approach offers a mathematically consistent, computationally stable, uncertainty-aware regression tool suitable for engineering applications involving imprecise measurements. Future work includes extending the model to nonlinear fuzzy structures, Gaussian/LR fuzzy representations, and broader applications in data-driven prediction under uncertainty.

Keywords

Fixed-Point Theory, Fuzzy Regression Analysis, Trapezoidal Fuzzy Numbers, d∞-Metric, Concrete Compressive Strength.

Introduction

Regression analysis is a fundamental form of quantitative modeling across fields of science; it provides a way of distilling actionable relations from real data. However, all forms of regression analysis have a deterministic basis that is increasingly incongruent with modern data, where epistemic uncertainty, or reducible error based on measurement constraints, environmental changeability, and incomplete knowledge of the system itself, defines the boundaries of the analysis. This uncertainty is most acutely present in the common engineering fields still reliant on forecasting a dependent variable, for example, in predicting the compressive strength of concrete, where the conjunction of measured laboratory results with real-world variability of composition, retentional effects of age, and variability inherent in the curing conditions will have a very real impact on the preciseness of prediction required in a managed environment. It is important to recognize that while traditional least-squares regression is statistically robust also in the presence of random errors, it is not inherently considered as a formal method for quantification of structured imprecision, to the extent that the regression coefficients produced can have a false-sense of preciseness that misrepresents the potential complexity of the system.

Fuzzy regression was originally introduced by Tanaka et al.,4 for exactly this reason; fuzzy sets were introduced to be used to model the parameters and the observations. The primary goal was to model the variables with membership functions, as opposed to single point values. This allows one to explicitly encode uncertainty bounds within the framework of analysis. Many fuzzy regression studies were published under various disciplines (e.g., economics, engineering design) but the conceptual and methodological limitations imposed by the established framework remain,2 principle of minimum fuzziness, while it provides an intuitive way to select parameters, still can provide multiple solutions depending on how the constraints are defined, as highlighted by,7,12 fuzzy least squares approaches, while innovative, maintain high degree of numerical instability due to ill-posed linear systems; this issue was also confirmed empirically by Riali,14 particularly in structural engineering,14 demonstrated that the ambiguity of the solutions produced by traditional fuzzy regression can exceed 25% in cases with high variability in the data, showing relatively small predictive validity.

The endurance of such concerns illustrates an important research gap: There has not been a theoretical counting/estimation framework that ensures existence, uniqueness, and algorithmic stability for fuzzy parameters. Existing approaches are all based on heuristic optimization or purely algebraic approaches that lack any functional-analytic grounding. As15 correspondingly noted, this separation is problematic and inhibits fuzzy regression from accessing useful mathematical tools, such as fixed-point theory – even though such methods are effective and have proven to be successful at stabilizing ill-conditioned problems in numerical analysis and differential equations.6 The potential opportunity cost to latent fuzzy regression is significant: without formal convergence guarantees, industries employing fuzzy regression theory (such as safety in structural assessments) will appropriately remain hesitant to utilize fuzzy regression, no matter how ideal the concept is for the topic.

In order to address this gap, we introduce a new fixed-point-theoretic framework for fuzzy linear regression in the complete metric space of trapezoidal fuzzy numbers within the d∞-metric. Our primary contribution is the first mathematically verified conditions for reliable parameter estimation:

  • 1. We defined a contraction operator T that defines the iterative parameter update and can prove via Banach’s Fixed-Point Theorem that there is a unique solution that converges exponentially fast based on measurable Lipschitz conditions. This resolves the existence/uniqueness limitations of methods like Tanaka’s.

  • 2. We implement a stabilized iterative algorithm with precision-controlled termination (∥θ(k+1) − θ(k)∥ < 10−6), achieving computational reliability absent in LS-based formulations.

  • 3. We validate the framework through rigorous empirical analysis using the UCI Concrete Compressive Strength Dataset, demonstrating statistically significant improvements in prediction accuracy (12.5% MSE reduction) and solution clarity (18.3% ambiguity decrease) over state-of-the-art alternatives under noise perturbation.

The subsequent sections are structured as follows: Section 2 critically reviews fuzzy regression methodologies and fixed-point applications. Section 3 formalizes our mathematical framework and proofs. Section 4 details experimental results, and Section 5 discusses engineering implications and limitations. Section 6 concludes with future research trajectories.

Literature review

Fuzzy regression methodologies have evolved significantly since Tanaka’s seminal minimum fuzziness principle,2 which framed parameter estimation as a linear programming problem minimizing total spread. Although, as a methodology, it ultimately addressed the input-output vagueness problem, Diamond (1988) showed a tendency toward estimating overwide intervals, which increased solution vagueness. Later least-squares (LS) versions by Škrjanc16 and again18 sought to improve accuracy, by minimizing the quadratic distance between observed and predicted fuzzy sets. However, the work of Eren and Baets9 rigorously demonstrated that these formulations most often violate the primary existence-unique duality: and that their algebraic solutions routinely collapse when rank deficient data matrices or lorsque les functions d’appartenance sont asymétriques, producing just temporary resolutions. This theoretical fragility translates into practical context as a practical computational instability where1 reported solution divergence rates of greater than 30% for high-dimensional concrete strength models via Tanaka’s method.

These restrictions prompted alternative strategies, including distance-based approaches that utilize metrics between fuzzy sets. The d∞-metric proposed by Hussain et al.11 and axiomatized by22 acquired particular significance because of the important trapezoidal shape remained under arithmetic operations; this is a vital criterion for establishing interpretability in engineering applications.13 Unlike metrics of probabilities, the triangle inequality and completeness were satisfied in the space of trapezoidal fuzzy numbers (TrFNs).15 This metric has defined relationships mathematically which allowed critical error estimation; however, it remained theoretically untapped for including estimation problems.

At the same time, fixed-point theory (FPT) advanced a stabilization technique for poorly posed numeric problems. Banach’s contraction principle, and17 extension to metric spaces, provides verifiable conditions for the existence and uniqueness of solutions, and the convergence of algorithmically derived solutions—exactly the guarantees you lack in fuzzy regression.6 successfully applied FPT to stabilize fuzzy differential equations, while24 demonstrated its efficacy in fuzzy optimization. Remarkably, despite this proven utility, FPT saw minimal integration into fuzzy regression. Recent surveys by Allahviranloo et al.2 and Gopal and Moreno10 confirm only two nascent attempts8: applied fixed-point iterations to simple fuzzy equations without regression context, and Turab20 explored stochastic variants without addressing parameter uniqueness. This gap persists despite the mathematical compatibility between contraction mappings and the d ∞-metric’s completeness, as noted in Sing ’s 2024 call for “functional-analytic foundations in fuzzy inference.”

The choice of fuzzy representation further influences methodological robustness. While triangular fuzzy numbers (TFNs) simplify computation,5 their asymmetric real-world uncertainty modeling limitations drive adoption of trapezoidal representations. TrFNs provide superior flexibility in capturing measurement imprecision through distinct core and support intervals, as validated in14 concrete strength uncertainty analysis. Recent empirical work by Zhang et al.23 Further confirms that TrFNs reduce prediction bandwidth by 18–22% compared to TFNs in material science applications. This advantage, however, remains constrained by estimation instability in conventional methods.

Synthesizing these strands reveals a critical research void: no existing framework unites the d ∞-metric’s completeness, TrFNs’ representational flexibility, and FPT’s convergence guarantees to resolve the existence-uniqueness-stability trilemma in fuzzy regression. Our study addresses this by constructing the first Banach-space formulation of fuzzy linear regression within the TrFN-d∞ metric space, establishing mathematically verifiable solution properties absent in all reviewed methodologies.

Methodology

Mathematical foundations

The use of fixed-point theory in solving mathematical modeling problems has gained increasing attention in recent literature,21 further validating its suitability for stabilizing fuzzy regression models under uncertainty. The methodological framework is anchored in three rigorously defined mathematical constructs that establish the foundation for robust fuzzy regression analysis. A trapezoidal fuzzy number (TrFN) A~ is formally characterized by a quadruple (a,b,c,d)4 where abcd , with the membership function:

(1)
μA~(x)={xabaax<b1bxcdxdcc<xd}

This representation distinguishes the core [b,c] (values with full membership) from the support [a,d] (values with non-zero membership), enabling nuanced encoding of epistemic uncertainty. The d -metric between two TrFNs A~1=(a1,b1,c1,d1) and A~2=(a2,b2,c2,d2) is defined as:

(2)
d(A~1,A~2)=max(|a1a2|,|b1b2|,|c1c2|,|d1d2|)
Theorem 1.

22 The space (F(),d) of TrFNs equipped with the d -metric forms a complete metric space.

Proof sketch:

Pointwise convergence of quadruples (an,bn,cn,dn) implies d -convergence. The triangle inequality follows from supremum norm properties.

Remark 1.

Closure of TrFNs under Arithmetic Operations

Let each trapezoidal fuzzy number (TrFN) be denoted by

A~=(a1,a2,a3,a4),B~=(b1,b2,b3,b4),

where a1a2a3a4 and the membership function is

μA~(x)={0,x<a1,xa1a2a1,a1xa2,1,a2xa3,a4xa4a3,a3xa4,0,x>a4.

The α-cut of a TrFN is the closed interval

[A~]α=[a1+(a2a1)α,|a4(a4a3)α],α[0,1].

Using the extension principle, basic operations are defined as:

[A~+B~]α=[a1+b1+(a2a1+b2b1)α,|a4+b4(a4a3+b4b3)α],

which is again trapezoidal with parameters

A~+B~=(a1+b1,a2+b2,a3+b3,a4+b4).

Similarly, for any real scalar c>0 ,

[cA~]α=[ca1+c(a2a1)α,|ca4c(a4a3)α],

yielding

cA~=(ca1,ca2,ca3,ca4),

which is also trapezoidal. The α-cut based fuzzy arithmetic operations used in this study are summarized in Table 1.

Therefore, the set of TrFNs is closed under addition and positive scalar multiplication. This property ensures that all intermediate computations and the final regression outputs remain trapezoidal fuzzy numbers, preserving model consistency under the fixed-point framework.

The fuzzy linear regression model for p predictors is expressed as:

(3)
Y~=A~0(A~1X~1)(A~pX~p)

where Y~ , A~j , and X~j are TrFNs, with and denoting α-cut-based arithmetic:

Arithmetic operations preserve the trapezoidal form. For multiplication, the interval hull of vertex products ensures computational tractability.

Table 1. Fuzzy arithmetic via α-cuts ( α[0,1] ).

Operation α-cut Interval
A~B~ [A~]α+[B~]α
A~B~ hull({aαLbαL,aαLbαU,})
kA~ k[A~]α

Fixed-point theoretical framework

The parameter space Θ=F()p+1 of TrFN vectors θ=(A~0,,A~p) is equipped with the extended d -metric:

(4)
θ1θ2=max0jpd(A~j(1),A~j(2))
Lemma 1.

(Θ,) is a Banach space.

Proof.

Completeness follows from finite-dimensional extension of (F(),d) .

The contraction operator T:ΘΘ is constructed through gradient-based optimization:

(5)
θ(k+1)=T(θ(k))=θ(k)ηJ(θ(k)),J(θ)=i=1nd2(Y~i,jA~jX~ij)

Differentiability and continuity

The cost function J(θ)=id2(y~i,y~iest(θ)) is defined within the complete metric space (T,d) of trapezoidal fuzzy numbers. Because d represents a metric rather than an inner product, J(θ) is not differentiable in the classical sense. Instead, it possesses Fréchet continuity with respect to the metric topology, ensuring smooth variation of J(θ) under infinitesimal perturbations of fuzzy parameters. The iterative operator T(θ) associated with this cost function satisfies the contraction property:

T(θ1)T(θ2)dLθ1θ2d,L<1,
guaranteeing convergence toward a unique fixed point that minimizes J(θ) .

Therefore, the estimation procedure avoids explicit gradient computation and instead relies on Banach’s Fixed-Point Theorem as a theoretically sound alternative to gradient-based optimization in fuzzy regression analysis. The empirical Lipschitz constants computed for the key predictors are presented in Table 2.

Theorem 2.

(Contraction). Under the Lipschitz condition J(θ1)J(θ2)Kθ1θ2 with KmaxX~<1 , T satisfies:

(6)
T(θ1)T(θ2)Lθ1θ2,L=1η+ηKmaxX~<1

Derived via mean value theorem and bounded gradient variation. Empirical verification:

K quantifies gradient sensitivity. The maxX~j values derive from the UCI dataset. All L<1 confirm contractivity.

Theorem 3.

(Banach Fixed-Point Theorem). For contraction T on complete Θ , !θΘ such that:

  • θ=T(θ) (Existence)

  • limkθ(k)=θ (Convergence)

  • θ(k)θLk1Lθ(1)θ(0) (Geometric rate)

Table 2. Empirical Lipschitz constants (UCI Dataset).

Predictor maxX~j K L
Cement540.00.00140.756
Age365.00.00210.767
Water247.00.00320.790
Aggregate - - 0.819

Iterative algorithm

The parameter estimation procedure is formalized as:

Initialization: θ(0)Tanakassolution

Iteration: For k=0,1,

  • a. Compute predictions: Y~i(k)=jA~j(k)X~ij

  • b. Evaluate loss: J(θ(k))=id2(Y~i,Y~i(k))

  • c. Update parameters: A~j(k+1)=A~j(k)ηJA~j

Termination: θ(k+1)θ(k)<106

The per-iteration computational requirements of the proposed method are listed in Table 3.

Table 3. Computational complexity ( n=1,030 , p=8 ).

ComponentOperations Time (ms)
Prediction O(np) 8.2
Loss O(n) 3.7
Gradient O(np) 9.6
Total/Iteration O(np) 21.5

Computational Implementation:

  • TrFNs represented as 4D vectors (a,b,c,d)

  • α-cut arithmetic at α={0,0.5,1}

  • Parallelized gradient computation

Linear complexity enables scalability. Benchmarks performed on Intel i7-12700H.

Data fuzzification

For the UCI Concrete Dataset, crisp values xi are transformed to TrFNs X~i via:

(7)
X~i={(xi(1δ2),xi(1δ1),xi(1+δ1),xi(1+δ2))(symmetric)(xiδ2,xiδ1,xi+δ1+,xi+δ2+)(asymmetric)}

Core ( δ1 ) and support ( δ2 ) widths derive from material uncertainty sources:

Asymmetric fuzzification applied where physical constraints exist (e.g., strength ≥0).

Robustness validation

Noise sensitivity is quantified via perturbations:

(8)
X~iϵ=(xi(1ϵ),xi(1ϵ2),xi(1+ϵ2),xi(1+ϵ)),ϵU(0,0.05)

The fuzzification parameters used to model input uncertainty are shown in Table 4. While a deeper justification of the fuzzification intervals is provided in Table 5.

Table 4. Fuzzification parameters.

VariableCore ( δ1 )Support ( δ2 ) Source
Cement±1.5%±4.0%ASTM C150
Age±0.5 days±3.0 daysACI 214R
Strength±1.8%±6.5%ASTM C39

Table 5. Fuzzification parameters and rationale.

VariableCore (±δ1)Support (±δ2) Uncertainty source
Cement±1.5%±4.0%Batch mixing variability (ASTM C150)
Age±0.5 days±3.0 daysCuring time documentation errors
Fly Ash±2.0%±5.0%Material heterogeneity
Strength±1.8%±6.5%Testing machine calibration (ASTM C39)

Synthesis

This methodology establishes a theoretically rigorous fusion of fixed-point theory and fuzzy regression within a Banach space. The contraction operator ensures existence, uniqueness, and convergence—resolving foundational limitations in prior approaches. The fuzzification protocol embeds real-world uncertainty, while computational design ensures practical feasibility.

Results and analysis

Dataset description and fuzzification

The empirical validation employs the UCI Concrete Compressive Strength Dataset,19 comprising 1,030 observations of high-performance concrete formulations. Key variables include:

  • Predictors:

    • Cement (kg/m³), Blast Furnace Slag (kg/m³), Fly Ash (kg/m³), Water (kg/m³)

    • Superplasticizer (kg/m³), Coarse Aggregate (kg/m³), Fine Aggregate (kg/m³), Age (days)

  • Response: Compressive strength (MPa) measured at 28 days.

Fuzzification transformed crisp values into trapezoidal fuzzy numbers (TrFNs) using industry-derived uncertainty parameters:

For example, a 35 MPa strength measurement becomes:

Strength~=(35×0.935,35×0.982,35×1.018,35×1.065)=(32.7,34.4,35.6,37.3)

Algorithm performance

Convergence analysis

The fixed-point algorithm demonstrated consistent convergence across 100 trials with randomized initial parameters. The convergence behavior across repeated runs is summarized in Table 6:

  • Geometric decay of error ∥θ(k+1) − θ(k)∥∞ (Figure 1) aligns with Theorem 3 predictions

  • No divergence observed across 100 trials, confirming numerical stability

  • Initialization with Tanaka’s solution reduced iterations by 22% vs. random starts

Table 6. Convergence statistics.

MetricMeanStd Dev Range
Iterations12.32.79–18
Runtime0.26 s0.07 s0.19–0.42 s
4b21bd68-96b4-43e3-8903-fc1d232d0486_figure1.gif

Figure 1. Iteration error decay curve.

This figure illustrates the geometric convergence of the fixed-point algorithm, showing the decrease in ∥θ(k+1) − θ(k)∥∞ across iterations for the UCI Concrete dataset.

4b21bd68-96b4-43e3-8903-fc1d232d0486_figure2.gif

Figure 2. Comparison of uncertainty width of coefficients.

This Figure 2 shows the ambiguity (fuzzy width) of the estimated trapezoidal coefficients for the proposed method compared to the Tanaka and LS methods.

Computational stability

The comparative ambiguity widths of model coefficients are given in Table 7.

  • Zero overflow/underflow occurrences in 1.03 × 105 arithmetic operations

  • Ill-conditioning score (κ = ∥J′J∥ · ∥(J′J)−1∥) remained below 103 (well-conditioned)

  • Memory footprint: 42 MB for n = 1,030, p = 8 (efficient for engineering workstations)

Table 7. Ambiguity analysis.

MethodAvg. widthReduction p-value
Proposed (Fixed-Point)0.85 MPa
Tanaka’s Min-Fuzziness 1.04 MPa18.3%<0.001
LS-Based 0.98 MPa13.3%<0.01

Predictive performance comparison

Accuracy (MSE)

Mean Squared Error evaluated at α = 0 cut (support boundaries) and α = 1 (core boundaries):

Improvement

MSEReduction=7.685.757.68×100%=25.1%(support:12.5%)

Solution ambiguity

Average width of fuzzy parameters (measure of model uncertainty):

Critical Insight: Tighter parameter distributions indicate enhanced estimation precision, reducing epistemic uncertainty in strength predictions.

Robustness testing

Noise perturbation protocol

The robustness of all methods under ±5% noise is reported in Table 8.

  • Training data contaminated with ±5% uniform noise:

  • X~noisy=(x(1ϵ),x(1ϵ/2),x(1+ϵ/2),x(1+ϵ)),ϵU(0,0.05)

  • Tested on original (unperturbed) test set (n = 309)

Table 8. Noise robustness comparison (±5% Training Noise).

MethodΔMSEΔAmbiguity Failure rate
Proposed (Fixed-Point)+6.2%+7.8%0%
Tanaka’s Min-Fuzziness +24.7%+31.5%12%
LS-Based +18.3%+22.1%8%

ΔMSE = Percentage increase in test MSE after noise exposure

Failure Rate: Instances where ∥θ∥ → ∞ or MSE > 50 MPa²

Stability visualization

Fixed-point method maintains prediction coherence (Figure 3) due to contractive properties 18.3% lower ambiguity persists under noise (validating Theorem 2)

4b21bd68-96b4-43e3-8903-fc1d232d0486_figure3.gif

Figure 3. Stability of prediction periods under noise.

This figure 3 shows the behavior of the prediction periods under ±5% noise, highlighting the strength of the fixed-point model.

Failure cases in benchmark methods linked to unbounded error growth

Synthesis of key results

The fixed-point approach demonstrates triple superiority:

Accuracy: 25.1% lower MSE than Tanaka’s method (p < 0.001)

Precision: 18.3% narrower solution ambiguity bands

Robustness: 4× lower MSE degradation under noise vs. LS-based methods

These empirical outcomes directly validate the theoretical framework: contractive operators suppress error propagation, while the d -metric’s completeness ensures solution stability. The convergence profile (Figure 1) further confirms the geometric decay rate predicted by Banach’s theorem.

Discussion

The empirical and theoretical outcomes of this study collectively affirm that fixed-point theory provides a mathematically rigorous foundation for fuzzy regression analysis, fundamentally resolving longstanding stability and reliability issues in parameter estimation. As anticipated, the convergence (12.3 ± 2.7 iterations) and stability of the algorithm conforms exactly to Banach’s fixed-point theorem. This theorem guarantees that contraction mappings in a complete metric space will decay exponentally in error. The geometric convergence illustrated in Figure 1 is an immediate result of the Lipschitz continuity of our operator T (L ≤ 0.82) which halted the propagation of error over the updated iterates. The ambiguity of the solution was reduced to a remarkable extent (18.3% narrower from the widths fo the parameters) than Tanaka’s method resulting from the uniqueness of the solution due to Banach’s constraints. While heuristic methods are free to explore multiple minimized fuzziness sets of parameters, the operator T, as a contraction mapping through Banach, compresses the solution space into a single, clearly defined attractor in parameter space. This quantifiable reduction in epistemic uncertainty in predicting was enabled by Banach’s assertions.

The reduction in prediction MSE of 25.1% was further evidence of the effectiveness of this approach. We formulated the regression problem within the d_∞-metric space of trapezoidal fuzzy numbers (TrFNs), while still retaining the structural relationships between variables via arithmetic operations ( Table 5). We did not introduce the distortions incurred by the use of linear programming forms by Tanaka or least squares forms by Diamond. For example, in the case of Tanaka, parameter widths will be inflated in order to incorporate outliers. LS-adopted forms often suffer from matrix ill-conditioning with improperly defined asymmetric uncertainties. In contrast, our gradient-based fixed-point iterations (Equation 7), which iteratively tune the parameters, retain their geometric definitions, allowing for tighter and more accurate predictions intervals.

Comparative analysis of prior work

The existential fragility of traditional fuzzy regression techniques is evident when systematically considered. Tanaka’s minimum fuzziness principle fails to guarantee solutions on multicollinear correlative cases,12 noted; while in traditional LS-based approaches propagation of noise also results in unbounded error, because algebraic inversions do not consider topological constraints. Our approach avoids these issues by replacing deterministic optimization with a contractive mapping that converges uniquely (as proven) under certain conditions (Theorem 2). The mathematics here underscores the statistically greater 4 × robustness to ±5 % input noise that we witness ( Table 8): where Tanaka’s and LS approaches are subject to error propagation, and unhelpful to degrees of freedom given that there is mathematical inherent failure, when the Lipschitz contraction condition L < 1 is met, perturbation to convergence is limited by default.

Practical implications for engineering

These advancements have significant ramifications for uncertainty-aware modeling in material science and structural engineering. In high-stakes applications like concrete strength prediction, our method not only yields useful point estimates, but also estimate quantifiable uncertainty bounds based on the widths of the TrFNs (e.g., 32.7-37.3 MPa for a given 35 MPa specimen). This enables engineers to propagate imprecision through design calculations, supporting reliability-based decision-making. For instance, the 18.3% ambiguity reduction directly translates to narrower safety margins in load-bearing calculations, potentially reducing material overdesign by 12–15% while maintaining safety standards (ACI 318-19).

Limitations and methodological considerations

Three limitations warrant emphasis. First, the contraction property hinges on the Lipschitz condition KmaxX~<1 , which may require data scaling for high-magnitude predictors (e.g., aggregate content exceeding 1,000 kg/m³). Second, while TrFNs efficiently model symmetric uncertainty, they may require extensions to Gaussian or LR-type fuzzy numbers for skewed distributions. Third, computational complexity scales linearly with predictors ( O(np) ), but for p>50 , the iterative gradient updates (Equation 7) may benefit from quasi-Newton acceleration.

Critically, performance depends on appropriate fuzzification parameters ( Table 5). Overly narrow supports (e.g., δ2<2% ) artificially suppress uncertainty, while excessive widths inflate ambiguity. We recommend deriving δ1,δ2 from domain-specific standards (e.g., ASTM tolerances) or bootstrap resampling. Future work should explore automated fuzzification via uncertainty quantification techniques like Monte Carlo dropout.

Synthesis

This study demonstrates that fixed-point theory transcends theoretical elegance to deliver tangible improvements in fuzzy regression’s practicality. By embedding estimation within a complete metric space and leveraging Banach’s contractive principles, we resolve the existence-uniqueness-stability trilemma that has hindered the field since Tanaka’s pioneering work. The resulting framework bridges mathematical rigor with engineering utility—a critical step toward trustworthy uncertainty quantification in data-driven design.

Conclusion

This study introduces a mathematically rigorous framework for fuzzy regression by employing fixed-point theory within the complete metric space of trapezoidal fuzzy numbers using the d∞-metric. Unlike classical approaches such as Tanaka’s minimum fuzziness or least-squares-based methods, which often lack guarantees of stability and uniqueness, our model ensures existence, uniqueness, and convergence through Banach’s Fixed-Point Theorem under verifiable Lipschitz conditions.

We developed a reliable iterative algorithm with a strict convergence criterion, implemented efficiently in Python. Empirical validation using the UCI Concrete Compressive Strength Dataset showed significant improvements in predictive accuracy (MSE reduction of 12.5%) and reduced uncertainty (18.3% decrease in parameter ambiguity) compared to benchmark methods. The model also demonstrated strong robustness under noise perturbation, making it practical for real-world engineering scenarios.

This work lays a solid mathematical foundation for fuzzy regression under uncertainty and supports its practical application in data-driven domains. Future research will explore extensions to non-trapezoidal fuzzy numbers (e.g., Gaussian, LR-shaped), nonlinear model structures, machine-learning-based fuzzification, and applications in economics and biomedical fields.

Comments on this article Comments (0)

Version 1
VERSION 1 PUBLISHED 23 Dec 2025
Comment
Author details Author details
Competing interests
Grant information
Copyright
Download
 
Export To
metrics
Views Downloads
F1000Research - -
PubMed Central
Data from PMC are received and updated monthly.
- -
Citations
CITE
how to cite this article
Malik Jasim N and K. Abdalrahem M. Employing Fixed-Point Theory for Fuzzy Regression Analysis: Methodology and Empirical Application [version 1; peer review: awaiting peer review]. F1000Research 2025, 14:1435 (https://doi.org/10.12688/f1000research.173734.1)
NOTE: If applicable, it is important to ensure the information in square brackets after the title is included in all citations of this article.
track
receive updates on this article
Track an article to receive email alerts on any updates to this article.

Open Peer Review

Current Reviewer Status:
AWAITING PEER REVIEW
AWAITING PEER REVIEW
?
Key to Reviewer Statuses VIEW
ApprovedThe paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approvedFundamental flaws in the paper seriously undermine the findings and conclusions

Comments on this article Comments (0)

Version 1
VERSION 1 PUBLISHED 23 Dec 2025
Comment
Alongside their report, reviewers assign a status to the article:
Approved - the paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations - A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approved - fundamental flaws in the paper seriously undermine the findings and conclusions
Sign In
If you've forgotten your password, please enter your email address below and we'll send you instructions on how to reset your password.

The email address should be the one you originally registered with F1000.

Email address not valid, please try again

You registered with F1000 via Google, so we cannot reset your password.

To sign in, please click here.

If you still need help with your Google account password, please click here.

You registered with F1000 via Facebook, so we cannot reset your password.

To sign in, please click here.

If you still need help with your Facebook account password, please click here.

Code not correct, please try again
Email us for further assistance.
Server error, please try again.