请写一段多元线性回归的python代码，中文注释： 采用多因变量线性回归， 以关怀总分，意愿总分为因变量，其他变量为自变量。

# 导入相关库
import pandas as pd
import numpy as np
import statsmodels.api as sm

# 读取数据
data = pd.read_csv('data.csv')

# 定义自变量和因变量
X = data[['变量1', '变量2', '变量3']]
y = data[['关怀总分', '意愿总分']]

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

# 构建模型
model = sm.OLS(y, X).fit()

# 输出模型结果
print(model.summary())