I ran a test case with a 4:1 xmfr after adding resistive loss to the xmfr.
A plot of return loss (from s11) and then transmission loss (from s21) is provided.
Looks like 1/(1-|s11|) where S11 is obtained from the return loss value matches up with the calculated transmission loss.