A Symbolic Algorithm for Checking the Identifiability of a Time-Series Model

Several authors have attempted to compute the asymptotic Fisher information matrix for a univariate or multivariate time-series model to check for its identifiability. This has the form of a contour integral of a matrix of rational functions. A recent paper has proposed a short Wolfram Mathematica notebook for VARMAX models that makes use of symbolic integration. It cannot be used in open-source symbolic computation software like GNU Octave and GNU Maxima. It was based on symbolic integration but the integrand lacked symmetry characteristics in the appearance of polynomial roots smaller or greater than 1 in modulus. A more symmetric form of the integrand is proposed for VARMA models that first allows a simpler approach to symbolic integration. Second, the computation of the integral through Cauchy residues is also possible. Third, an old numerical algorithm by Söderström is used symbolically. These three approaches are investigated and compared on a pair of examples, not only for the Wolfram Language in Mathematica but also for GNU Octave and GNU Maxima. As a consequence, there are now sufficient conditions for exact model identifiability with fast procedures.