Curve fitting -- including convexity constraint

You’re not wrong. In a 2010 paper (linked above), @avehtari proposes using a probit likelihood to enforce the constraint and that’s what I have tried to replicate here, without much success.