Recording the activity of large populations of neurons requires new methods to analyze and use the large volumes of time series data thus created. Fast and clear methods for finding functional connectivity are an important step toward the goal of understanding neural processing. This problem presents itself readily in somatosensory neuroprosthesis (SSNP) research, which uses microstimulation (MiSt) to activate neural tissue to mimic natural stimuli, and has the capacity to potentiate, depotentiate, or even destroy functional connections. As the aim of SSNP engineering is artificially creating neural responses that resemble those observed during natural inputs, a central goal is describing the influence of MiSt on activity structure among groups of neurons, and how this structure may be altered to affect perception or behavior. In this paper, we demonstrate the concept of Granger causality, combined with maximum likelihood methods, applied to neural signals recorded before, during, and after natural and electrical stimulation. We show how these analyses can be used to evaluate the changing interactions in the thalamocortical somatosensory system in response to repeated perturbation. Using LFPs recorded from the ventral posterolateral thalamus (VPL) and somatosensory cortex (S1) in anesthetized rats, we estimated pair-wise functional interactions between functional microdomains. The preliminary results demonstrate input-dependent modulations in the direction and strength of information flow during and after application of MiSt. Cortico-cortical interactions during cortical MiSt and baseline conditions showed the largest causal influence differences, while there was no statistically significant difference between pre- and post-stimulation baseline causal activities. These functional connectivity changes agree with physiologically accepted communication patterns through the network, and their particular parameters have implications for both rehabilitation and brain—machine interface SSNP applications.