Conpensate for delay from a sgolay filter

I'm creating an sgolay filter using the sgolay function. And applying it afterwards using the https://www.mathworks.com/help/signal/ref/sgolay.html algorithm that allows for calculating Savitzky-Golay Differentiation.. I'm trying to calculate smoothed velocity/acceleration given a set of positions over time (coordinates duplets) ... although I now need to compensate for delay introduced by the sgolay filter applied using the conv function.

I've seen some methods using the method: delay = mean(grpdelay(firf)) --> where firf is the filter matrix sos.

But I'm not able to do the same for the sgolay filter because I don't have that matrix.. what do I need to do in order to compensate delay introduced by the filter?