12月1日-12月7日地铁站点日客流量预测结果：

由十二月份一号到七号的预测数据画出的图表观察得出，站点135与137的客流量普遍较高，地铁工作人员应该在人流量多的站点多安排值班人员进行巡逻，维护站点的现场秩序，以免发生人员拥挤和安全事故的发生。特别的，图中12月7号147站点出现不同寻常的高峰客流量，考虑意外因素等情况造成的客流量增多，这一天我们可以提前预警并采取相应的安保措施，保证当天地铁站点的秩序。预测结果仅作为参考一个权重值，还需要专家意见，按照一定的权重来计算。

首先，我将加载并查看提供的CSV文件以了解其结构和内容。然后，我将使用折线图来展示24小时的交通流量变化趋势。让我们开始吧。

python复制import pandas as pd# 加载CSV文件file_path = '/mnt/data/.csv'df = pd.read_csv(file_path)# 显示前几行数据以了解其结构df.head()

数据集包含多个字段，其中date_time和traffic_volume是我们关注的主要字段。date_time字段表示日期和时间，而traffic_volume字段表示交通流量。
接下来，我将使用折线图来展示24小时内的交通流量变化趋势。python复制import matplotlib.pyplot as plt# 将date_time列转换为datetime类型df['date_time'] = pd.to_datetime(df['date_time'])# 对数据进行排序以确保时间顺序df = df.sort_values('date_time')# 绘制折线图plt.figure(figsize=(15, 6))plt.plot(df['date_time'], df['traffic_volume'], marker='o', linestyle='-')plt.title('24小时内交通流量变化趋势')plt.xlabel('时间')plt.ylabel('交通流量')plt.grid(True)plt.xticks(rotation=45)plt.tight_layout()plt.show()

这是一个展示24小时内交通流量变化趋势的折线图。从图中可以看出，交通流量在不同时间段内有所波动。

为了更准确地展示24小时内的交通流量变化趋势，我需要先筛选出特定日期的数据。然后，我将再次绘制折线图。python复制# 筛选出特定日期的数据（例如最早的日期）specific_date = df['date_time'].dt.date.iloc[0]df_specific_date = df[df['date_time'].dt.date == specific_date]# 对数据进行排序以确保时间顺序df_specific_date = df_specific_date.sort_values('date_time')# 绘制折线图plt.figure(figsize=(15, 6))plt.plot(df_specific_date['date_time'], df_specific_date['traffic_volume'], marker='o', linestyle='-')plt.title(f'{specific_date} 交通流量变化趋势')plt.xlabel('时间')plt.ylabel('交通流量')plt.grid(True)plt.xticks(rotation=45)plt.tight_layout()plt.show()

这是一个展示特定日期（最早的日期）内24小时交通流量变化趋势的折线图。从图中可以看出，交通流量在不同时间段内有所波动。

我将按照以下步骤进行操作：

1. 使用散点图矩阵来探索traffic_volume与weather_main、temp、rain_1h、snow_1h和clouds_all等变量之间的关系。
2. 使用柱状图或箱线图来比较holiday为“None”与特定假日（如Columbus Day, Veterans Day, Thanksgiving Day）的traffic_volume差异。

首先，我将进行第一步的分析。

python复制import seaborn as sns# 选择相关的列columns_to_explore = ['traffic_volume', 'weather_main', 'temp', 'rain_1h', 'snow_1h', 'clouds_all']# 由于weather_main是分类变量，我们需要将其转换为数值编码df['weather_main_code'] = df['weather_main'].astype('category').cat.codes# 创建散点图矩阵sns.pairplot(df[columns_to_explore], diag_kind='kde', markers=["o", "s", "D", "v", "<", ">"])plt.suptitle('交通流量与天气条件的关系', y=1.02)plt.show()

这是一个散点图矩阵，用于探索traffic_volume与weather_main、temp、rain_1h、snow_1h和clouds_all等变量之间的关系。
接下来，我将进行第二步的分析，即比较holiday为“None”与特定假日（如Columbus Day, Veterans Day, Thanksgiving Day）的traffic_volume差异。python复制# 筛选出特定的假日specific_holidays = ['Columbus Day', 'Veterans Day', 'Thanksgiving Day']df_holidays = df[df['holiday'].isin(specific_holidays) | (df['holiday'] == 'None')]# 创建箱线图plt.figure(figsize=(12, 6))sns.boxplot(x='holiday', y='traffic_volume', data=df_holidays, palette=["rgb(202,206,210, max =255)", "rgb(230,232,234, max =255)"])plt.title('假日与非假日交通流量对比')plt.xlabel('假日类型')plt.ylabel('交通流量')plt.show()python复制# 修正颜色格式并重新创建箱线图plt.figure(figsize=(12, 6))sns.boxplot(x='holiday', y='traffic_volume', data=df_holidays, palette=["#C0C4C8", "#E6E8EA"])plt.title('假日与非假日交通流量对比')plt.xlabel('假日类型')plt.ylabel('交通流量')plt.show()

这是一个箱线图，用于比较holiday为“None”与特定假日（如Columbus Day, Veterans Day, Thanksgiving Day）的traffic_volume差异。
从图中可以看出，不同类型的假日和非假日的交通流量有所不同。

以下是我们实验过程的详细描述：

首先，我们对特征集X和目标变量y进行了训练集和测试集的划分，其中测试集占比为25%，并确保了数据的随机性，以避免任何潜在的偏差。具体实现如下：

X_train, X_test, y_train, y_test = train_test_split(X,                                                     y,                                                     test_size=0.25,                                                     shuffle=True,                                                     random_state=42)tpot = TPOTRegressor(max_time_mins=60,                     verbosity=2,                     n_jobs=-1)tpot.fit(X_train,y_train)tpot.export('metro_traffic.py')

接着，我们初始化并训练了一个TPOTRegressor模型，设置最大运行时间为60分钟，以确保模型有足够的时间探索可能的解决方案。此外，我们将verbosity设置为2，以便在训练过程中获得详细的输出，并将n_jobs设置为-1，以利用所有可用的处理器核心。

在模型训练完成后，我们使用TPOT生成的模型对测试集进行了预测，并通过R²、平均绝对误差（MAE）、均方误差（MSE）和均方根误差（RMSE）等指标对模型性能进行了评估。

y_predictions = tpot.predict(X_test)r2 = sklearn.metrics.r2_score(y_test, y_predictions)mae = sklearn.metrics.mean_absolute_error(y_test, y_predictions)mse = sklearn.metrics.mean_squared_error(y_test, y_predictions)rmse = np.sqrt(mse)

此外，我们还尝试了其他几种机器学习模型，包括随机森林回归器和极限梯度提升回归器（XGBRegressor），并对其进行了参数调优以优化性能。以下是随机森林回归器的一个示例：

training_features, testing_features, training_target, testing_target = \            train_test_split(features, tpot_data['traffic_volume'].values, random_state=None)exported_pipeline = make_pipeline(    StandardScaler(),    RandomForestRegressor(bootstrap=False, max_features=0.45, min_samples_leaf=2, min_samples_split=2, n_estimators=100))

exported_pipeline = RandomForestRegressor(bootstrap=True, max_features=0.9000000000000001, min_samples_leaf=9, min_samples_split=19, n_estimators=100)exported_pipeline.fit(training_features, training_target)results = exported_pipeline.predict(testing_features)

features = tpot_data.drop('traffic_volume', axis=1).valuestraining_features, testing_features, training_target, testing_target = \            train_test_split(features, tpot_data['traffic_volume'].values, random_state=None)exported_pipeline = ExtraTreesRegressor(bootstrap=False, max_features=0.8, min_samples_leaf=8, min_samples_split=18, n_estimators=100)exported_pipeline.fit(training_features, training_target)results = exported_pipeline.predict(testing_features)

我们还尝试了极限梯度提升回归器，并对其进行了参数调整：

a=tpot_data['traffic_volume'].valuesfeatures = tpot_data.drop('traffic_volume', axis=1).valuestraining_features, testing_features, training_target, testing_target = \            train_test_split(features, tpot_data['traffic_volume'].values, random_state=None)# Average CV score on the training set was:-259017.37025310827exported_pipeline = XGBRegressor(learning_rate=0.1, max_depth=3, min_child_weight=15, n_estimators=100, nthread=1, subsample=0.6000000000000001,objective = 'reg:squarederror')exported_pipeline.fit(training_features, training_target)results = exported_pipeline.predict(testing_features)

在每种情况下，我们都记录了模型在训练集上的平均交叉验证分数，并使用训练好的模型对测试集进行了预测。

通过这些实验，我们能够识别出最适合我们数据集的模型，并对其进行了详细的性能评估。这些结果为我们提供了有价值的见解，有助于进一步优化交通流量预测模型，并为城市规划和交通管理部门提供数据支持。

关于作者

Xinyi He

✉

     联系我们

我们对Xinyi He表示诚挚的感谢。她专注于数据科学与大数据技术。Xinyi He在数据分析与数据挖掘、Hadoop大数据平台部署、Spark基础编程等领域具有扎实的理论基础和实践经验。她对机器学习、深度学习以及大数据平台搭建等前沿技术有着深入的理解和应用能力。