import pandas as pd
import numpy as np
import statsmodels.formula.api as smf
import matplotlib.pyplot as plt
import yfinance as yf
stocks="KRZ.IR", "^ISEQ"
data = yf.download(stocks, start="2020-01-09", end="2022-01-09")
data = data[['Adj Close']]
data.columns = ["KRZ", "ISEQ"]
data.head()
returns = data.pct_change()
returns = returns.dropna()
returns.columns = ['rKRZ', "rISEQ"]
returns
# OLS regression find the relationship between the daily return of KRZ.IR and ISEQ
formula = "rKRZ ~ rISEQ"
results = smf.ols(formula, returns).fit()
print(results.summary())
# Null hypothesis test
# H 0 : β 1 = 1
# H 1 : β 1 ≠ 1
formula = "rKRZ ~ rISEQ"
hypotheses = "rISEQ = 1"
results = smf.ols(for