Y_average=[222.0, 222.667, 226.0, 227.0, 223.333, 223.333, 226.667, 224.333]
x0_factor=[1, 1, 1, 1, 1, 1, 1, 1]
x1_factor=[-1, -1, 1, 1, -1, -1, 1, 1]
x2_factor=[-1, 1, -1, 1, -1, 1, -1, 1]
x3_factor=[-1, 1, 1, -1, 1, -1, -1, 1]
x1x2_factor=[1, -1, -1, 1, 1, -1, -1, 1]
x2x3_factor=[1, 1, -1, -1, -1, -1, 1, 1]
x1x3_factor=[1, -1, 1, -1, -1, 1, -1, 1]
x1x2x3_factor=[-1, -1, -1, -1, 1, 1, 1, 1]
N=8

beta_list_test,beta_list=[],[]
def beta_calc():
    for i in range(len(x0_factor)):
        beta_list_test.append((Y_average[i]*x0_factor[i])/N)
    beta_list.append(sum(beta_list_test))
beta_calc()

print(f"beta_list_test:{beta_list_test}")
print(f"beta_list:{beta_list}")