用python写一个能输出含有16个自变量的回归方程的的多元线性回归程序

可以使用Python的statsmodels库来实现多元线性回归。以下是一个可以输出含有16个自变量的回归方程的程序示例：

import numpy as np
import pandas as pd
import statsmodels.api as sm

# 读入数据，其中y是因变量，X是自变量
data = pd.read_csv('data.csv')
y = data['y']
X = data.iloc[:,1:]

# 添加常数列
X = sm.add_constant(X)

# 构建多元线性回归模型
model = sm.OLS(y, X).fit()

# 输出回归系数和截距
print(model.params)

请注意，以上代码仅仅是输出了回归系数和截距，如果需要输出完整的回归方程，需要将上述代码进行修改。