Fuzzy Logistic Ridge Estimators (FLRE)

• FLRE is a fuzzy logistic regression approach that models both responses and coefficients as triangular fuzzy numbers, integrating classical ridge penalization.
• It stabilizes parameter estimates under severe multicollinearity by adding an ℓ2 penalty to the fuzzy log-likelihood, effectively reducing variance inflation.
• Empirical studies show FLRE reduces mean squared error compared to FMLE, with further improvements by Liu-type extensions in practical applications.

Fuzzy Logistic Ridge Estimators (FLRE) are regularization-based estimators specifically developed for fuzzy logistic regression models exhibiting significant multicollinearity. They extend classical ridge penalization into the fuzzy setting, in which both the response variable and the parameter vector are represented as triangular fuzzy numbers. FLREs aim to stabilize parameter estimation by counteracting variance inflation induced by collinear predictors, while incorporating fuzziness intrinsic to both outcomes and coefficients, as formulated in recent work by Shemail et al. (Shemail et al., 27 Dec 2025).

1. Fuzzy Logistic Regression Model Formulation

In the fuzzy logistic regression framework, every observation and regression coefficient is modeled as a triangular fuzzy number:

• Each response, $\tilde y_i = (\ell_{y_i}, m_{y_i}, u_{y_i})$, possesses a membership function

$$\mu_{\tilde y_i}(y) = \begin{cases} \frac{y - \ell_{y_i}}{m_{y_i} - \ell_{y_i}}, & \ell_{y_i} \leq y \leq m_{y_i} \ \frac{u_{y_i} - y}{u_{y_i} - m_{y_i}}, & m_{y_i} < y \leq u_{y_i} \ 0, & \text{otherwise} \end{cases}$$

• The parameter vector is also fuzzy: $\tilde\beta = (\ell_\beta, m_\beta, u_\beta)$, each entry being triangular.

The success probability for observation $i$ is given by the fuzzy logistic function,

$$\tilde\pi_i = \Pr(\tilde y_i=1|x_i) = \frac{\exp(\tilde\beta^{\top}x_i)}{1+\exp(\tilde\beta^{\top}x_i)},$$

with the logit transform,

$$\operatorname{logit}(\tilde\pi_i) = x_i^\top\tilde\beta.$$

While covariates $x_i\in\mathbb{R}^p$ are crisp, separate equations for lower, center, and upper levels of fuzziness are applied throughout.

2. Construction of the Fuzzy Logistic Ridge Estimator

FLREs are derived by augmenting the fuzzy log-likelihood with a ridge penalty. The log-likelihood for the crisp case is adapted to the fuzzy setting—for each $\alpha$-cut, the form persists:

$$\ell(\beta) = \sum_{i=1}^n \left[ y_i\ln\pi_i + (1 - y_i)\ln(1 - \pi_i) \right], \qquad \pi_i = \frac{e^{x_i^\top\beta}}{1+e^{x_i^\top\beta}}.$$

The penalized criterion is

$$\ell_R(\beta; k) = \ell(\beta) - k\,\beta^\top\beta, \qquad k > 0.$$

The estimating equations become:

$$\mu_{\tilde y_i}(y) = \begin{cases} \frac{y - \ell_{y_i}}{m_{y_i} - \ell_{y_i}}, & \ell_{y_i} \leq y \leq m_{y_i} \ \frac{u_{y_i} - y}{u_{y_i} - m_{y_i}}, & m_{y_i} < y \leq u_{y_i} \ 0, & \text{otherwise} \end{cases}$$0

where $$\mu_{\tilde y_i}(y) = \begin{cases} \frac{y - \ell_{y_i}}{m_{y_i} - \ell_{y_i}}, & \ell_{y_i} \leq y \leq m_{y_i} \ \frac{u_{y_i} - y}{u_{y_i} - m_{y_i}}, & m_{y_i} < y \leq u_{y_i} \ 0, & \text{otherwise} \end{cases}$$1 is the model matrix and $$\mu_{\tilde y_i}(y) = \begin{cases} \frac{y - \ell_{y_i}}{m_{y_i} - \ell_{y_i}}, & \ell_{y_i} \leq y \leq m_{y_i} \ \frac{u_{y_i} - y}{u_{y_i} - m_{y_i}}, & m_{y_i} < y \leq u_{y_i} \ 0, & \text{otherwise} \end{cases}$$2 is the vector of fitted probabilities for each cut.

A closed-form FLRE solution, akin to a one-step Newton–Raphson update around the Fuzzy Maximum Likelihood Estimator (FMLE), is:

$$\mu_{\tilde y_i}(y) = \begin{cases} \frac{y - \ell_{y_i}}{m_{y_i} - \ell_{y_i}}, & \ell_{y_i} \leq y \leq m_{y_i} \ \frac{u_{y_i} - y}{u_{y_i} - m_{y_i}}, & m_{y_i} < y \leq u_{y_i} \ 0, & \text{otherwise} \end{cases}$$3

with

$$\mu_{\tilde y_i}(y) = \begin{cases} \frac{y - \ell_{y_i}}{m_{y_i} - \ell_{y_i}}, & \ell_{y_i} \leq y \leq m_{y_i} \ \frac{u_{y_i} - y}{u_{y_i} - m_{y_i}}, & m_{y_i} < y \leq u_{y_i} \ 0, & \text{otherwise} \end{cases}$$4

A compact “shrinkage-follow-MLE” representation (Eq. 13 of (Shemail et al., 27 Dec 2025)):

$$\mu_{\tilde y_i}(y) = \begin{cases} \frac{y - \ell_{y_i}}{m_{y_i} - \ell_{y_i}}, & \ell_{y_i} \leq y \leq m_{y_i} \ \frac{u_{y_i} - y}{u_{y_i} - m_{y_i}}, & m_{y_i} < y \leq u_{y_i} \ 0, & \text{otherwise} \end{cases}$$5

3. Theoretical Properties under Multicollinearity

Let $$\mu_{\tilde y_i}(y) = \begin{cases} \frac{y - \ell_{y_i}}{m_{y_i} - \ell_{y_i}}, & \ell_{y_i} \leq y \leq m_{y_i} \ \frac{u_{y_i} - y}{u_{y_i} - m_{y_i}}, & m_{y_i} < y \leq u_{y_i} \ 0, & \text{otherwise} \end{cases}$$6. The theoretical risk properties under multicollinearity are characterized as follows:

• Bias:

$$\mu_{\tilde y_i}(y) = \begin{cases} \frac{y - \ell_{y_i}}{m_{y_i} - \ell_{y_i}}, & \ell_{y_i} \leq y \leq m_{y_i} \ \frac{u_{y_i} - y}{u_{y_i} - m_{y_i}}, & m_{y_i} < y \leq u_{y_i} \ 0, & \text{otherwise} \end{cases}$$7

• Variance:

$$\mu_{\tilde y_i}(y) = \begin{cases} \frac{y - \ell_{y_i}}{m_{y_i} - \ell_{y_i}}, & \ell_{y_i} \leq y \leq m_{y_i} \ \frac{u_{y_i} - y}{u_{y_i} - m_{y_i}}, & m_{y_i} < y \leq u_{y_i} \ 0, & \text{otherwise} \end{cases}$$8

• Mean Squared Error (MSE) matrix:

$$\mu_{\tilde y_i}(y) = \begin{cases} \frac{y - \ell_{y_i}}{m_{y_i} - \ell_{y_i}}, & \ell_{y_i} \leq y \leq m_{y_i} \ \frac{u_{y_i} - y}{u_{y_i} - m_{y_i}}, & m_{y_i} < y \leq u_{y_i} \ 0, & \text{otherwise} \end{cases}$$9

In scalar form,

$\tilde\beta = (\ell_\beta, m_\beta, u_\beta)$0

• Asymptotic Properties: As $\tilde\beta = (\ell_\beta, m_\beta, u_\beta)$1 and $\tilde\beta = (\ell_\beta, m_\beta, u_\beta)$2, bias vanishes and variance approaches the fuzzy Fisher information bound.

This framework controls the inflation of estimator variance associated with severe collinearity while introducing a bias that can be analytically characterized and, if $\tilde\beta = (\ell_\beta, m_\beta, u_\beta)$3 is suitably chosen, minimized in predictive risk.

4. Choice of Ridge Penalty Parameter

Several methods are established for selecting the optimal ridge penalty $\tilde\beta = (\ell_\beta, m_\beta, u_\beta)$4 in the fuzzy context:

• Cross-Validation minimizes prediction error on held-out folds.
• Generalized Cross-Validation (GCV) provides a computationally efficient approximation to leave-one-out error.
• Information Criteria: $\tilde\beta = (\ell_\beta, m_\beta, u_\beta)$5 can be chosen by minimizing $\tilde\beta = (\ell_\beta, m_\beta, u_\beta)$6, where the degrees-of-freedom depend on the effective number of parameters.
• Analytical Shrinkage Formula: Asar & Genç (2016) propose using the eigen-decomposition $\tilde\beta = (\ell_\beta, m_\beta, u_\beta)$7 with $\tilde\beta = (\ell_\beta, m_\beta, u_\beta)$8 eigenvalues and $\tilde\beta = (\ell_\beta, m_\beta, u_\beta)$9, yielding a median-based shrinkage factor

$i$0

and penalty

$i$1

These approaches adapt classical criterion to the fuzzy parameter setting, allowing $i$2 to be robustly tuned to the magnitude and structure of multicollinearity.

5. Empirical Performance: Simulation Results

Shemail et al. report simulation studies comparing FMLE, FLRE, FLLE, FLLTE, and FLLTPE at high collinearity ($i$3), $i$4, with $i$5 determined via median-shrinkage:

Estimator	Goodness-of-Fit $i$6	MSE
FMLE	1.9456	2.6909
FLRE	1.5345	2.2798
FLLE	1.1345	1.8798
FLLTE	0.5345	1.2798
FLLTPE	0.4123	1.1576

The FLRE outperforms the fuzzy maximum likelihood estimator (FMLE) in both mean squared error and goodness-of-fit, but FLLTE and FLLTPE show further marked improvements in both measures. This suggests that FLRE's penalization is effective in regularizing fuzzy logistic regression estimates under collinearity, with greater gains realized through Liu-type regularization strategies.

6. Real Data Evaluation: Kidney Failure Study

FLRE was evaluated on a real dataset comprising $i$7 kidney failure patients with fuzzy outcome labels ("renal failure rate") and crisp covariates: blood sugar ($i$8), hemoglobin ($i$9), and age ($$\tilde\pi_i = \Pr(\tilde y_i=1|x_i) = \frac{\exp(\tilde\beta^{\top}x_i)}{1+\exp(\tilde\beta^{\top}x_i)},$$0). Severe multicollinearity was diagnosed (condition number $$\tilde\pi_i = \Pr(\tilde y_i=1|x_i) = \frac{\exp(\tilde\beta^{\top}x_i)}{1+\exp(\tilde\beta^{\top}x_i)},$$1).

Estimator performance is as follows (from Table 9 of (Shemail et al., 27 Dec 2025)):

Method	$$\tilde\pi_i = \Pr(\tilde y_i=1|x_i) = \frac{\exp(\tilde\beta^{\top}x_i)}{1+\exp(\tilde\beta^{\top}x_i)},$$2	$$\tilde\pi_i = \Pr(\tilde y_i=1|x_i) = \frac{\exp(\tilde\beta^{\top}x_i)}{1+\exp(\tilde\beta^{\top}x_i)},$$3	$$\tilde\pi_i = \Pr(\tilde y_i=1|x_i) = \frac{\exp(\tilde\beta^{\top}x_i)}{1+\exp(\tilde\beta^{\top}x_i)},$$4	MSE	$$\tilde\pi_i = \Pr(\tilde y_i=1|x_i) = \frac{\exp(\tilde\beta^{\top}x_i)}{1+\exp(\tilde\beta^{\top}x_i)},$$5
FMLE	3.2632	4.5128	0.4649	1.8073	0.5422
FLRE	2.0395	2.8205	0.2906	1.4152	0.4846
FLLE	2.2661	3.1339	0.3228	1.2882	0.3865
FLLTE	1.8129	2.5071	0.2583	0.5672	0.1702
FLLTPE	1.5863	2.1937	0.2260	0.5342	0.1603

FLRE reduced MSE by approximately 22% relative to FMLE. FLLTE and FLLTPE provided further reductions, nearly halving MSE and doubling the goodness-of-fit compared to ridge penalization alone. This confirms the robust bias–variance trade-off offered by ridge-type and related Liu-type penalizations in fuzzy models under multicollinearity.

7. Summary and Comparative Assessment

FLREs generalize classical ridge regularization to the fuzzy logistic regression domain by incorporating an $$\tilde\pi_i = \Pr(\tilde y_i=1|x_i) = \frac{\exp(\tilde\beta^{\top}x_i)}{1+\exp(\tilde\beta^{\top}x_i)},$$6 penalty into the fuzzy log-likelihood, offering closed-form shrinkage estimators centered at the FMLE. Explicit bias, variance, and risk expressions can be derived, and $$\tilde\pi_i = \Pr(\tilde y_i=1|x_i) = \frac{\exp(\tilde\beta^{\top}x_i)}{1+\exp(\tilde\beta^{\top}x_i)},$$7 selection is adaptable to the fuzzy setting using cross-validation, information criteria, or median-based shrinkage approaches. Empirical evaluations, both in simulation and application to kidney failure data, indicate substantive improvements over naive fuzzy maximum likelihood estimation. However, Liu-type extensions (FLLTE, FLLTPE) further lower MSE by introducing additional bias structure, suggesting their superiority in scenarios where two-parameter regularization is warranted (Shemail et al., 27 Dec 2025).