A methodology for modeling the Far Field (FF) radiated by antennas subject to random variabilities with surrogate models of high efficiency is presented herein. The FF is expanded onto Vector Spherical Harmonics, then the stochastic variability of their coefficients is quantitatively modeled with the Polynomial Chaos expansion. The accuracy of the surrogate model is guaranteed by a robust and efficient adaptive procedure; “extreme parsimony” of the model is achieved thanks to various “compression” techniques; moreover, the characterization of the antenna’s frequency behavior over a band of interest has been taken into account in the modeling for the first time. The methodology is applied to a practical wearable textile patch antenna designed to cover the 2.4 – 2.5 GHz band and which is subject to five geometric and material random parameters. Specifically, the “Weierstrass” curves crumpling effect is modeled for the first time. Comparison to the full-wave simulations shows that the derived surrogate model predicts the FF with good accuracy and with a speed up factor of 105. Such type of surrogate models could be beneficial not only for antenna design and optimization purpose, but also, e.g., for joint antenna-channel stochastic analyses.