By using the Kronecker product,An identification model for multivariable ARX-like stochastic systems is derived and developed a hierarchical least squares parameter estimation algorithm by the hierarchical identification principle.The proposed algorithm has less computational eorts than the recursive least squares algorithm.A simulation example is included.