To understand the complex tectonic interactions between subduction margins and transcurrent faults, we investigated the spatial distribution of stress orientations in the Philippines by using focal mechanisms derived from the waveform data of regional broadband seismic stations and Global Centroid Moment Tensor solutions. We investigated the spatial distribution of stress orientations by dividing the region into a central region containing the Philippine fault (PF) and regions to the east and west of the PF. Our stress tensor inversion results show that the sigma(1) axes in the central and eastern regions are WNW-ESE oriented and parallel to the relative plate motion orientations. However, in the western region, the orientations of sigma(1) axes were different. In particular, the orientations of the sigma(1) axes in the southern part (Bohol region) differ substantially (differences exceed approximately 60 degrees) from those of the relative plate motion. Furthermore, the orientations of the sigma(1) axes in the northern part (Mindoro and southwestern Luzon regions) differ by approximately 30 degrees from those found throughout the Philippines. The sigma(1) axes in the Bohol region and the sigma(3) axes in the Mindoro and southwestern Luzon regions are parallel to the strike of the PF. It was demonstrated that our numerical model incorporating slip on the PF and interplate coupling is able to successfully reproduce the features of the observed stress pattern.