A structured matrix factorization method for computational modeling of hierarchical polarization in social interactions