Abstract:The inference performance of the pseudolikelihood method is discussed in the framework of the inverse Ising problem when the l (2)-regularized (ridge) linear regression is adopted. This setup is introduced for theoretically investigating the situation where the data generation model is different from the inference one, namely the model mismatch situation. In the teacher-student scenario under the assumption that the teacher couplings are sparse, the analysis is conducted using the replica and cavity methods, with a special focus on whether the presence/absence of teacher couplings is correctly inferred or not. The result indicates that despite the model mismatch, one can perfectly identify the network structure using naive linear regression without regularization when the number of spins N is smaller than the dataset size M, in the thermodynamic limit N -> infinity. Further, to access the underdetermined region M < N, we examine the effect of the l (2) regularization, and find that biases appear in all the coupling estimates, preventing the perfect identification of the network structure. We, however, find that the biases are shown to decay exponentially fast as the distance from the center spin chosen in the pseudolikelihood method grows. Based on this finding, we propose a two-stage estimator: in the first stage, the ridge regression is used and the estimates are pruned by a relatively small threshold; in the second stage the naive linear regression is conducted only on the remaining couplings, and the resultant estimates are again pruned by another relatively large threshold. This estimator with the appropriate regularization coefficient and thresholds is shown to achieve the perfect identification of the network structure even in 0 < M/N < 1. Results of extensive numerical experiments support these findings.

Representative of L1/2 Regularization among Lq (0 < Q ≤ 1) Regularizations: an Experimental Study Based on Phase Diagram

Combining Regularization Frameworks for Solving the Electrocardiography Inverse Problem.

Choosing Near-Optimal Regularization Parameter for the Inverse Problem of Electrocardiography

Structure Learning in Inverse Ising Problems Using ℓ_2-Regularized Linear Estimator

Sparse Estimation Via ℓ_q Optimization Method in High-Dimensional Linear Regression

A Sharp Nonasymptotic Bound and Phase Diagram of L 1/2 Regularization

Sparse Estimation Via Lower-Order Penalty Optimization Methods in High-Dimensional Linear Regression.

Learning Rates of L^q Coefficient Regularization Learning with Gaussian Kernel

ℓ _1/2,1 group sparse regularization for compressive sensing

Group Sparse Optimization Via Lp, Q Regularization.

Regularization for Compressive Sensing

L-1/2 NORM REGULARIZATION FOR 3D SEISMIC DATA INTERPOLATION

Revisiting $L_q(0\leq q<1)$ Norm Regularized Optimization

Group sparse optimization via l

A Cyclic Coordinate Descent Algorithm for lq Regularization

$L_{1/2}$ Regularization: Convergence of Iterative Half Thresholding Algorithm

Linear Convergence of Inexact Descent Method and Inexact Proximal Gradient Algorithms for Lower-Order Regularization Problems

AnLq–Lpoptimization Framework for Image Reconstruction of Electrical Resistance Tomography

&Lt;formula Formulatype="inline"> <tex Notation="TeX">$L_{1/2}$</tex></formula> Regularization: A Thresholding Representation Theory and a Fast Solver

Kernel Regression with Coefficient-based l(q)-regularization

The Essential Ability of Sparse Reconstruction of Different Compressive Sensing Strategies