from functools import reduce import pandas as pd import statsmodels.api as sm import numpy as np from scipy.optimize import minimize df = pd.read_csv("pheno.asd", sep="\s+", header=0) gmat = np.loadtx