Thanks! So using the data file I found online https://github.com/jodyabney/DBDA2Eprograms/blob/master/OrdinalProbitData-1grp-1.csv, I would start your code with:
# Load data
df = pd.read_csv('OrdinalProbitData-1grp-1.csv').rename(columns={'Y': 'y'})-1
df['intercept'] = 1