The Kuramoto-Sakaguchi model for coupled phase oscillators with phase frustration is often studied in the thermodynamic limit of infinitely many oscillators. Here we extend a model reduction method based on collective coordinates to capture the collective dynamics of finite-size Kuramoto-Sakaguchi models. We find that the inclusion of the effects of rogue oscillators is essential to obtain an accurate description, in contrast to the original Kuramoto model, where we show that their effects can be ignored. We further introduce a more accurate ansatz function to describe the shape of synchronized oscillators. Our results from this extended collective coordinate approach reduce in the thermodynamic limit to the well-known mean-field consistency relations. For finite networks we show that our model reduction describes the collective behavior accurately, reproducing the order parameter, the mean frequency of the synchronized cluster, and the size of the cluster at a given coupling strength, as well as the critical coupling strength for partial and for global synchronization.