#建模分析 import pandas as pd import numpy as np import statsmodels.api as sm from sklearn.model_selection import train_test_split from sklearn.metrics import mean_squared_error from sklearn.ensemble import GradientBoostingRegressor # 去除异常值 diabetes = diabetes[(diabetes['bmi'] > 10) & (diabetes['HbA1c_level'] < 15)] # 划分训练集和测试集 train, test = train_test_split(diabetes, test_size=0.3, random_state=42) # 构建线性回归模型 xtrain, ytrain = train.drop('diabetes', axis=1), train['diabetes'] xtest, ytest = test.drop('diabetes', axis=1), test['diabetes'] Xtrain = sm.add_constant(xtrain) Xtest = sm.add_constant(xtest) print(diabetes.info()) reg = sm.OLS(ytrain, Xtrain).fit() print(reg.summary()) # 计算线性回归的预测误差 ypred = reg.predict(Xtest) mse = mean_squared_error(ytest, ypred) rmse = np.sqrt(mse) print('Linear Regression RMSE:', rmse) # 构建GBDT模型 gbdt = GradientBoostingRegressor(learning_rate=0.3).fit(xtrain, ytrain) print('GBDT R^2:', gbdt.score(xtrain, ytrain)) # 计算GBDT的预测误差 ypred = gbdt.predict(xtest) mse = mean_squared_error(ytest, ypred) rmse = np.sqrt(mse) print('GBDT RMSE:', rmse)

这段代码是一个用于糖尿病数据集的建模分析，主要使用了线性回归和GBDT（梯度提升决策树）两种模型进行预测。在代码中，首先通过去除异常值的方法对数据进行预处理，然后将数据集划分为训练集和测试集。接着，使用Statsmodels库中的OLS函数构建线性回归模型，并计算线性回归的预测误差。同时，使用sklearn库中的GradientBoostingRegressor函数构建GBDT模型，并计算GBDT的预测误差。最后，输出线性回归和GBDT模型的预测误差。