矩阵OLS估计验证

import pandas as pd
import numpy as np

def solve():
    url = '____'
    df = pd.read_excel(____)
    n = len(____)
    X = np.column_stack([np.ones(____), df['____'].values])
    y = df['____'].values
    beta = np.linalg.inv(____) @ X.T @ y
    print(____)