Electroencephalography (EEG) and simultaneously-recorded electromyography (EMG) data are a means to assess integrity of the functional connection between the cortex and the muscle during movement. EEG-EMG coupling is typically assessed with pair-wise squared coherence, resulting in a small, but statistically-significant coherence between a single EEG and a single EMG channel. However, a means to combine results across subjects is not straightforward with this approach because the exact frequency of maximal EEG-EMG coupling may vary between individuals, and it emphasizes the role of an individual locus in the brain in driving the muscle activity, when interactions between head regions may in fact be more influential on ongoing EMG activity. To deal with these issues, we implemented a multiblock Partial Least Squares (mbPLS) procedure, previously proposed in chemical applications, which incorporates a hierarchical structure into the ordinary two-block PLS often used in neuroimaging studies. In the current implementation, each subject's data features are collected in individual data blocks on a sub-level, while simultaneously aggregating the sub-level information to obtain a super-level group "consensus". We further extended the mbPLS model to include 3-dimensional matrices: time-frequency-EEG channel and a time-frequency-connection utilizing Partial Directed Coherence (PDC). We applied the proposed method to concurrent EEG and EMG data collected from ten normal subjects and nine patients with mild-moderate Parkinson's disease (PD) performing a dynamic motor task-that of sinusoidal squeezing. The results demonstrate that connections between EEG electrodes, rather than activity at individual electrodes, correspond more closely to ongoing EMG activity. In PD subjects, there was enhanced connectivity to and from occipital regions, likely related to the previously-described enhanced use of visual information during motor performance in this group. The proposed mbPLS framework is a promising technique for performing multi-subject, multi-modal data analysis and it allows for robust group inferences even in the face of large inter-subject variability.