We propose a novel semi-automatic approach to design biomarkers for capturing pharmacodynamic effects induced by pharmacological agents on the spectral power of electroencephalography (EEG) recordings. We apply this methodology to investigate the pharmacodynamic effects of methylphenidate (MPH) and atomoxetine (ATX) on attention deficit/hyperactivity disorder (ADHD), using rodent models. We inject the two agents into the spontaneously hypertensive rat (SHR) model of ADHD, the Wistar-Kyoto rat (WKY), and the Wistar rat (WIS), and record their EEG patterns. To assess individual EEG patterns quantitatively, we use an integrated methodological approach, which consists of calculating the mean, slope and intercept parameters of temporal records of EEG spectral power using a smoothing filter, outlier truncation, and linear regression. We apply Fisher discriminant analysis (FDA) to identify dominant discriminants to be heuristically consolidated into several new composite biomarkers. Results of the analysis of variance (ANOVA) and t-test show benefits in pharmacodynamic parameters, especially the slope parameter. Composite biomarker evaluation confirms their validity for genetic model stratification and the effects of the pharmacological agents used. The methodology proposed is of generic use as an approach to investigating thoroughly the dynamics of the EEG spectral power.