Estimating neural connection strengths from firing intervals

Abstract: We propose and analyse a procedure for using a standard activity-based neuron network model and firing data to compute the effective connection strengths between neurons in a network. We assume a Heaviside response function, that the external inputs are given and that the initial state of the neural activity is known. The associated forward operator for this problem, which maps given connection strengths to the time intervals of firing, is highly nonlinear. Nevertheless, it turns out that the inverse problem of determining the connection strengths can be solved in a rather transparent manner, only employing standard mathematical tools. In fact, it is sufficient to solve a system of decoupled ODEs, which yields a linear system of algebraic equations for determining the connection strengths. The nature of the inverse problem is investigated by studying some mathematical properties of the aforementioned linear system and by a series of numerical experiments. Finally, under an assumption preventing the effective contribution of the network to each neuron from staying at zero, we prove that the involved forward operator is continuous. Sufficient criteria on the external input ensuring that the needed assumption holds are also provided.