plt.plot(x, y)

log_x = [np.log(xi) for xi in x]
log_y = [np.log(yi) for yi in y]

plt.plot(log_x, log_y)
plt.grid()

np.polyfit(log_x, log_y, 1)

array([-1.00082759,  0.32878786])