Currently, research predominantly focuses on evaluating clinical effects at specific time points while neglecting underlying patterns within the treatment process. This study aims to analyze the dynamic alterations in PANSS total scores and prolactin levels in patients with schizophrenia treated with risperidone, along with the influencing covariates. Using data from an 8-week randomized, double-blind, multicenter clinical trial, a population pharmacodynamic model was established for the PANSS total scores of and prolactin levels in patients treated with risperidone. The base model employed was the E_max model. Covariate selection was conducted using a stepwise forward inclusion and backward elimination approach. A total of 144 patients were included in this analysis, with 807 PANSS total scores and 531 prolactin concentration values. The PANSS total scores of the patients treated with risperidone decreased over time, fitting a proportionally parameterized sigmoid E_max model with covariates including baseline score, course of the disease, gender, plasma calcium ions, and lactate dehydrogenase levels. The increase in prolactin levels conformed to the ordinary E_max model, with covariates encompassing course of the disease, gender, weight, red blood cell count, and triglyceride levels. The impacts of the baseline scores and the course of the disease on the reduction of the PANSS scores, as well as the influence of gender on the elevation of prolactin levels, each exceeded 20%. This study provides valuable quantitative data regarding PANSS total scores and prolactin levels among patients undergoing risperidone treatment across various physiological conditions.