Simples regressão linear múltipla

regression.py

# -*- coding: utf-8 -*-

import numpy as np

def mmq(x, y):
    u"""
    x é uma matriz com as variáveis independentes
    y é um array com a variável dependente 
    O retorno é um array contendo os valores dos beta's
    """
    x = np.insert(x, 0, 1, axis=1)
    x_t = np.transpose(x)
    xt_x = np.dot(x_t, x)
    inverse_xt_x = np.linalg.inv(xt_x)
    xt_y = np.dot(x_t, y)
    return np.dot(inverse_xt_x, xt_y)