df['date_delta'] = (df['date'] - df['date'].min()) / np.timedelta64(1, 'D')
city_data = df[df['city'] == 'London']
result = sm.ols(formula='sales ~ date_delta', data=city_data).fit()
Преобразование даты в число с плавающей запятой для линейной регрессии во фрейме данных Pandas