On PC solutions of PARMA (<i>p</i>, <i>q</i>) models