How to get all S-parameters from a circuit?

How can one get all four S-parameters of an active two-port circuit?
I routinely do calculations like (S21(V1)-S12(V1))/(1-S11(V1) * S22(V1)+S21(V1) * S12(V1)-2 * S12(V1)) (Randall/Hock mismatched gain correction) using .net directive in LTSpice.
Is there a similar trick in QSpice?
I am aware of simple S11 and S21 solution.
Leo