We present a novel technique to perform the model-order reduction (MOR) of multiport systems under the effect of statistical variability of geometrical or electrical parameters. The proposed approach combines a deterministic MOR phase with the use of the Polynomial Chaos (PC) expansion to perform the variability analysis of the system under study very efficiently. The combination of MOR and PC techniques generates a final reduced-order model able to accurately perform stochastic computations and variability analysis. The novel proposed method guarantees a high-degree of flexibility, since different MOR schemes can be used and different types of modern electrical systems (e. g., filters and connectors) can be modeled. The accuracy and efficiency of the proposed approach is verified by means of two numerical examples and compared with other existing variability analysis techniques.