Traditional theoretical models, while conservative, overlook the complex, nonlinear interactions between material and geometric parameters that influence wall flexural strength. Additionally, the extent to which these models maintain their conservative assumptions across varying parameters remains uncertain. This study aims to develop advanced predictive models for the flexural strength of planar composite plate shear walls—concrete filled (C-PSW/CF) by employing Linear Regression (LR) and Response Surface Methodology (RSM). Central Composite Design (CCD) was utilized to design numerical experiments to analyze the effects of varying parameters, including steel plate yield strength, concrete compressive strength, wall depth, depth-to-thickness ratio, and reinforcement ratio. A validated complex finite element modeling (FEM) procedure was used to analyze numerical experiments involving varying parameters of the selected variables. Both LR and RSM were applied to the generated dataset to formulate predictive equations for the maximum moment capacity of walls. Interaction plots were developed to reveal the nonlinear relationships between the variables and identify the most effective parameters influencing the wall strength. Comparisons between theoretical predictions, RSM, and LR models revealed that RSM provided the most accurate predictions, closely aligning with FEM results and effectively capturing the nonlinear interactions between wall parameters. The findings of this paper contribute to the development of more practical, accurate and optimized design methods for C-PSW/CF, ultimately improving both structural safety and efficiency.