Even if it has more parameters (4 x 4 matrix) insted of 6 scalars, it leads to much more cleaner code, a little faster as no conversions to homogeneous matrix are needed, and much more stable as no euler angles are used, avoiding conventions/representations/singularities.