In a recent paper, we proposed a system identification method for constructing reduced-order models for the electromechanical dynamics of large power systems, divided into multiple coherent clusters, using Synchrophasors. Every cluster in the actual model was represented as an aggregate generator in the reduced-order model. An aggregate network graph connected one aggregate generator to another. In this paper, we extend this identification approach to differential-algebraic (DAE) models. First, every cluster is associated with a unique terminal bus, referred to as the pilot bus, that couples its internal network to the rest of the system. The proposed algorithm uses Synchrophasor measurements from the pilot buses to identify the dynamic model of the aggregate generator for each cluster using nonlinear least squares while retaining the identity of all the pilot buses. The resulting reduced-order model is in the form of a nonlinear electric circuit described by aggregate differential and algebraic equations. We illustrate our results using two case studies, one for the IEEE 9-bus power system and another for the IEEE 39-bus power system. We also discuss how these reduced-order DAE models may be useful for designing shunt controllers at the pilot buses by using Synchrophasor feedback.