Background Non-steroidal anti-inflammatory drugs (NSAIDs) are considered first-line medications for acute migraine attacks. However, the response exhibits considerable variability among individuals. Thus, this study aimed to explore a machine learning model based on the percentage of amplitude oscillations (PerAF) and gray matter volume (GMV) to predict the response to NSAIDs in migraine treatment.Methods Propensity score matching was adopted to match patients having migraine with response and nonresponse to NSAIDs, ensuring consistency in clinical characteristics and migraine-related features. Multimodal magnetic resonance imaging was employed to extract PerAF and GMV, followed by feature selection using the least absolute shrinkage and selection operator regression and recursive feature elimination algorithms. Multiple predictive models were constructed and the final model with the smallest predictive residuals was chosen. The model performance was evaluated using the area under the receiver operating characteristic (ROCAUC) curve, area under the precision-recall curve (PRAUC), balance accuracy (BACC), sensitivity, F1 score, positive predictive value (PPV), and negative predictive value (NPV). External validation was performed using a public database. Then, correlation analysis was performed between the neuroimaging predictors and clinical features in migraine.Results One hundred eighteen patients with migraine (59 responders and 59 non-responders) were enrolled. Six features (PerAF of left insula and left transverse temporal gyrus; and GMV of right superior frontal gyrus, left postcentral gyrus, right postcentral gyrus, and left precuneus) were observed. The random forest model with the lowest predictive residuals was selected and model metrics (ROCAUC, PRAUC, BACC, sensitivity, F1 score, PPV, and NPV) in the training and testing groups were 0.982, 0.983, 0.927, 0.976, 0.930, 0.889, and 0.973; and 0.711, 0.648, 0.639, 0.667,0.649, 0.632, and 0.647, respectively. The model metrics of external validation were 0.631, 0.651, 0.611, 0.808, 0.656, 0.553, and 0.706. Additionally, a significant positive correlation was found between the GMV of the left precuneus and attack time in non-responders.Conclusions Our findings suggest the potential of multimodal neuroimaging features in predicting the efficacy of NSAIDs in migraine treatment and provide novel insights into the neural mechanisms underlying migraine and its optimized treatment strategy.
