High-Dimensional Inference with the generalized Hopfield Model: Principal Component Analysis and Corrections

Simona Cocco (1, 2) , Remi Monasson (1, 3) , Vitor Sessak (3)

1 Simons Center for Systems Biology
2 LPS - Laboratoire de Physique Statistique de l'ENS
3 LPTENS - Laboratoire de Physique Théorique de l'ENS [École Normale Supérieure]

We consider the problem of inferring the interactions between a set of N binary variables from the knowledge of their frequencies and pairwise correlations. The inference framework is based on the Hopfield model, a special case of the Ising model where the interaction matrix is defined through a set of patterns in the variable space, and is of rank much smaller than N. We show that Maximum Lik elihood inference is deeply related to Principal Component Analysis when the amp litude of the pattern components, xi, is negligible compared to N^1/2. Using techniques from statistical mechanics, we calculate the corrections to the patterns to the first order in xi/N^1/2. We stress that it is important to generalize the Hopfield model and include both attractive and repulsive patterns, to correctly infer networks with sparse and strong interactions. We present a simple geometrical criterion to decide how many attractive and repulsive patterns should be considered as a function of the sampling noise. We moreover discuss how many sampled configurations are required for a good inference, as a function of the system size, N and of the amplitude, xi. The inference approach is illustrated on synthetic and biological data.

Domains: Statistical Mechanics [cond-mat.stat-mech]; Machine Learning [stat.ML]

Simona Cocco, Remi Monasson, Vitor Sessak. High-Dimensional Inference with the generalized Hopfield Model: Principal Component Analysis and Corrections. Physical Review E : Statistical, Nonlinear, and Soft Matter Physics, 2011. ⟨hal-00586950⟩