There are well known identities relating elementary symmetric polynomials in $k$ variables and complete homogeneous symmetric polynomials in the same set of variables. For example,
$$e_{2,k}(x) = h_{1,k}(x)^2 - h_{2,k}(x).$$
I'd like to do the same but for the factorial versions that also have $y$ variables. There is some subtlety that needs to be worked out, in particular the $y$ variables need to be altered to make the same identities work. For example, $$E_{2,4}(x_1,x_2,x_3,x_4;y_1,y_2,y_3)=H_{1, 4}(x_1, x_2, x_3, x_4; y_1, y_2, y_3, y_5)H_{1, 4}( x_1, x_2, x_3, x_4; y_1, y_2, y_3, y_4) - H_{2,4}(x_1,x_2,x_3,x_4;y_1,y_2,y_3,y_4,y_5).$$
I have a systematic way of computing this programmatically, I'm just looking for a formulaic way rather than algorithmic, which would very likely be much faster.