We present a new method for the construction of parametric macromodels for admittance and impedance input- output representations starting from multivariate data samples that depend on frequency and additional design variables such as geometric and material parameters. Poles and residues are parameterized indirectly, while stability and passivity of the parametric macromodel are guaranteed over a user defined range of design parameter values. Pertinent numerical examples validate the proposed parametric macromodeling approach.