Abstract
Insomnia, as a common sleep disorder, is the most common complaints in medical practice affecting a large proportion of the population on a situational, recurrent or chronic basis. It has been demonstrated that, during wakefulness, patients with insomnia exhibit increased EEG power in theta, beta, and gamma band. However, the relevant mechanisms underlying such power changes are still lack of understanding. In this paper, by combining the neural computational model with the real EEG data, we focus on exploring what's behind the EEG power changes for insomniac. We first develop a modified Liley model, named FSR-Liley, by respectively considering the fast and slow synaptic responses in inhibitory neurons along with the one-way projection between them. Then we introduce a parameter selection and evaluation method based on Markov chain Monte Carlo algorithm and Wasserstein distance, by which the sensitive parameters are selected automatically, and meanwhile, the optimal values of selected parameters are evaluated. Finally, through combining with EEG data, we determine the sensitive parameters in FSR-Liley and accordingly provide the mechanistic hypotheses: (1) decrease in Peif , corresponding to the input from the thalamus to cortical inhibitory population with fast synaptic response, leads to the increased theta and beta power; (2) decrease in Neif , corresponding to the projection from cortical excitatory population to inhibitory population with fast synaptic response, causes the increased gamma power. The results in this paper provide insights into the mechanisms of EEG power changes in insomnia and establish a theoretical foundation to support further experimental research.