Electroencephalography (EEG) recordings are often contaminated by muscle artifacts. In the literature, a number of methods have been proposed to deal with this problem. Yet most denoising muscle artifact methods are designed for either single-channel EEG or hospital-based, high-density multichannel recordings, not the few-channel scenario seen in most ambulatory EEG instruments. In this paper, we propose utilizing interchannel dependence information seen in the few-channel situation by combining multivariate empirical mode decomposition and canonical correlation analysis (MEMD-CCA). The proposed method, called MEMD-CCA, first utilizes MEMD to jointly decompose the few-channel EEG recordings into multivariate intrinsic mode functions (IMFs). Then, CCA is applied to further decompose the reorganized multivariate IMFs into the underlying sources. Reconstructing the data using only artifact-free sources leads to artifact-attenuated EEG. We evaluated the performance of the proposed method through simulated, semisimulated, and real data. The results demonstrate that the proposed method is a promising tool for muscle artifact removal in the few-channel setting.