# -*- coding: utf-8 -*-
"""
Created on Tue Nov 30 13:13:35 2021

@author: Sim
"""
import numpy as np
import pandas as pd
import statsmodels.api as sm
from statsmodels.formula.api import ols

iq = np.array([100 ,123 ,132,  98 ,105 ,110 ,117 , 95, 127 ,135])
vt =np.array([89 , 93 , 97,  80 , 83 , 88 , 90 , 70, 93 , 97])

df = pd.DataFrame({'iq': iq, 'vt': vt})
lm1 = ols('iq~vt', data=df).fit()
print(lm1.conf_int())

vt1 = sm.add_constant(vt)
lm2 = sm.OLS(iq, vt1).fit()
print(lm2.conf_int())