Corticomuscular coupling analysis, i.e., examining the relations between simultaneously recorded brain (e.g., electroencephalography--EEG) and muscle (e.g., electro-myography-EMG) signals, is a useful tool for understanding aspects of human motor control. Traditionally, the most popular method to assess corticomuscular coupling has been the pairwise magnitude-squared coherence (MSC) between EEG and concomitant EMG recordings. In this paper, we propose assessing corticomuscular coupling by combining partial least squares (PLS) and independent component analysis (ICA), which addresses many of the limitations of MSC, such as difficulty in robustly assessing group inference and relying on the biologically implausible assumption of pairwise interaction between brain and muscle recordings. In the proposed framework, response relevance and statistical independence are jointly incorporated into a multiobjective optimization function to meaningfully combine the goals of PLS and ICA under the same mathematical umbrella. Simulations, performed under realistic assumptions, illustrated the utility of such an approach. The method was extended to address intersubject variability to robustly discover common corticomuscular coupling patterns across subjects. We then applied the proposed framework to concurrent EEG and EMG data collected in a Parkinson's disease (PD) study. The results from applying the proposed technique revealed temporal components in the EEG and EMG that were significantly correlated with one another. In addition to the expected motor areas, the corresponding spatial activation patterns demonstrated enhanced occipital connectivity in PD subjects, consistent with previous studies suggesting that PD subjects rely excessively on visual information to counteract the deficiency in being able to generate internal commands from their affected basal ganglia.