This function implements the Phenomenological Renormalization Group (PRG) method, introduced by Meshulam et al. (2018, 2019). It also can implement the extended version of PRG based on connectivity.