R语言SVM、决策树与因子分析对城市空气质量分类与影响因素可视化研究

数据处理和分析在数据科学领域中扮演着至关重要的角色。

由Kaizong Ye,Liao Bao撰写

确保数据的准确性和完整性是数据处理的首要任务。在本研究中,我们以空气质量数据为例,帮助客户进行了数据处理和分析。

根据空气质量的指标(x1~x7),我们将30个城市分成两类,并使用Y1来评估分类的效果。

为了便于分类,在本研究中,我们将使用决策树方法和支持向量机(SVM)方法来进行分类分析。

这两种方法在机器学习领域被广泛应用,能够有效地处理分类问题,并提高模型的准确性和泛化能力。

读取数据

首先,我们使用R语言中的read.csv函数来导入名为”air.csv”的数据集,并通过skip参数跳过第一行进行读取。

空气质量数据:

image.png

随后,使用head函数来查看数据的前几行,以初步了解数据的结构和内容。

data=read.csv("air.csv",skip = 1) head(data)
image.png

第一部分:

SVM模型可视化解释 | 支持向量机交互式演示

支持向量机 (SVM) 可视化解释

通过交互式动画理解SVM如何找到最佳分隔超平面,以及它在机器学习中的应用

开始探索

什么是支持向量机?

支持向量机(Support Vector Machine,简称SVM)是一种强大的监督学习模型,主要用于分类任务,也可用于回归分析。

SVM的核心思想是找到一个最佳分隔超平面,将不同类别的数据点分开,同时使两类数据点到超平面的最小距离最大化

距离超平面最近的那些数据点被称为支持向量,它们决定了超平面的位置和方向,是SVM模型的关键。

核心优势

  • 在高维空间中表现良好
  • 对小样本数据集效果显著
  • 通过核函数处理非线性问题

最佳分隔超平面

交互式SVM演示

点击图表添加数据点,观察SVM如何找到最佳分隔超平面。您也可以拖动已有点来改变它们的位置。

控制选项

较小的C值表示更强的正则化

演示说明

  • 点击图表区域添加数据点,先选择左侧的类别A或类别B
  • 拖动已有的数据点可以改变其位置,观察超平面的变化
  • 支持向量会以加粗的边框显示,它们是决定超平面的关键
  • 尝试不同的核函数,观察它们如何处理线性不可分的数据

核函数的作用

核函数是SVM处理非线性问题的关键,它能将低维空间中的非线性数据映射到高维空间,使其线性可分。

线性核

最简单的核函数,适用于线性可分的数据,计算效率高。

K(x, y) = x · y
适用场景: 线性可分数据,文本分类

RBF核

径向基函数核,适用于大多数非线性场景,是默认选择。

K(x, y) = exp(-γ||x-y||²)
适用场景: 非线性数据,图像识别

多项式核

将数据映射到多项式特征空间,适用于特定非线性问题。

K(x, y) = (x·y + c)^d
适用场景: 图像处理,文本分类

不同核函数处理非线性数据的效果

线性核

RBF核

多项式核

SVM的应用场景

支持向量机在多个领域都有广泛应用,尤其在小样本、高维特征的场景中表现出色。

图像识别

SVM在图像分类、人脸识别等领域表现优异,能够处理高维图像特征并实现精确分类。

文本分类

在垃圾邮件检测、情感分析和新闻分类等任务中,SVM能够有效处理高维文本特征向量。

生物信息学

SVM被用于基因分类、蛋白质结构预测等,能从复杂的生物数据中提取有效模式。

手写识别

在OCR技术中,SVM能有效识别手写字符,在数字识别等任务中准确率高。

异常检测

SVM可用于检测欺诈交易、网络入侵等异常行为,尤其适用于不平衡数据集。

回归分析

SVM不仅用于分类,也可扩展为支持向量回归(SVR),用于预测连续值输出。

数据处理 : 我们将Y1的指标转化为0和1,分别表示该城市的空气质量好和差。这种分类方法有助于后续的分析和建模。


 data$Y1<-ifelse(data$Y1>365*0.8,'空气好','空气差') 

构建因变量

data$Y1<-as.factor(data$Y1)  ## 将因变量格式转为因子型

视频

从决策树到随机森林:R语言信用卡违约分析信贷数据实例

探索见解

去bilibili观看

探索更多视频


视频

决策树模型原理和R语言预测心脏病实例

探索见解

去bilibili观看

探索更多视频


视频

因子分析简介及R语言应用实例

探索见解

去bilibili观看

探索更多视频


视频

支持向量机SVM、支持向量回归SVR和R语言网格搜索超参数优化实例

探索见解

去bilibili观看

探索更多视频


视频

支持向量机算法原理

探索见解

去bilibili观看

探索更多视频

构建训练集、测试集

在编程的过程中,我们需要确保每个步骤都得到充分的考虑和完善。从数据预处理、特征工程、数据可视化到建模、验证和优化,每个环节都至关重要,不能有任何遗漏。此外,构建训练集和测试集也是非常关键的一步。我们采用了分层抽样的方法,将数据集分为70%的训练集和30%的测试集。通过这种划分方式,我们可以在训练集上建立模型,并在测试集上验证模型的准确性和泛化能力,从而评估模型的有效性和可靠性。


train<-data[trainindex, ]  ## 去除price变量的训练集 test<-data[-trainindex, ]  ## 去除price变量的测试集

训练svm模型

一、使用线性核函数去拟合SVM模型

在训练SVM模型的过程中,我们首先使用线性核函数进行拟合。

1)模型拟合


,data=train,kernel='linear',              cost=10,scale=F) #kernel='linear'  ## 选择线性核函数 #scale=F  ## 对数据不进行标准化处理,支持向量机的的损失函数为凸函数,是否标准化不影响最优解,但标准化之后可以使求解速度变快 #cost=10  ## 参数代表犯错的成本,越大模型对误差的惩罚越大,生成的分类边界越复杂

在上述代码中,我们使用svm函数拟合了一个SVM模型,其中指定了使用线性核函数(kernel='linear')进行分类。参数cost=10表示对误差的惩罚程度,这个值越大,模型对误差的惩罚越大,生成的分类边界也会更复杂。同时,我们选择不对数据进行标准化处理(scale=FALSE),因为SVM的损失函数是凸函数,标准化数据不会影响最优解,但可以加快求解速度。


 summary ( svmfit1 )

image.png 在上述代码中,我们使用summary函数对拟合的SVM模型进行了摘要。该摘要包含了模型的关键参数和性能指标,如下所示:


  • SVM-Type: C-classification
  • SVM-Kernel: linear
  • cost: 10
  • gamma: 0.02173913
  • Support Vectors数量: 8

从摘要中可以看出,我们拟合的SVM模型是一个C-classification类型的模型,使用了线性核函数,cost参数为10。此外,模型中有8个支持向量,这些支持向量在决定分类边界时起到关键作用。最后,模型中共有2个类别,分别为"空气差"和"空气好"。

2)对svmfit1模型进行改进,选择最优的cost值

在对已拟合的svmfit1模型进行改进时,我们选择了最优的cost值。以下是具体步骤和结果的解释:


    ranges =list(cost=c(0.001 ,                                    0.01, 0.1, 1,5,10,100) )) summary (tune.out )
image.png

在上述代码中,我们使用tune函数对svm模型进行参数调优,通过10折交叉验证的方式选择最优的cost值。summary函数用于查看调优结果摘要,包括最佳参数和性能指标。


图片

数据分享|R语言分析上海空气质量指数数据:kmean聚类、层次聚类、时间序列分析:arima模型、指数平滑法

阅读文章


从上述结果中可以看出,经过参数调优后,最佳的cost值为0.1,对应的误差率为0.05。在调优过程中,我们对不同的cost值进行了评估,并选择了性能最优的参数值。


随时关注您喜欢的主题



summary ( best.mode1 )
image.png

在最优模型摘要中,我们可以看到调优后的最佳模型参数为C-classification类型的线性SVM模型,使用线性核函数,cost参数为0.1。模型中共有11个支持向量,共包含2个类别:"空气差"和"空气好"。这些结果表明通过参数调优,我们成功选择了最优的cost值,优化了SVM模型的性能。

3)模型评估

在对模型进行评估的过程中,我们首先对SVM模型进行评估,然后训练决策树模型。


ypred<-predict(
image.png

在上述代码中,我们使用predict函数对训练好的SVM模型进行预测,得到了预测结果。这些结果可以与真实值进行对比,从而评估模型的性能。以上结果展示了SVM模型的预测情况,包括真实值和预测值的对比,可以进一步分析模型的准确性和性能。

训练决策树模型

接下来,我们训练决策树模型,以下是具体步骤和结果的解释:


.-label-Y1,              method="class", dat
image.png

想了解更多关于模型定制、咨询辅导的信息?


在上述代码中,我们使用rpart函数训练了一个决策树模型,其中method="class"表示进行分类分析。模型使用了数据集中除了label和Y1以外的其他变量作为预测因子。 通过plotcp函数可视化交叉验证结果,进一步评估模型的性能。以下是可视化结果:


plotcp(fit) # visualize cross-validation results
image.png
summary(fit) # 


image.png 最后,利用summary函数对训练的决策树模型进行摘要,包括模型的性能指标、变量重要性以及节点信息。通过摘要结果,我们可以深入了解模型的构建过程和性能表现,为进一步的模型评估和优化提供参考。

2.png

第二部分:

根据第一部分的结果分成的两类城市来, 对第一类城市(空气好的城市):先根据城市发展指标(x8~x15)做因子分析,再将分出的因子和AQI值(y2)做对应分析,来分析它们的相关关系;

data1=data[data$Y1=="空气好",]

对第一类城市(空气好的城市)进行因子分析和相关关系分析的步骤如下:

因子分析

画出协方差阵和相关系数矩阵

cov(data1[,-c(1,9:10)]  )
image.png

协方差阵展示了不同变量之间的协方差关系,可以帮助我们了解变量之间的线性关系。

cor(data1[,-c(1,9:10)]  )
image.png

相关系数矩阵显示了各个变量之间的相关性程度,可以帮助我们理解变量之间的相关关系。

采用斜交旋转提取因子

a1[ ,-c(1,9,10,11)]), nfactors=8, rotate=  cor(data1$Y2 ,fm$scores)
image.png

在这一步中,我们采用斜交旋转提取因子的方法,将城市发展指标(x8~x15)进行因子分析,并提取8个因子。接着,我们将提取的因子与AQI值(Y2)进行对应分析,以探究它们之间的相关关系。

以上代码展示了AQI值(Y2)与提取的因子之间的相关系数,帮助我们分析城市发展指标与空气质量之间的关联关系。

对第二类城市(即空气质量较差的城市)的研究中,我们首先针对城市发展指标(x8~x15)进行了因子分析,以探究这些指标之间的潜在结构。随后,我们将因子分析得到的因子与空气质量指数(AQI值,即y2)进行了对应分析,以揭示它们之间的相关关系。

首先,我们从原始数据集中筛选出空气质量较差的城市数据

  data1=data[data$Y1=="空气差",]

因子分析

接着,我们对筛选出的城市发展指标(x8~x15)进行了因子分析。为了初步了解这些指标之间的关联性,我们计算了它们的协方差矩阵。协方差矩阵能够展示各指标之间的变异程度以及它们之间的线性相关程度。

cov(data1[,-c(1,9:10)]  ) cor(data1[,-c(1,9:10)]  )
image.png
image.png

采用斜交旋转提取因子

对第二类城市(即空气质量较差的城市)进行深入研究时,我们首先采用因子分析方法来探究城市发展指标(x8~x15)之间的潜在结构。为了更清晰地解释因子并使其具有实际意义,我们使用了斜交旋转(特别是“varimax”旋转)来提取因子。这种方法有助于我们理解各个因子所代表的原始指标组合,以及这些因子在解释城市发展特征时的相对重要性。

在因子分析过程中,我们设定了提取8个因子的目标,这是因为我们希望找到能够代表原始指标中大部分信息的少数几个因子。通过斜交旋转,我们得到了旋转后的因子载荷矩阵,这有助于我们解释每个因子所代表的含义。

随后,为了探究这些因子与空气质量指数(AQI值,即y2)之间的相关关系,我们计算了因子得分与AQI值之间的相关系数。相关系数矩阵显示了每个因子与AQI值之间的线性相关程度。

ctors=8, rotate="varimax cor(data1$Y2 ,fm$scores)
image.png

从相关系数矩阵中,我们可以看到MR5因子与AQI值之间存在较高的正相关关系(相关系数为0.8802801),这意味着该因子所代表的城市发展特征与空气质量较差有显著的正向关联。同时,MR6因子也与AQI值呈现出一定的正相关(相关系数为0.3007104),表明该因子同样与较差的空气质量有关。

其他因子与AQI值之间的相关系数虽然较低,但也可能存在某种关联。这些结果为我们提供了关于城市发展指标与空气质量之间关系的初步线索,有助于我们进一步理解和分析这些城市在发展过程中所面临的空气质量挑战。


可下载资源

关于作者

Kaizong Ye拓端研究室(TRL)的研究员。在此对他对本文所作的贡献表示诚挚感谢,他在上海财经大学完成了统计学专业的硕士学位,专注人工智能领域。擅长Python.Matlab仿真、视觉处理、神经网络、数据分析。

本文借鉴了作者最近为《R语言数据分析挖掘必知必会 》课堂做的准备。

​非常感谢您阅读本文,如需帮助请联系我们!

 
QQ在线咨询
售前咨询热线
15121130882
售后咨询热线
0571-63341498

关注我们,永远不要错过任何见解。


技术干货二维码

技术干货

最新洞察二维码

最新洞察

视频号二维码

视频号

This will close in 0 seconds