A patient's responses to Transcranial Magnetic Stimulation (TMS) pulses on the motor cortex have a complex spatial pattern, making it challenging to understand the response patterns across multiple patients. We developed a mixture of Dirichlet process models to distinguish between patient-specific and shared spatial patterns across multiple patients to provide insight into consistent response patterns essential for developing personalized treatment procedures. The Metropolis-Hastings within Gibbs sampler of the Markov Chain Monte Carlo algorithm was developed for estimation. The model was used to analyse the TMS data of 3 healthy subjects. The study revealed that the primary motor cortex of the hand consistently emerges as a promising region for eliciting optimal responses. This area serves as a key target for brain mapping using TMS to identify cortical hotspots. However, the excitability patterns in this region can vary significantly among patients. C1 [Egbon, Osafu Augustine] Univ Fed Sao Carlos, Dept Stat, Sao Carlos, Brazil. [Egbon, Osafu Augustine; Louzada, Francisco] Univ Sao Paulo, Inst Math & Comp Sci, 400 Ave Trabalhador Sao carlense, BR-13566590 Sao Carlos, SP, Brazil. [Egbon, Osafu Augustine; Heumann, Christian] Ludwig Maximilians Univ Munchen, Dept Stat, Munich, Germany. [Nascimento, Diego Carvalho] Univ Atacama, Dept Math, Copiapo, Chile. C3 Universidade Federal de Sao Carlos; Universidade de Sao Paulo; University of Munich; Universidad de Atacama