Я хочу запустить множественную линейную регрессию на python, используя statsmodels.Мне удалось научиться делать простую линейную регрессию, используя следующий код.
import pandas as pd
import numpy as np
from scipy import stats
import statsmodels.api as sm
import matplotlib.pyplot as plt
from IPython.html.widgets import interact
from IPython.html import widgets
import matplotlib
import statsmodels.api as sm
import statsmodels.formula.api as smf
Y = [1,3,4,5,2,3,4]
X = range(1,8)
Z=np.array([10.0, 8.0, 13.0, 9.0, 11.0, 14.0,1.0])+X
X = sm.add_constant(X)
model = sm.OLS(Y,X)
results = model.fit()
results.params
results.tvalues
print(results.summary())
Итак, теперь я хочу запустить модель, например: Y = a1X + a2Z + константа.Как я могу это сделать, следуя логике, которая у меня есть до сих пор в моем коде?