This paper presents an improved gearbox fault diagnosis approach by integrating complementary ensemble empirical mode decomposition (CEEMD) with permutation entropy (PE). The presented approach identifies faults appearing in a gearbox system based on PE values calculated from selected intrinsic mode functions (IMFs) of vibration signals decomposed by CEEMD. Specifically, CEEMD is first used to decompose vibration signals characterizing various defect severities into a series of IMFs. Then, filtered vibration signals are obtained from appropriate selection of IMFs, and correlation coefficients between the filtered signal and each IMF are used as the basis for useful IMFs selection. Subsequently, PE values of those selected IMFs are utilized as input features to a support vector machine (SVM) classifier for characterizing the defect severity of a gearbox. Case study conducted on a gearbox system indicates the effectiveness of the proposed approach for identifying the gearbox faults.