marginaleffects
margins
mfx
emmeans
probemod
interactions
interactionR
sjPlot
τ(x,c)=E[Y1i−Y0i|Xi=x,Ci=c]