文章目录

• 📚引言
• 📖数据加载以及基本观察
• • 📑缺失值观察及处理
  • • 🔖缺失值观察以及可视化
    • 🔖缺失值处理
• 📖用户数据探索
• • 📑什么时间预定酒店将会更经济实惠？
  • 📑哪个月份的酒店预订是最繁忙的？
• 📖商家数据探索
• • 📑按市场细分的不同预定情况是怎样的？
  • 📑什么样的人更容易取消预订？
  • • 🔖数据编码
    • 🔖特征筛选
    • 🔖构建模型并预测
    • 🔖根据特征重要性得出结论
• 📍总结

📚引言

🙋‍♂️作者简介：生鱼同学，大数据科学与技术专业硕士在读👨‍🎓，曾获得华为杯数学建模国家二等奖🏆，MathorCup 数学建模竞赛国家二等奖🏅，亚太数学建模国家二等奖🏅。

✍️研究方向：复杂网络科学

🏆兴趣方向：利用python进行数据分析与机器学习，数学建模竞赛经验交流，网络爬虫等。

在有预定酒店的需求时，你是否考虑过以下的问题：

• 你是否曾想过每年什么时候是预订酒店房间的最佳时间？
• 为了获得最佳的每日价格，最佳的逗留时间？
• 如果你想预测一家酒店是否有可能收到过多的特殊要求呢？

在本文中我们就对Kaggle平台上的酒店需求数据集进行分析，这个数据集包含了一家城市酒店和一家度假酒店的预订信息，并包括诸如预订时间、逗留时间、成人、儿童和/或婴儿的数量，以及可用的停车位数量等信息。话不多说，我们开始吧。

本项目中的数据来源于Kaggle开放数据Hotel booking demand链接如下：
Hotel booking demand
需要的小伙伴可以自行下载获取。

📖数据加载以及基本观察

在进行数据加载之前，我们首先对数据的各个列进行解释，具体情况如下表所示：

可以看到，数据所给出的特征还是比较多的。接下来我们就基于上述数据来解决一些商业角度可能关心的问题，本文解决的问题如下：

首先是从用户的角度来看关心的问题：

• 什么时间预定酒店将会更经济实惠？
• 哪个月份的酒店预订是最繁忙的？

其次是商家更容易关心的问题：

• 客人来自哪里？
• 按市场细分的不同预定情况是怎样的？
• 对酒店提出特殊要求的人有什么共同点？
• 什么样的人更容易取消预订？

首先我们调用**info()**函数来对数据进行初步的观察。代码和结果如下：

我们可以看到，公司这列有明显的缺失情况。与此同时，在这里我们也能了解到不同的数据类型。

📑缺失值观察及处理

🔖缺失值观察以及可视化

首先提取所有缺失的列以及它们缺失的个数情况，代码如下：

data_missing = data.isnull().sum()
data_missing = data_missing[data_missing > 0]
data_missing

结果如下：

接下来，我们对数据的缺失值进行一些简单的柱状图可视化。代码如下：

data_missing.plot.bar()

🔖缺失值处理

我们首先来看一下缺失值的实际意义代表情况，如下：

基于上述特征的实际意义，我们对其进行如下的处理：

• children ：缺失数据相对较少，我们选择删除携带孩子缺失的数据行。
• country：因为后续我们可能会用到该数据，所以设置为‘Unknown’。
• agent：预定的ID在本次数据分析中，不太重要，删除该特征。
• company：公司的数据在本次数据中，不太重要，而且确实程度过大，直接删除。

相关代码如下：

# inplace表示是否在原数据进行填充
data.dropna(axis=0, how='any',subset=['children'], inplace=True)
# 用Unknown填充'country'为Nan的数据
data.fillna({'country':'Unknown'}, inplace=True)
# 删除'agent','company'的列
data.drop(['agent','company'],axis=1, inplace=True)

📖用户数据探索

在进行了用户数据的基础探索后，接下来就开始对其进行一些分析工作。

📑什么时间预定酒店将会更经济实惠？

作为游客或者用户，我们较为关注的通常是酒店是否经济实惠。在下面的工作中，我们首先提取没有取消的那些订单作为我们分析的数据，然后对其进行一些处理最后进行可视化，代码如下：

import seaborn as sns
# 提取没有取消订单的数据
data_no_canceled = data[data['is_canceled']==0]
# 算出人均价格
data_no_canceled['adr_deal'] = data_no_canceled['adr']/ (data_no_canceled['adults'] + data_no_canceled['children'])
# 对结果进行可视化
ax = sns.lineplot(data=data_no_canceled, x = 'arrival_date_month',y='adr_deal',hue='hotel')
# 设置图片尺寸
ax.figure.set_size_inches(12,6)

结果如下：

在图中我们可以发现，在七月Resort Hotel价格最高，而在其他时间City Hote的价格会稍高一些。另外，在十一月以及十二月两种酒店的价格是全年最低的，也就是最经济实惠的。

另外，我还探究了不同房型的价格分布情况，代码如下：

# 提取数据并进行排序
roomtype = data_no_canceled[["hotel", "reserved_room_type", "adr_deal"]].sort_values("reserved_room_type")
# 绘制箱线图
ax_box = sns.boxplot(data = roomtype, x = 'reserved_room_type',y='adr_deal',hue='hotel')
ax_box.figure.set_size_inches(12,6)

结果如下：

可以看到，在图中City Hotel的E房型价格偏高，而Resort Hotel的F房型价格偏高。

📑哪个月份的酒店预订是最繁忙的？

基于上述的分析，我们进一步探索哪个月份的酒店预订是最繁忙的，代码如下：

# 提取用户们到达的日期，并对其进行排序可视化
data_no_canceled['arrival_date_month'].value_counts().sort_values(ascending=False).plot.bar()

结果如下：

在图中可以看到，5，7，8月是最繁忙的时间，而秋冬季节的酒店人数相对来说较少。

📖商家数据探索

作为商家，较为关心的问题即用户的组成，以便更好的提供服务。另外，做为商家了解不同用户的订购渠道对于进一步的进行市场营销将会有积极的作用。

📑按市场细分的不同预定情况是怎样的？

我们首先基于数据来探索预订酒店的不同国家的人数占比，并进行可视化，代码如下：

import plotly.express as px
# 计算不同国家的百分比
country_data = pd.DataFrame(data['country'].value_counts())
country_data.columns = ['guest_num']
country_data['guest_persent'] = round(country_data['guest_num'] / country_data['guest_num'].sum() * 100,2)

# 对小于一定比例的国家归类
country_data.loc['OTHER','guest_num'] = country_data[country_data['guest_persent'] < 2]['guest_num'].sum()
country_data.loc['OTHER','guest_persent'] = round(country_data.loc['OTHER','guest_num'] / country_data['guest_num'].sum(),2)
# 可视化
country_data.drop(country_data[country_data['guest_persent'] < 0.5].index, inplace=True)
fig = px.pie(country_data,values = 'guest_persent',names = country_data.index)
fig.update_traces(textposition="inside", textinfo="value+percent+label")

结果如下：

进一步的，我们对不同预订酒店的渠道进行了统计，代码如下：

data_no_canceled['market_segment'].value_counts().plot.bar()

结果如下：

从图中可以看到，预订的主力军都是旅行社居多。另外，在线预定远比线下预订的数目要多很多。针对商家，可以对线上的广告推荐等进一步的进行策划营销。

📑什么样的人更容易取消预订？

最后，我们希望了解什么样的人更容易取消预订。商家可以针对这部分用户进行服务优化，并进一步的改善酒店的经营模式。在这里，我们准备基于数据进行机器学习模型的构建。

🔖数据编码

在这一步中，我们将对字符型的数据做编码处理，从而更好的适应不同模型，代码如下：

from sklearn import preprocessing
data_to_ml = data.copy()
# 因为我们预测的目标是是否取消了预订，所以删除这两列
data_to_ml.drop(['reservation_status','reservation_status_date'], axis=1, inplace=True)
# 对大部分特征进行编码处理
for col in ['hotel','arrival_date_month','meal','country','market_segment','distribution_channel','customer_type','deposit_type']:
    encoder = preprocessing.LabelEncoder()
    encoder.fit(data_to_ml[col])
    data_to_ml[f'{col}_labeled']  = encoder.transform(data_to_ml[col])
    data_to_ml.drop([col], axis=1, inplace=True)

在这里，我们有一个特殊的操作，即我们根据现有的特征提取了一个新的特征，即酒店预留的房间以及用户预定的是否是相同的。因为如果作为用户发现预定的房间和实际房间不同很可能会取消预订。具体代码如下：

import numpy as np
data_to_ml['is_reserved_assigned_equal'] = np.where(data_to_ml['reserved_room_type']==data_to_ml['assigned_room_type'],1,0)
data_to_ml.drop(['reserved_room_type','assigned_room_type'], axis=1, inplace=True)

🔖特征筛选

在这一步中，我们对所有处理后的数据进行了皮尔逊相关分析并做了可视化，代码如下所示：

import matplotlib.pylab as plt
correlation = data_to_ml.corr('pearson')

f, ax = plt.subplots(figsize = (9, 9))
plt.title('Correlation of Numeric Features with',y=1,size=16)
sns.heatmap(correlation,square = True,  vmax=0.8)

结果如下：

可以看到，我们所挑选的特征自相关性不算太高可以接受。与此同时，其针对我们所选定的目标is_canceled来说相关性也是可以接受的，所以我们暂时不对其进行处理。

🔖构建模型并预测

在这一部分中，我们构建了主流的模型并使用10折交叉认证对其进行验证，代码如下：

from sklearn import svm
from sklearn.ensemble import RandomForestClassifier
from sklearn.neighbors import KNeighborsClassifier
from sklearn import model_selection

x_all = data_to_ml.drop('is_canceled',axis=1)
y_all = data_to_ml['is_canceled']

# 为了避免训练时间过长，这里我选择了部分数据集
x_train = x_all.loc[:19999,:]
y_train = y_all[:20000]

model_svm = svm.SVC()
model_knn = KNeighborsClassifier()
model_rf = RandomForestClassifier()

model_dict = {
    'SVM':model_svm,
    'KNeighborsClassifier':model_knn,
    'RandomForestClassifier':model_rf
}
# 训练模型
for model in model_dict:
    model_dict[model].fit(x_train, y_train)
    scores = model_selection.cross_val_score(model_dict[model], X=x_train, y=y_train, verbose=1, cv = 10, scoring='f1')
    print(model, scores.mean())

>>>SVM 0.6480724503120129
>>>KNeighborsClassifier 0.6577070795293697
>>>RandomForestClassifier 0.7363861778481174

可以看到，随机森林的效果是可以接受的，接下来我们对其特征重要性进行可视化。

🔖根据特征重要性得出结论

调用随机森林中自带的特征重要性函数，并对其进行了可视化，代码如下：

# 构建Series让特征重要性和特征名称一一对应
feature_importances_series = pd.Series(list(model_rf.feature_importances_), index = x_all.columns, )
# 对其进行排序
feature_importances_series =feature_importances_series.sort_values(ascending=False)
# 进行可视化
sns.barplot(x = feature_importances_series.values, y = feature_importances_series.index, orient='h')

结果如下：

如图所示，用户取消预订和逗留时长、用户国籍、到达年份都有很强的相关性。另外，我们提取的特征也是非常重要的，即用户预定的房间和实际的房间不同时，用户将会更容易取消预订。

📍总结

在本文中，我们基于python对酒店预订需求进行分析，并从多种角度对其展开了探索性的工作。这对于养成数据分析习惯有很大的帮助，在实际工作或者学习中还需要不断练习。

感兴趣的朋友们可以自己按照上述步骤进行操作，或在评论区与我讨论。

需要源码的朋友可以私信我进行索取，我们下次再见。