d = {'a': [1, 2, 3, 4, 3, 3, 3], 'b': [5, 6, 7, 8, 4, 4, 4], 'c': [9, 10, 
    11, 12, 3, 3, 3], 'd': ['red', 'blue', 'green', 'red', 'orange', 'blue',
    'red']}
df = pd.DataFrame(d)
model = smf.ols('a ~ b + c + C(d, Treatment(reference="green"))', df).fit()
print(model.summary())
Relevel() и факторные переменные R в линейной регрессии в пандах