SKiNN’s role in a joint modeling framework. SKiNN takes, as input, a parametric description of the lens mass profile, lens light profile, anisotropy of the stellar orbits in the lens galaxy, and lens inclination and outputs a υ_rms map used to evaluate a kinematic likelihood (). This likelihood is then combined with the lensing likelihood () to give a total likelihood (), which is then sampled. In the case of a variable lensed source with measured time delays, a time-delay likelihood (, not shown in the figure) can be added to the inference (see Sect. 5).