Deep brain stimulation (DBS) provides dramatic tremor relief in patients with severe essential tremor (ET). Typically, the VIM nucleus is the most effective brain area to target for high-frequency electrical stimulation in these patients. Correlation analysis between electrical local field potential (LFP) recordings from the thalamic DBS leads and electrical muscle activity from the contralateral tremulous limb has become an attractive practical tool to interpret the LFPs and their association with the tremulous clinical manifestations. Although functional connectivity analysis between brain electrical recordings and electromyographic (EMG) signals from the tremor has been of interest to an increasing number of engineering researchers, there is no well-accepted tailored framework to consistently characterise the association between thalamic electrical recordings and the tremorogenic EMG activity. Methods. This paper proposes a novel framework to address this challenge, including an estimation of the interaction strength using wavelet cross-spectrum and phase lag index while demonstrating the statistical significance of the findings. Results. Consistent results were estimated for single and multiple trials of consecutive or partially overlapping epochs of data. The latter approach reveals a substantial increase on the range of statistically significant dynamic low-frequency interrelationships while decreasing the dynamic range of high-frequency interactions. Conclusion. Results from both simulation and real data demonstrate the feasibility and robustness of the proposed framework. Significance. This study offers the proof of principle required to implement this methodology to uncover VIM thalamic LFP-EMG interactions for (i) better understanding of the pathophysiology of tremor; (ii) objective selection of the DBS electrode contacts with the highest strength of association with the tremorogenic EMG, a particularly useful feature for the implementation of novel multicontact directional leads in clinical practice; and (iii) future research on DBS closed-loop devices
