the quasistationary (QS) probability distribution of small systems. The QS distribution (or [*Yaglom limit*]{}, as it is known in the probability literature), is the probability distribution at long times, conditioned on survival of the process [@QSS]. For the one-dimensional CP and allied models [@exact], and an activated random walker model [@sleepy], finite-size scaling analysis of numerically exact QS results on a sequence of lattice sizes yields good estimates for the critical point, exponents and moment ratios. In the present case, with four states per site, attaining the sizes required for a precise analysis appears to be very costly, computationally, and we shall merely attempt to obtain reasonable estimates for the phase boundary $\lambda_{c}(\mu)$.

As described in detail in [@exact], obtaining the QS distribution numerically requires (1) enumerating all configurations on a lattice of a given size; (2) enumerating all transitions between configurations, and their associated rates; and (3) using this information in an iterative procedure to generate the QS distribution. Once the latter is known, one may calculate properties such as the order parameter or lifetime. For small systems these quantities are smooth functions of the control parameter and show no hint of the critical singularity. It is known, however, 