We study the impact of including the baryonic decay Λb→Λ(→pπ-)μ+μ- in a Bayesian analysis of |ΔB|=|ΔS|=1 transitions. We perform fits of the Wilson coefficients C9, C9′, C10 and C10′, in addition to the relevant nuisance parameters. Our analysis combines data for the differential branching fraction and three angular observables of Λb→Λ(→pπ-)μ+μ- with data for the branching ratios of Bs→μ+μ- and inclusive b→sℓ+ℓ- decays. Newly available precise lattice QCD results for the full set of Λb→Λ form factors are used to evaluate the observables of the baryonic decay. Our fits prefer shifts to C9 that are opposite in sign compared to those found in global fits of only mesonic decays, and the posterior odds show no evidence of physics beyond the Standard Model. We investigate a possible hadronic origin of the observed tensions between theory and experiment.