https://stats.stackexchange.com/a/591651/69508

591651_69508.py

import numpy as np
from scipy.stats import t
import statsmodels.api as sm

# Prepare data
x = np.array([1,2,3,3,4,5,5])
y = np.array([3,7,5,8,11,14,12])
n = y.size
X = sm.add_constant(x)

# Define and train
model = sm.OLS(y, X)
results = model.fit()
print(results.summary())

# t-test
t_obs = (results.params[1] - 2) / results.bse[1]
print(1 - t.sf(t_obs, n-2))