We simulated the spatiotemporal modelling of 3D stress and strain distributions during the triaxial compression laboratory test on a westerly granite sample using finite-difference numerical modelling implemented with FLAC3D software. The modelling was performed using a ubiquitous joint constitutive law with strain softening. The applied procedure is capable of reproducing the macroscopic stress and strain evolution in the sample during triaxial deformation until a failure process occurs. In addition, we calculated focal mechanisms of acoustic emission (AE) events and resolved local stress field orientations. This detailed stress information was compared with that from numerical modelling. The comparison was made based on the 3D rotation angle between the cardinal axes of the two stress tensors. To infer the differences in rotation, we applied ANOVA. We identified the two time levels as the plastic deformation phase and the after-failure phase. Additionally, we introduced the bin factor, which describes the location of the rotation scores in the rock sample. The p values of the test statistics F for the bin and phase effects are statistically significant. However, the interaction between them is insignificant. We can, therefore, conclude that there was a significant difference in the time between the rotation means in the particular bins, and we ran post hoc tests to obtain more information where the differences between the groups lie. The largest rotation of the stress field provided by the focal mechanisms of AE events from the numerically calculated stress field is observed in the edge bins, which do not frame the damage zone of the sample.