Linear predictive models are applied to functional MRI (fMRI) data to estimate boundaries that predict experimental task states for scans. These boundaries are visualized as statistical parametric maps (SPMs) and range from low to high spatial reproducibility across subjects (e.g., Strother , 2004; LaConte , 2003). Such inter-subject pattern reproducibility is an essential characteristic of interpretable SPMs that generalize across subjects. Therefore, we introduce a flexible hybrid model that optimizes reproducibility by simultaneously enhancing the prediction power and reproducibility. This hybrid model is formed by a weighted summation of the optimization functions of a linear discriminate analysis (LDA) model and a generalized canonical correlation (gCCA) model (Afshin-Pour , 2012). LDA preserves the model's ability to discriminate the fMRI scans of multiple brain states while gCCA finds a linear combination for each subject's scans such that the estimated boundary map is reproducible. The hybrid model is implemented in a split-half resampling framework (Strother , 2010) which provides reproducibility (r) and prediction (p) quality metrics. Then the model was compared with LDA, and Gaussian Naive Bayes (GNB). For simulated fMRI data, the hybrid model outperforms the other two techniques in terms of receiver operating characteristic (ROC) curves, particularly for detecting less predictable but spatially reproducible networks. These techniques were applied to real fMRI data to estimate the maps for two task contrasts. Our results indicate that compared to LDA and GNB, the hybrid model can provide maps with large increases in reproducibility for small reductions in prediction, which are jointly closer to the ideal performance point of (p=1, r=1).