import statsmodels.formula.api as sm
model = sm.ols(formula = 'output ~ input1 + input2', data=df).fit()
print(model.summary())
corr = df.corr()
plt.figure(figsize=(10,5))
sns.heatmap(corr, annot=True, fmt='.1f', cmap='Blues')
import scipy.stats as stats
stats.pearsonr(df['a], df['b'])
stats.spearmanr(df['a'], df['b'])
stats.kendalltau(df['a'], df['b'])