How to use fit_mle for Gamma

Hi all,
I have a series of deterioration data and the corresponding time for these measurements:

# deterioration
r = [0.1803, 0.1978, 0.2033, 0.2092, 0.2157, 0.2254]
# time
t = [0, 185, 348, 515, 658, 830]

The increments are:

t_inc = [185, 163, 167, 143, 172]
r_inc = [0.0175, 0.0055, 0.0059, 0.0065, 0.0097]

I need to estimate the shape and scale parameter of a Gamma process for the above-measured data.
I see there is a fit_mle function defined in Distributions that support Gamma, but how should I define the date for:
fit_mle function?
Thank you very much for checking!