We analyze the phase diagram of a topological insulator model including antiferromagnetic interactions in the form of an extended Su-Schrieffer Heeger model. To this end, we employ a recently introduced operational definition of topological order based on the ability of a system to perform topological error correction. We show that the necessary error correction statistics can be obtained efficiently using a Monte-Carlo sampling of a matrix product state representation of the ground state wave function. Specifically, we identify two distinct symmetry-protected topological phases corresponding to two different fully dimerized reference states. Finally, we extend the notion of error correction to classify thermodynamic phases to those exhibiting local order parameters, finding a topologically trivial antiferromagnetic phase for sufficiently strong interactions.