I made this nls model
formula <- nls(y ~ a+(x/500)^b, data=my_data, start =list(a=1, b=1))
then I got the result below,
summary(formula)
Param