本文将使用Python的绘图库Matplotlib,通过数据可视化分析影响房价的主要因素。
波士顿房价数据集
波士顿房价数据集来源于1978年美国某经济学杂志上。该数据集包含若干波士顿房屋的价格及其各项数据,每个数据项包含14个数据,分别是犯罪率、是否在河边和平均房间数等相关信息,其中最后一个数据是房屋中间价。
变量名称解释
CRIM: per capita crime rate by town 每个城镇人均犯罪率
ZN: proportion of residential land zoned for lots over 25,000 sq.ft. 超过25000平方英尺用地划为居住用地的百分比
INDUS: proportion of non-retail business acres per town 非零售商用地百分比
CHAS: Charles River dummy variable (= 1 if tract bounds river; 0 otherwise) 是否靠近查尔斯河
NOX: nitric oxides concentration (parts per 10 million) 氮氧化物浓度
RM: average number of rooms per dwelling 住宅平均房间数目
AGE: proportion of owner-occupied units built prior to 1940 1940年前建成自用单位比例
DIS: weighted distances to five Boston employment centres 到5个波士顿就业服务中心的加权距离
RAD: index of accessibility to radial highways 无障碍径向高速公路指数
TAX: full-value property-tax rate per $10,000 每万元物业税率
PTRATIO: pupil-teacher ratio by town 小学师生比例
B: 1000(Bk - 0.63)^2 where Bk is the proportion of blacks by town 黑人比例指数
LSTAT: % lower status of the population 下层经济阶层比例
MEDV: Median value of owner-occupied homes in $1000's 业主自住房屋中值
数据导入和查看
波士顿房价数据集在机器学习的sklearn.datasets包中,用datasets.load_boston()导入数据。接着使用.data和.target分别加载房屋属性参数值和房屋中间价。
import pandas as pd
from sklearn import datasets
boston_dataset = datasets.load_boston()
X_full = boston_dataset.data
Y = boston_dataset.target
用房屋属性参数值创建DataFrame,使用.columns修改行标签为房屋属性名称,并添加一列房屋中间价数据。分别使用.head()和.info()查看数据,可以看到数据导入和整理成功,便于下一步的可视化分析。
boston = pd.DataFrame(X_full)
boston.columns = boston_dataset.feature_names
boston['PRICE'] = Y
boston.head()
boston.info()
RangeIndex: 506 entries, 0 to 505
Data columns (total 14 columns):
CRIM 506 non-null float64
ZN 506 non-null float64
INDUS 506 non-null float64
CHAS 506 non-null float64
NOX 506 non-null float64
RM 506 non-null float64
AGE 506 non-null float64
DIS 506 non-null float64
RAD 506 non-null float64
TAX 506 non-null float64
PTRATIO 506 non-null float64
B 506 non-null float64
LSTAT 506 non-null float64
PRICE 506 non-null float64
dtypes: float64(14)
memory usage: 55.4 KB
数据可视化分析
从上文的数据可以看出,接下来要分析13个房屋属性参数与房屋中间价的关系。房屋属性参数中有一个数据是离散数据(CHAS,0 或者 1),其余的均为连续数据。
首先,使用matplotlib绘图分析河边CHAS(离散数据)和价格PRICE是否有相关性,可以看到二者没有明显的相关性,接下来的数据分析可以剔除属性CHAS。
import matplotlib.pyplot as plt
%matplotlib inline
%config InlineBackend.figure_format = 'retina'
plt.style.use('seaborn-whitegrid')
plt.scatter(boston.CHAS, boston.PRICE)
plt.xlabel('CHAS')
plt.ylabel('PRICE')
plt.show()
接下来,需要分析12个连续数据与房屋中间价的相关性。因为变量比较多,为了迅速找出影响房屋中间价的因素,可以使用matplotlib中的Seaborn模块绘制多变量图。多变量图的对角线上是变量分布的直方图,非对角线上是两个变量的散点图。
把12个房屋属性分成两组与中间价绘图,第一组:CRIM, ZN, INDUS, NOX, RM 和 AGE。
从图中可以看出,CRIM(犯罪率)和 RM(平均房间数目)与PRICE(房屋中间价)有明显相关性,RM和 PRICE呈现正态分布。
import seaborn as sns
sns.pairplot(boston, vars=['CRIM', 'ZN', 'INDUS', 'NOX', 'RM', 'AGE', 'PRICE'])
第二组:DIS, RAD, TAX, PTRATIO, B 和 LSTAT。从图中可以看出,LSTAT(下层经济阶层比例)与PRICE有明显相关性。
sns.pairplot(boston, vars=['DIS', 'RAD', 'TAX', 'PTRATIO', 'B', 'LSTAT', 'PRICE'])
将CRIM,RM,LSTAT 与 PRICE绘图,可以更清楚的观察和分析数据的相关性。从图中看出 PRICE 与 LSTAT 及 CRIM 呈非线性关系,而与 RM 更呈线性关系。
sns.pairplot(boston, vars=['CRIM', 'RM', 'LSTAT', 'PRICE'])
总结:
通过数据可视化分析,获得波士顿地区自住房屋中间价与下层经济阶层比例、平均房间数目和犯罪率有明显的相关性。具体符合哪种相关性,进一步分析将会涉及到线性回归等模型的应用。