simple_regression.py

simple_regression.py

import numpy as np

def beta(x, y):
    mean_x = np.mean(x)
    mean_y = np.mean(y)
    numerator = 0
    for pair in zip(x, y):
        numerator += (pair[0] - mean_x) * (pair[1] - mean_y)
    denominator = np.sum([np.square(item - mean_x) for item in x])
    return numerator / denominator